Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,061,157	G→A	*237* (TAG→TAA)	yihL →	putative DNA‑binding transcriptional regulator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,061,157	0	G	A	100.0%	190.5 / NA	60	*237* (TAG→TAA)	yihL	putative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (31/29);  total (31/29)

GGCGGTTCTGCGGGTATGCCAGCGATCTTAATTGAAAAGCATCGCGCCGACCAGCAAGGCAATATTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTATGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGAT  >  NC_000913/4061032‑4061281                                                                                                                              |                                                                                                                             ggCGGTTCTGCGGGTATGCCAGCGATCTTAATTGAAAAGCATCGCGCCGACCAGCAAGGCAATATTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATgg                                                                                                                   >  1:477990/1‑137 (MQ=255) ggCGGTTCTGCGGGTATGCCAGCGATCTTAATTGAAAAGCATCGCGCCGACCAGCAAGGCAATATTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATgg                                                                                                                   <  2:477990/137‑1 (MQ=255)     gTTCTGCGGGTATGCCAGCGATCTTAATTGAAAAGCATCGCGCCGACCAGCAAGGCAATATTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATgg                                                                                                                   >  2:133317/1‑133 (MQ=255)     gTTCTGCGGGTATGCCAGCGATCTTAATTGAAAAGCATCGCGCCGACCAGCAAGGCAATATTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATgg                                                                                                                   <  1:133317/133‑1 (MQ=255)       tCTGCGGGTATGCCAGCGATCTTAATTGAAAAGCATCGCGCCGACCAGCAAGGCAATATTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACaata                                                                                                           >  2:1425153/1‑139 (MQ=255)                aTGCCAGCGATCTTAATTGAAAAGCATCGCGCCGACCAGCAAGGCAATATTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCa                                                                                                  >  2:356433/1‑139 (MQ=255)                  gCCAGCGATCTTAATTGAAAAGCATCGCGCCGACCAGCAAGGCAATATTGTCCAGATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCaag                                                                                                >  2:633218/1‑139 (MQ=255)                    cAGCGATCTTAATTGAAAAGCATCGCGCCGACCAGCAAGGCAATATTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCaagaa                                                                                              <  1:356433/139‑1 (MQ=255)                                     aaGCATCGCGCCGACCAGCAAGGCAATATTGTCCAAATAGCTATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACgtgt                                                                              <  1:1001370/138‑1 (MQ=255)                                     aaGCATCGCGCCGACCAGCAAGGCAATATTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACgtgt                                                                              >  2:1001370/1‑138 (MQ=255)                                     aaGCATCGCGCCGACCAGCAAGGCAATATTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTc                                                                             >  1:351352/1‑139 (MQ=255)                                           cgcgCCGACCAGCAAGGCAATATTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCaagaa                                                                                              <  1:1589996/116‑1 (MQ=255)                                           cgcgCCGACCAGCAAGGCAATATTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCaagaa                                                                                              >  2:1589996/1‑116 (MQ=255)                                           cgcgCCGACCAGCAAGGCAATATTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACAct                                                                       >  2:667322/1‑139 (MQ=255)                                               ccGACCAGCAAGGCAATATTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTc                                                                   <  1:1425153/139‑1 (MQ=255)                                                          ggCAATATTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCtttattta                                                        <  2:947207/139‑1 (MQ=255)                                                          ggCAATATTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCtttattta                                                        <  1:667322/139‑1 (MQ=255)                                                                    tCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTAACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTATg                                              >  1:589461/1‑139 (MQ=255)                                                                     ccAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGaca                                                                         >  2:469715/1‑111 (MQ=255)                                                                     ccAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGaca                                                                         <  1:469715/111‑1 (MQ=255)                                                                      cAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACgtgt                                                                              <  2:805954/105‑1 (MQ=255)                                                                      cAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACgtgt                                                                              >  1:805954/1‑105 (MQ=255)                                                                       aaaTAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTATGCCt                                           >  2:262372/1‑139 (MQ=255)                                                                       aaaTAGATATTGAATATTGGCGATTTGAGGCCGTAGACATCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTATGCCt                                           <  2:589461/139‑1 (MQ=255)                                                                              tatTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATgg                                                                                                                   >  1:1313022/1‑60 (MQ=255)                                                                              tatTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATgg                                                                                                                   <  2:1313022/60‑1 (MQ=255)                                                                              tatTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCaagaa                                                                                              <  1:1112951/81‑1 (MQ=255)                                                                              tatTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCaagaa                                                                                              >  2:1112951/1‑81 (MQ=255)                                                                                 tGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGa                                                                                                                >  2:962439/1‑60 (MQ=255)                                                                                 tGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGa                                                                                                                <  1:962439/60‑1 (MQ=255)                                                                                  gAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCaagaa                                                                                              <  1:276285/77‑1 (MQ=255)                                                                                  gAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCaagaa                                                                                              >  2:276285/1‑77 (MQ=255)                                                                                   aaTATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGc                                                                                                   >  2:142048/1‑71 (MQ=255)                                                                                   aaTATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGc                                                                                                   <  1:142048/71‑1 (MQ=255)                                                                                   aaTATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTATGCCTTTCATTTAAata                               <  1:633218/139‑1 (MQ=255)                                                                                       ttGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATgg                                                                                                                   >  2:1469900/1‑51 (MQ=255)                                                                                       ttGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATgg                                                                                                                   <  1:1469900/51‑1 (MQ=255)                                                                                       ttGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTATGCCTTTCATTTAAATATGCg                           >  2:720910/1‑139 (MQ=255)                                                                                       ttGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTATGCCTTTCATTTAAATATGCg                           <  2:351352/139‑1 (MQ=255)                                                                                         ggCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACgt                                                                                >  2:876349/1‑84 (MQ=255)                                                                                         ggCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACgt                                                                                <  1:876349/84‑1 (MQ=255)                                                                                         ggCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTATGCCTTTCATTTAAATATGCGGc                         >  2:703648/1‑139 (MQ=255)                                                                                          gCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTATGCCTTTCATTTAAATATGCGGCt                        >  2:1579608/1‑139 (MQ=255)                                                                                            gATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTATGCCTTTCATTTAAATATGCGGCTgg                      <  1:720910/139‑1 (MQ=255)                                                                                             aTTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTATGCCTTTCATTTAAATATGCGGCTGGa                     >  2:719983/1‑139 (MQ=255)                                                                                                        tAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTATGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGc          >  2:1495733/1‑139 (MQ=255)                                                                                                          gACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTATGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCct        <  1:719983/139‑1 (MQ=255)                                                                                                             cTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGa                                                                           >  2:400122/1‑69 (MQ=255)                                                                                                             cTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGa                                                                           <  1:400122/69‑1 (MQ=255)                                                                                                               catcctTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGa                                                                           <  1:1092519/67‑1 (MQ=255)                                                                                                               catcatTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCaagaaa                                                                                             >  2:413303/1‑49 (MQ=255)                                                                                                               catcatTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCaagaaa                                                                                             <  1:413303/49‑1 (MQ=255)                                                                                                               catcatTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGa                                                                           >  2:1092519/1‑67 (MQ=255)                                                                                                               catcatTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTATGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGcc         >  1:203397/1‑133 (MQ=255)                                                                                                               catcatTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTATGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGcc         <  2:203397/133‑1 (MQ=255)                                                                                                                atcatTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGac                                                                          <  2:628191/67‑1 (MQ=255)                                                                                                                atcatTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGac                                                                          >  1:628191/1‑67 (MQ=255)                                                                                                                atcatTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTATGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGAt  <  1:1579608/139‑1 (MQ=255)                                                                                                                          tGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTATGCCTTTCATTTAAATATGCGGCTGGaa                    >  1:1420628/1‑111 (MQ=255)                                                                                                                          tGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTATGCCTTTCATTTAAATATGCGGCTGGaa                    <  2:1420628/111‑1 (MQ=255)                                                                                                                              |                                                                                                                             GGCGGTTCTGCGGGTATGCCAGCGATCTTAATTGAAAAGCATCGCGCCGACCAGCAAGGCAATATTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTATGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGAT  >  NC_000913/4061032‑4061281

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