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%	128.8 / NA	35	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 (18/17);  total (18/17)

CTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTGCAGTTTATGGTCTATTGCATAGACTGAGGGGGCAGCACGAAACGTTAAAGTGCGATAGCCGCCACAGTGCGACGTAATGATGATGTTATTTTCAGCAAGGGATCGCTATAACTCATCAGGCCAACACCCACGATCAACAGATGAATTTATTATGATAAGAAATGTGTTGTTAGCGTTCATGATATGCAGCGGAATGACATTACTCGGAGGATGCTCCAGCGTGATGTCCCATACC  >  NC_000913/3867353‑3867811                                                                                                                                                                                                                                |                                                                                                                                                                                                                                            cTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTATTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGt                                                                                                                                                                                                                      >  1:161197‑M1/1‑246 (MQ=255)   cagTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatt                                                                                                                                                                                                                 <  1:83153‑M1/249‑7 (MQ=255)   aaGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagattt                                                                                                                                                                                                                >  2:382741‑M1/1‑244 (MQ=255)       ggTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgg                                                                                                                                                                                                            >  2:28902‑M1/1‑240 (MQ=255)         ttAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCCGTTTATGGTCTATTGCATAGtagatttgcggg                                                                                                                                                                                                           >  2:266052‑M1/1‑238 (MQ=255)           aaacaaacGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGAGAAATCAAAGTTACGCATAATCAATATCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggat                                                                                                                                                                                                         <  2:108010‑M1/250‑15 (MQ=255)               aaacGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataa                                                                                                                                                                                                    <  1:210261‑M1/251‑20 (MQ=255)                       aaaTCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggat                                                                                                                                                                                            <  2:481598‑M1/251‑28 (MQ=255)                                    gAACGGTAAAGCGGGGATAAATCAAAGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaa                                                                                                                                                                                >  1:482716‑M1/1‑211 (MQ=255)                                         gTAAAGCGGGGATAAATCAAAGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggc                                                                                                                                                                          >  1:17272‑M1/1‑206 (MQ=255)                                                   gATAAATCAAAGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatc                                                                                                                                                                >  1:20212‑M1/1‑196 (MQ=255)                                                    aTAAATCAAAGTTACGCATAATCAATAGCGCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatcg                                                                                                                                                               <  2:36805‑M1/251‑57 (MQ=255)                                                                tACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCGCAATGGACAGGTCAGGTAGCCAAAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTGTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATTGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaggggcggcgatgatctccggccttttc                                                                                                                                                    <  2:412604‑M1/250‑68 (MQ=255)                                                                        ctcagtagctcCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGACAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACATTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggctgataattgcggattgcagaactaaaaatggcggagatgatgtccgaccttttcttatagact                                                                                                                                           >  2:385037‑M1/6‑175 (MQ=255)                                                                                               agaATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgttaattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtggacgc                                                                                                                    >  2:523242‑M1/1‑152 (MQ=255)                                                                                                    gTAAGCCCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGAGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGCACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcggggtgataattgcggattgcagcaataaaaagggcggagatgatctcc                                                                                                                                                             >  2:389009‑M1/1‑147 (MQ=255)                                                                                                    gTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctcc                                                                                                                                                             <  1:389009‑M1/205‑59 (MQ=255)                                                                                                      aaGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaacgccttatc                                                                                                              <  1:358109‑M1/250‑106 (MQ=255)                                                                                                       aGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaacgccttatccc                                                                                                            >  2:293636‑M1/1‑144 (MQ=255)                                                                                                             tCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCCTCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgctcttttcttatagcttcttgccggatgcggcgggaacgccttatccggcct                                                                                                        <  2:161197‑M1/249‑112 (MQ=255)                                                                                                              ccACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaacgccttatccggcctaca                                                                                                     <  1:382741‑M1/251‑115 (MQ=255)                                                                                                                                                                 ggCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaacgccttatccggcctacaaaatcatgaaaattcaatacattgcaagattttcgtaggcctgataagcgt                                                  >  2:482601‑M1/1‑86 (MQ=255)                                                                                                                                                                       gATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgttaattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaacgccttatccggcctacaaaatcatgaaaattcaatacattgcaagattttcgtaggcctgataagcgtagcgca                                            <  1:523242‑M1/251‑172 (MQ=255)                                                                                                                                                                          gTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCTGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaacgccttatccggcctacaaaatcatgaaaattcaatacattgcaagattttcgtaggcctgataagcgtagcgcatc                                          <  2:203087‑M1/250‑174 (MQ=255)                                                                                                                                                                            aaaaaTGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaacgccttatccggcctacaaaatcatgaaaattcaatacattgcaagattttcgtaggcctgataagcgtagcgcatcag                                        <  1:293636‑M1/250‑176 (MQ=255)                                                                                                                                                                            aaaaaTGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaacgccttatccggcctacaaaatcatgaaaattcaatacattgcaagattttcgtaggcctgataagcgtagcgcatcag                                        <  1:239843‑M1/250‑176 (MQ=255)                                                                                                                                                                              aaaTGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaacgccttatccggcctacaaaatcatgaaaattcaatacattgcaagattttcgtaggcctgataagcgtagcgcatcaggca                                     <  2:245797‑M1/251‑179 (MQ=255)                                                                                                                                                                                  gggTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaacgccttatccggcctacaaaatcatgaaaattcaatacattgcaagattttcgtaggcctgataagcgtagcgcatcaggcacgct                                 >  1:187941‑M1/1‑69 (MQ=255)                                                                                                                                                                                             tGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaacgccttatccggcctacaaaatcatgaaaattcaatacattgcaagattttcgtaggcctgataagcgtagcgcatcaggcacgctcgcatggtt                        >  1:91725‑M1/1‑58 (MQ=255)                                                                                                                                                                                               aaCTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaacgccttatccggcctacaaaatcatgaaaattcaatacattgcaagattttcgtaggcctgataagcgtagcgcatcaggcacgctcgcatggttagcg                    >  1:234297‑M1/1‑56 (MQ=255)                                                                                                                                                                                                          ccTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaacgccttatccggcctacaaaatcatgaaaattcaatacattgcaagattttcgtaggcctgataagcgtagcgcatcaggcacgcccgcatggttagcgcccttaaatat         >  2:69646‑M1/1‑45 (MQ=255)                                                                                                                                                                                                             tATCAGCGGTCTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaacgccttatccggcctacaaaatcatgaaaattcaatacattgcaagattttcgtaggcctgataagcgtagcgcatcaggcacgctcgcatggttagcgccattaaatatcg       <  2:175229‑M1/250‑209 (MQ=255)                                                                                                                                                                                                                cAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaacgccttatccggcctacaaaatcatgaaaattcaatacattgcaagattttcgtaggcctgataagcgtagcgcatcaggcacgctcgcatggttagcgccattaaatatcgatat   >  1:140919‑M1/1‑39 (MQ=255)                                                                                                                                                                                                                 aGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaacgccttatccggcctacaaaatcatgaaaattcaatacattgcaagattttcgtaggcctgataagcgtagcgcatcaggcacgctcgcatggttagcgccattaaatatcgatatt  <  2:17272‑M1/251‑214 (MQ=255)                                                                                                                                                                                                                 aGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaacgccttatccggcctacaaaatcatgaaaattcaatacattgcaagattttcgtaggcctgataagcgtagcgcatcaggcacgctcgcatggttagcgccattaaatatcgatatt  >  1:357081‑M1/1‑38 (MQ=255)                                                                                                                                                                                                                                |                                                                                                                                                                                                                                            CTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTGCAGTTTATGGTCTATTGCATAGACTGAGGGGGCAGCACGAAACGTTAAAGTGCGATAGCCGCCACAGTGCGACGTAATGATGATGTTATTTTCAGCAAGGGATCGCTATAACTCATCAGGCCAACACCCACGATCAACAGATGAATTTATTATGATAAGAAATGTGTTGTTAGCGTTCATGATATGCAGCGGAATGACATTACTCGGAGGATGCTCCAGCGTGATGTCCCATACC  >  NC_000913/3867353‑3867811

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