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%	130.5 / NA	38	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 (16/22);  total (16/22)

GGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTGCAGTTTATGGTCTATTGCATAGACTGAGGGGGCAGCACGAAACGTTAAAGTGCGATAGCCGCCACAGTGCGACGTAATGATGATGTTATTTTCAGCAAGGGATCGCTATAACTCATCAGGCCAACACCCACGATCAACAGATGAATTTATTATGATAAGAAATGTGTTGTTAGCGTTCATGATATGCAGCGGAATGACATTACTCGGAGGATGCTCCAGCGTGATGTCCCATACCGGAGGG  >  NC_000913/3867338‑3867817                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                  ggctctggTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGt                                                                                                                                                                                                                                        <  2:35140/250‑1 (MQ=255)         ttgttTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGAATAAATCAAAGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGCCCTTCCACAATGGAAAGGTCAGGTAGCAAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCAt                                                                                                                                                                                                                               >  1:237811/1‑251 (MQ=255)          tgttTTCTAAGTGGTTAAACAAACGGGCAAATCCAAAAGCAGAACGGTAAAGCGGGGATAAAACAAAGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCTTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTTCAt                                                                                                                                                                                                                               <  2:359701/250‑1 (MQ=255)          tggtttCTAAGTGGTTAAACAAATGGTCAAATCCAATAGAAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATAATCAATAGCTCCTGAAATCATCGAGAATGTAAGACCTGCCACAATGTACAGGTCAGGTAGCCAGATCACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTTGAAATTAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGTTCTATTGCAt                                                                                                                                                                                                                               <  2:339135/247‑1 (MQ=255)                    gTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAAGGGACAGGTCAGGTAGCCAGAACACCCCTCAGTCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgc                                                                                                                                                                                                                    <  2:144227‑M1/251‑10 (MQ=255)                      ggTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgg                                                                                                                                                                                                                  <  1:116568‑M1/251‑12 (MQ=255)                        ttAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcggg                                                                                                                                                                                                                 >  1:252137‑M1/1‑238 (MQ=255)                               aacGGGCAAAACCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATAATCAATAGCTCTTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACGCCCATCAGGCGTGTCCTCATCGGCTGCGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCGTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcggtatgataat                                                                                                                                                                                                         <  2:44097‑M1/251‑21 (MQ=255)                                 cGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATAATCAAGAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGTACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattg                                                                                                                                                                                                       <  2:219857‑M1/251‑23 (MQ=255)                                      aaaaccaataGCAGAACGGGAAAGCGGGGATAAATCAAAGTTACGCATAATCAATAGCTCCGGAAATCAGCGATAATGAAAGACCTTCCACAATGGACAGGACAGGTTGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTTGAAATGATCTTTTCAAGCCTTATCATCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggat                                                                                                                                                                                                  <  2:455442‑M1/247‑28 (MQ=255)                                        aTCCAATAGCAGAACGGTAAAGCGGGGATAANACAAAGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattg                                                                                                                                                                                                <  2:248499‑M1/251‑30 (MQ=255)                                                      cGGTAAAGCGGGGATAAATCAAAGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagg                                                                                                                                                                                  <  1:399786‑M1/251‑44 (MQ=255)                                                         tAAAGCGGGGATAAATCAAAGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGACCTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggaggataattgaggattgcagaaataaaaaggg                                                                                                                                                                                 >  2:365619‑M1/1‑205 (MQ=255)                                                                gggATAAATCAAAGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatg                                                                                                                                                                         >  1:436637‑M1/1‑198 (MQ=255)                                                                       aTCAAAGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgc                                                                                                                                                                 >  1:430256‑M1/1‑191 (MQ=255)                                                                            aGTTACGCATAAGCAAGAGCTCCTGAAATCAGCGAGAATGTAAGACCTGCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGTCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTGTATGGTCTATTGCATAGtagatttgcgtgatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttt                                                                                                                                                            <  2:317907‑M1/251‑66 (MQ=255)                                                                                          aaTAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttctt                                                                                                                                              >  1:404948‑M1/1‑172 (MQ=255)                                                                                          aaTAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttctt                                                                                                                                              >  1:416534‑M1/1‑172 (MQ=255)                                                                                             aGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGAGCGTCCTCATCGGCTACGATGTAAAAATGGGACTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaacgggcggagatgatctacgcccttttcttaaagcttcttgcc                                                                                                                                           >  1:17602‑M1/1‑169 (MQ=255)                                                                                              gCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagctccttgccg                                                                                                                                          <  2:229214‑M1/251‑84 (MQ=255)                                                                                                  cTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatg                                                                                                                                      >  1:131546‑M1/1‑164 (MQ=255)                                                                                                                aaTGTAAGACCGCCCACAATGTACAGGTCAGGTAGCCAGAACACCCATAAGGCGCGTCCTCATCGGCTAAGATGTAAAAGTGGGGCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaacgcct                                                                                                                        <  2:1249‑M1/251‑102 (MQ=255)                                                                                                                   gTAAGCCCTTCCACAATAGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaacgccttat                                                                                                                     <  1:409216‑M1/251‑105 (MQ=255)                                                                                                                     aaGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaacgccttatcc                                                                                                                   >  1:468897‑M1/1‑145 (MQ=255)                                                                                                                            tcctcaatgtACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagattcttgccggatgcggcgtgaacgccttatccggcct                                                                                                              <  2:345004‑M1/245‑112 (MQ=255)                                                                                                                            tCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagaggatctccgcccttttcttatagcttcttgcnggatgcggcgtgaacgccttatccggcctac                                                                                                            >  2:235813‑M1/1‑138 (MQ=255)                                                                                                                                          gtcaggtAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaacgccttatccggcctacaaaatcatgaaaat                                                                                              <  1:517295‑M1/251‑128 (MQ=255)                                                                                                                                            caggtAGCCAGAACACCCAGCAGGCGCGTCCTCAGCGGCTACTATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGTACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgccattttcttatagcttcttgccgtatgcggcgttaacgccttatccgtcctacaaaatcatgaaaattc                                                                                            <  1:365654‑M1/251‑130 (MQ=255)                                                                                                                                                                 aGGCGCGTCCTCATCGGCTACGAGTTAAAAATGGTTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaacgccttatccggcctacaaaatcatgaaaattcaatacattgcaagattttcgt                                                                       <  2:63399‑M1/251‑151 (MQ=255)                                                                                                                                                                                                                     cAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaaccccttatccggcctacaaaatcatgaaaattcaatacattgcaagattttcgtaggcctgataagcgtagcgcatcaggcacgctcgcatggttagcgccatta                    >  2:13517‑M1/1‑49 (MQ=255)                                                                                                                                                                                                                      aaGCCTTATCAGCGGACTTTTTTTTTCAGGTTATGGTCTATTGCATAGtagatttgcgggatgataatggcggattgcagaaataaaaagggcggagatgatcgccgcccttttcttatagcttcttgccggatgcggcgtgaacgccttatccggcctacaaaatcatgaaaatgcaatacattgcaatattttcgtaggcctgataagcgtagctcatcaggcacgctcgcatggttagcgccattaaa                  <  2:414687‑M1/251‑204 (MQ=255)                                                                                                                                                                                                                             aTCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaacgccttatccgg                                                                                                                 <  1:122392‑M1/149‑109 (MQ=255)                                                                                                                                                                                                                             aTCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaacgccttatccgg                                                                                                                 >  2:122392‑M1/1‑41 (MQ=255)                                                                                                                                                                                                                                agcgtgcttttttTTTCAGTTTATGGTCTATTGAATAGgagatttgtgggatgataattgcggattgcagaaataaaaatggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaacgccttatccggcctacaaaatcatgaaaattcaatacattgctagattttcgtaggcctgataagcgtagcgcatcaggcccgctcgcatggttagcgccattaaatatcgatatt        <  2:519924‑M1/245‑214 (MQ=255)                                                                                                                                                                                                                                aGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaacgccttatccggcctacaaaatcatgaaaattcaatacattgcaagattttcgtaggcctgataagcgtagcgcatcaggcacgctcgcatggttagcgccattaaatatcgatatt        >  1:27026‑M1/1‑38 (MQ=255)                                                                                                                                                                                                                                  cGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaacgccttatccggcctacaaaatcatgaaaattcaatacattgcaagattttcgtaggcctgataagcggagcgcatcagggacgctcgcatggttagcgccattaaatatcgatattcg      >  1:246321‑M1/1‑36 (MQ=255)                                                                                                                                                                                                                                   ggACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaacgccttatccggcctacaaaatcctgaaaattcaatacattgcaagattttcgtaggcctgataagagtagcgcatcaagcccgctcgcatggttatcgccaataaatatcgatattcgc     >  2:95995‑M1/1‑35 (MQ=255)                                                                                                                                                                                                                                    gACTTTTTTTTTCAGTTTATTGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaacgccttatccggcctacaaaatcatgaaaattcaatacattgcaagattttcgtaggcctgataagcgtagcgcatcaggcacgctcgcatggttagcgccattaaatatcgatattcgcc    <  2:105459‑M1/251‑218 (MQ=255)                                                                                                                                                                                                                                        ttttttttCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaacgccttatccggcctacaaaatcatgaaaattcaatacattgcaagattttcgtaggcctgataagcgtagcgcatcaggcacgctcgcatggttagcgccattaaatatcgatattcgccgc  <  1:169675‑M1/249‑220 (MQ=255)                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                  GGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGGGTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTGCAGTTTATGGTCTATTGCATAGACTGAGGGGGCAGCACGAAACGTTAAAGTGCGATAGCCGCCACAGTGCGACGTAATGATGATGTTATTTTCAGCAAGGGATCGCTATAACTCATCAGGCCAACACCCACGATCAACAGATGAATTTATTATGATAAGAAATGTGTTGTTAGCGTTCATGATATGCAGCGGAATGACATTACTCGGAGGATGCTCCAGCGTGATGTCCCATACCGGAGGG  >  NC_000913/3867338‑3867817

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