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%	134.0 / NA	36	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 (14/22);  total (14/22)

TTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGAC‑TTTTTTTTGCAGTTTATGGTCTATTGCATAGACTGAGGGGGCAGCACGAAACGTTAAAGTGCGATAGCCGCCACAGTGCGACGTAATGATGATGTTATTTTCAGCAAGGGATCGCTATAACTCATCAGGCCAACACCCACGATCAACAGATGAATTTATTATGATAAGAAATGTGTTGTTAGCGTTCATGATATGCAGCGGAATGACATTACTCGGAGGATGCTCCAGCGTGATGTCCCATACCGGAGGGAAAG  >  NC_000913/3867331‑3867821                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                      ttACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCAACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGAC‑TTTTTTTTTCAGtt                                                                                                                                                                                                                                                  >  1:413220/1‑251 (MQ=255)   aCTCTGGCTTTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGTGATAAATCAAAGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAATGCCTTCCCCAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGAC‑TTTTTTTTTCAGTTTa                                                                                                                                                                                                                                                <  2:335013/251‑1 (MQ=255)                              gTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGTGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGAC‑TTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcggg                                                                                                                                                                                                                     <  1:19414‑M1/251‑13 (MQ=255)                                           acaAATCCAAAAGCAGAACGGTAAAGCGGCGATAAATCAAATCTACGCATAAGCAATAGCACCTGAACTCAGCGAGAATGTAAGACCTTCCCCAATGTACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGAC‑TTTGTTTTTCAGCTTATGGTATATTGCATAGgagaggtgcgggaggattattgcgc                                                                                                                                                                                                        <  2:159297‑M1/250‑26 (MQ=255)                                                           aaCGGTAAAGCGGGGATAAATCAAAGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGAC‑TTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattg                                                                                                                                                                                                           >  2:464195‑M1/1‑210 (MQ=255)                                                           aaCGGTAAAGCGGGGATAAATCAAAGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGAC‑TTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattg                                                                                                                                                                                                           <  1:464195‑M1/232‑23 (MQ=255)                                                            aCGGTAAAGCGGGGATAAATCAAAGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGAC‑TTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaag                                                                                                                                                                                       <  1:356437‑M1/251‑43 (MQ=255)                                                              ggTAAAGCGGGGATAAATCAAAGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGAC‑TTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaaggg                                                                                                                                                                                     >  1:475695‑M1/1‑207 (MQ=255)                                                                     cGGGGATAAATCAAAGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGAC‑TTTTTTTTTCAGTTTATGGTCTATTGCATAGta                                                                                                                                                                                                                               >  1:338891‑M1/1‑200 (MQ=255)                                                                     cGGGGATAAATCAAAGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGAC‑TTTTTTTTTCAGTTTATGGTCTATTGCATAGta                                                                                                                                                                                                                               <  2:338891‑M1/202‑3 (MQ=255)                                                                       gggATAAATCAAAGTTACGCATAATCAATAGCTCCGGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCAGCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCCGGAAATGAACTTTTCAAGACTTATCAGCGGAC‑TTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaagaaaaagggcggagat                                                                                                                                                                              <  2:218666‑M1/249‑52 (MQ=255)                                                                          aTAAATCAAAGTTACGCATCATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGAC‑TTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagg                                                                                                                                                                                      >  2:425535‑M1/1‑195 (MQ=255)                                                                          aTAAATCAAAGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGAC‑TTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatct                                                                                                                                                                         >  1:27114‑M1/1‑195 (MQ=255)                                                                          aTAAATCAAAGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGAC‑TTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagg                                                                                                                                                                                      <  1:425535‑M1/238‑44 (MQ=255)                                                                           tAAATCAAAGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGCATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGAC‑TTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctc                                                                                                                                                                        >  2:340955‑M1/1‑194 (MQ=255)                                                                                cAAAGTTACGCATAAACAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGAC‑TTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggctgagatgagctccgcc                                                                                                                                                                    <  2:164623‑M1/250‑62 (MQ=255)                                                                                        cgaataaTCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGAC‑TTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttctta                                                                                                                                                           <  2:248184‑M1/248‑71 (MQ=255)                                                                                                   tAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGAC‑TTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggatt                                                                                                                                                                                                     <  2:458830‑M1/198‑29 (MQ=255)                                                                                                   tAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGAC‑TTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggatt                                                                                                                                                                                                     >  1:458830‑M1/1‑170 (MQ=255)                                                                                                              aTCAGCGAGAATGTAAGACCTTCCACAATGGACAGGACAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGAC‑TTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggatttcagaaataaaaagggcggagatgatctccgcccttttcttagagcttcttgccggatgcggcg                                                                                                                                     >  2:336754‑M1/1‑159 (MQ=255)                                                                                                                cAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGAC‑TTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtg                                                                                                                                   <  1:115546‑M1/251‑95 (MQ=255)                                                                                                                          gTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGAC‑TTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaacgccttat                                                                                                                         <  1:94659‑M1/251‑105 (MQ=255)                                                                                                                            aaGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGAC‑TTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaacgccttatcc                                                                                                                       <  1:336754‑M1/251‑107 (MQ=255)                                                                                                                                             aCAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGAC‑TTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaacgccttatccggcctacaaaatcatg                                                                                                       >  1:229665‑M1/1‑128 (MQ=255)                                                                                                                                               aggtcaggtAGCCAGAACACCCATCAGGTGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGAC‑TTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaacgccttatccggcctacaaaatcatgaaa                                                                                                    <  1:219918‑M1/251‑126 (MQ=255)                                                                                                                                                     ggtAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGAC‑TTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaacgccttatccggcctacaaaatcatgaaaattcaa                                                                                              <  1:240953‑M1/251‑132 (MQ=255)                                                                                                                                                        aGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGAC‑TTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaacgccttatccggcctacaaaatcatgaaaattcaata                                                                                            <  1:218276‑M1/250‑134 (MQ=255)                                                                                                                                                              aaCACCCAGCAGTCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaacgccttatccggcctacaaaatcatgaaaattcaatacattgc                                                                                      <  2:234236‑M1/251‑140 (MQ=255)                                                                                                                                                                       cAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGAC‑TTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaacgccttatccggcctacaaaatcatgaaaattcaatacattgcaagattttcg                                                                            >  1:499694‑M1/1‑102 (MQ=255)                                                                                                                                                                       cAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGAC‑TTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaacgccttatccggcctacaaaatcatgaaaattcaatacattgcaagattttcg                                                                            <  1:486568‑M1/251‑150 (MQ=255)                                                                                                                                                                                 cTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGAC‑TTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaacgccttatccggcctacaaaatcctgaacattcaataaattgcaagattttcgcaggcctgat                                                                  >  2:24510‑M1/1‑92 (MQ=255)                                                                                                                                                                                         cTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGAC‑TTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaacgccttatccggcctacaaaatcatgaaaattcaatacattgcaagattttcgtaggcctgataagcgt                                                            >  1:269061‑M1/1‑84 (MQ=255)                                                                                                                                                                                                              ggAAATGAACTTTTCAAGCCTTATCAGCGGAC‑TTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaacgccttatccggcctacaaaatcatgaaaattcaatacattgcaagattttcgtaggcctgataagcgtagcgcatcaggcacgctcgc                                        <  2:475695‑M1/248‑186 (MQ=255)                                                                                                                                                                                                                       cTTTTCAAGCCTTATCAGCGGAC‑TTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaacgccttatccggcctacaaaatcatgaaaattcaatacattgcaagattttcgtaggcctgataagcgtagcgcatcaggcacgctcgcatggttagcgcc                            >  1:369112‑M1/1‑54 (MQ=255)                                                                                                                                                                                                                                            aC‑TTTTTTTTTCAGGTGATTATCTATTGCATAGtatatttgcgggatgataattgctgattgcataaataaaaacgtcgtatatgatctccgcccttttcttatagctacttgccggatgcgccgtgaacgccttatccgtcctacaaaatcattaaaattcaatacattgcaagattttcgtaggcctgataagcgtagcgcatcaggcacgctcgcatggttatcgccattaaatatcgatattctccg       <  1:221199‑M1/251‑219 (MQ=255)                                                                                                                                                                                                                                                ttttttttCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaacgccttatccggcctacaaaatcatgaaaattcaatacattgcaagattttcgtaggcctgataagcgtagcgcatcaggcacgctcgcatggttagcgccattaaatatcgatattcgccgc      <  1:261535‑M1/249‑220 (MQ=255)                                                                                                                                                                                                                                                   ttttTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggatcgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaacgccttatccggcctacaaaatcatgaaaattcaatacattgcaagattttcgtaggcctgataagcgtagcgcatcaggcacgctcgcatggttagcgccattaaatagcgatattcgccgctttc  <  1:54176‑M1/250‑224 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                      TTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGAC‑TTTTTTTTGCAGTTTATGGTCTATTGCATAGACTGAGGGGGCAGCACGAAACGTTAAAGTGCGATAGCCGCCACAGTGCGACGTAATGATGATGTTATTTTCAGCAAGGGATCGCTATAACTCATCAGGCCAACACCCACGATCAACAGATGAATTTATTATGATAAGAAATGTGTTGTTAGCGTTCATGATATGCAGCGGAATGACATTACTCGGAGGATGCTCCAGCGTGATGTCCCATACCGGAGGGAAAG  >  NC_000913/3867331‑3867821

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