Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913357,8800TC43.5% 104.1 / 17.1 70K271E (AAA→GAA) cynRtranscriptional activator of cyn operon; autorepressor
Reads supporting (aligned to +/- strand):  ref base T (31/8);  new base C (0/30);  total (31/39)
Fisher's exact test for biased strand distribution p-value = 1.65e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGA‑TTCATTTCCGCC‑AA‑CAACCGCGCATTTATCCAACGCCATGTGCAAAAAT‑GCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGA  >  NC_000913/357732‑358024
                                                                                                                                                         |                                                                                                                                                  
tCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGA‑TTCATTTCCGCC‑AA‑CAACCGCGCATTTATCCAACGCCATGTGCAAAAAT‑GCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTTcc                                                                                                                                                  >  2:1083342/1‑151 (MQ=255)
 cGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGA‑TTCATTTCCGCC‑AA‑CAACCGCGCATTTATCCAACGCCATGTGCAAAAAT‑GCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTTcc                                                                                                                                                  <  3:1044853/150‑1 (MQ=255)
   tAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGA‑TTCATTTCCGCC‑AA‑CAACCGCGCATTTATCCAACGCCATGTGCAAAAAT‑GCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTTccgcc                                                                                                                                               >  3:330078/1‑151 (MQ=255)
   tAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGA‑TTCATTTCCGCC‑AA‑CAACCGCGCATTTATCCAACGCCATGTGCAAAAAT‑GCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTTccgcc                                                                                                                                               >  2:1072839/1‑151 (MQ=255)
    aGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGA‑TTCATTTCCGCC‑AA‑CAACCGCGCATTTATCCAACGCCATGTGCAAAAAT‑GCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTTccgccg                                                                                                                                              >  2:264657/1‑151 (MQ=255)
    aGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGA‑TTCATTTCCGCC‑AA‑CAACCGCGCATTTATCCAACGCCATGTGCAAAAAT‑GCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTTccgccg                                                                                                                                              >  2:577806/1‑151 (MQ=255)
       ccTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGA‑TTCATTTCCGCC‑AA‑CAACCGCGCATTTATCCAACGCCATGTGCAAAAAT‑GCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGca                                                                                                                                            >  2:991986/1‑150 (MQ=255)
         tGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGA‑TTCATTTCCGCC‑AA‑CAACCGCGCATTTATCCAACGCCATGTGCAAAAAT‑GCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGcaac                                                                                                                                          >  3:490141/1‑150 (MQ=255)
             aaGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGA‑TTCATTTCCGCC‑AA‑CAACCGCGCATTTATCCAACGCCATGTGCAAAAAT‑GCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAc                                                                                                                                     <  3:1425133/151‑1 (MQ=255)
              aGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGA‑TTCATTTCCGCC‑AA‑CAACCGCGCATTTATCCAACGCCATGTGCAAAAAT‑GCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAcc                                                                                                                                    <  2:896948/151‑1 (MQ=255)
              aGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGA‑TTCATTTCCGCC‑AA‑CAACCGCGCATTTATCCAACGCCATGTGCAAAAAT‑GCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAcc                                                                                                                                    <  2:1446198/151‑1 (MQ=255)
               gCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGA‑TTCATTTCCGCC‑AA‑CAACCGCGCATTTATCCAACGCCATGTGCAAAAAT‑GCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAcc                                                                                                                                    <  2:1232859/150‑1 (MQ=255)
                 gTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGA‑TTCATTTCCGCC‑AA‑CAACCGCGCATTTATCCAACGCCATGTGCAAAAAT‑GCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAccgcc                                                                                                                                 <  3:24959/151‑1 (MQ=255)
                     ccaataaGGCATTCCCAGCCGCAGCCTTGTGTCAGCGGCTACCGTGA‑TTCATTCCCGCC‑AA‑CACCCGCGCATTTATCCAACGCCTTGTGAAAAAAT‑CCCTCC‑GGGGCGGCTGTCTCCCAGCTTTTTTTCCGCCGCAACAAAAccgccg                                                                                                                                <  3:976774/142‑1 (MQ=18)
                         tCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGA‑TTCATTTCCGCC‑AA‑CAACCGCGCATTTATCCAACGCCATGTGCAAAAAT‑GCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTtctctc                                                                                                                         >  2:1422117/1‑151 (MQ=255)
                           aGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGA‑TTCATTTCCGCC‑AA‑CAACCGCGCATTTATCCAACGCCATGTGCAAAAAT‑GCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTcc                                                                                                                        >  3:303283/1‑150 (MQ=255)
                           aGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGA‑TTCATTTCCGCC‑AA‑CAACCGCGCATTTATCCAACGCCATGTGCAAAAAT‑GCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTcc                                                                                                                        >  2:1591011/1‑150 (MQ=255)
                             gCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGA‑TTCATTTCCGCC‑AA‑CAACCGCGCATTTATCCAACGCCATGTGCAAAAAT‑GCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCagt                                                                                                                     >  2:115315/1‑151 (MQ=255)
                             gCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGA‑TTCATTTCCGCC‑AA‑CAACCGCGCATTTATCCAACGCCATGTGCAAAAAT‑GCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCagt                                                                                                                     >  3:678947/1‑151 (MQ=255)
                               aaTTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGA‑TTCATTTCCGCC‑AA‑CAACCGCGCATTTATCCAACGCCATGTGCAAAAAT‑GCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCagt                                                                                                                     >  3:1626660/1‑149 (MQ=255)
                                aTTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGA‑TTCATTTCCGCC‑AA‑CAACCGCGCATTTATCCAACGCCATGTGCAAAAAT‑GCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCagtagt                                                                                                                  >  2:619217/1‑151 (MQ=255)
                                aTTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGA‑TTCATTTCCGCC‑AA‑CAACCGCGCATTTATCCAACGCCATGTGCAAAAAT‑GCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCagtagt                                                                                                                  >  2:752390/1‑151 (MQ=255)
                                 ttCCAGCCGCAGCCCTGTGTCAGCGGCTCCCGTGA‑TTCATTCCCGCC‑AA‑CACCCGCGCTTTTACCCAACGCCATGTGAAAAAAT‑GCCTTC‑GCGGCGGCTGTCTCCCATTTATTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCagtagt                                                                                                                  <  2:394864/150‑1 (MQ=18)
                                 ttCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGA‑TTCATTTCCGCC‑AA‑CAACCGCGCATTTATCCAACGCCATGTGCAAAAAT‑GCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTg                                                                                                                 >  3:1270600/1‑151 (MQ=255)
                                 ttCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGA‑TTCATTTCCGCC‑AA‑CAACCGCGCATTTATCCAACGCCATGTGCAAAAAT‑GCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTg                                                                                                                 >  3:536020/1‑151 (MQ=255)
                                      gCCGCAGACCTGTGTCAGCGGCTACCGTGA‑TTCATTTCCGCC‑AA‑CAACCGCGCATTTATCCAACGCCATGTGCAAAAAT‑GCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCggg                                                                                                            >  3:636234/1‑151 (MQ=255)
                                       ccGCAGCCTTGTTCCAGCGCCTCCCGTG‑‑TTCCTTTCCCCCAAA‑CACCCGCGCTTTTACCCACCGCAGTGTGAAAAAAT‑GCCTTC‑GGGGGGGCTGTCTGCCGGTTTTTTTCCCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCgggg                                                                                                           <  3:70444/151‑1 (MQ=255)
                                       ccGCAGCCTTGTGTCGGGGCCTCCCGTTA‑TTCTTCCCCCCC‑AA‑CACCCCCTCTTTTCCCCACCGCCTTGCAAAAAAAT‑GCCTTC‑GGGGGGTCTGTCTGCCTGCTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCgggg                                                                                                           <  3:165028/151‑1 (MQ=255)
                                       ccGCAGACCTGTGTCAGCGGCTACCGTGA‑TTCATTTCCGCC‑AA‑CAACCGCGCATTTATCCAACGCCATGTGCAAAAAT‑GCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggg                                                                                                           >  2:94718/1‑151 (MQ=255)
                                        cGCAGCCTTGTGCCGGCGGCT‑CCGTTA‑TTTCTTTCCCCC‑AACCACCCGCGCTTTTCCCAACCCCCTTTTAAAAAAAT‑GCTTCC‑GCGGCGGTTGTCGCCCACTTTTTTTCCCGCCCCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGc                                                                                                          <  2:475221/151‑1 (MQ=255)
                                        cGCAAACCTGTGGCAGCGCCT‑CCGTTATTTCTTCCCCGCA‑AA‑CAACCGCCCTTTTTCCCACCGCCTTGTCCAAAAGT‑GCTTCG‑GGGGCGGCTGTGTGCCGGTTTTTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGc                                                                                                          <  3:71883/151‑1 (MQ=255)
                                         gCAGACCTGTGTCAGCGGCTACCGTGA‑TTCATTTCCGCC‑AA‑CAACCGCGCATTTATCCAACGCCATGTGCAAAAAT‑GCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggg                                                                                                           >  2:784155/1‑149 (MQ=255)
                                         gCAGACCTGTGTCAGCGGCTACCGTGA‑TTCATTTCCGCC‑AA‑CAACCGCGCATTTATCCAACGCCATGTGCAAAAAT‑GCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggg                                                                                                           >  2:995855/1‑149 (MQ=255)
                                            gACCTGTGTCGGCGCCTCCCGTGT‑TTCATTCCCGCC‑AA‑CACCCGCGCATTTACCCACCGCCTTGTGAAAAAAT‑GCCTTC‑GCGGGGGTTGTCGCCCAGTTTTTTTCCCCCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCaa                                                                                                        <  2:903562/149‑1 (MQ=255)
                                                tgtgTCAGCGGCTACCGTGA‑TTCATTTCCGCC‑AA‑CAACCGCGCATTTATCCAACGCCATGTGCAAAAAT‑GCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggg                                                                                                           >  2:396177/1‑142 (MQ=255)
                                                tgtgTCAGCGGCTACCGTGA‑TTCATTTCCGCC‑AA‑CAACCGCGCATTTATCCAACGCCATGTGCAAAAAT‑GCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGA‑Gaa                                                                                                  >  2:322155/1‑150 (MQ=255)
                                                tgggccgggggCTCCCGTGT‑TTCTTTCCCCCC‑AA‑CACCCGCGCTTTTACCCACCCCCATGTGAAAAATT‑CCCTCC‑GCGGCGGTTGTCTGCAAGTTTTTTTCCCGCCACAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑Gaa                                                                                                  <  3:148192/141‑1 (MQ=255)
                                                     cAGCGGCTACCGTGA‑TTCATTTCCGCC‑AA‑CAACCGCGCATTTATCCAACGCCATGTGCAAAAAT‑GCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGCAGA‑Gaaaaaaa                                                                                             >  2:669123/1‑146 (MQ=255)
                                                       gCGGCTACCGTGT‑TTCATTCCCCCC‑AA‑CCCCCGCGCTTTTCCCCACCGCCATGGAAAAAAAT‑GCCTGC‑GGGGGGGCTGCCCGCCATTTTTTTTCCCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGc                                                                                             <  3:1105348/148‑1 (MQ=255)
                                                         ggCCACCGTAA‑TTTTTTCCCGCC‑AA‑CACCCGCGCTTTTATCAACCGACTTGTGAAAAAAT‑GCCTTC‑GCGGGGGTTGTCTGCCATTTTTTTTGCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGC‑tt                                                                                          <  2:1397977/148‑1 (MQ=255)
                                                          gCTACCGTGA‑TTCATTTCCGCC‑AA‑CAACCGCGCATTTATCCAACGCCATGTGCAAAAAT‑GCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGCAGG‑Gaaaaaaaattt                                                                                         >  2:1147582/1‑141 (MQ=25)
                                                              ccgtgttttCTTTTCCGCC‑AA‑CACCCGCGCATTTGCCAACCGCCATATACAAAATT‑GCCTCC‑GCGGCGACTCCCTCCAACCTTTTTTCCCGCCACAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGC‑TTTaa                                                                                       <  2:1067470/140‑1 (MQ=255)
                                                              ccGTTA‑TTCTTTCCCGCC‑AA‑CACCCCCCCATTTATCCACCCCATTGGGCAAAAAC‑GCCTCC‑GCGGCGGTTGTCTCCCAGCTTTTTTCCCGCCGCAATAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGC‑t                                                                                           <  2:1634742/142‑1 (MQ=255)
                                                               cGTGA‑TTCATTTCCGCC‑AA‑CAACCGCGCATTTATCCAACGCCATGTGCAAAAAT‑GCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGA‑Ggaaaaagatatagca                                                                                     >  3:991055/1‑136 (MQ=25)
                                                                     ttCATTCCCGCC‑AA‑CACCCGCGCTTTTACCCACCGCCATGTGAAAAAAT‑GCCTCC‑GCGGCGGCTGTCTGCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGC‑TTTAAGCCCGTCAtg                                                                             <  2:884911/150‑1 (MQ=255)
                                                                               cc‑AA‑CACCCGCGTTTTTACCAACCGCCATGTAAAAAATT‑GCCTGC‑GGGGGGGCTTTCTGCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCaa                                                                  <  2:1595156/151‑1 (MQ=255)
                                                                                       cccGCGCTTTTACCCACCCCAATGTGAAAAAAT‑GCCTCC‑GCGGGGGCTTTCTCCCAGTTTTTTTCCCGCCGAAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATc                                                                <  3:1335253/146‑1 (MQ=255)
                                                                                          gcgcTTTTACCAACCGCAATGTG‑AAAAAA‑GCCTTCGGCGGCGGCTGTCTCCCAGTTTTTTTCCCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCg                                                               <  2:377633/145‑1 (MQ=255)
                                                                                           cgtaTTTACCCACCCCCTTGT‑CAAAAATGGCCTTC‑GGGGCGGCTTTCTGCCAGTTTTTTTCCCCCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTgg                                                        <  3:1023012/148‑1 (MQ=255)
                                                                                                    ccacccccATGTAAAAAAAT‑GCCTCC‑GGGGCGGCTGTCGCCCAGTTTTTTTCCCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAAtg                                                <  3:1157216/144‑1 (MQ=255)
                                                                                                     cAACGCCATGTGCAAAAAT‑GCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGA‑GAAAAAAA‑TTTTAACCCCCCCTGTTTGGTGGGGAACGCTGCTgggtaaaa                                                  >  1:10347/1‑141 (MQ=25)
                                                                                                        cGCCATGTGAAAAAAT‑GCCTCC‑GCGGCGGCTGTCTCCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGaa                                           <  2:52696/151‑1 (MQ=255)
                                                                                                         gCCATGTGCAAAAAT‑GCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGCAGG‑GAAAAAACTTTTAAACCCG‑CCTATTTTGTGGGGAACGCTGCTGGttgcaa                                                  >  2:730861/1‑138 (MQ=21)
                                                                                                               tgCAAAAAT‑CCCTTG‑GCGGCGGCTGTCTGCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGaa                                     <  2:155345/150‑1 (MQ=255)
                                                                                                               gtttaataaggCCCTC‑GCGGCGGCGGTCTGCCAGCTATTTTTCCGCCATAAGAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAAGCGCTGCTGGTAACAATGTGGAAAGGGaa                                     <  2:1379316/141‑1 (MQ=255)
                                                                                                                  aaaatT‑CCCTCC‑GCGGCGGCTGTCTGCCAGTTTTTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTg                                  <  2:922209/150‑1 (MQ=255)
                                                                                                                  aaaagT‑GCCTTC‑GCGGCGTCTGCCGGCCGGTTTTTTTCCCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTg                                  <  3:1534811/150‑1 (MQ=255)
                                                                                                                  aaaaaT‑CCCTTC‑GGGGGGGTTGTCTGCCAGTTTTTTTCCCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAggg                                       <  2:660253/145‑1 (MQ=255)
                                                                                                                    aatt‑GCCTCC‑GCGGGGGTTGTCTCCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTgcgg                               <  2:506257/148‑1 (MQ=255)
                                                                                                                    aact‑GCCTCG‑GCGGCGGCTTCCTGCAAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAg                                    <  3:1343315/143‑1 (MQ=255)
                                                                                                                    aaaaccccTCC‑GGGGCGCCTGTTTGCCATCTTTTTTCCCGCCCAAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTgcg                                <  2:1071955/145‑1 (MQ=255)
                                                                                                                     aattcccTTC‑GGGGCGGTTTTCTCCAAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTgcgg                               <  3:1253449/146‑1 (MQ=255)
                                                                                                                      aT‑GCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAAA‑GAAAAAAA‑TTTTAACCCCCCCTGTTTTGTTGGGAACGCCGCTGGGTACAAAGTGGGAAGGGGAgggggggg                             >  2:813301/1‑144 (MQ=255)
                                                                                                                      aT‑GCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCGAGAGGAAAAAAT‑TTTAACCCCCCCAT‑TTTTGTGGGGAACCCTGCTGGGAAAAAAATGGGAAggggaaggggg                               >  2:1181885/1‑141 (MQ=255)
                                                                                                                          ccTTC‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGA‑GAAAAAAC‑ATTAAATCCGCCATGTTTTGTGGGAAATCCTGCCGGTGAAAAAATGGGAAgggggggggggggg                            >  2:554405/1‑138 (MQ=255)
                                                                                                                                              cctttttttttcCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCgcgc         <  3:847256/142‑1 (MQ=255)
                                                                                                                                              ccAGTTTTTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTa       <  2:927245/151‑1 (MQ=255)
                                                                                                                                              ccAGTTTTTTTCCCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAcc             <  3:1442934/145‑1 (MQ=255)
                                                                                                                                              ccAGTTTTTTTCCCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACcgcg          <  2:1199181/148‑1 (MQ=255)
                                                                                                                                              ccAGTTTTTTTCCCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCGGCTCCAGAACCGCGCTa       <  2:166166/151‑1 (MQ=255)
                                                                                                                                              ccAGTTTTTTTCCCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACcg            <  2:490176/146‑1 (MQ=255)
                                                                                                                                              caattttttttcCCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGaa               <  3:1144159/136‑1 (MQ=255)
                                                                                                                                              caagtttttttcCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACcgc           <  3:881320/140‑1 (MQ=255)
                                                                                                                                                     ttttcCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGAGCGGCGAATCAGCTCCAGAACCGCGCTAATTGa  <  2:457546/149‑1 (MQ=255)
                                                                                                                                                     ttttcCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACcgcg          <  3:371708/141‑1 (MQ=255)
                                                                                                                                                         |                                                                                                                                                  
TCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGA‑TTCATTTCCGCC‑AA‑CAACCGCGCATTTATCCAACGCCATGTGCAAAAAT‑GCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGA  >  NC_000913/357732‑358024

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