Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913357,8800TC50.0% 91.4 / 37.5 84K271E (AAA→GAA) cynRtranscriptional activator of cyn operon; autorepressor
Reads supporting (aligned to +/- strand):  ref base T (35/7);  new base C (0/42);  total (35/49)
Fisher's exact test for biased strand distribution p-value = 1.02e-16
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAA‑CAACCGCGCATTTA‑TCCAACGCCAT‑GTGCAAAAA‑T‑G‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATGTTGT‑GTGGCAATCGCTGCT‑GGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTC  >  NC_000913/357732‑358028
                                                                                                                                                           |                                                                                                                                                          
tCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAA‑CAACCGCGCATTTA‑TCCAACGCCAT‑GTGCAAAAA‑T‑G‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑cc                                                                                                                                                         >  2:849935/1‑151 (MQ=255)
 cGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAA‑CAACCGCGCATTTA‑TCCAACGCCAT‑GTGCAAAAA‑T‑G‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑ccg                                                                                                                                                        >  2:1078544/1‑151 (MQ=255)
 cGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAA‑CAACCGCGCATTTA‑TCCAACGCCAT‑GTGCAAAAA‑T‑G‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑cc                                                                                                                                                         <  2:957394/150‑1 (MQ=255)
    aGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAA‑CAACCGCGCATTTA‑TCCAACGCCAT‑GTGCAAAAA‑T‑G‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑ccgccg                                                                                                                                                     >  2:938952/1‑151 (MQ=255)
     ggCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAA‑CAACCGCGCATTTA‑TCCAACGCCAT‑GTGCAAAAA‑T‑G‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑Ccgccgc                                                                                                                                                    >  1:336945/1‑151 (MQ=255)
         tGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAA‑CAACCGCGCATTTA‑TCCAACGCCAT‑GTGCAAAAA‑T‑G‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGcaac                                                                                                                                                 >  2:1673857/1‑150 (MQ=255)
          gATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAA‑CAACCGCGCATTTA‑TCCAACGCCAT‑GTGCAAAAA‑T‑G‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGcaaca                                                                                                                                                <  1:485522/150‑1 (MQ=255)
          gATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAA‑CAACCGCGCATTTA‑TCCAACGCCAT‑GTGCAAAAA‑T‑G‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGcaaca                                                                                                                                                >  2:1141075/1‑150 (MQ=255)
           aTAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAA‑CAACCGCGCATTTA‑TCCAACGCCAT‑GTGCAAAAA‑T‑G‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGcaacaa                                                                                                                                               >  2:1496254/1‑150 (MQ=255)
             actcgcaGCGCACCAGGCAATTCCAGCCGCAGACCTGTTTCTTCGGCTACCGTGATTCATTTCCGCCAA‑CACCCGCGCATTTA‑TCCAACGCCAT‑GTGCAAACA‑T‑G‑CATT‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACaaaa                                                                                                                                             <  2:444192/144‑1 (MQ=255)
               gCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAA‑CAACCGCGCATTTA‑TCCAACGCCAT‑GTGCAAAAA‑T‑G‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAA‑‑ccg                                                                                                                                        <  2:686945/151‑1 (MQ=255)
                cGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAA‑CAACCGCGCATTTA‑TCCAACGCCAT‑GTGCAAAAA‑T‑G‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAA‑‑ccgc                                                                                                                                       <  2:271784/151‑1 (MQ=255)
                      gcATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAA‑CAACCGCGCATTTA‑TCCAACGCCAT‑GTGCAAAAA‑T‑G‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAA‑‑CCGCCGTtct                                                                                                                                 >  1:300182/1‑151 (MQ=255)
                       cATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAA‑CAACCGCGCATTTA‑TCCAACGCCAT‑GTGCAAAAA‑T‑G‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAA‑‑CCGCCGTtctc                                                                                                                                >  2:1175717/1‑151 (MQ=255)
                        aTCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAA‑CAACCGCGCATTTA‑TCCAACGCCAT‑GTGCAAAAA‑T‑G‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAA‑‑CCGCCGTtctct                                                                                                                               >  1:1761495/1‑151 (MQ=255)
                          cAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAA‑CAACCGCGCATTTA‑TCCAACGCCAT‑GTGCAAAAA‑T‑G‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAA‑‑CCGCCGTTCTCTcc                                                                                                                             >  1:501054/1‑151 (MQ=255)
                              cAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAA‑CAACCGCGCATTTA‑TCCAACGCCAT‑GTGCAAAAA‑T‑G‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAA‑‑CCGCCGTTCTCTCCagt                                                                                                                          >  1:221927/1‑150 (MQ=255)
                              cAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAA‑CAACCGCGCATTTA‑TCCAACGCCAT‑GTGCAAAAA‑T‑G‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAA‑‑CCGCCGTTCTCTCCagt                                                                                                                          >  2:364711/1‑150 (MQ=255)
                              cAATTCCAGCCGCAGACCTGTGTCAGCGGCTACAGTGATTCATTTCCGCCAA‑CAACCGCGCATTTA‑TCCAACGCCAT‑GTGCAAAAA‑T‑G‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAC‑‑CCGCCGAtctct                                                                                                                               >  2:886532/1‑145 (MQ=255)
                               aaTTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAA‑CAACCGCGCATTTA‑TCCAACGCCAT‑GTGCAAAAA‑T‑G‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAA‑‑CCGCCGTTCTCTCCagt                                                                                                                          >  2:282525/1‑149 (MQ=255)
                                   ccgcccGCGGCCCTGTGTCAGCGGCTCCCGTGTTTCATTCCCGCCAC‑CACCCGCGCATTTA‑CCCACCGCCTT‑GTGAAAAAA‑T‑G‑CCTC‑CGCGGCGGCTTTCTGCCAGTTTTTTTC‑CCGCCGCAACAAAA‑‑CCGCCGTTCTCTCCAGTAGTgg                                                                                                                     <  1:1360504/146‑1 (MQ=255)
                                   ccgcccGCAGCCCTGTGTCAGCGCCTCCCGTATTTCATTCCCCCCAC‑CACCCGCGCTTTTA‑CCCACCCCCAT‑GTGAAAAAA‑T‑C‑CTTC‑CGCGGGGGCTGTCTGCCAGTTTTTTTC‑CCGCCGCAAAAAAA‑‑CCGCCGTTCTCTCCAGTAGTgg                                                                                                                     <  1:621951/146‑1 (MQ=255)
                                   ccAACCGCAACCCTGTCTCGCCGCCTACGGTGATTCATTCCCTCCAA‑CAGCCGCGCTTTTC‑CCCAACCCTAT‑GTGCAAAAA‑T‑C‑CCTC‑CGGGGGGGCTTCCTCCCAGTTTTTTTC‑CCCCCCCAAAAAAG‑‑CCGCCGTTCTCTCCAGTAGTGGc                                                                                                                    <  1:1208048/151‑1 (MQ=255)
                                       ccGCAGCCCTGTGTCAGGGGCTCCCGTGATTCTTTCCCCCCAA‑CACCCGCTCTTTTA‑CCCACCGCCTT‑GTGAAAAAC‑T‑G‑CCTC‑CGCGGCGGCTGTCTGCCAGTTTTTTTC‑CCGCCGCAAAAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCggg                                                                                                                 <  1:1184444/150‑1 (MQ=255)
                                       ccGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAA‑CAACCGCGCATTTA‑TCCAACGCCAT‑GTGCAAAAA‑T‑G‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCggg                                                                                                                 >  1:1618270/1‑150 (MQ=255)
                                         gCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAA‑CAACCGCGCATTTA‑TCCAACGCCAT‑GTGCAAAAA‑T‑G‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCgggg                                                                                                                >  1:355433/1‑149 (MQ=255)
                                          cAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAA‑CAACCGCGCATTTA‑TCCAACGCCAT‑GTGCAAAAA‑T‑G‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCgggg                                                                                                                >  1:1358088/1‑148 (MQ=255)
                                           accccTGTGTCAGCGTCTCCCGTTTTTCTTTCCCCCCAA‑CACCCGCGCTTTTC‑CCCACCCCCAT‑GTAAAAAAA‑T‑G‑CTTT‑CGGGGCGGTTGTCTGCCAGTTTTTTTC‑CCCCCGCAAAAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAag                                                                                                            <  1:391625/148‑1 (MQ=255)
                                            gACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAA‑CAACCGCGCATTTA‑TCCAACGCCAT‑GTGCAAAAA‑T‑G‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCggggg                                                                                                               >  2:1459979/1‑146 (MQ=255)
                                             aCCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAA‑CAACCGCGCATTTA‑TCCAACGCCAT‑GTGCAAAAA‑T‑G‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCgggg                                                                                                                >  1:690899/1‑145 (MQ=255)
                                              ccTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAA‑CAACCGCGCATTTA‑TCCAACGCCAT‑GTGCAAAAA‑T‑G‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCgggggggg                                                                                                            >  2:243700/1‑144 (MQ=255)
                                                  ttccAGCGCCTCCCGTATTTCATTCCCCCAAA‑CACCCGCGCATTTT‑CCCAACGCCAT‑TTGCAAAAT‑T‑G‑CCTT‑CGCGGCGGCTGTCTCCCAGCTTTTTTC‑CCGCCGCAAAAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAAt                                                                                                      <  1:1132613/147‑1 (MQ=255)
                                                  tgtCAGCGGCTACCGTGATTCATTTCCGCCAA‑CAACCGCGCATTTA‑TCCAACGCCAT‑GTGCAAAAA‑T‑G‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGGAGGAGaaaa                                                                                                      >  1:258826/1‑149 (MQ=255)
                                                   gccggggcctcccGTTATTCTTTCCCCCAAA‑CACCCGCCCTTTTA‑CCCACCCCCAG‑GTGAAAAAA‑G‑G‑CCTC‑CGGGGCGGCTGTCGCCCAGTTTTTTCC‑CCGCCGCAAAAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAgaga                                                                                                         <  2:87153/135‑1 (MQ=255)
                                                     cggcgcCTCCCGTGTTTCTTTCCCCCAAA‑CACCCCCGCTTTTA‑CCCACCGCCGT‑GTAAAAAAA‑T‑C‑CCTC‑CGCGGCGGCTTTCCGCCATTTTTTTTC‑CCGCCGAAAAAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCaa                                                                                                             <  2:1342987/138‑1 (MQ=255)
                                                     cAGCGGCTACCGTGATTCATTTCCGCCAA‑CAACCGCGCATTTA‑TCCAACGCCAT‑GTGCAAAAA‑T‑G‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGGgaaagaaaaa                                                                                                     >  1:150684/1‑146 (MQ=255)
                                                         ggCTACCGTGATTCTTTCCCCCCAA‑CACCCCCGCTTTTA‑CCCACCCCCAT‑G‑GCAAAAATT‑C‑CCTC‑CGCGGCGGTTGTCGCCCAGCTTTTTTC‑CCCCCGAAAAAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGC‑TTTaa                                                                                             <  1:1674945/151‑1 (MQ=255)
                                                         gcttcccGTGTTTCTTTCCCCCAAA‑CACCCGCGCTTTTA‑CCCACCGCCAT‑GTCAAAAAT‑T‑G‑CTTT‑CGGGGGGGTTGTCTGCCAGTTTTTTTC‑CCGCCGCAACAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGC‑TTTa                                                                                              <  1:121229/145‑1 (MQ=255)
                                                         gcctgccGGGATTTTTTTCCCCCAA‑CACCCGCG‑ATTTACTCCACCGCTAT‑GTGAAAAAG‑T‑C‑CCTT‑CGCGGCGTCTGT‑TGGCAGGTTTTTTTCCCGCCACAAAAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGc                                                                                                   <  1:1726483/141‑1 (MQ=255)
                                                         gcctcccGTGTTTCATTCCCCCAAC‑CCCCCGCCCTTTTA‑CCCAACCCCTT‑GTAAAAAAA‑T‑C‑CCTG‑CGGGGCGGCTGTCGGCAACTTTTTTTC‑CCGCCGAAAAAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGaaa                                                                                                       <  2:1551206/137‑1 (MQ=255)
                                                          ggcaccgttaattctttCCCCCAACCACCCGCCCTTTTC‑CCCACCCCCAT‑GTGAAAAAA‑G‑C‑CCTC‑CGGGGCGGCT‑TCCGGCCGCTTTTTTTCCCCCCACAAAAAAG‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGC‑ttt                                                                                               <  1:145877/135‑1 (MQ=255)
                                                              ccgtttttCTTTCCCCCCAA‑CACCCGCGCTTTTA‑TCCACCCCATT‑GTGCAAAAA‑T‑G‑CCTT‑CGCGGCGGTTTTCTGCCAGTTTTTTTC‑CCGCCCAAAAAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGC‑TTTAAGCCCg                                                                                        <  1:396902/145‑1 (MQ=255)
                                                              ccgtttttCATTCCCCCAAA‑CACCCGC‑CATTTC‑CCCCACCCCATGGTAAAAAAT‑T‑G‑CCTG‑CGGGGCGGTTTTCTCCCAGTTTTTTTC‑CCGCCGAAAAAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGC‑TTTAAGcc                                                                                          <  2:1490047/143‑1 (MQ=255)
                                                              ccGTGTTTCTTTCCCCCCAA‑CACCCGCGCTTTTA‑CCCACCCCCAT‑GTGAAAAAA‑T‑G‑CCTG‑CGGGGGGGTTGTCTGCCAGTTTTTTTC‑CCGCCGAAAAAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGC‑TTTAAGCCCg                                                                                        <  1:895558/151‑1 (MQ=255)
                                                              ccGTGATTCATTTCCGCCAA‑CAACCGCGCATTTA‑TCCAACGCCAT‑GTGCAAAAA‑T‑G‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGGAGGCGaaaaaaactttaaaacc                                                                                         >  1:1112961/1‑137 (MQ=25)
                                                              ccGTGATTCATTCCCCCCAC‑CACCCGCGCTTTTA‑CCCACCGCCAT‑GTGCAAAAA‑T‑G‑CTTT‑CGGGGCGGCTGTCTCCCAGCTTTTTTC‑CCGCCGAAAAAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGC‑TTTAAGCCCg                                                                                        <  1:1570986/151‑1 (MQ=255)
                                                                gTGATTCTTTTCCCCAAA‑CACCCGCGCTTTTA‑CCCAACGCCAT‑TTGAAAAAA‑T‑C‑CCTT‑CGCGGCGGTTTTCTGCCAGTTTTTTTC‑CCGCCGCAAAAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGC‑TTTAAGCCCg                                                                                        <  1:1500348/149‑1 (MQ=255)
                                                                gTGATTCATTTCCGCCAA‑CAACCGCGCATTTA‑TCCAACGCCAT‑GTGCAAAAA‑T‑G‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGGGGGAGaaaaaaattttaaccc                                                                                          >  1:833361/1‑135 (MQ=25)
                                                                   tttCATTCCCCCAAA‑CACCCGCTCTTTTA‑CCCACCCCCAT‑GT‑AAAAAA‑T‑T‑CCTTGCGCGGGGGCTGTCCGCCAGTTTTTTTC‑CCGCCGCAAAAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGC‑TTTAAGCCCGTCa                                                                                     <  2:1278433/148‑1 (MQ=255)
                                                                    ttCTTTCCCCCCAA‑CACCCGCCCTTTTA‑CCCACCGCCAT‑GTGAAAAAA‑T‑C‑CCTC‑CGGGGCGGCTGTCTGCAAGTTTTTTTC‑CCGCCGAAAAAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGC‑TTTAAGCCCGTCa                                                                                     <  1:1754055/148‑1 (MQ=255)
                                                                    ttCATTCCCCCAAA‑CACCCGCGCTTTTA‑CCCACCGCCTT‑GTAAAAAAA‑T‑C‑CCTG‑CGGGGCGGCTTTCTCCCAGTTTTTTTC‑CCGCCGAAAAAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGC‑TTTAAGCCCGTCAt                                                                                    <  1:814042/149‑1 (MQ=255)
                                                                        tttCCGCCAA‑CAACCGCGCATTTA‑TCCAACGCCAT‑GTGCAAAAA‑T‑G‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCGCGAGGAAAAGG‑ttttaatcccccccggt                                                                                 >  2:92223/1‑134 (MQ=25)
                                                                        tccccGCCAA‑CACCCTCCAATTTA‑GCCAGCGCCAT‑GTCAAAAAC‑T‑G‑CCTC‑CGCGGTGGTTGTCTCCACGTTATTTCC‑CCGCGCAAAAAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGC‑TTTAAGCCCGTCAtgttgt                                                                               <  1:1213426/147‑1 (MQ=255)
                                                                          ccccccAA‑CACCCGCGCTTTTA‑CCCAACGCCTT‑TTGAAAAAA‑G‑C‑CCTC‑CGCGGCGGCGTTCTCCCTGTTTTTTTC‑CCGCCACAACAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATGTtgt‑gtg                                                                           <  1:938787/147‑1 (MQ=255)
                                                                           ccccaaa‑caccCGCGTTTTTA‑CCCACCGCCAG‑GTCAAAAAA‑T‑C‑CCTC‑CGGGGCGGCTGTCTGCCAGTTTTTTTC‑CCGCCAAAAAAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATGTTGT‑GTgg                                                                          <  2:609523/146‑1 (MQ=255)
                                                                           ccccaaa‑caccCGCGCTTTTA‑CCCACCGCCAT‑GTAAAAAAA‑G‑G‑CCTT‑CGCGGGGGTTTTCTGCCGGTTTTTTTC‑CCGCCGCAAAAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATGTTGT‑GTgg                                                                          <  1:197563/146‑1 (MQ=255)
                                                                                      ccGCGCATTTA‑TCCAACGCCAT‑GTGCAAAAA‑T‑G‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCCGAAGAAAAAACATTTAAACCCG‑CAtgttgttgtg                                                                           >  2:1530859/1‑137 (MQ=25)
                                                                                       cgcataTTTA‑CCCACCGCC‑T‑TTGAAAAAA‑T‑GCCCTT‑GGGGGGGGTTGTCGCCCAGTTTTTTTC‑CCGCCTAAAAAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATGTTGT‑GTGGCAATCgctgct                                                               <  2:840735/146‑1 (MQ=255)
                                                                                         cgcATTTA‑TCCAACGCCAT‑GTGCAAAAA‑T‑G‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGGGAAAGAAAAAAA‑TTTAAAACCGCCCTGTTTT‑GTGGGGGTAGctgct‑g                                                             >  2:1356024/1‑150 (MQ=25)
                                                                                             tttA‑CCCACCGCCAT‑GTAAAAAAT‑T‑G‑CTTC‑CGGGGGGGCTTTCCCCAAGTTTTTTTC‑CCGCCGCAACAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATGTTGT‑GTGGCAATCGCTGCT‑GGTa                                                          <  1:1385519/149‑1 (MQ=255)
                                                                                                    caacccccttgccaaaaag‑T‑C‑CCTC‑CGGGGCGCTTGTCCCCCAGTTTTTTTC‑CCCCCGAAAAAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATGTTGT‑GTGGCAATCGCTGCT‑GGTAAc                                                        <  2:833322/132‑1 (MQ=255)
                                                                                                     accGCCAT‑GGGAAAAAA‑G‑G‑CCTT‑CGGGGCGGCTGTCTGCCAGTTTTTTTC‑CCGCCGAAAAAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATGTTGT‑GTGGCAATCGCTGCT‑GGTAACAAtgt                                                   <  1:849974/147‑1 (MQ=255)
                                                                                                         ccAT‑GTGAAAAAG‑T‑C‑CCTC‑CGCGGCGGCTGTCTGCCAGTTTTTTTC‑CCGCCGCAAAAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATGTTGT‑GTGGCAATCGCTGCT‑GGTAACAATGTGGaa                                               <  1:1175796/149‑1 (MQ=255)
                                                                                                         caaT‑GTGAAAAAA‑C‑G‑CTTC‑CGGGGGGGCTTTCCCCCAGTTTTTTTC‑CCCCCGAAAAAAAC‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATGTTGT‑GTGGCAATCGCTGCT‑GGTAACAATGTGGaaa                                              <  2:1244408/148‑1 (MQ=255)
                                                                                                          cAT‑GTGCAAAAA‑T‑G‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGGCAGAGAAAAAAA‑TTTAAAACCCCCCTGTTGT‑GTGGGGATCGCTGCT‑GGGTACAATGTgg                                                 >  1:155441/1‑146 (MQ=37)
                                                                                                           aT‑GTGCAAAAA‑T‑G‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGGCAGAGAAAAAAC‑TTTTAAACCCTCCTGTTGT‑GTGGGGAACGCTGCT‑GGGGACAAAGTGGGAAgg                                            >  1:151766/1‑150 (MQ=255)
                                                                                                            t‑gtgCAAAAA‑T‑G‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAAAAAAAAA‑ATTTAAACCGCCCTGCTGT‑TTGGGAAACGCCGCT‑GGGAACAAtgatg                                                 >  2:1421831/1‑141 (MQ=255)
                                                                                                            t‑gtgAAAAAA‑T‑G‑CCTC‑CGGGGCGGTTTCCTGCCAGTTATTTTC‑CCGCCGAAAAAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATGTTGT‑GTGGCAATCGCTGCT‑GGTAACAATGTGGAAAggg                                           <  2:1037665/150‑1 (MQ=255)
                                                                                                               taaaaaaa‑T‑C‑CCTC‑CGGGGGGGCTGTCGCCAAGTTTTTTTC‑CCGCCAAAAAAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATGTTGT‑GTGGCAATCGCTGCT‑GGTAACAATGTGGaaa                                              <  1:774948/142‑1 (MQ=255)
                                                                                                               gtaaaaaa‑TCG‑CCTT‑CGGGGCGGCTGTCTCCCATTTTTTTTC‑CCCCCAAAAAAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATGTTGT‑GTGGCAATCGCTGCT‑GGTAACAATGTGGAAAGGGaa                                         <  2:1626142/148‑1 (MQ=255)
                                                                                                                 aaaaat‑T‑G‑CTTC‑CGGGGGGGCTGTCTCCCAGTTTTTTTC‑CCGCCGAAAAAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATGTTGT‑GTGGCAATCGCTGCT‑GGTAACAATGTGGAAAGGGAAGTg                                      <  1:1460675/150‑1 (MQ=255)
                                                                                                                  aaaaa‑T‑G‑CCTC‑CGCGGCGGTTTTCTCCCAGTTTTTTTC‑CCGCCGAAAAAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATGTTGT‑GTGGCAATCGCTGCT‑GGTAACAATGTGGAAAGGGaa                                         <  2:524827/147‑1 (MQ=255)
                                                                                                                  aaaaa‑G‑C‑CCTG‑CGCGGCGGCTGTCGCCCAGTTTTTTTC‑CCGCCGAAAAAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATGTTGT‑GTGGCAATCGCTGCT‑GGTAACAATGTGGAAAGGGa                                          <  2:1367030/146‑1 (MQ=255)
                                                                                                                      aaaagtccttccggggcggcTGTCGCCCAGTTTTTTTC‑CCGCCGCAAAAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGGAATGGC‑TTTAAGCCCGTCATGTTGT‑GTGGCAATCGCTGCT‑GGTAACAATGTGGAAAGGGAAGt                                       <  1:1674015/135‑1 (MQ=255)
                                                                                                                              ctcggggggggtggTCTCCCAGTTTTTTTC‑CCCCCGAAAAAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATGTTGT‑GTGGCAATCGCTGCT‑GGTAACAATGTGGAAAGGGAAGTGCGGCGa                                <  2:336928/135‑1 (MQ=255)
                                                                                                                                  gggggggcTGTCTCCCAGTTTTTTTC‑CCGCCGAAAAAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATGTTGT‑GTGGCAATCGCTGCT‑GGTAACAATGTGGAAAGGGAAGTGCGGCGAAt                              <  1:1360944/141‑1 (MQ=255)
                                                                                                                                    ggcggcTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGAGAAAAAAA‑TTTAAACCCCCCCTGTTGT‑GTGGGGAAAGCCGCT‑TGGGACAAAGTGGGAAGGGGAGTGGGGGGGATCAACTcc                       >  1:925086/1‑151 (MQ=255)
                                                                                                                                    ggcggcTGTCCCCCATCTTTTTTC‑CCGCCGAAAAAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATGTTGT‑GTGGCAATCGCTGCT‑GGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTcc                       <  1:174616/151‑1 (MQ=255)
                                                                                                                                           gTCTGCCAGCTATTTTT‑CCGCCGCAACAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGGAAAAGAAAAAACATTTAAGACC‑CCATGTTGT‑GGGGGAACAGCTGCC‑GGGAACAAATTGGGAAGGGGAGGGGGGGGGAACAGCCCCCGaaac                 >  1:813919/1‑148 (MQ=255)
                                                                                                                                            tCTGCCAGCTATTTTT‑CCGCCGCAACAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGGGAAAAAAA‑TTTTAAACCCCCCTTTTTT‑GTGGGGAACGCTGCTGGGTAAAAAGGGGGAAGGGGAGGGGGGGGGGA‑CAGCTCCAGaa                   >  1:1226859/1‑147 (MQ=255)
                                                                                                                                               gCCAGCTATTTTT‑CCGCCGCAACAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCGAGAGAAAAAAC‑TTTTAAACCCTCAATTTGT‑GTGGCAAACCCCGCC‑GGGTACAATTTGGGGAAGGGAGGGGGGGGGATCAGCACCCGAAAcgcccc            >  1:1765955/1‑148 (MQ=255)
                                                                                                                                                ccattttttttc‑CCGCCACAAAAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATGTTGT‑GTGGCAATCGCTGCT‑GGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTa           <  2:1504798/144‑1 (MQ=255)
                                                                                                                                                ccattttttttc‑CCCCCACAAAAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATGTTGT‑GTGGCAATCGCTGCT‑GGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTa           <  2:1761331/144‑1 (MQ=255)
                                                                                                                                                ccAGCTTTTTTC‑CCCCCGCAAAAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATGTTGT‑GTGGCAATCGCTGCT‑GGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACcgc               <  1:243709/147‑1 (MQ=255)
                                                                                                                                                ccAGCTATTTTC‑CCCCCCAAAAAAAC‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATGTTGT‑GTGGCAATCGCTGCT‑GGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAg                     <  1:1011785/141‑1 (MQ=255)
                                                                                                                                                 ccgCTTTTTTC‑CCGCCGCAAAAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATGTTGT‑GTGGCAATCGCTGCT‑GGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTaa          <  1:481697/149‑1 (MQ=255)
                                                                                                                                                     ttttttc‑CCGCCGCAACAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATGTTGT‑GTGGCAATCGCTGCT‑GGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAAt         <  1:25767/146‑1 (MQ=255)
                                                                                                                                                     ttttttc‑CCGCCGCAAAAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATGTTGT‑GTGGCAATCGCTGCT‑GGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTg       <  1:504264/148‑1 (MQ=255)
                                                                                                                                                     ttttttc‑CCGCCACAAAAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATGTTGT‑GTGGCAATCGCTGCT‑GGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAAtt        <  1:26326/147‑1 (MQ=255)
                                                                                                                                                      tttttt‑CCCCCCCAAAAAAAACCCGCCGTTCCCCCAGGTATGGGCGGGGAAAGAGAAATAGT‑TTTAAGCCC‑TTAAGTTGTGGGGGCATCCCCTGCT‑GGT‑ACAATGTGGGAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTg       <  2:211673/149‑1 (MQ=255)
                                                                                                                                                       ttttc‑CCGCCGCAAAAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATGTTGT‑GTGGCAATCGCTGCT‑GGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAg     <  1:830430/150‑1 (MQ=255)
                                                                                                                                                        tttc‑CCGCCGAAAAAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATGTTGT‑GTGGCAATCGCTGCT‑GGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGt    <  1:775460/150‑1 (MQ=255)
                                                                                                                                                        tttc‑CCCCCGCAACAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATGTTGT‑GTGGCAATCGCTGCT‑GGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGt    <  2:1350731/150‑1 (MQ=255)
                                                                                                                                                         ttt‑CCGCCGCAACAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGAGAAAAAAC‑TTTTAGACCGCCCTGTTTT‑GTGGGGAACGCTGCT‑GGGAAAAAAATGGGAAGGGGAGGGCGGGGGATCAACTCCCCAACCGCGCCAATTTAGttt  >  1:693299/1‑150 (MQ=255)
                                                                                                                                                           |                                                                                                                                                          
TCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAA‑CAACCGCGCATTTA‑TCCAACGCCAT‑GTGCAAAAA‑T‑G‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATGTTGT‑GTGGCAATCGCTGCT‑GGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTC  >  NC_000913/357732‑358028

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