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%	209.2 / NA	65	*237* (TAG→TAA)	yihL	putative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (32/33);  total (32/33)

TATCTTGGCGGTTCTGCGGGTATGCCAGCGATCTTAATTGAAAAGCATCGCGCCGACCAGCAAGGCAATATTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTATGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGC  >  NC_000913/4061026‑4061277                                                                                                                                    |                                                                                                                         tATCTTGGCGGTTCTGCGGGTATGCCAGCGATCTTAATTGAAAAGCATCGCGCCGACCAGCAAGGCAATATTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGtttt                                                                                                                   >  2:1379025/1‑139 (MQ=255) tATCTTGGCGGTTCTGCGGGTATGCCAGCGATCTTAATTGAAAAGCATCGCGCCGACCAGCAAGGCAATATTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGtttt                                                                                                                   >  2:550109/1‑139 (MQ=255)  aTCTTGGCGGTTCTGCGGGTATGCCAGCGATCTTAATTGAAAAGCATCGCGCCGACCAGCAAGGCAATATTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTACACCTCATCATTAATCTGTAAGTGTTTTa                                                                                                                  >  1:1517452/1‑139 (MQ=255)   tCTTGGCGGTTCTGCGGGTATGCCAGCGATCTTAATTGAAAAGCATCGCGCCGACCAGCAAGGCAATATTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTAt                                                                                                                 >  2:1210094/1‑139 (MQ=255)     ttGGCGGTTCTGCGGGTATGCCAGCGATCTTAATTGAAAAGCATCGCGCCGACCAGCAAGGCAATATTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATgg                                                                                                               <  2:337503/139‑1 (MQ=255)       ggCGGTTCTGCGGGTATGCCAGCGATCTTAATTGAAAAGCATCGCGCCGACCAGGAAGGCAATATTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTGATCATTAAACTGTAAGTGTTTTATGGTg                                                                                                             >  2:1106491/1‑139 (MQ=255)           gTTCTGCGGGTATGCCAGCGATCTTAATTGAAAAGCATCGCGCCGACCAGCAAGGCAATATTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATgg                                                                                                               <  1:1309475/133‑1 (MQ=255)           gTTCTGCGGGTATGCCAGCGATCTTAATTGAAAAGCATCGCGCCGACCAGCAAGGCAATATTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATgg                                                                                                               >  2:1309475/1‑133 (MQ=255)                  gggTATGCCAGCGATCTTAATTGAAAAGCATCGCGCCGACCAGCAAGGCAATATTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAataaataa                                                                                                  >  2:572024/1‑139 (MQ=255)                     tATGCCAGCGATCTTAATTGAAAAGCATCGCGCCGACCAGCAAGGCAATATTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGc                                                                                               <  1:1379025/139‑1 (MQ=255)                        gCCAGTGATCTTAATTGAAAAGCATCGCGCCGACCAGCAAGGCAATATTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGTCCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCaag                                                                                            <  1:558757/139‑1 (MQ=255)                        gCCAGCGATCTTAATTGAAAAGCATCGCGCCGACCAGCAAGGCAATATTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCaag                                                                                            <  1:1498177/139‑1 (MQ=255)                          cAGCGATCTTAATTGAAAAGCATCGCGCCGACCAGCAAGGCAATATTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCaagaa                                                                                          <  1:572024/139‑1 (MQ=255)                           aGCGATCTTAATTGAAAAGCATCGCGCCGACCAGCAAGGCAATATTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCaagaaa                                                                                         <  2:1517452/139‑1 (MQ=255)                                          aaaGCATCGCGCCGACCAGCAAGGCAATATTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTAc                                                                                 >  1:11534/1‑132 (MQ=255)                                          aaaGCATCGCGCCGACCAGCAAGGCAATATTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTAc                                                                                 <  2:11534/132‑1 (MQ=255)                                             gCATCGCGCCGACCAGCAGGGCAATATTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGa                                                                       <  2:201025/139‑1 (MQ=255)                                               aTCGCGCCGACCAGCAAGGCAATATTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGaca                                                                     <  1:1210094/139‑1 (MQ=255)                                                 cgcgCCGACCAGCAAGGCAATATTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACa                                                                                                          <  1:784971/100‑1 (MQ=255)                                                 cgcgCCGACCAGCAAGGCAATATTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACa                                                                                                          >  2:784971/1‑100 (MQ=255)                                                       gACCAGCAAGGCAATATTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCct                                                             <  1:316440/139‑1 (MQ=255)                                                             cAAGGCAATATTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTAc                                                                                 >  2:1194770/1‑113 (MQ=255)                                                             cAAGGCAATATTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTAc                                                                                 <  1:1194770/113‑1 (MQ=255)                                                                          tCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCGTCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTAc                                                   >  2:1086495/1‑130 (MQ=255)                                                                          tCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATgg                                                                                                               <  1:1067051/70‑1 (MQ=255)                                                                          tCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATgg                                                                                                               >  2:1067051/1‑70 (MQ=255)                                                                          tCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTAc                                                   <  1:1086495/130‑1 (MQ=255)                                                                          tCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTATg                                          >  1:581212/1‑139 (MQ=255)                                                                               aTAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTc                                                               >  1:852231/1‑113 (MQ=255)                                                                               aTAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTc                                                               <  2:852231/113‑1 (MQ=255)                                                                               aTAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTATGCCttt                                     >  2:894285/1‑139 (MQ=255)                                                                                 aGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACACaagaa                                                                                          <  2:1236917/84‑1 (MQ=255)                                                                                 aGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACACaagaa                                                                                          >  1:1236917/1‑84 (MQ=255)                                                                                   atatTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGc                                               >  2:308846/1‑125 (MQ=255)                                                                                   atatTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGc                                               <  1:308846/125‑1 (MQ=255)                                                                                    tatTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGaca                                                                     <  1:1088188/102‑1 (MQ=255)                                                                                    tatTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCAGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGaca                                                                     >  2:1088188/1‑102 (MQ=255)                                                                                        gAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTATGCCTTTCATTTAAat                            >  2:1508396/1‑139 (MQ=255)                                                                                          atatTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTATGCCTTTCATTTAAatat                          >  2:70439/1‑139 (MQ=255)                                                                                           tatTGGCGATTTGAGGCTGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCCACACTCTTCCTCTTTATTTACATGCTTATGCCTTTCATTTAAATATg                         <  1:1106491/139‑1 (MQ=255)                                                                                           tatTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTATGCCTTTCATTTAAATATg                         <  1:894285/139‑1 (MQ=255)                                                                                              tGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACg                                                                             <  2:127426/84‑1 (MQ=255)                                                                                              tGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACg                                                                             >  1:127426/1‑84 (MQ=255)                                                                                              tGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTATGCCTTTCATTTAAATATg                         >  1:288893/1‑136 (MQ=255)                                                                                              tGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTATGCCTTTCATTTAAATATg                         <  2:288893/136‑1 (MQ=255)                                                                                                   aTTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTATGCCTTTCATTTAAATATGCGGCTGGa                 >  1:907704/1‑139 (MQ=255)                                                                                                             gTAGACCTCATCATTAATCTGTAAGTGTTTTATgg                                                                                                               >  1:37668/1‑35 (MQ=255)                                                                                                             gTAGACCTCATCATTAATCTGTAAGTGTTTTATgg                                                                                                               <  2:37668/35‑1 (MQ=255)                                                                                                             gTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCctct                                                           <  2:777723/87‑1 (MQ=255)                                                                                                             gTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCctct                                                           >  1:777723/1‑87 (MQ=255)                                                                                                              tAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCaagaa                                                                                          <  1:180624/55‑1 (MQ=255)                                                                                                              tAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCaagaa                                                                                          >  2:180624/1‑55 (MQ=255)                                                                                                              tAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTc                                                               >  1:931222/1‑82 (MQ=255)                                                                                                              tAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTc                                                               <  2:931222/82‑1 (MQ=255)                                                                                                               aGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAATAAAGATTCTACAACGTGTCGacac                                                                    <  2:346252/76‑1 (MQ=255)                                                                                                               aGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGacac                                                                    >  1:346252/1‑76 (MQ=255)                                                                                                                  ccTCATCATCAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTATGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCctgc  <  2:905317/139‑1 (MQ=255)                                                                                                                    tcatcaTTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTAc                                                   <  2:435594/88‑1 (MQ=255)                                                                                                                    tcatcaTTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTAc                                                   >  1:435594/1‑88 (MQ=255)                                                                                                                     catcatTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACa                                                  >  2:3655/1‑88 (MQ=255)                                                                                                                     catcatTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACa                                                  <  1:3655/88‑1 (MQ=255)                                                                                                                     catcatTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTATGCCTTTCATTTAAATATg                         >  1:851442/1‑113 (MQ=255)                                                                                                                     catcatTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCGCTTTATTTACATGCTTATGCCTTTCATTTAAATATg                         <  2:851442/113‑1 (MQ=255)                                                                                                                      atcatTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCaagaa                                                                                          <  2:489932/47‑1 (MQ=255)                                                                                                                      atcatTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCaagaa                                                                                          >  1:489932/1‑47 (MQ=255)                                                                                                                                    |                                                                                                                         TATCTTGGCGGTTCTGCGGGTATGCCAGCGATCTTAATTGAAAAGCATCGCGCCGACCAGCAAGGCAATATTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTATGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGC  >  NC_000913/4061026‑4061277

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