Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,867,576	G→T	intergenic (‑154/‑152)	ibpA ← / → yidQ	small heat shock protein IbpA/DUF1375 domain‑containing putative lipoprotein YidQ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,867,576	0	G	T	100.0%	99.1 / NA	27	intergenic (‑154/‑152)	ibpA/yidQ	small heat shock protein IbpA/DUF1375 domain‑containing putative lipoprotein YidQ
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (15/12);  total (15/12)

CTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTGCAGTTTATGGTCTATTGCATAGACTGAGGGGGCAGCACGAAACGTTAAAGTGCGATAGCCGCCACAGTGCGACGTAATGATGATGTTATTTTCAGCAAGGGATCGCTATAACTCATCAGGCCAACACCCACGATCAACAGATGAATTTATTATGATAAGAAATGTGTTGTTAGCGTTCATGATATGCAGCGGAATGACATTACTC  >  NC_000913/3867334‑3867781                                                                                                                                                                                                                                                   |                                                                                                                                                                                                              ctctggctctggTTGTTTTCTAAGTGGTTAAACAAACGGTCACATCCCATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATAATCAATAGCTCCTGCAATCAGCGAGAATGTAAGACCTTCCAAAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTAt                                                                                                                                                                                                       >  1:329869/1‑251 (MQ=255)     ggctctggTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGt                                                                                                                                                                                                    <  1:423015/250‑1 (MQ=255)                      aaGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTAGTGCATAGtagat                                                                                                                                                                                    <  1:402705‑M1/249‑6 (MQ=255)                                acaaacGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTATCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGAGGTAAAAAGGTGTCTGGAAATTAACTTTTCAAGCCTTATCAGCGGGCTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgat                                                                                                                                                                        <  1:244905‑M1/251‑18 (MQ=255)                                  aaacGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgata                                                                                                                                                                       >  1:362968‑M1/1‑232 (MQ=255)                                          aaaTCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggat                                                                                                                                                              >  1:393432‑M1/1‑224 (MQ=255)                                             tCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCCGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgc                                                                                                                                                           >  2:453801‑M1/1‑221 (MQ=255)                                                                  cGGGGATAAATCAAAGTTACGCATAATCAATAGCTCCTGAAAGCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCGGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGGCTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagattggcgggatgataattgcggcttgcagaaataaaaagggcggagat                                                                                                                                      <  1:279427‑M1/251‑52 (MQ=255)                                                                    gggATAACTCAAAGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatg                                                                                                                                     >  1:431960‑M1/1‑198 (MQ=255)                                                                         aaaTCAAAGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctcc                                                                                                                               >  1:341512‑M1/1‑193 (MQ=255)                                                                                  ttACGCATAATCAATAGCACCTGAACTCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagctttgcgggatgataattgcggattgcagaaataaaaagggaggagatgatctccgcccttttc                                                                                                                      >  2:80055‑M1/1‑184 (MQ=255)                                                                                             cAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCAACAATGGCCAGGTCAGGTAGCCAGAACACACCTAAGGCGAGTACTCATCGGCTACGATGTAAAAATGGGTCTGGAAATCAACTTTTCAAGCCTTATCAGCGGCCTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttct                                                                                                           >  1:36235‑M1/1‑173 (MQ=255)                                                                                               aTAGCTCCTGAAATCAGCGAGCATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAACATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatcgcttcttg                                                                                                         >  2:260312‑M1/1‑171 (MQ=255)                                                                                                     ccggaaaTCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggat                                                                                                   <  1:80055‑M1/248‑87 (MQ=255)                                                                                                          aaTCAGCGTGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCCGTTTATgg                                                                                                                                                                                                     >  2:175171/1‑148 (MQ=255)                                                                                                          aaTCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTGTTTTTCAGTTTATgg                                                                                                                                                                                                     <  1:175171/148‑1 (MQ=255)                                                                                                              aGCGAGAATGTAAGACCTTCCACACTGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTACTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCCTAGtagatttgcgggatgataattgcggattgcagaaataaacagggcggagatgctctcagcccttttcttatagcttcttgccggatgcggcgtgg                                                                                          >  2:443987‑M1/1‑156 (MQ=255)                                                                                                                cGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGtt                                                                                                                                                                                                          <  2:477959/137‑1 (MQ=255)                                                                                                                cGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGtt                                                                                                                                                                                                          >  1:477959/1‑137 (MQ=255)                                                                                                                         aacaccgtcaacaatgaacaggtCAGGTAGACAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCAGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATTGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctcctcccttttcttatagcttcttgccggatgcggcgtgaacgccttatcc                                                                               <  2:21287‑M1/241‑107 (MQ=255)                                                                                                                                   acaATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaacgccttatccggcctacaa                                                                      >  1:479652‑M1/1‑135 (MQ=255)                                                                                                                                                caggtAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaacgccttatccggcctacaaaatcatgaaaatt                                                         <  2:362968‑M1/250‑129 (MQ=255)                                                                                                                                                 aggtAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgggaaccccttctccggcctacaaaatcatggaaattc                                                        >  2:434173‑M1/1‑121 (MQ=255)                                                                                                                                                                  aTCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaacgccttatccggcctacaaaatcatgaaaattcaatacattgcaagatttt                                      <  1:438694‑M1/251‑148 (MQ=255)                                                                                                                                                                    cAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaacgccttatccggcctacaaaatcatgaaaattcaatacattgcaagattttc                                     <  2:89542‑M1/250‑149 (MQ=255)                                                                                                                                                                                   cGGCTACGATGTAAAACTGGGTATGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCTGTTTATGGTCTATTGCATAGtagatttgcgggatgataatttcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgacagccttatccggcctacaaaatcaagacaattcaattcattgcacgattttctaaggcctgtaaagcg                     >  2:344077‑M1/1‑87 (MQ=255)                                                                                                                                                                                                      ggTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaacgccttatccggcctacaaaatcatgaaaattcaatacattgcaagattttcgtaggcctgataagcgtagcgcatcaggcacgctc  <  1:223472‑M1/251‑184 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                              CTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTGCAGTTTATGGTCTATTGCATAGACTGAGGGGGCAGCACGAAACGTTAAAGTGCGATAGCCGCCACAGTGCGACGTAATGATGATGTTATTTTCAGCAAGGGATCGCTATAACTCATCAGGCCAACACCCACGATCAACAGATGAATTTATTATGATAAGAAATGTGTTGTTAGCGTTCATGATATGCAGCGGAATGACATTACTC  >  NC_000913/3867334‑3867781

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