Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913357,8750AT47.1% 106.0 / 53.7 87N272K (AAT→AAA) cynRtranscriptional activator of cyn operon; autorepressor
Reads supporting (aligned to +/- strand):  ref base A (37/9);  new base T (0/41);  total (37/50)
Fisher's exact test for biased strand distribution p-value = 2.26e-16
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAA‑CAACCGCGCATTTA‑TCCAACGCCAT‑GTGCAAAAA‑T‑G‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAA  >  NC_000913/357726‑358020
                                                                                                                                                            |                                                                                                                                                    
cTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAA‑CAACCGCGCATTTA‑TCCAACGCCAT‑GTGCAAAAA‑T‑G‑CCTT‑CGCGGCGGCTGTCTGCCAGCTAt                                                                                                                                                     >  1:1242707/1‑151 (MQ=255)
  gAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAA‑CAACCGCGCATTTA‑TCCAACGCCAT‑GTGCAAAAA‑T‑G‑CCTT‑CGCGGCGGCTGTCTGCCAGCTAttt                                                                                                                                                   >  1:780739/1‑151 (MQ=255)
  gAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAA‑CAACCGCGCATTTA‑TCCAACGCCAT‑GTGCAAAAA‑T‑G‑CCTT‑CGCGGCGGCTGTCTGCCAGCTAt                                                                                                                                                     <  1:1319754/149‑1 (MQ=255)
      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‑CCGCCGCAACAAAAccg                                                                                                                               <  2:686945/151‑1 (MQ=255)
                      cGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAA‑CAACCGCGCATTTA‑TCCAACGCCAT‑GTGCAAAAA‑T‑G‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAAccgc                                                                                                                              <  2:271784/151‑1 (MQ=255)
                            gcATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAA‑CAACCGCGCATTTA‑TCCAACGCCAT‑GTGCAAAAA‑T‑G‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTtct                                                                                                                        >  1:300182/1‑151 (MQ=255)
                             cATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAA‑CAACCGCGCATTTA‑TCCAACGCCAT‑GTGCAAAAA‑T‑G‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTtctc                                                                                                                       >  2:1175717/1‑151 (MQ=255)
                              aTCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAA‑CAACCGCGCATTTA‑TCCAACGCCAT‑GTGCAAAAA‑T‑G‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTtctct                                                                                                                      >  1:1761495/1‑151 (MQ=255)
                                cAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAA‑CAACCGCGCATTTA‑TCCAACGCCAT‑GTGCAAAAA‑T‑G‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTcc                                                                                                                    >  1:501054/1‑151 (MQ=255)
                                    cAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAA‑CAACCGCGCATTTA‑TCCAACGCCAT‑GTGCAAAAA‑T‑G‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCagt                                                                                                                 >  1:221927/1‑150 (MQ=255)
                                    cAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAA‑CAACCGCGCATTTA‑TCCAACGCCAT‑GTGCAAAAA‑T‑G‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCagt                                                                                                                 >  2:364711/1‑150 (MQ=255)
                                    cAATTCCAGCCGCAGACCTGTGTCAGCGGCTACAGTGATTCATTTCCGCCAA‑CAACCGCGCATTTA‑TCCAACGCCAT‑GTGCAAAAA‑T‑G‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAACCCGCCGAtctct                                                                                                                      >  2:886532/1‑145 (MQ=255)
                                     aaTTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAA‑CAACCGCGCATTTA‑TCCAACGCCAT‑GTGCAAAAA‑T‑G‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCagt                                                                                                                 >  2:282525/1‑149 (MQ=255)
                                         ccgcccGCGGCCCTGTGTCAGCGGCTCCCGTGTTTCATTCCCGCCAC‑CACCCGCGCATTTA‑CCCACCGCCTT‑GTGAAAAAA‑T‑G‑CCTC‑CGCGGCGGCTTTCTGCCAGTTTTTTTC‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTgg                                                                                                            <  1:1360504/146‑1 (MQ=255)
                                         ccgcccGCAGCCCTGTGTCAGCGCCTCCCGTATTTCATTCCCCCCAC‑CACCCGCGCTTTTA‑CCCACCCCCAT‑GTGAAAAAA‑T‑C‑CTTC‑CGCGGGGGCTGTCTGCCAGTTTTTTTC‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTgg                                                                                                            <  1:621951/146‑1 (MQ=255)
                                         ccAACCGCAACCCTGTCTCGCCGCCTACGGTGATTCATTCCCTCCAA‑CAGCCGCGCTTTTC‑CCCAACCCTAT‑GTGCAAAAA‑T‑C‑CCTC‑CGGGGGGGCTTCCTCCCAGTTTTTTTC‑CCCCCCCAAAAAAGCCGCCGTTCTCTCCAGTAGTGGc                                                                                                           <  1:1208048/151‑1 (MQ=255)
                                             ccGCAGCCCTGTGTCAGGGGCTCCCGTGATTCTTTCCCCCCAA‑CACCCGCTCTTTTA‑CCCACCGCCTT‑GTGAAAAAC‑T‑G‑CCTC‑CGCGGCGGCTGTCTGCCAGTTTTTTTC‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCggg                                                                                                        <  1:1184444/150‑1 (MQ=255)
                                             ccGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAA‑CAACCGCGCATTTA‑TCCAACGCCAT‑GTGCAAAAA‑T‑G‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCggg                                                                                                        >  1:1618270/1‑150 (MQ=255)
                                               gCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAA‑CAACCGCGCATTTA‑TCCAACGCCAT‑GTGCAAAAA‑T‑G‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggg                                                                                                       >  1:355433/1‑149 (MQ=255)
                                                cAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAA‑CAACCGCGCATTTA‑TCCAACGCCAT‑GTGCAAAAA‑T‑G‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggg                                                                                                       >  1:1358088/1‑148 (MQ=255)
                                                 accccTGTGTCAGCGTCTCCCGTTTTTCTTTCCCCCCAA‑CACCCGCGCTTTTC‑CCCACCCCCAT‑GTAAAAAAA‑T‑G‑CTTT‑CGGGGCGGTTGTCTGCCAGTTTTTTTC‑CCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAag                                                                                                   <  1:391625/148‑1 (MQ=255)
                                                  gACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAA‑CAACCGCGCATTTA‑TCCAACGCCAT‑GTGCAAAAA‑T‑G‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCggggg                                                                                                      >  2:1459979/1‑146 (MQ=255)
                                                   aCCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAA‑CAACCGCGCATTTA‑TCCAACGCCAT‑GTGCAAAAA‑T‑G‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggg                                                                                                       >  1:690899/1‑145 (MQ=255)
                                                    ccTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAA‑CAACCGCGCATTTA‑TCCAACGCCAT‑GTGCAAAAA‑T‑G‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggggggg                                                                                                   >  2:243700/1‑144 (MQ=255)
                                                        ttccAGCGCCTCCCGTATTTCATTCCCCCAAA‑CACCCGCGCATTTT‑CCCAACGCCAT‑TTGCAAAAT‑T‑G‑CCTT‑CGCGGCGGCTGTCTCCCAGCTTTTTTC‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAAt                                                                                             <  1:1132613/147‑1 (MQ=255)
                                                        tgtCAGCGGCTACCGTGATTCATTTCCGCCAA‑CAACCGCGCATTTA‑TCCAACGCCAT‑GTGCAAAAA‑T‑G‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAGGAGaaaa                                                                                             >  1:258826/1‑149 (MQ=255)
                                                         gccggggcctcccGTTATTCTTTCCCCCAAA‑CACCCGCCCTTTTA‑CCCACCCCCAG‑GTGAAAAAA‑G‑G‑CCTC‑CGGGGCGGCTGTCGCCCAGTTTTTTCC‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAgaga                                                                                                <  2:87153/135‑1 (MQ=255)
                                                           cggcgcCTCCCGTGTTTCTTTCCCCCAAA‑CACCCCCGCTTTTA‑CCCACCGCCGT‑GTAAAAAAA‑T‑C‑CCTC‑CGCGGCGGCTTTCCGCCATTTTTTTTC‑CCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCaa                                                                                                    <  2:1342987/138‑1 (MQ=255)
                                                           cAGCGGCTACCGTGATTCATTTCCGCCAA‑CAACCGCGCATTTA‑TCCAACGCCAT‑GTGCAAAAA‑T‑G‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGgaaagaaaaa                                                                                            >  1:150684/1‑146 (MQ=255)
                                                               ggCTACCGTGATTCTTTCCCCCCAA‑CACCCCCGCTTTTA‑CCCACCCCCAT‑G‑GCAAAAATT‑C‑CCTC‑CGCGGCGGTTGTCGCCCAGCTTTTTTC‑CCCCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGC‑TTTaa                                                                                    <  1:1674945/151‑1 (MQ=255)
                                                               gcttcccGTGTTTCTTTCCCCCAAA‑CACCCGCGCTTTTA‑CCCACCGCCAT‑GTCAAAAAT‑T‑G‑CTTT‑CGGGGGGGTTGTCTGCCAGTTTTTTTC‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGC‑TTTa                                                                                     <  1:121229/145‑1 (MQ=255)
                                                               gcctgccGGGATTTTTTTCCCCCAA‑CACCCGCG‑ATTTACTCCACCGCTAT‑GTGAAAAAG‑T‑C‑CCTT‑CGCGGCGTCTGT‑TGGCAGGTTTTTTTCCCGCCACAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGc                                                                                          <  1:1726483/141‑1 (MQ=255)
                                                               gcctcccGTGTTTCATTCCCCCAAC‑CCCCCGCCCTTTTA‑CCCAACCCCTT‑GTAAAAAAA‑T‑C‑CCTG‑CGGGGCGGCTGTCGGCAACTTTTTTTC‑CCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGaaa                                                                                              <  2:1551206/137‑1 (MQ=255)
                                                                ggcaccgttaattctttCCCCCAACCACCCGCCCTTTTC‑CCCACCCCCAT‑GTGAAAAAA‑G‑C‑CCTC‑CGGGGCGGCT‑TCCGGCCGCTTTTTTTCCCCCCACAAAAAAGCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGC‑ttt                                                                                      <  1:145877/135‑1 (MQ=255)
                                                                    ccgtttttCTTTCCCCCCAA‑CACCCGCGCTTTTA‑TCCACCCCATT‑GTGCAAAAA‑T‑G‑CCTT‑CGCGGCGGTTTTCTGCCAGTTTTTTTC‑CCGCCCAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGC‑TTTAAGCCCg                                                                               <  1:396902/145‑1 (MQ=255)
                                                                    ccgtttttCATTCCCCCAAA‑CACCCGC‑CATTTC‑CCCCACCCCATGGTAAAAAAT‑T‑G‑CCTG‑CGGGGCGGTTTTCTCCCAGTTTTTTTC‑CCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGC‑TTTAAGcc                                                                                 <  2:1490047/143‑1 (MQ=255)
                                                                    ccGTGTTTCTTTCCCCCCAA‑CACCCGCGCTTTTA‑CCCACCCCCAT‑GTGAAAAAA‑T‑G‑CCTG‑CGGGGGGGTTGTCTGCCAGTTTTTTTC‑CCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGC‑TTTAAGCCCg                                                                               <  1:895558/151‑1 (MQ=255)
                                                                    ccGTGATTCATTTCCGCCAA‑CAACCGCGCATTTA‑TCCAACGCCAT‑GTGCAAAAA‑T‑G‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAGGCGaaaaaaactttaaaacc                                                                                >  1:1112961/1‑137 (MQ=25)
                                                                    ccGTGATTCATTCCCCCCAC‑CACCCGCGCTTTTA‑CCCACCGCCAT‑GTGCAAAAA‑T‑G‑CTTT‑CGGGGCGGCTGTCTCCCAGCTTTTTTC‑CCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGC‑TTTAAGCCCg                                                                               <  1:1570986/151‑1 (MQ=255)
                                                                      gTGATTCTTTTCCCCAAA‑CACCCGCGCTTTTA‑CCCAACGCCAT‑TTGAAAAAA‑T‑C‑CCTT‑CGCGGCGGTTTTCTGCCAGTTTTTTTC‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGC‑TTTAAGCCCg                                                                               <  1:1500348/149‑1 (MQ=255)
                                                                      gTGATTCATTTCCGCCAA‑CAACCGCGCATTTA‑TCCAACGCCAT‑GTGCAAAAA‑T‑G‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGGGAGaaaaaaattttaaccc                                                                                 >  1:833361/1‑135 (MQ=25)
                                                                         tttCATTCCCCCAAA‑CACCCGCTCTTTTA‑CCCACCCCCAT‑GT‑AAAAAA‑T‑T‑CCTTGCGCGGGGGCTGTCCGCCAGTTTTTTTC‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGC‑TTTAAGCCCGTCa                                                                            <  2:1278433/148‑1 (MQ=255)
                                                                          ttCTTTCCCCCCAA‑CACCCGCCCTTTTA‑CCCACCGCCAT‑GTGAAAAAA‑T‑C‑CCTC‑CGGGGCGGCTGTCTGCAAGTTTTTTTC‑CCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGC‑TTTAAGCCCGTCa                                                                            <  1:1754055/148‑1 (MQ=255)
                                                                          ttCATTCCCCCAAA‑CACCCGCGCTTTTA‑CCCACCGCCTT‑GTAAAAAAA‑T‑C‑CCTG‑CGGGGCGGCTTTCTCCCAGTTTTTTTC‑CCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGC‑TTTAAGCCCGTCAt                                                                           <  1:814042/149‑1 (MQ=255)
                                                                              tttCCGCCAA‑CAACCGCGCATTTA‑TCCAACGCCAT‑GTGCAAAAA‑T‑G‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCGCGAGGAAAAGG‑ttttaatcccccccggt                                                                        >  2:92223/1‑134 (MQ=25)
                                                                              tccccGCCAA‑CACCCTCCAATTTA‑GCCAGCGCCAT‑GTCAAAAAC‑T‑G‑CCTC‑CGCGGTGGTTGTCTCCACGTTATTTCC‑CCGCGCAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGC‑TTTAAGCCCGTCAtgttgt                                                                      <  1:1213426/147‑1 (MQ=255)
                                                                                ccccccAA‑CACCCGCGCTTTTA‑CCCAACGCCTT‑TTGAAAAAA‑G‑C‑CCTC‑CGCGGCGGCGTTCTCCCTGTTTTTTTC‑CCGCCACAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATGTtgtgtg                                                                   <  1:938787/147‑1 (MQ=255)
                                                                                 ccccaaa‑caccCGCGTTTTTA‑CCCACCGCCAG‑GTCAAAAAA‑T‑C‑CCTC‑CGGGGCGGCTGTCTGCCAGTTTTTTTC‑CCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATGTTGTGTgg                                                                  <  2:609523/146‑1 (MQ=255)
                                                                                 ccccaaa‑caccCGCGCTTTTA‑CCCACCGCCAT‑GTAAAAAAA‑G‑G‑CCTT‑CGCGGGGGTTTTCTGCCGGTTTTTTTC‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATGTTGTGTgg                                                                  <  1:197563/146‑1 (MQ=255)
                                                                                            ccGCGCATTTA‑TCCAACGCCAT‑GTGCAAAAA‑T‑G‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCCGAAGAAAAAACATTTAAACCCG‑CAtgttgttgtg                                                                  >  2:1530859/1‑137 (MQ=25)
                                                                                             cgcataTTTA‑CCCACCGCC‑T‑TTGAAAAAA‑T‑GCCCTT‑GGGGGGGGTTGTCGCCCAGTTTTTTTC‑CCGCCTAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCgctgct                                                       <  2:840735/146‑1 (MQ=255)
                                                                                               cgcATTTA‑TCCAACGCCAT‑GTGCAAAAA‑T‑G‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAAAGAAAAAAA‑TTTAAAACCGCCCTGTTTTGTGGGGGTAGctgct‑g                                                     >  2:1356024/1‑150 (MQ=25)
                                                                                                   tttA‑CCCACCGCCAT‑GTAAAAAAT‑T‑G‑CTTC‑CGGGGGGGCTTTCCCCAAGTTTTTTTC‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTa                                                  <  1:1385519/149‑1 (MQ=255)
                                                                                                          caacccccttgccaaaaag‑T‑C‑CCTC‑CGGGGCGCTTGTCCCCCAGTTTTTTTC‑CCCCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAAc                                                <  2:833322/132‑1 (MQ=255)
                                                                                                           accGCCAT‑GGGAAAAAA‑G‑G‑CCTT‑CGGGGCGGCTGTCTGCCAGTTTTTTTC‑CCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAAtgt                                           <  1:849974/147‑1 (MQ=255)
                                                                                                               ccAT‑GTGAAAAAG‑T‑C‑CCTC‑CGCGGCGGCTGTCTGCCAGTTTTTTTC‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAATGTGGaa                                       <  1:1175796/149‑1 (MQ=255)
                                                                                                               caaT‑GTGAAAAAA‑C‑G‑CTTC‑CGGGGGGGCTTTCCCCCAGTTTTTTTC‑CCCCCGAAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAATGTGGaaa                                      <  2:1244408/148‑1 (MQ=255)
                                                                                                                cAT‑GTGCAAAAA‑T‑G‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGCAGAGAAAAAAA‑TTTAAAACCCCCCTGTTGTGTGGGGATCGCTGCT‑GGGTACAATGTgg                                         >  1:155441/1‑146 (MQ=37)
                                                                                                                 aT‑GTGCAAAAA‑T‑G‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGCAGAGAAAAAAC‑TTTTAAACCCTCCTGTTGTGTGGGGAACGCTGCT‑GGGGACAAAGTGGGAAgg                                    >  1:151766/1‑150 (MQ=255)
                                                                                                                  t‑gtgCAAAAA‑T‑G‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAAAAAAAAA‑ATTTAAACCGCCCTGCTGTTTGGGAAACGCCGCT‑GGGAACAAtgatg                                         >  2:1421831/1‑141 (MQ=255)
                                                                                                                  t‑gtgAAAAAA‑T‑G‑CCTC‑CGGGGCGGTTTCCTGCCAGTTATTTTC‑CCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAATGTGGAAAggg                                   <  2:1037665/150‑1 (MQ=255)
                                                                                                                     taaaaaaa‑T‑C‑CCTC‑CGGGGGGGCTGTCGCCAAGTTTTTTTC‑CCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAATGTGGaaa                                      <  1:774948/142‑1 (MQ=255)
                                                                                                                     gtaaaaaa‑TCG‑CCTT‑CGGGGCGGCTGTCTCCCATTTTTTTTC‑CCCCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAATGTGGAAAGGGaa                                 <  2:1626142/148‑1 (MQ=255)
                                                                                                                       aaaaat‑T‑G‑CTTC‑CGGGGGGGCTGTCTCCCAGTTTTTTTC‑CCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAATGTGGAAAGGGAAGTg                              <  1:1460675/150‑1 (MQ=255)
                                                                                                                        aaaaa‑T‑G‑CCTC‑CGCGGCGGTTTTCTCCCAGTTTTTTTC‑CCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAATGTGGAAAGGGaa                                 <  2:524827/147‑1 (MQ=255)
                                                                                                                        aaaaa‑G‑C‑CCTG‑CGCGGCGGCTGTCGCCCAGTTTTTTTC‑CCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAATGTGGAAAGGGa                                  <  2:1367030/146‑1 (MQ=255)
                                                                                                                            aaaagtccttccggggcggcTGTCGCCCAGTTTTTTTC‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGGAATGGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAATGTGGAAAGGGAAGt                               <  1:1674015/135‑1 (MQ=255)
                                                                                                                                    ctcggggggggtggTCTCCCAGTTTTTTTC‑CCCCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAATGTGGAAAGGGAAGTGCGGCGa                        <  2:336928/135‑1 (MQ=255)
                                                                                                                                        gggggggcTGTCTCCCAGTTTTTTTC‑CCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAATGTGGAAAGGGAAGTGCGGCGAAt                      <  1:1360944/141‑1 (MQ=255)
                                                                                                                                          ggcggcTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGAGAAAAAAA‑TTTAAACCCCCCCTGTTGTGTGGGGAAAGCCGCT‑TGGGACAAAGTGGGAAGGGGAGTGGGGGGGATCAACTcc               >  1:925086/1‑151 (MQ=255)
                                                                                                                                          ggcggcTGTCCCCCATCTTTTTTC‑CCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTcc               <  1:174616/151‑1 (MQ=255)
                                                                                                                                                 gTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAAAGAAAAAACATTTAAGACC‑CCATGTTGTGGGGGAACAGCTGCC‑GGGAACAAATTGGGAAGGGGAGGGGGGGGGAACAGCCCCCGaaac         >  1:813919/1‑148 (MQ=255)
                                                                                                                                                  tCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGGGAAAAAAA‑TTTTAAACCCCCCTTTTTTGTGGGGAACGCTGCTGGGTAAAAAGGGGGAAGGGGAGGGGGGGGGGA‑CAGCTCCAGaa           >  1:1226859/1‑147 (MQ=255)
                                                                                                                                                     gCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCGAGAGAAAAAAC‑TTTTAAACCCTCAATTTGTGTGGCAAACCCCGCC‑GGGTACAATTTGGGGAAGGGAGGGGGGGGGATCAGCACCCGAAAcgcccc    >  1:1765955/1‑148 (MQ=255)
                                                                                                                                                      ccAGCTTTTTTC‑CCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACcgc       <  1:243709/147‑1 (MQ=255)
                                                                                                                                                      ccAGCTATTTTC‑CCCCCCAAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAg             <  1:1011785/141‑1 (MQ=255)
                                                                                                                                                       ccgCTTTTTTC‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTaa  <  1:481697/149‑1 (MQ=255)
                                                                                                                                                            |                                                                                                                                                    
CTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAA‑CAACCGCGCATTTA‑TCCAACGCCAT‑GTGCAAAAA‑T‑G‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑GGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAA  >  NC_000913/357726‑358020

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