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

TATCTTGGCGGTTCTGCGGGTATGCCAGCGATCTTAATTGAAAAGCATCGCGCCGACCAGCAAGGCAATATTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTATGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATA  >  NC_000913/4061026‑4061291                                                                                                                                    |                                                                                                                                       tATCTTGGCGGTTCTGCGGGTATGCCAGCGATCTTAATTGAAAAGCATCGCGCCGACCAGCAAGGCAATATTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGtttt                                                                                                                                 >  2:1493806/1‑139 (MQ=255) tATCTTGGCGGTTCTGCGGGTATGCCAGCGATCTTAATTGAAAAGCATCGCGCCGACCAGCAAGGCAATATTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGtttt                                                                                                                                 >  2:239804/1‑139 (MQ=255)     ttGGCGGTTCTGCGGGTATGCCAGCGATCTTAATTGAAAAGCATCGCGCCGACCAGCAAGGCAATATTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATgg                                                                                                                             >  2:1132310/1‑139 (MQ=255)           gTTCTGCGGGTATGCCAGCGATCTTAATTGAAAAGCATCGCGCCGACCAGCAAGGCAATATTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACaa                                                                                                                       >  1:1259968/1‑139 (MQ=255)               tGCGGGTATGCCAGCGATCTTAATTGAAAAGCATCGCGCCGACCAGCAAGGCAATATTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACaataaa                                                                                                                   <  1:1132310/139‑1 (MQ=255)                          cAGCGATCTTAATTGAAAAGCATCGCGCCGACCAGCAAGGCAATATTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCaagaa                                                                                                        <  1:774292/139‑1 (MQ=255)                           aGCGATCTTAATTGAAAAGCATCGCGCCGACCAGCAAGGCAATATTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCaagaaa                                                                                                       <  1:1493806/139‑1 (MQ=255)                            gCGATCTTAATTGAAAAGCATCGCGCCGACCAGCAAGGCAATATTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCaagaaag                                                                                                      <  1:239804/139‑1 (MQ=255)                             cGATCTTAATTGAAAAGCATCGCGCCGACCAGCAAGGCAATATTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCaagaaaga                                                                                                     >  2:1032880/1‑139 (MQ=255)                                              cATCGCGCCGACCAGCAAGGCAATATTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGac                                                                                    >  2:161565/1‑139 (MQ=255)                                               aTCGCGCCGACCAGCAAGGCAATATTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGaca                                                                                   <  1:960637/139‑1 (MQ=255)                                                     ccGACCAGCAAGGCAATATTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTc                                                                             <  2:1259968/139‑1 (MQ=255)                                                             cAAGGCGATATTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCtttat                                                                     >  1:535789/1‑139 (MQ=255)                                                                  cAATATTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACaata                                                                                                                     <  2:440158/86‑1 (MQ=255)                                                                  cAATATTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACaata                                                                                                                     >  1:440158/1‑86 (MQ=255)                                                                        tGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTa                                                          <  1:1032880/139‑1 (MQ=255)                                                                          tCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATgg                                                                                                                             >  1:1226737/1‑70 (MQ=255)                                                                          tCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATgg                                                                                                                             <  2:1226737/70‑1 (MQ=255)                                                                                 aGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGc                                                                                                             <  1:1112881/79‑1 (MQ=255)                                                                                 aGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGc                                                                                                             >  2:1112881/1‑79 (MQ=255)                                                                                   atatTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAATTGTTTTATGGTGACa                                                                                                                        <  2:1008937/66‑1 (MQ=255)                                                                                   atatTGAATATTGGCGATTTGAGGCCGTAGACCTCATAATTAAACTGTAAGTGTTTTATGGTGACa                                                                                                                        >  1:1008937/1‑66 (MQ=255)                                                                                    tatTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTAc                                                                 <  1:1024424/120‑1 (MQ=255)                                                                                    tatTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTAc                                                                 >  2:1024424/1‑120 (MQ=255)                                                                                     atTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTATGCCTTTCATTTa                                             >  2:370635/1‑139 (MQ=255)                                                                                         aaTATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACa                                                                <  2:315409/116‑1 (MQ=255)                                                                                         aaTATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACa                                                                >  1:315409/1‑116 (MQ=255)                                                                                         aaTATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTATGCCTTTCATTTAAata                                         <  1:161565/139‑1 (MQ=255)                                                                                           tatTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCaagaa                                                                                                        <  1:795201/74‑1 (MQ=255)                                                                                           tatTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCaagaa                                                                                                        >  2:795201/1‑74 (MQ=255)                                                                                           tatTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTATGCCTTTCATTTAAATATg                                       <  1:75123/139‑1 (MQ=255)                                                                                           tatTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTATGCCTTTCATTTAAATATg                                       <  1:370635/139‑1 (MQ=255)                                                                                              tGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATgg                                                                                                                             >  1:1072368/1‑50 (MQ=255)                                                                                              tGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATgg                                                                                                                             <  2:1072368/50‑1 (MQ=255)                                                                                              tGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCaag                                                                                                          <  2:797442/69‑1 (MQ=255)                                                                                              tGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCaag                                                                                                          >  1:797442/1‑69 (MQ=255)                                                                                              tGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTATGCCTTTCATTTAAATATGCgg                                    >  2:675066/1‑139 (MQ=255)                                                                                                gCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCtt                                                                              <  1:761866/95‑1 (MQ=255)                                                                                                gCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCtt                                                                              >  2:761866/1‑95 (MQ=255)                                                                                                gCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTATGCCTTTCATTTAAATATGCGGCt                                  >  2:503181/1‑139 (MQ=255)                                                                                                   aTTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTATGCCTTTCATTTAAata                                         >  2:1389671/1‑129 (MQ=255)                                                                                                   aTTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTATGCCTTTCATTTAAata                                         <  1:1389671/129‑1 (MQ=255)                                                                                                   aTTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTATGCCTTTCATTTAAATATg                                       <  1:492611/131‑1 (MQ=255)                                                                                                   aTTTGAGGCCGTAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTATGCCTTTCATTTAAATATg                                       >  2:492611/1‑131 (MQ=255)                                                                                                              tAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCtttattta                                                                  <  1:671432/93‑1 (MQ=255)                                                                                                              tAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCtttattta                                                                  >  2:671432/1‑93 (MQ=255)                                                                                                              tAGACCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTATGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGc                    >  1:485689/1‑139 (MQ=255)                                                                                                                 aCCTCATCATTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTATGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCctg                 <  2:485689/139‑1 (MQ=255)                                                                                                                     catcatTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTATGCCTTTCATTTAAATATGCGGCTGGaa                              >  2:469374/1‑122 (MQ=255)                                                                                                                     catcatTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTATGCCTTTCATTTAAATATGCGGCTGGaa                              <  1:469374/122‑1 (MQ=255)                                                                                                                      atcatTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTAc                                                                                               <  1:148962/56‑1 (MQ=255)                                                                                                                      atcatTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTAc                                                                                               >  2:148962/1‑56 (MQ=255)                                                                                                                      atcatTAATCTGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTATGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGAt            <  1:503181/139‑1 (MQ=255)                                                                                                                                tGTAAGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTATGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGata  >  2:636628/1‑139 (MQ=255)                                                                                                                                    |                                                                                                                                       TATCTTGGCGGTTCTGCGGGTATGCCAGCGATCTTAATTGAAAAGCATCGCGCCGACCAGCAAGGCAATATTGTCCAAATAGATATTGAATATTGGCGATTTGAGGCCGTAGACCTCATCATTAATCTGTAGGTGTTTTATGGTGACAATAAATAACGCAAGAAAGATTCTACAACGTGTCGACACTCTTCCTCTTTATTTACATGCTTATGCCTTTCATTTAAATATGCGGCTGGAAAGAGTGTTGCCTGCTGATTTACTTGATA  >  NC_000913/4061026‑4061291

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: ‑