Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,487,163	0	C	T	42.1%	81.8 / 19.7	57	intergenic (+200/‑72)	hrpA/ydcF	putative ATP‑dependent helicase/DUF218 superfamily protein, SAM‑binding
Reads supporting (aligned to +/- strand):  ref base C (5/28);  new base T (24/0);  total (29/28)
Fisher's exact test for biased strand distribution p-value = 1.85e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01

GTTTTTTATTAATTCTAAAACGGTAAGGGTAAAAATTCAGGAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGG‑TTAAGGACATACGATTTCCTCCTTTCAGA‑GTGCTCCGCTTCTCACTATTATCTCACGCAGTA‑TTC‑TTAAGGGAAC‑GATAAGGAGGAACCATGAACATTA‑CCCCGTTTCCGACGCTTTCGCCGGCAACTATAGATGCCATAAATGTTATCGGACAGTGGCTGGCGCAGG  >  NC_000913/1487014‑1487313                                                                                                                                                       |                                                                                                                                                            gTTTTTTATTAATTCTAAAACGGTAAGGGTAAAAATTCAGGAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGG‑TTAAGGACATACGATTtcc                                                                                                                                                            <  1:775747/151‑1 (MQ=255)  ttttttATTAATTCTAAAACGGTAAGGGTAAAAATTCAGGAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGG‑ATTATGACATAAAATTtctt                                                                                                                                                           >  2:193013/1‑149 (MQ=255)       tattaATTCTAAAACGGTAAGGGTAAAAATTCAGGAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGG‑TTAAGGACATACGATTTCCTCCtt                                                                                                                                                       <  1:264695/150‑1 (MQ=255)        attaATTCTAAAACGGTAAGGGTAAAAATTCAGGAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGG‑TTAAGGACATACGATTTCCTCCTTTc                                                                                                                                                     <  2:815593/151‑1 (MQ=255)                    aCGGTAAGGGTAAAAATTCAGGAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGG‑TTTAGGAAAAAAGATTTTCTCCttt                                                                                                                                                      >  1:655045/1‑138 (MQ=255)                     cGGTAAAGGTAAAAATTCAGGAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGG‑TTAAGGACATACGATTTCCTCCTTTCAGA‑GTGCTCCGCt                                                                                                                                       <  2:421190/151‑1 (MQ=255)                         aaGGGTAAAAATTCAGGAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGG‑TTAAGGACATACGATTTCCTCCTTTCAGA‑GTGCTCCGCTtctc                                                                                                                                   <  1:931956/151‑1 (MQ=255)                         aaGGGTAAAAATTCAGGAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGG‑TTAAGGAAAAAACATTTTCTCCTTTTTGA‑GTGGTCCCCctct                                                                                                                                    >  1:678740/1‑150 (MQ=255)                          aGGGTAAAAATTCAGGAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGG‑TTAAGGACATACGATTTCCTCCTTTCAGA‑GTGCTCCGCTTCTCa                                                                                                                                  <  2:196241/151‑1 (MQ=255)                           gggTAAAAATTCAGGAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGG‑TTAAAGAAAAAAGATTTTCTCCTTTTTGA‑GTGCGcc                                                                                                                                          >  1:15237/1‑142 (MQ=255)                               aaaaaTTCAGGAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGG‑TTAAAGAAAAAAGATTTTCTCCCTTCAGA‑GGGATCCGc                                                                                                                                        >  1:1470/1‑140 (MQ=255)                                aaaaTTCAGGAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGG‑TTAAGGACATAAAATTTTCTTCTTTCTTA‑GGGGTTCGCTTCTTCCtatt                                                                                                                             >  2:571958/1‑150 (MQ=255)                                    ttCAGGAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGG‑TTAAAGACATACCATTTTCTACTTTTTTA‑GTGGTCCGCTTCTCTCTATTTtct                                                                                                                         >  1:796973/1‑150 (MQ=255)                                          aaTTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGG‑TTTAAGAAAAACGATTTTCTCCTTTCAGA‑GTGCTCCGCCTCTCCCCATTATTTCACGc                                                                                                                    >  1:751871/1‑149 (MQ=255)                                           aTTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGG‑TTAAGGACATACGATTTCCTCCTTTCAGA‑GTGCTCCGCTTCTCACTATTATCTCACGCAGt                                                                                                                 <  2:461361/151‑1 (MQ=255)                                               aGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGG‑GTAAAGACAAAACAATTTCTCCTTTCTTA‑GGGGTCCCCCTCTCACTATTTTTTCACCCAGTA‑tt                                                                                                             >  2:56197/1‑150 (MQ=255)                                                    aaTACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGG‑TTAAAGAAATAAAATTTTTTTCTTTTTTA‑GGGGTCCGCCTCTCTCTATTTTCTCCCCCAGAA‑Ttt‑ttt                                                                                                        >  2:109963/1‑148 (MQ=255)                                                     aTACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGG‑TTAAGGACATACGATTTCCTCCTTTCAGA‑GTGCTCCGCTTCTCACTATTATCTCACGCAGTA‑TTC‑TTAAgg                                                                                                     <  2:73298/151‑1 (MQ=255)                                                     aTACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGG‑TTAAAGACAAAAAATTTTTTCCCTTTTGA‑GGGCTCCGCTTCTCCCTATTTTTTCACGCCGTA‑Ttt‑tt                                                                                                         >  2:189090/1‑147 (MQ=255)                                                        cAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGG‑TTAAGGACATACGATTTCCTCCTTTCAGA‑GTGCTCCGCTTCTCACTATTATCTCACGCAGTA‑TTC‑TTAAGGGaa                                                                                                  <  1:879443/151‑1 (MQ=255)                                                          aTTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGG‑TTAAAGAAATAAAATTTTCTTCTTTTATA‑GGGGTTCGCTTCTTCCTATTATTTCACCCCGGA‑TTCTTTAAGGGa                                                                                                   >  2:138094/1‑149 (MQ=255)                                                           ttCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGG‑ATTAAGAAATAAAAATTTCTCCTTTCTTA‑GGGGTTCGCTTCTCTCTATTTTCTTACCCAGGA‑TTT‑TTTAggg                                                                                                    >  2:112347/1‑146 (MQ=255)                                                            tctcTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGG‑TTAAGGACATACGATTTCCTCCTTTCAGA‑GTGCTCCGCTTCTCACTATTATCTCACGCAGTA‑TTC‑TTAAGGGAAC‑Ga                                                                                              <  1:595158/150‑1 (MQ=255)                                                            tctcTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGG‑TTAAAGAAAAAAGATTTTCTCCCTTTTGA‑GTGCTCCGCCTCTCCCTATTTTTTCCCCCAGGA‑TTT‑TTTATGGaa                                                                                                  >  1:772402/1‑147 (MQ=255)                                                               ctgctgCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGG‑TTAAGGACATACGATTTCCTCCTTTCAGA‑GTGCTCCGCTTCTCACTATTATCTCACGCAGTA‑TTC‑TTAAGGGAAC‑GAt                                                                                             <  1:228627/148‑1 (MQ=255)                                                                 gctgcAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGG‑TTAAGGACATACGATTTCCTCCTTTCAGA‑GTGCTCCGCTTCTCACTATTATCTCACGCAGTA‑TTC‑TTAAGGGAAC‑GATaa                                                                                           <  1:361750/148‑1 (MQ=255)                                                                    gcAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGG‑TTAAGGACATACGATTTCCTCCTTTCAGA‑GTGCTCCGCTTCTCACTATTATCTCACGCAGTA‑TTC‑TTAAGGGAAC‑GATAaggagg                                                                                      <  1:939755/150‑1 (MQ=255)                                                                     cAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGG‑TTAAGGACATACGATTTCCTCCTTTCAGA‑GTGCTCCGCTTCTCACTATTATCTCACGCAGTA‑TTC‑TTAAGGGAAC‑GATAAggagga                                                                                     <  2:220790/150‑1 (MQ=255)                                                                     cAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGG‑TTAAAGAAAAACGATTTTCTCCTTTTATA‑GTGGTCCCCCTCTCCCTATTATTTCACCCAGCA‑TTT‑TTTAGGGGAA‑GATAaagagg                                                                                      >  1:220804/1‑149 (MQ=255)                                                                      aaGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGG‑TTAAGGACATACGATTTCCTCCTTTCAGA‑GTGCTCCGCTTCTCACTATTATCTCACGCAGTA‑TTC‑TTAAGGGAAC‑GATAAGGAGGAAc                                                                                   <  2:212335/151‑1 (MQ=255)                                                                       aGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGG‑TTAAGGACATACGATTTCCTCCTTTCAGA‑GTGCTCCGCTTCTCACTATTATCTCACGCAGTA‑TTC‑TTAAGGGAAC‑GATAAGGAGGAAcc                                                                                  <  2:863773/151‑1 (MQ=255)                                                                         aTGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGG‑TTAAAGAAAAAAGAATTTCTCCTTTTTGA‑GTGCTCCGCCTTTCCCTATTTTTTCCCCCCGTA‑TTT‑TTTAAGGGAC‑GATAaagagggaa                                                                                   >  1:299192/1‑145 (MQ=255)                                                                                atGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGG‑TTAAGGACATACGATTTCCTCCTTTCAGA‑GTGCTCCGCTTCTCACTATTATCTCACGCAGTA‑TTC‑TTAAGGGAAC‑GATAAGGAGGAACCATGAACAtt                                                                         <  1:249282/151‑1 (MQ=255)                                                                                  gTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGG‑TTTAAGAGAAAAGAATTTCTCCCTTTTGA‑GGGGTCCCCTTCTCACCATTTTTTCACCCAGTA‑TTCTTTAAGGGAA‑‑GATAAAGAGGGAACATGAACa                                                                           >  1:732260/1‑147 (MQ=255)                                                                                         tACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGG‑TTTAAGAAAAAAAATTTTCTCCTTTCCGT‑TGGGTCCGCGTCTCCCTATTAATTCACCCAGAA‑TTT‑TTTAAGGGGC‑GATAAGAAGGAACCACAAAAATTA‑cccccctt                                                               >  1:332396/1‑147 (MQ=255)                                                                                          aCAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGG‑TTTAGGACATAAAAATTTCTTCTTTTTGA‑GGGGTCCGCTTCTCACTTTTATTTCAAGCCGTA‑TTT‑TTTAGGGGGA‑GATAAAGAGGGAACATGAACATTA‑ctcccctt                                                               >  2:149591/1‑146 (MQ=255)                                                                                          aCAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGG‑TTAAGGACATACGATTTCCTCCTTTCAGA‑GTGCTCCGCTTCTCACTATTATCTCACGCAGTA‑TTC‑TTAAGGGAAC‑GATAAGGAGGAACCATGAACATTA‑CCCCg                                                                  <  1:867398/147‑1 (MQ=255)                                                                                            aGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGG‑TTAAGGAAAAACGATTTTCTCCTTTTCGA‑GGGGTCCGCCTCTCCCCATTATTTCACCCAGGA‑TTT‑TTTAGGGGAA‑GATAAAGAGGGAACAAGAAAAATACCCCCGTTTc                                                              >  1:234500/1‑150 (MQ=255)                                                                                                 ttCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGG‑GTTATGGCATAAAAGTTTTTCCTTTCTGAGGGGGTCCCCTT‑TCCCCATTTTTTTACGCCGGA‑TTT‑TTTAGGGGACGGAAAAGGGGGAACCAGAAAAATAA‑CCCCCTTTc                                                              >  2:282986/1‑145 (MQ=255)                                                                                                   ccTTAAAAGATATGTCAGGCTTGCGGAGTGGCGG‑TTTAGTACATAACAATTTCTCCTTTTAGA‑GGGGTCCGCTTCTCACTATTATCTCACCCCCGA‑TTT‑TTTATGGAAA‑GATAAGGGGGAAACATAAAAAATA‑CCCCCCTTTCCACGCttt                                                     >  2:393088/1‑151 (MQ=255)                                                                                                      tAAAAGATATGTCAGGCTTGCGGAGTGGCGG‑TTAAGGACATACGATTTCCTCCTTTCAGA‑GTGCTCCGCTTCTCACTATTATCTCACGCAGTA‑TTC‑TTAAGGGAAC‑GATAAGGAGGAACCATGAACATTA‑CCCCGTTTCCGACGCTTTc                                                    <  1:935800/149‑1 (MQ=255)                                                                                                        aaaGATATGTCAGGCTTGCGGAGTGGCGG‑TTAAAGACATAAGATTTTCTCCTTTCTGA‑GTGCTCCGCCTCTCCCTATTTTTTCACCCCGGATTTT‑TTAAGGGAA‑‑GATAAGGAGGGAACAAGAAAAATA‑CCCCCTTTCCCACGCtttttcc                                                 >  1:111858/1‑146 (MQ=255)                                                                                                            atatGTCAGGCTTGCGGAGTGGCGGCTTAAAAAC‑TACGAATTTCTCCTTTCTCA‑GTGGTCGGCCTCTCCCTATTTTATCTCGCAGTA‑TTT‑TTTAGGGGAC‑GGGTAAGAGGAACCTTGAAAATTA‑AACCCTTTTCCACGCTTTCTcc                                                 >  2:213773/1‑146 (MQ=255)                                                                                                                gTCAGGCTTGCGGAGTGGCGG‑TTAAGGACATACGATTTCCTCCTTTCAGA‑GTGCTCCGCTTCTCACTATTATCTCACGCAGTA‑TTC‑TTAAGGGAAC‑GATAAGGAGGAACCATGAACATTA‑CCCCGTTTCCGACGCTTTCGCCGGCAACta                                         <  1:726921/150‑1 (MQ=255)                                                                                                                gTCAGGCTTGCGGAGTGGCGG‑TTAAGGACATACGATTTCCTCCTTTCAGA‑GTGCTCCGCTTCTCACTATTATCTCACGCACTA‑TTC‑TTAAGGGAAC‑GATAAGGAGGAACCATGAACATTA‑CCCCGTTTCCGACGCTTTCGCCGGCAGCta                                         <  2:584702/150‑1 (MQ=255)                                                                                                                   aGGCTTGCGGAGTGGCGG‑TTTAAGAAATAAGATTTTCTCCCTTTTGA‑GGGGTCCGCCTCTCCCTTTTATTTTCCGCAGTATTTT‑TTAAGGGAA‑‑GAAAAGGAGGGAACATGAAAAATA‑ACCCCTTTTCCACGCCTTTCCCgg                                               >  2:928407/1‑141 (MQ=255)                                                                                                                   aGGCTTGCGGAGTGGCGG‑TTAAGGACATACGATTTCCTCCTTTCAGA‑GTGCTCCGCTTCTCACTATTATCTCACGCAGTA‑TTC‑TTAAGGGAAC‑GATAAGGAGGAACCATGAACATTA‑CCCCGTTTCCGACGCTTTCGCCGGCAACTATAGa                                     <  1:649632/151‑1 (MQ=255)                                                                                                                            gAGTGGCGG‑TTAAGGACATACGATTTCCTCCTTTCAGA‑GTGCTCCGCTTCTCACTATTATCTCACGCAGTA‑TTC‑TTAAGGGAAC‑GATAAGGAGGAACCATGAACATTA‑CCCCGTTTCCGACGCTTTCGCCGGCAACTATAGATGCCATaaa                            <  2:163267/151‑1 (MQ=255)                                                                                                                            gAGTGGCGG‑TTAAGGACATACGATTTCCTCCTTTCAGA‑GTGCTCCGCTTCTCACTATTATCTCACGCAGTA‑TTC‑TTAAGGGAAC‑GATAAGGAGGAACCATGAACATTA‑CCCCGTTTCCGACGCTTTCGCCGGCAACTATAGATGCCATaaa                            <  2:192217/151‑1 (MQ=255)                                                                                                                              gTGGCGG‑TTAAGGACATACGATTTCCTCCTTTCAGA‑GTGCTCCGCTTCTCACTATTATCTCACGCAGTA‑TTC‑TTAAGGGAAC‑GATAAGGAGGAACCATGAACATTA‑CCCCGTTTCCGACGCTTTCGCCGGCAACTATAGATGCCATAAATg                          <  1:967108/151‑1 (MQ=255)                                                                                                                                ggCGG‑TTAAGGACATACGATTTCCTCCTTTCAGA‑GTGCTCCGCTTCTCACTATTATCTCACGCAGTA‑TTC‑TTAAGGGAAC‑GATAAGGAGGAACCATGAACATTA‑CCCCGTTTCCGACGCTTTCGCCGGCAACTATAGATGCCATAAATGtt                        >  1:554061/1‑151 (MQ=255)                                                                                                                                 gCGG‑TTAAGGACATACGATTTCCTCCTTTCAGA‑GTGCTCCGCTTCTCACTATTATCTCACGCAGTA‑TTC‑TTAAGGGAAC‑GATAAGGAGGAACCATGAACATTA‑CCCCGTTTCCGACGCTTTCGCCGGCAACTATAGATGCCATAAATGTTa                       <  2:710412/151‑1 (MQ=255)                                                                                                                                           gACATACGATTTCCTCCTTTCAGA‑GTGCTCCGCTTCTCACTATTATCTCACGCAGTA‑TTC‑TTAAGGGAAC‑GATAAGGAGGAACCATGAACATTA‑CCCCGTTTCCGACGCTTTCGCCGGCAACTATAGATGCCATAAATGTTATCGGACAGt              <  1:189103/151‑1 (MQ=255)                                                                                                                                            aCATACGATTTCCTCCTTTCAGA‑GTGCTCCGCTTCTCACTATTATCTCACGCAGTA‑TTC‑TTAAGGGAAC‑GATAAGGAGGAACCATGAACATTA‑CCCCGTTTCCGACGCTTTCGCCGGCAACTATAGATGCCATAAATGTTATCGGACAGtg             >  1:685142/1‑151 (MQ=255)                                                                                                                                               tACGATTTCCTCCTTTCAGA‑GTGCTCCGCTTCTCACTATTATCTCACGCAGTA‑TTC‑TTAAGGGAAC‑GATAAGGAGGAACCATGAACATTA‑CCCCGTTTCCGACGCTTTCGCCGGCAACTATAGATGCCATAAATGTTATCGGACAGtgg            >  1:138973/1‑149 (MQ=255)                                                                                                                                                aCGATTTCCTCCTTTCAGA‑GTGCTCCGCTTCTCACTATTATCTCACGCAGTA‑TTC‑TTAAGGGAAC‑GATAAGGAGGAACCATGAACATTA‑CCCCGTTTCCGACGCTTTCGCCGGCAACTATAGATGCCATAAATGTTATCGGACAGtggct          <  2:855336/150‑1 (MQ=255)                                                                                                                                                 cGATTTCCTCCTTTCAGA‑GTGCTCCGCTTCTCACTATTATCTCACGCAGTA‑TTC‑TTAAGGGAAC‑GATAAGGAGGAACCATGAACATTA‑CCCCGTTTCCGACGCTTTCGCCGGCAACTATAGATGCCATAAATGTTATCGGACAGtggctgg        >  2:125433/1‑151 (MQ=255)                                                                                                                                                  gATTTCCTCCTTTCAGA‑GTGCTCCGCTTCTCACTATTATCTCACGCAGTA‑TTC‑TTAAGGGAAC‑GATAAGGAGGAACCATGAACATTA‑CCCCGTTTCCGACGCTTTCGCCGGCAACTATAGATGCCATAAATGTTATCGGACAGtggctggc       <  2:350231/151‑1 (MQ=255)                                                                                                                                                   aTTTCCTCCTTTCAGA‑GTGCTCCGCTTCTCACTATTATCTCACGCAGTA‑TTC‑TTAAGGGAAC‑GATAAGGAGGAACCATGAACATTA‑CCCCGTTTCCGACGCTTTCGCCGGCAACTATAGATGCCATAAATGTTATCGGACAGTGGCTGgcg      >  2:939464/1‑151 (MQ=255)                                                                                                                                                       cctcctTTCAGA‑GTGCTCCGCTTCTCACTATTATCTCACGCAGTA‑TTC‑TTAAGGGAAC‑GATAAAGAGGAACCATGAACATTA‑CCCCGTTTCCGACGCTTTCGCCGGCAACTATAGATGCCATAAATGTTATCGGACAGTGGCTGGCGCAgg  <  1:112357/151‑1 (MQ=255)                                                                                                                                                       |                                                                                                                                                            GTTTTTTATTAATTCTAAAACGGTAAGGGTAAAAATTCAGGAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGG‑TTAAGGACATACGATTTCCTCCTTTCAGA‑GTGCTCCGCTTCTCACTATTATCTCACGCAGTA‑TTC‑TTAAGGGAAC‑GATAAGGAGGAACCATGAACATTA‑CCCCGTTTCCGACGCTTTCGCCGGCAACTATAGATGCCATAAATGTTATCGGACAGTGGCTGGCGCAGG  >  NC_000913/1487014‑1487313

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