Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	345,649	G→A	*82* (TAG→TAA)	yahM →	uncharacterized protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	345,649	0	G	A	100.0%	105.5 / NA	33	*82* (TAG→TAA)	yahM	uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (17/16);  total (17/16)

AGCCCACACAGGGGAGAGTGGGCTGAAATGGGAAGCTAAAGACTCAAGTAAACTTATCGGAAATAAGGACCACGCATTACGGGGGCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAGAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAATCAGGCGTAGCGCGAATACACCAATAATTGCACCAATAACCCGACTGGCAACGCGTTGCATACGCCCATAAGCACGACGCACAGCGGG  >  NC_000913/345517‑345779                                                                                                                                     |                                                                                                                                   aGCCCACACAGGGGAGAGTGGGCTGAAATGGGAAGCTAAAGACTCAAGTAAACTTATCGGAAATAAGGACCACGCATTACGGGGGCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTa                                                                                                                              >  1:770336/1‑139 (MQ=255)   cccACACAGGGGAGAGTGGGCTGAAATGGGAAGCTAAAGACTCAAGTAAACTTATCGGAAATAAGGACCACGCATTACGGGGGCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGt                                                                                                                               <  2:549601/136‑1 (MQ=255)   cccACACAGGGGAGAGTGGGCTGAAATGGGAAGCTAAAGACTCAAGTAAACTTATCGGAAATAAGGACCACGCATTACGGGGGCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGt                                                                                                                               >  1:549601/1‑136 (MQ=255)                 agTGGGCTGAAATGGGAAGCTAAAGACTCAAGTACACTTATCGGAAATAAGGACCCCGCATTACGGGGGCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGATGATTACGTAAAGAGTAAATTATCTGCTCACCg                                                                                                              >  1:15336/1‑139 (MQ=255)                  gTGGGCTGAAATGGGAAGCTAAAGACTCAAGTAAACTTATCGGAAATAAGGACCACGCATTACGGGGGCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGc                                                                                                             >  2:49828/1‑139 (MQ=255)                      gCTGAAATGGGAAGCTAAAGACTCAAGTAAACTTATCGGAAATAAGGACCACGCATTACGGGGGCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCg                                                                                                         >  1:448759/1‑139 (MQ=255)                                gAAGCTAAAGACTCAAGTAAACTTATCGGAAATAAGGACCACGCATTACGGGGGCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTc                                                                                               >  1:65736/1‑139 (MQ=255)                                      aaaGACTCAAGTAAACTTATCGGAAATAAGTACCACGCATTACGGGGGCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAAt                                                                                         <  2:1038771/139‑1 (MQ=255)                                           cTCAAGTAAACTTATCGGAAATAAGGACCACGCATTACGGGGGCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAATCAGGc                                                                                    <  1:1190430/139‑1 (MQ=255)                                                      ttATCGGAAATAAGGACCACGCATTACGGGGGCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAATCAGGCGTAGCGCGAAt                                                                         <  1:49828/139‑1 (MQ=255)                                                      ttATCGGAAATAAGGACCACGCATTACGGGGGCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAATCAGGCGTAGCGCGAAt                                                                         <  1:1228383/139‑1 (MQ=255)                                                      ttATCGGAAATAAGGACCACGCATTACGGGGGCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAATCAGGCGTAGCGCGAAt                                                                         <  2:448759/139‑1 (MQ=255)                                                       tATCGGAAATAAGGACCACGCATTACGAGGGCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAATCAGGCGTAGCGCGAATa                                                                        <  1:786993/139‑1 (MQ=255)                                                        aTCGGAAATAAGGACCACGCATTACGGGGGCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATT‑TCTGCTCACCGCTGCGTCACCCCTTCGTAAATCAGGCGTAGCGCGAATaca                                                                      >  2:534830/1‑139 (MQ=255)                                                        aTCGGAAATAAGGACCACGCATTACGGGGGCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATT‑TCTGCTCACCGCTGCGTCACCCCTTCGTAAATCAGGCGTAGCGCGAATaca                                                                      >  2:450477/1‑139 (MQ=255)                                                        aTCGGAAATAAGGACCACGCATTACGGGGGCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATT‑TCTGCTCACCGCTGCGTCACCCCTTCGTAAATCAGGCGTAGCGCGAATaca                                                                      >  2:1342622/1‑139 (MQ=255)                                                          cGGAAATAAGGACCACGCATTACGGGGGCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATT‑TCTGCTCACCGCTGCGTCACCCCTTCGTAAATCAGGCGTAGCGCGAATACAcc                                                                    <  1:534830/139‑1 (MQ=255)                                                          cGGAAATAAGGACCACGCATTACGGGGGCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATT‑TCTGCTCACCGCTGCGTCACCCCTTCGTAAATCAGGCGTAGCGCGAATACAcc                                                                    <  1:1342622/139‑1 (MQ=255)                                                                  aGGACCACGCATTACGGGGGCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGt                                                                                                        <  1:414036/96‑1 (MQ=255)                                                                  aGGACCACGCATTACGGGGGCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGt                                                                                                        >  2:414036/1‑96 (MQ=255)                                                                        aCGCATTACGGGGGCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTAt                                                                                                                        >  2:452699/1‑74 (MQ=255)                                                                        aCGCATTACGGGGGCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTAt                                                                                                                        <  1:452699/74‑1 (MQ=255)                                                                         cGCATTACGGGGGCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGACACCCCTTCGTAAATCAGGCGTAGCGCGAATACACCAATAATTGCACCaa                                                      >  2:867397/1‑139 (MQ=255)                                                                          gCATTACGGGGGCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAATCAGGCGTAGCGCGAATACACCAATAATTGCACCAAt                                                     >  2:762727/1‑139 (MQ=255)                                                                                  ggggCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAATCAGGCGTAGCGCGAATACACCAATAATTGCACCAATAACCCGAc                                             <  2:770336/139‑1 (MQ=255)                                                                                   gggCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAATCAGGCGTAGCGCGAATACACCAATAATTGCACCAATAACCCGACt                                            >  2:1134080/1‑139 (MQ=255)                                                                                     gCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATT‑TCTGCTCACCGCTGCGTCACCCCTTCGTAAATCAGGCGTAGCGCGAATACACCAATAATTGCACCAATAACCCGACTGGc                                         <  1:450477/139‑1 (MQ=255)                                                                                            tcGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAATCAGGCGTAGCGCGAATACACCAATAATTGCACCAATAACCCGACTGGCAAcgcg                                   <  1:603613/139‑1 (MQ=255)                                                                                              gCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAATCAGGCGTAGCGCGAATACACCaataa                                                               >  1:627794/1‑109 (MQ=255)                                                                                              gCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAATCAGGCGTAGCGCGAATACACCaataa                                                               <  2:627794/109‑1 (MQ=255)                                                                                                   tGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAATCAGGCGTAGCGCGAATACACCAATAATTGCACCAATAACCCGACTGGCAACGCGTTGCATa                            <  1:762727/139‑1 (MQ=255)                                                                                                               ggCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAATCAGGCGTAGCGCGAATACACCAATAATTGCACCAATAACCCGACTGGCAACGCGTTGCATACGCCCATAAGCa                >  1:1078208/1‑139 (MQ=255)                                                                                                                             aTTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAATCAGGCGTAGCGCGAATACACCAATAATTGCACCAATAACCCGACTGGCAACGCGTTGCATACGCCCATAAGCACGACGCACAGCggg  >  1:781116/1‑139 (MQ=255)                                                                                                                                     |                                                                                                                                   AGCCCACACAGGGGAGAGTGGGCTGAAATGGGAAGCTAAAGACTCAAGTAAACTTATCGGAAATAAGGACCACGCATTACGGGGGCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAGAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAATCAGGCGTAGCGCGAATACACCAATAATTGCACCAATAACCCGACTGGCAACGCGTTGCATACGCCCATAAGCACGACGCACAGCGGG  >  NC_000913/345517‑345779

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