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%	149.3 / NA	46	*82* (TAG→TAA)	yahM	uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (23/23);  total (23/23)

AGCCCACACAGGGGAGAGTGGGCTGAAATGGGAAGCTAAAGACTCAAGTAAACTTATCGGAAATAAGGACCACGCATTACGGGGGCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAGAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAATCAGGCGTAGCGCGAATACACCAATAATTGCACCAATAACCCGACTGGCAACGCGTTGCATACGCCCATAAGCACGACGCACAGCGGG  >  NC_000913/345517‑345779                                                                                                                                     |                                                                                                                                   aGCCCACACAGGGGAGAGTGGGCTGAAATGGGAAGCTAAAGACTCAAGTAAACTTATCGGAAATAAGGACCACGCATTACGGGGGCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTa                                                                                                                              >  1:1380114/1‑139 (MQ=255)            gggAGAGTGGGCTGAAATGGGAAGCTAAAGACTCAAGTAAACTTATCGGAAATAAGGACCACGCATTACGGGGGCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCt                                                                                                                   >  1:626143/1‑139 (MQ=255)                  gTGGGCTGAAATGGGAAGCTAAAGACTCAAGTAAACTTATCGGAAATAAGGACCACGCATTACGGGGGCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTaaag                                                                                                                                  >  1:250949/1‑118 (MQ=255)                                gAAGCTAAAGACTCAAGTAAACTTATCGGAAATAAGGACCACGCATTACGGGGGCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGt                                                                                                                               <  1:957418/107‑1 (MQ=255)                                gAAGCTAAAGACTCAAGTAAACTTATCGGAAATAAGGACCACGCATTACGGGGGCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGt                                                                                                                               >  2:957418/1‑107 (MQ=255)                                              aaGTAAACTTATCGGAAATAAGGACCACGCATTACGGGGGCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAATCAGGCGTa                                                                                 >  2:1473898/1‑139 (MQ=255)                                                    aCTTATCGGAAATAAGGACCACGCATTACGGGGGCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAATCAGGCGTAGCGCGa                                                                           >  2:34668/1‑139 (MQ=255)                                                     cTTATCGGAAATAAGGACCACGCATTACGGGGGCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAATCAGGCGTAGCGCGaa                                                                          >  1:1141955/1‑139 (MQ=255)                                                     cTTATCGGAAATAAGGACCACGCATTACGGGGGCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAATCAGGCGTAGCGCGaa                                                                          <  2:858851/139‑1 (MQ=255)                                                        aTCGGAAATAAGGACCACGCATTACGGGGGCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTg                                                                                                           <  2:322899/103‑1 (MQ=255)                                                        aTCGGAAATAAGGACCACGCATTACGGGGGCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTg                                                                                                           >  1:322899/1‑103 (MQ=255)                                                        aTCGGAAATAAGGACCACGCATTACGGGGGCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCAccc                                                                                                   <  2:1049982/111‑1 (MQ=255)                                                        aTCGGAAATAAGGACCACGCATTACGGGGGCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAATCAGGCGTAGCGCGAATac                                                                       >  1:219575/1‑139 (MQ=255)                                                        aTCGGAAATAAGGACCACGCATTACGGGCGCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTAATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCAccc                                                                                                   >  1:1049982/1‑111 (MQ=255)                                                           ggAAATAAGGACCACGCATTACGGGGGCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAATCAGGCGTAGCGCGAATACAcc                                                                    <  2:1380114/139‑1 (MQ=255)                                                                   ggACCACGCATTACGGGGGCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAATCAGGCGTAGCGCGAATACACCAATAATTg                                                            <  2:219575/139‑1 (MQ=255)                                                                       cACGCATTACGGGGTCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAATCAGGCGTAGCGCGAAt                                                                         <  2:661168/122‑1 (MQ=255)                                                                       cACGCATTACGGGGTCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAATCAGGCGTAGCGCGAAt                                                                         >  1:661168/1‑122 (MQ=255)                                                                         cGCATTACGGGGGCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCt                                                                                                                   <  2:467758/78‑1 (MQ=255)                                                                         cGCATTACGGGGGCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCt                                                                                                                   >  1:467758/1‑78 (MQ=255)                                                                                   gggCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAATCAGGCGTAGCGCGAATACACCAATAATTGCACCAATAACCCGACt                                            >  1:1403047/1‑139 (MQ=255)                                                                                     gCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAATCAGGCGTAGCGCGAATACACCAATAATTGCACCAATAACCCGACTgg                                          >  1:1277630/1‑139 (MQ=255)                                                                                       tATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAATCAGGCGTAGCGCGAATACACCAATAATTGCACCAATAACCCGACTGGCa                                        <  2:1141955/139‑1 (MQ=255)                                                                                       tATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAATCAGGCGTAGCGCGAATACACCAATAATTGCACCAATAACCCGACTGGCa                                        <  1:34668/139‑1 (MQ=255)                                                                                        atcatcGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAATCAGGCGTAGCGCGAATACACCAATAATTGCACCAATAACCCGACTGGCaa                                       >  1:441960/1‑139 (MQ=255)                                                                                             cGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCAcc                                                                                                    <  1:854479/73‑1 (MQ=255)                                                                                             cGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCAcc                                                                                                    >  2:854479/1‑73 (MQ=255)                                                                                                 gATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAATCAGGCGTAGCGCGAATACACCAATAATTGCACCAATAACCCGACTGGCAACGCGTTGCa                              <  2:626143/139‑1 (MQ=255)                                                                                                   tGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCAcc                                                                                                    <  2:213298/67‑1 (MQ=255)                                                                                                   tGGCCGTGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCAcc                                                                                                    >  1:213298/1‑67 (MQ=255)                                                                                                         tGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAATCAGGCGTAGCGCGaa                                                                          >  1:1590581/1‑87 (MQ=255)                                                                                                         tGATTCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAATCAGGCGTAGCGCGaa                                                                          <  2:1590581/87‑1 (MQ=255)                                                                                                            ttCGGCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAATCAGGCGTAGCGCGAATACACCAATAATTGCACCAATAACCCGACTGGCAACGCGTTGCATACGCCCATaa                   <  1:1473898/139‑1 (MQ=255)                                                                                                               ggCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAATCAGGCGTAGCGCGAATACACCAATAATTGCACCAATAACCCGACTGGCAACGCGTTGCATACGCCCATAAGCa                <  2:1403047/139‑1 (MQ=255)                                                                                                                gCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAATCAGGCGTAGCGCGAATACACCAATAATTGCACCAATAACCCGACTGGCAACGCGTTGCATa                            >  1:362222/1‑126 (MQ=255)                                                                                                                gCAAATCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAATCAGGCGTAGCGCGAATACACCAATAATTGCACCAATAACCCGACTGGCAACGCGTTGCATa                            <  2:362222/126‑1 (MQ=255)                                                                                                                    aTCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAATCAGGCGTAGCGCGAATACACCAATAATTGCACCAATAACCCGACTGGCAACGCGTTGCATa                            >  1:716873/1‑122 (MQ=255)                                                                                                                    aTCAGGTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAATCAGGCGTAGCGCGAATACACCAATAATTGCACCAATAACCCGACTGGCAACGCGTTGCATa                            <  2:716873/122‑1 (MQ=255)                                                                                                                         gTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAATCAGGCGTAGCGCGAATACACCAATAATTGCACCAATAACCCGACTg                                           <  1:1205416/102‑1 (MQ=255)                                                                                                                         gTTGATTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAATCAGGCGTAGCGCGAATACACCAATAATTGCACCAATAACCCGACTg                                           >  2:1205416/1‑102 (MQ=255)                                                                                                                             aTTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCAcc                                                                                                    <  1:1173447/41‑1 (MQ=255)                                                                                                                             aTTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCAcc                                                                                                    >  2:1173447/1‑41 (MQ=255)                                                                                                                             aTTACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAATCAGGCGTAGCGCGAATACACCAATAATTGCACCAATAACCCGACTGGCAACGCGTTGCATACGCCCATAAGCACGACGCACAGCggg  <  2:1277630/139‑1 (MQ=255)                                                                                                                              ttACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGt                                                                                             >  2:446854/1‑47 (MQ=255)                                                                                                                              ttACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGt                                                                                             <  1:446854/47‑1 (MQ=255)                                                                                                                              ttACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAATCAGGCGTAGCGCGAATACACCAATAATTGCACCAATAACCGGACTGGCAACGCGTTGCATACGCCCATAAGCa                >  1:1010762/1‑124 (MQ=255)                                                                                                                              ttACGTAAAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAATCAGGCGTAGCGCGAATACACCAATAATTGCACCAATAACCCGACTCGCAACGCGTTGCATACGCCCATAAGCa                <  2:1010762/124‑1 (MQ=255)                                                                                                                                     |                                                                                                                                   AGCCCACACAGGGGAGAGTGGGCTGAAATGGGAAGCTAAAGACTCAAGTAAACTTATCGGAAATAAGGACCACGCATTACGGGGGCTATCATCGCCGATGGCCGTGATTCGGCAAATCAGGTTGATTACGTAGAGAGTAAATTATCTGCTCACCGCTGCGTCACCCCTTCGTAAATCAGGCGTAGCGCGAATACACCAATAATTGCACCAATAACCCGACTGGCAACGCGTTGCATACGCCCATAAGCACGACGCACAGCGGG  >  NC_000913/345517‑345779

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