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%	194.4 / NA	53	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 (23/30);  total (23/30)

CATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGG‑GTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTGCAGTTTATGGTCTATTGCATAGACTGAGGGGGCAGCACGAAACGTTAAAGTGCGATAGCCGCCACAGTGCGACGTAATGATGATGTTATTTTCAGCAAGGGATCGCTATAACTCATCAGGCCAACACCCACGATCAACAGATGAATTTATTATGATAAGAAATGTGTTGTTAGCGTTCATGATATGCAGCGGAATGACATTACTCGGAGGATGCTCCAGCGTGATGTCCCATACC  >  NC_000913/3867329‑3867811                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                            cATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGG‑GTCTGGAAATGAACTTTTAAAGCCTTATCAGCGGACTTTTTTTTTCAg                                                                                                                                                                                                                                          <  1:322358/251‑1 (MQ=255)    tACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGGCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAAGCAAAGTGACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGGCCTCATCGGCTACGATGTAAAAATGGTGT‑TGGAAATGAACTTTTCAAGCCTTATCAGCGGACGTTTTTTTTCAGtt                                                                                                                                                                                                                                        <  1:410854/250‑1 (MQ=255)     aCTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGG‑GTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGttt                                                                                                                                                                                                                                       <  1:360774/250‑1 (MQ=255)           gctctggTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGG‑GTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCt                                                                                                                                                                                                                                <  1:270586/251‑1 (MQ=255)                ggTTGTTTTCTAAGTGGTTAACCAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGG‑GTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTAtt                                                                                                                                                                                                                             <  1:244217/249‑1 (MQ=255)                ggTTGTTTTCTAAGTGGTTAAACAAACGGTCAAAGCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGTTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGG‑GTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGGTTATGGTCTATTg                                                                                                                                                                                                                            <  2:35418/250‑1 (MQ=255)                ggTTGTGTTCTAAGGGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAACCGGGGATAAATCAAAGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGG‑GTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGc                                                                                                                                                                                                                           <  1:36769/251‑1 (MQ=255)                           aaGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTGCGCATAATCAATAGCCCCTGGAAGCAGGGAGAATGGAAGACCTTCCACAAGGGACAAGTCAGGGAGGCAGAACCCCCATTAGGCTCGGCCGCATTGGCTACGATGTAAAAATGG‑GTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagattt                                                                                                                                                                                                                <  1:32997‑M1/251‑8 (MQ=255)                                     acaaacGGTCAAAACCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGG‑GTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgat                                                                                                                                                                                                      <  1:314567‑M1/251‑18 (MQ=255)                                      caaacGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCCAAGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGG‑GTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgat                                                                                                                                                                                                      >  2:476400‑M1/1‑233 (MQ=255)                                         acGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGG‑GTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataatt                                                                                                                                                                                                  <  1:61440‑M1/251‑22 (MQ=255)                                          cGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGG‑GTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattg                                                                                                                                                                                                 >  1:357694‑M1/1‑229 (MQ=255)                                           ggTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGG‑GTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgc                                                                                                                                                                                                >  1:151980‑M1/1‑228 (MQ=255)                                                        aGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGG‑GTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaat                                                                                                                                                                                   <  1:470746‑M1/251‑37 (MQ=255)                                                          caggacgggAAAGCGGGGATAAATAAAAGTTNCGCATAATCAATATCGCCTGAAATCAGGGAGAATGTAAGGCCCTCCACAATGGACAGGTCAGGTAGCCAGAACACCCAGCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGG‑GTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTGTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataa                                                                                                                                                                                 <  2:228905‑M1/247‑39 (MQ=255)                                                               cGGTAAAGCGGGGATAAATCAAAGTTCCGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGG‑GTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagg                                                                                                                                                                            <  2:92646‑M1/251‑44 (MQ=255)                                                               cGGTAAAGCGGGGATAAATCAAAGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGG‑GTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagg                                                                                                                                                                            <  1:104909‑M1/251‑44 (MQ=255)                                                                 gTAAAGCTGGGATAAATCAAAGTTACGCAAAATCAATAGCTCCTGAAAACAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGG‑GTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggc                                                                                                                                                                          <  2:92074‑M1/251‑46 (MQ=255)                                                                        ggggATAAATCAAAGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGG‑GTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatg                                                                                                                                                                   <  1:4672‑M1/251‑53 (MQ=255)                                                                        ggggATAAATCAAAGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGG‑GTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaaaaaaaagggcggagatg                                                                                                                                                                   <  1:97367‑M1/251‑53 (MQ=255)                                                                         gggATAAATCAAAGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGG‑GTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggaaatg                                                                                                                                                                   >  2:381420‑M1/1‑198 (MQ=255)                                                                            aTAAATCAAAGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGG‑GTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatct                                                                                                                                                               >  1:456203‑M1/1‑195 (MQ=255)                                                                              aaaTCAAAGTTACGTATAAGCAAGAGCTCCTGAAATCAGCGAGAATGTAAGACCGTCCCCAATGGCCAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGG‑GTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctcc                                                                                                                                                             <  2:436847‑M1/251‑59 (MQ=255)                                                                                  caaagtatcgcataatcAATAGCGCCTGAAANCAGCGAGAATGTAAGACCGTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTTATCGGCTACGATGTAAAAATGG‑GTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataatttcggattgcagaaataaaaagggcggagatgatctccgccc                                                                                                                                                         <  2:242606‑M1/243‑63 (MQ=255)                                                                                          cGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGG‑GTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttctt                                                                                                                                                  >  1:11393‑M1/1‑181 (MQ=255)                                                                                               aaTCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGG‑GTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagctgatctccgcccttttcttatagct                                                                                                                                            >  1:203525‑M1/1‑176 (MQ=255)                                                                                                      aGCTCCTGCAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGAAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGG‑GTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCCTAGtagatttgcgggatgataattgaggattgcagaaataaaaagggcggcgatgatccccgcccttttcttatcgcttcttgcc                                                                                                                                     >  2:418982‑M1/1‑169 (MQ=255)                                                                                                             gAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGG‑GTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcg                                                                                                                              >  1:459268‑M1/1‑162 (MQ=255)                                                                                                                  cAGCGAGAATGTAAGACCTTCCACAATGGACGGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGG‑GTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgt                                                                                                                          <  1:219342‑M1/250‑94 (MQ=255)                                                                                                                  cAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGG‑GTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtg                                                                                                                         >  1:43728‑M1/1‑157 (MQ=255)                                                                                                                              aaGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGG‑GTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaacgccttatcc                                                                                                             <  1:250257‑M1/251‑107 (MQ=255)                                                                                                                                     tccaaaagtgacaggtcagGTAACCAGACCTCCCAACAGGCGCTTCCTCAGCGGCTACGATGGAAAAATGG‑GTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatcgccgccctttgcttatagcttcttgccggaggcggcgttaacgccttatccggcct                                                                                                        <  2:193234‑M1/240‑112 (MQ=255)                                                                                                                                     tCCACAATGGACAGGTCAGGTAGCCAGAANCCCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGG‑GTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaacgccttatccggcct                                                                                                        <  2:219956‑M1/249‑112 (MQ=255)                                                                                                                                          aaTGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGCTGTAAAAATGG‑GTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATAGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttctgatagcttcttgccagtagctgcgtgcacgccttctccggcccaaaaaat                                                                                                 >  2:144358‑M1/1‑133 (MQ=255)                                                                                                                                             ggACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGG‑GTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaacgccttatccggcctacaaaatcat                                                                                              <  1:476400‑M1/251‑122 (MQ=255)                                                                                                                                             ggACAGGTAAGGCAGCCAGAACACCAATCTGGCGCGTCCTCATCGGCTACGATGTAAAAATGG‑GTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaacgccttatccggcctacaaaatcat                                                                                              >  2:451869‑M1/1‑130 (MQ=255)                                                                                                                                              gACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCGTCGGCTACGCTGTAAAAATGG‑GTCTGGAAAGGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGgagatttgcgggatgcgaattgcggattgcagaaataaaaagggcggggctgatctacgcccttttcttatagatcattgcgggttgcggcgagaacgccttacccggcctacaacataatg                                                                                             >  2:138398‑M1/1‑129 (MQ=255)                                                                                                                                                caggtcaggtAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGG‑GTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaacgccttatccggcctacaaaatcatgaa                                                                                           <  1:473727‑M1/251‑125 (MQ=255)                                                                                                                                                     caggtAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGG‑GTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaacgccttatccggcccacaaaatcatgaaaattc                                                                                      >  2:396181‑M1/1‑122 (MQ=255)                                                                                                                                                       ggtAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGG‑GTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaaaggcgtggacgatctccgcccttttcttatagcttcttgccggatgcgtcttggaccccttttccgcccttcaaaatcattacaattccc                                                                                    >  2:14991‑M1/1‑120 (MQ=255)                                                                                                                                                             aataacacCCATCAGGCGCGTCCTCATCGGGTAAGTTGTAAAAATGG‑GTCTTGAAATGAACGTTTCAAGCCTTTTCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGgagagttgcgtgatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaacgccttatccggcctacaaaatcatgaatattcaatacatt                                                                              <  2:75825‑M1/248‑138 (MQ=255)                                                                                                                                                                         cAGGCGCGTCCTCATCGGCTACGATGTAAAAATGG‑GTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaacgccttatccggcctacaaaatcatgaaaattcaatacattgcaagattttc                                                                   <  1:14545‑M1/250‑149 (MQ=255)                                                                                                                                                                                     cATCGGCTACGATGTAAAAATGG‑GTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttctttccggaggcggcgtgaacgccttatccggcctacaaaaacatgaaaattcaatacattgcaagattttcgtaagcctgata                                                       >  2:69896‑M1/1‑90 (MQ=255)                                                                                                                                                                                         ggCTACGATGTAAAAATGG‑GTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaacgccttatccggcctacaaaatcatgaaaattcaatacattgcaagattttcgtaggcctgataagcgt                                                  >  2:417243‑M1/1‑86 (MQ=255)                                                                                                                                                                                             aCGATGTGCAAATGG‑GTCGGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaacgccttatccggcctacaaaatcatgcaaattaaatacattgcacgattttcgtaggcctgttaagcgtag                                                >  2:411716‑M1/1‑82 (MQ=255)                                                                                                                                                                                                  gTAAAAATGG‑GTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaacgccttatccggcctacaaaatcatgaaaattcaatacattgcaagattttcgtaggcctgataagcgtagcgcatca                                         <  1:451869‑M1/251‑175 (MQ=255)                                                                                                                                                                                                       aaTGG‑GTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaacgccttatccggcctacaaaatcatgaaaattcaatacattgcaagattttcgtaggcctgataagcgtagcgcatcaggcac                                    <  1:324758‑M1/251‑180 (MQ=255)                                                                                                                                                                                                             gTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaacgccttatccggcctacaaaatcatgaaaattcaatacattgcaagattttcgtaggcctgataagcgtagcgcatcaggcacgctcg                               >  1:364054‑M1/1‑67 (MQ=255)                                                                                                                                                                                                                       gAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaacgccttatccggcctacaaaatcatgaaaattcaatacattgcaagattttcgtaggcctgataagcgtagcgcatcaggcacgctcgcatggttagc                     >  1:139797‑M1/1‑57 (MQ=255)                                                                                                                                                                                                                              gcaAGCCTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaacgccttatccggcctacaaaatcatgaaaattcaatacattgcaagattttcgtaggcctgataagcgtagcgcatcaggcacgctcgcatggttagcgccatt               <  1:282622‑M1/249‑201 (MQ=255)                                                                                                                                                                                                                                    cTTATCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaacgccttatccggcctacaaaatcatgaaaattcaataccttgcaagattttcgtaagcctgataagcgtagcgcatcaggcacgctcgcctggttagcgccattaaatatc        >  1:295130‑M1/1‑44 (MQ=255)                                                                                                                                                                                                                                       aTCAGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaaagccttatccggcctacaaaatcatgaaaattcaatacattgcaagattttcgtaggcctgataagcgtagcgcctcaggaccgctcgcctggttaacgcccttaaatatcgat     >  2:409397‑M1/1‑41 (MQ=255)                                                                                                                                                                                                                                          aGCGGACTTTTTTTTTCAGTTTATGGTCTATTGCATAGtagatttgcgggatgataattgcggattgcagaaataaaaagggcggagatgatctccgcccttttcttatagcttcttgccggatgcggcgtgaacgccttatccggcctacaaaatcatgaaaattcaatacattgcaagattttcgtaggcctgataagcgtagcgcatcaggcacgctcgcatggttagcgccattaaatatcgatatt  >  1:79675‑M1/1‑38 (MQ=255)                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                            CATTACTCTGGCTCTGGTTGTTTTCTAAGTGGTTAAACAAACGGTCAAATCCAATAGCAGAACGGTAAAGCGGGGATAAATCAAAGTTACGCATAATCAATAGCTCCTGAAATCAGCGAGAATGTAAGACCTTCCACAATGGACAGGTCAGGTAGCCAGAACACCCATCAGGCGCGTCCTCATCGGCTACGATGTAAAAATGG‑GTCTGGAAATGAACTTTTCAAGCCTTATCAGCGGACTTTTTTTTGCAGTTTATGGTCTATTGCATAGACTGAGGGGGCAGCACGAAACGTTAAAGTGCGATAGCCGCCACAGTGCGACGTAATGATGATGTTATTTTCAGCAAGGGATCGCTATAACTCATCAGGCCAACACCCACGATCAACAGATGAATTTATTATGATAAGAAATGTGTTGTTAGCGTTCATGATATGCAGCGGAATGACATTACTCGGAGGATGCTCCAGCGTGATGTCCCATACC  >  NC_000913/3867329‑3867811

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