Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	402,167	G→A	S136N (AGT→AAT)	SAUSA300_RS01860 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	402,167	0	G	A	100.0%	114.1 / NA	42	S136N (AGT→AAT)	SAUSA300_RS01860	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (40/2);  total (40/2)

CTGCAGACAGTTATTGCTTTTTTAGTATTGCTTATTTTTGTGGTAGGGCACGCTGCAGCAAATGCAATACTATTCTTTTTAATACCATTTTTTGCTAGTGCTATTTTCAATACACAATTTACACTGTTTAATAGAAGGTTTGATGACAGAATACCAAAAAT  >  NC_007793/402070‑402230                                                                                                  |                                                                cTGCAGACAGTTATTGCTTTTTTAGTATTGCTTATTTTTGTGGTAGGGCACGCTGCAGCAAATGCAATACTATTCTTTTTAATACCATTTTTTGCTAATg                                                               >  1:2898594/1‑100 (MQ=255)   gCAGACAGTTATTGCTTTTTTAGTATTGCTTATTTTTGTGGTAGGGCACGCTGCAGCAAATGCAATACTATTCTTTTTAATACCATTTTTTGCTAATGCt                                                             >  1:3632481/1‑100 (MQ=255)   gCAGACAGTTATTGCTTTTTTAGTATTGCTTATTTTTGTGGTAGGGCACGCTGCAGCAAATGCAATACTATTCTTTTTAATACCATTTTTTGCTAATGCt                                                             >  1:1332804/1‑100 (MQ=255)   gCAGACAGTTATTGCTTTTTTAGTATTGCTTATTTTTGTGGTAGGGCACGCTGCAGCAAATGCAATACTATTCTTTTTAATACCATTTTTTGCTAATGCt                                                             >  1:5058416/1‑100 (MQ=255)   gCAGACAGTTATTGCTTTTTTAGTATTGCTTATTTTTGCGGTAGGGCACGCTGCAGCAAATGCAATACTATTCTTTTTAATACCATTTTTTGCTAATGCt                                                             >  1:3908331/1‑100 (MQ=255)     aGACAGTTATTGCTTTTTTAGTATTGCTTATTTTTGTGGTAGGGCACGCTGCAGCAAATGCAATACTATTCTTTTTAATACCATTTTTTGCTAATGCTAt                                                           <  1:3003744/100‑1 (MQ=255)      gACAGTTATTGCTTTTTTAGTATTGCTTATTTTTGTGGTAGGGCACGCTGCAGCAAATGCAATACTATTCTTTTTAATACCATTTTTTGCTAATGCTAtg                                                          <  1:159825/100‑2 (MQ=255)           tgaTTGCTTTTTTAGTATTGCTTATTTTTGTGGTAGGGCACGCTGCAGCAAATGCAATACTATTCTTTTTAATACCATTTTTTGCTAATGCTATTTTCaa                                                     >  1:4053337/3‑100 (MQ=255)             aTTGCTTTTTTAGTATTGCTTATTTTTGTGGTAGGGCACGCTGCAGCAAATGCAATACTATTCTTTTTAATACCATTTTTTGCTAATGCTATTTTCAATa                                                   >  1:2235655/1‑100 (MQ=255)                             tGCTTATTTTTGTGGTAGGGCACGCTGCAGCAAATGCAATACTATTCTTTTTAATACCATTTTTTGCTAATGCTATTTTCAATACACAATTTACACTGtt                                   >  1:3615231/1‑100 (MQ=255)                             tGCTTATTTTTGTGGTAGGGCACGCTGCAGCAAATGCAATACTATTCTTTTTAATACCATTTTTTGCTAATGCTATTTTCAATACACAATTTACACTGtt                                   >  1:27886/1‑100 (MQ=255)                              gCTTATTTTTGTGGTAGGGCACGCTGCAGCAAATGCAATACTATTCTTTTTAATACCATTTTTTGCTAATGCTATTTTCAATACACAATTTACACTGttt                                  >  1:520855/1‑100 (MQ=255)                                     tttGTGGTAGGGCACGCTGCAGCAAATGCAATACTATTCTTTTTAATACCATTTTTTGCTAATGCTATTTTCAATACACAATTTACACTGTTTAATAGaa                           >  1:2526638/1‑100 (MQ=255)                                     tatgtgGTAGGGCACGCTGCAGCAAATGCAATACTATTCTTTTTAATACCATTTTTTGCTAATGCTATTTTCAATACACAATTTACACTGTTTAATAGaa                           >  1:1630004/3‑100 (MQ=255)                                      atgtgGTAGGGCACGCTGCAGCAAATGCAATACTATTCTTTTTAATACCATTTTTTGCTAATGCTATTTTCAATACACAATTTACACTGTTTAATAGAAg                          >  1:5459754/2‑100 (MQ=255)                                       tgtgGTAGGGCACGCTGCAGCAAATGCAATACTATTCTTTTTAATACCATTTTTTGCTAATGCTATTTTCAATACACAATTTACACTGTTTAATAGAAgg                         >  1:1221060/1‑100 (MQ=255)                                       tgtgGTAGGGCACGCTGCAGCAAATGCAATACTATTCTTTTTAATACCATTTTTTGCTAATGCTATTTTCAATACACAATTTACACTGTTTAATAGAAgg                         >  1:3662403/1‑100 (MQ=255)                                       tgtgGTAGGGCACGCTGCAGCAAATGCAATACTATTCTTTTTAATACCATTTTTTGCTAATGCTATTTTCAATACACAATTTACACTGTTTAATAGAAgg                         >  1:3076149/1‑100 (MQ=255)                                       tgtgGTAGGGCACGCTGCAGCAAATGCAATACTATTCTTTTTAATACCATTTTTTGCTAATGCTATTTTCAATACACAATTTACACTGTTTAATAGAAgg                         >  1:5578988/1‑100 (MQ=255)                                       tgtgGTAGGGCACGCTGCAGCAAATGCAATACTATTCTTTTTAATACCATTTTTTGCTAATGCTATTTTCAATACACAATTTACACTGTTTAATAGAAgg                         >  1:2179629/1‑100 (MQ=255)                                       tctgGTAGGGCACGCTGCAGCAAATGCAATACTATTCTTTTTAATACCATTTTTTGCTAATGCTATTTTCAATACACAATTTACACTGTTTAATAGAAgg                         >  1:558921/3‑100 (MQ=255)                                        gtgGTAGGGCACGCTGCAGCAAATGCAATACTATTCTTTTTAATACCATTTTTTGCTAATGCTATTTTCAATACACAATTTACACTGTTTAATAGAAGGt                        >  1:1469182/1‑100 (MQ=255)                                        gtgGTAGGGCACGCTGCAGCAAATGCAATACTATTCTTTTTAATACCATTTTTTGCTAATGCTATTTTCAATACACAATTTACACTGTTTAATAGAAGGt                        >  1:1201961/1‑100 (MQ=255)                                           ctAGGGCACGCTGCAGCAAATGCAATACTATTCTTTTTAATACCATTTTTTGCTAATGCTATTTTCAATACACAATTTACACTGTTTAATAGAAGGTTtg                     >  1:1332797/2‑100 (MQ=255)                                            tAGGGCACGCTGCAGCAAATGCAATACTATTCTTTTTAATACCATTTTTTGCTAATGCTATTTTCAATACACAATTTACACTGTTTAATAGAAGGTTtga                    >  1:1905990/1‑100 (MQ=255)                                            caGGGCACGCTGCAGCAAATGCAATACTATTCTTTTTAATACCATTTTTTGCTAATGCTATTTTCAATACACAATTTACACTGTTTAATAGAAGGTTtga                    >  1:2805329/2‑100 (MQ=255)                                             aGGGCACGCTGCAGCAAATGCAATACTATTCTTTTTAATACCATTTTTTGCTAATGCTATTTTCAATACACAATTTACACTGTTTAATAGAAGGTTtgat                   >  1:3187588/1‑100 (MQ=255)                                             aGGGCACGCTGCAGCAAATGCAATACTATTCTTTTTAATACCATTTTTTGCTAATGCTATTTTCAATACACAATTTACACTGTTTAATAGAAGGTTtgat                   >  1:4563294/1‑100 (MQ=255)                                               ggCACGCTGCAGCAAATGCAATACTATTCTTTTTAATACCATTTTTTGCTAATGCTATTTTCAATACACAATTTACACTGTTTAGTAGAAgatcggaaga                 >  1:1868490/1‑91 (MQ=255)                                               ggCACGCTGCAGCAAATGCAATACTATTCTTTTTAATACCATTTTTTGCTAATGCTATTTTCAATACACAATTTACACTGTTTACTAGAAGGTTtgatga                 >  1:1462009/1‑100 (MQ=255)                                                gCACGCTGCAGCAAATGCAATACTATTC‑TTTTAATACCATTTTTTGCTAATGCTATTTTCAATACACAATTTACACTGTTTAATAGAAGGTTTGATGACa               >  1:163756/1‑100 (MQ=255)                                                 tgcGCTGCAGCAAATGCAATACTATTCTTTTTAATACCATTTTTTGCTAATGCTATTTTCAATACACAATTTACACTGTTTAATAGAAGGTTTGATGACa               >  1:2004080/3‑100 (MQ=255)                                                  ccGCTGCAGCAAATGCAATACTATTCTTTTTAATACCATTTTTTGCTAATGCTATTGTCAATACACAATTTACACTATTTAATAGAGGGATTGATGACAg              >  1:1669112/2‑100 (MQ=255)                                                  aCGCTGCAGCAAATGCAATACTATTCTTTTTAATACCATTTTTTGCTAATGCTATTTTCAATACACAATTTACACTGTTTAATAGAAGGTTTGATGACAg              >  1:4574646/1‑100 (MQ=255)                                                   cGCTGCAGCAAATGCAATACTATTCTTTTTAATACCATTTTTTGCTAATGCTATTTTCAATACACAATTTACACTGTTTAATAGAAGGTTTGATGACAGa             >  1:4739250/1‑100 (MQ=255)                                                   cGCTGCAGCAAATGCAATACTATTCTTTTTAATACCATTTTTTGCTAATGCTATTTTCAATACACAATTTACACTGTTTAATAGAAGGTTTGATGACAGa             >  1:4808882/1‑100 (MQ=255)                                                    ccTGCAGCAAATGCAATACTATTCTTTTTAATACCATTTTTTGCTAATGCTATTTTCAATACACAATTTACACTGTTTAATAGAAGGTTTGATGACAGaa            >  1:2440500/2‑100 (MQ=255)                                                      tGCAGCAAATGCAATACTATTCTTTTTAATACCATTTTTTGCTAATGCTATTTTCAATACACAATTTACACTGTTTAATAGAAGGTTTGATGACAGAATa          >  1:2965159/1‑100 (MQ=255)                                                       gcagcaAATGCAATACTATTCTTTTTAATACCATTTTTTGCTAATGCTATTTTCAATACACAATTTACACTGTTTAATAGAAGGTTTGATGACAGAATAc         >  1:959029/1‑100 (MQ=255)                                                       gcagcaAATGCAATACTATTCTTTTTAATACCATTTTTTGCTAATGCTATTTTCAATACACAATTTACACTGTTTAATAGAAGGTTTGATGACAGAATAc         >  1:4378472/1‑100 (MQ=255)                                                        cagcaAATGCAATACTATTCTTTTTAATACCATTTTTTGCTAATGCTATTTTCAATACACAATTTACACTGTTTAATAGAAGGTTTGATGACAGAATAcc        >  1:4699143/1‑100 (MQ=255)                                                              aTGCAATACTATTCTTTTTAATACCATTTTTTGCTAATGCTATTTTCAATACACAATTTACACTGTTTAATAGACGGTTTGATGACAGAATACCAAAAAt  >  1:3102210/1‑100 (MQ=255)                                                                                                  |                                                                CTGCAGACAGTTATTGCTTTTTTAGTATTGCTTATTTTTGTGGTAGGGCACGCTGCAGCAAATGCAATACTATTCTTTTTAATACCATTTTTTGCTAGTGCTATTTTCAATACACAATTTACACTGTTTAATAGAAGGTTTGATGACAGAATACCAAAAAT  >  NC_007793/402070‑402230

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑