Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913357,8900CA55.7% 30.1 / 112.4 88L267F (TTG→TTT) cynRtranscriptional activator of cyn operon; autorepressor
Reads supporting (aligned to +/- strand):  ref base C (29/10);  new base A (0/49);  total (29/59)
Fisher's exact test for biased strand distribution p-value = 4.20e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.04e-11

TGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAA‑CAACCGCGCATTTA‑TCCAACGCCAT‑GTGCAAAAA‑T‑G‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATGTTGT‑GTGGCAATCGCTGCT‑GGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCAATGACC  >  NC_000913/357741‑358040
                                                                                                                                                             |                                                                                                                                                           
tGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAA‑CAACCGCGCATTTA‑TCCAACGCCAT‑GTGCAAAAA‑T‑G‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGcaac                                                                                                                                                             >  2:1673857/1‑150 (MQ=255)
 gATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAA‑CAACCGCGCATTTA‑TCCAACGCCAT‑GTGCAAAAA‑T‑G‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGcaaca                                                                                                                                                            >  2:1141075/1‑150 (MQ=255)
 gATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAA‑CAACCGCGCATTTA‑TCCAACGCCAT‑GTGCAAAAA‑T‑G‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGcaaca                                                                                                                                                            <  1:485522/150‑1 (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)
                                                                                                                                           tttttttcccccccgaAAAAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATGTTGT‑GTGGCAATCGCTGCT‑GGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAAt                     <  2:300166/138‑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)
                                                                                                                                                tttccccccaaaacaaAA‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATGTTGT‑GTGGCAATCGCTGCT‑GGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGt                <  2:1765756/139‑1 (MQ=255)
                                                                                                                                                ttt‑CCGCCGCAACAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGAGAAAAAAC‑TTTTAGACCGCCCTGTTTT‑GTGGGGAACGCTGCT‑GGGAAAAAAATGGGAAGGGGAGGGCGGGGGATCAACTCCCCAACCGCGCCAATTTAGttt              >  1:693299/1‑150 (MQ=255)
                                                                                                                                                   ccccccgcAACAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATGTTGT‑GTGGCAATCGCTGCT‑GGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCg             <  2:501033/146‑1 (MQ=255)
                                                                                                                                                    ccgccgCAAAAAAC‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATGTTGT‑GTGGCAATCGCTGCT‑GGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCg             <  2:355416/149‑1 (MQ=255)
                                                                                                                                                             caaAA‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGAGAAAAAAA‑TTTTAACCCCCCCTTTTTT‑GTGGGCATCGCTGCC‑GGGTAAAAAGTGGGAAGGGGAGGGGGGGGGATCAACTCCCCAAACCCGCCAATTTAGTTTGCCCCCATGacc  >  1:1612348/1‑151 (MQ=255)
                                                                                                                                                             |                                                                                                                                                           
TGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAA‑CAACCGCGCATTTA‑TCCAACGCCAT‑GTGCAAAAA‑T‑G‑CCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATGTTGT‑GTGGCAATCGCTGCT‑GGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCAATGACC  >  NC_000913/357741‑358040

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