Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913357,8800TC41.6% 112.2 / 22.2 78K271E (AAA→GAA) cynRtranscriptional activator of cyn operon; autorepressor
Reads supporting (aligned to +/- strand):  ref base T (34/11);  new base C (0/32);  total (34/44)
Fisher's exact test for biased strand distribution p-value = 2.00e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGT  >  NC_000913/357733‑358026
                                                                                                                                                      |                                                                                                                                                    
cGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTccg                                                                                                                                                   <  1:1955930/151‑1 (MQ=255)
cGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTccg                                                                                                                                                   <  1:287073/151‑1 (MQ=255)
cGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTcc                                                                                                                                                    <  1:987141/150‑1 (MQ=255)
cGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTcc                                                                                                                                                    <  2:627353/150‑1 (MQ=255)
  tAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTccgcc                                                                                                                                                 >  1:361430/1‑151 (MQ=255)
  tAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTccgc                                                                                                                                                  >  2:1743184/1‑150 (MQ=255)
   aGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTccgccg                                                                                                                                                <  2:878369/151‑1 (MQ=255)
   aGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTccgcc                                                                                                                                                 >  2:1140213/1‑150 (MQ=255)
         gATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGcaaca                                                                                                                                           >  2:2026022/1‑150 (MQ=255)
          aTAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACaaa                                                                                                                                         <  1:1243860/151‑1 (MQ=255)
          aTAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACaaa                                                                                                                                         <  1:1848082/151‑1 (MQ=255)
           tAAGCGTAGCGCTTCAGGCAATTCCAGCCGCGGACCTGTGTCAGCGGCTACCGGGATTCATTTCCGCCAACACCCGCGCATTCATCCAACGCCATGTGCAAAAATGCATTCGCGTCGGCTGTCTGCCAGCTATTTTCCGCCGcaacaa                                                                                                                                          <  2:1237731/148‑1 (MQ=14)
           tAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACaaa                                                                                                                                         >  2:768295/1‑150 (MQ=255)
            aaGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACaaaa                                                                                                                                        <  2:1258916/150‑1 (MQ=255)
            aaGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAc                                                                                                                                       >  1:1792281/1‑151 (MQ=255)
             aGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAcc                                                                                                                                      >  2:352727/1‑151 (MQ=255)
                gTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAccgcc                                                                                                                                   >  1:727587/1‑151 (MQ=255)
                gTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAccgac                                                                                                                                   >  1:1616099/1‑149 (MQ=255)
                   gcgcATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGtt                                                                                                                                >  1:1635819/1‑151 (MQ=255)
                     gcATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTtct                                                                                                                              >  1:668219/1‑151 (MQ=255)
                      cATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTtctc                                                                                                                             >  1:1114480/1‑151 (MQ=255)
                      cATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTtctc                                                                                                                             >  1:1683744/1‑151 (MQ=255)
                         cAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTcc                                                                                                                          >  2:1068227/1‑151 (MQ=255)
                               aTTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCagtagt                                                                                                                    >  1:281082/1‑151 (MQ=255)
                               aTTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCagtag                                                                                                                     >  2:1875817/1‑150 (MQ=255)
                                tccccgCCGCAGCCCTGTGTCGGCGGCTCCCGTGTTTCTTTCCCCCCAACACCCGCGCTTTTACCCAACGCCTTGTGAAAAAATCCCTTCGCGGCGGCTGTCTGCCAGCTATTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCagta                                                                                                                      <  1:1319676/143‑1 (MQ=255)
                                   cAGCCGCAGACCTGTGTCAGAGGCTAACGTGATTCATTTCCGAAAACAACCGCGCATTTATCCAACGCCATGTGCAAAACTGCCTTCGCGGACGCTGTCTGCCAGCTATTTTTCAGCAGCAACAACACCGCCGTTCACTCCAGTAGTGGc                                                                                                                 >  2:2139920/1‑150 (MQ=255)
                                              cTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCggggg                                                                                                           >  2:1700009/1‑143 (MQ=255)
                                               tgtgCCGGCGCCTCCCGTGATTATTTCCCCCCAACACCCGCGCTTTTACCCACCGCCTTGTGAAAAAATGCCTCCGCGGCGGTGGTCGGCCGGTTTTTTTCCCGCCAAAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCaa                                                                                                         <  2:1704582/145‑1 (MQ=255)
                                                    cggcgcCTCCCGTCTTTCATTCCCCCCAACACCCCCGCATTTACCCACCCCGATGGGAAAAAGGGCCTCGGCGCCGGCTCTCTCCCGTTTATTTCCCGGCCGCACAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAg                                                                                               <  1:1490716/147‑1 (MQ=255)
                                                     aGCGGCTACCGTGAATCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGAGGCTGTCTGCCAGCTCTTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCggggggggaagca                                                                                                   >  2:1358711/1‑136 (MQ=25)
                                                      gcgcCTCCCGTGTTTCATTCCCCCCAACACCCGCCCTTTTACCCACCGCCATGTGAAAAAAGGCCTCCGCGGCGGCTGTCTCCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCttt                                                                                           <  1:36197/151‑1 (MQ=255)
                                                      gCGGCTCCCGTGTTTTTTTCCCCCAAACACCCGCGCTTTTACCCAACGCCATGTAAAAAAATCCCTTCGGGGGGGCTGTCTCCCAGTTTTTTTCCCCCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCt                                                                                             <  2:2131483/149‑1 (MQ=255)
                                                        ggCTCCCTTGATTCATTCCCGCCAACACCCGCGTTTTTACCCACCCCATTGTGAAAAAATGCCTCCGCGGCGGCTGTCGGCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCt                                                                                             <  1:420232/147‑1 (MQ=255)
                                                        ggCTCCCGTATTTCTTTCCCGCCAACACCCGCGCTTTTATCCACCCCCATGTGAAAAAATGCCTCCGGGGCGGCTGTCTGCCAGTTTTTTTTCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTa                                                                                          <  1:640800/150‑1 (MQ=255)
                                                         gCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggtgccggggaata                                                                                                 >  2:1180879/1‑131 (MQ=255)
                                                          atACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGTAcgagggaaagatct                                                                                            >  2:1299634/2‑137 (MQ=25)
                                                             cggTGATTCATCCCCCCCAACACCCGCGAATTAATCCACCGCATTTAGAAAAATTGCCTCCGCGGCGGCTGTCCCCCAGCTTTTTTCCCGCGGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGccc                                                                                     <  1:1910652/148‑1 (MQ=255)
                                                             ccGTGTTTCATTCCCGCCAACAACCGCCATTTTACCCAACCCCTTGTAAAAATTTGCCTCCGGGGCGCCTGTTTCCCAGCTTTTTCCCCCCGCAAAAAAGCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGccc                                                                                     <  1:2177338/149‑1 (MQ=255)
                                                             ccGTGTTTCTTTCCCCCAAAAACCCGCCCTTTTACCCACCGCCATGTCAAAAAATCCCTTCGCGGCGGTTTTCTCCAAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAg                                                                                        <  1:972508/147‑1 (MQ=255)
                                                               gtttttCTTTCCCCCCACCACCCGCGCTTTTATCCACCCCCTTGTGAAAAAAGGCCTTGGGGGCGGTTGTCTGCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGt                                                                                   <  1:482270/146‑1 (MQ=255)
                                                                tgttcctttccccccAACACCCGCCCTTTCCCCCACCGCCTTGTGAAAAAGCGCCTGCGCGGCGGTTTTCGCCCAGTTTTTTTCCCGCCACAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGt                                                                                   <  2:1097151/136‑1 (MQ=255)
                                                                 gATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGAAAAAACTTTaatacc                                                                                     >  1:146954/1‑142 (MQ=25)
                                                                       tttCCCCCAACACCCGCCCTTTTCACCACCCCCATGTGAAAAAATCCCTTCGGGGCGGTTTTCTGCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCAtgttg                                                                            <  2:560526/149‑1 (MQ=255)
                                                                       ttccccccAACACCCCCCCTTTTACCCACCGCCATGCGAAAAAATCCCTTGGCGGCGGCTGTCGGCCAGTTATTTTGCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGt                                                                                   <  2:1653482/136‑1 (MQ=255)
                                                                        ttCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAAAAACTTTTAACCCCCCCTGTtttgt                                                                         >  1:267943/1‑151 (MQ=25)
                                                                          cccccaccaCCCGCGCATTTATCCACCGCCTTGTGAAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTtgtgt                                                                         <  1:464253/146‑1 (MQ=255)
                                                                               aacaccCGCCCTTTTACCCACCGCCTTGTGAAAAAAGCCCTTCGCGGCGGTTGTCTGCCAGCTTTTTTCCCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCa                                                                     <  1:1032545/148‑1 (MQ=255)
                                                                                  aacCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAAAAACTTTAAACCCGTCATGTTTTGTGGCaa                                                                    >  1:382756/1‑146 (MQ=255)
                                                                                      cgcgCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGCAAGGGAAAAGACTTTAAACCCTCATGGTGTGTggggttag                                                                 >  2:333337/1‑140 (MQ=25)
                                                                                       gccaaTTTACCCAACCCCTTGTCAAAAATCCCCTTCGGGGCGGTTTTCCGCCGTCTTTTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgctgct                                                            <  1:1508304/146‑1 (MQ=255)
                                                                                            tttATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGCAAAGAAAAAAATTTTAAACCCCCCTGTTTTGTGGGGAACGCTGCTgg                                                          >  1:787669/1‑147 (MQ=21)
                                                                                            tttATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCGAGAGAAAAAAATTTTAAACCCCCATGTTTTGTGGGGAAAGCTGCTgg                                                          >  1:465384/1‑147 (MQ=25)
                                                                                            tttATCCAACCCCATGTAAAAAAATCCCTTCGCGGCGGCTGTCTCCCAGTTATTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTgg                                                          <  1:342028/147‑1 (MQ=255)
                                                                                            tttACCCACCGCCATGTCAAAAAATCCCTCCGGGGCGGTTGTCTGCCAGTTATTTTCCCGCCCCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGctgctg                                                           <  1:904170/146‑1 (MQ=255)
                                                                                            tttACCCACCCCCATGTGAAAAAATGCTTCCGCGGCGGCTTTCTCCCATTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTa                                                        <  1:1969844/149‑1 (MQ=255)
                                                                                            tttACCCAACGCCATGTGAAAAATTCCCTCGGGGGCGGTTGTCTCCCAGCTTTTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgctgct                                                            <  1:962394/145‑1 (MQ=255)
                                                                                                 ccAACCCCATGTGCAAAAATCCCTCCGCGGCGGTTGTCTGCCAGCTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAAtgt                                                 <  1:1429417/151‑1 (MQ=255)
                                                                                                  caccccaagggcaaaaaaTCCCTCCGGGGGGGCTTTCTCCCAGGTTTTTTCCCCCACCAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACa                                                     <  2:721553/133‑1 (MQ=255)
                                                                                                  cAACCCCATGTGAAAAAAGTCCTTCCGGGCGTGTTGCCCCCATGTTTTTTCCCCCCACAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACa                                                     <  2:2177442/146‑1 (MQ=255)
                                                                                                    aCGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGAAAAAAATTTAAACCCCGCCTGTTTTGTGGGGAACGCTGCTGGGAaaaaaaggg                                               >  1:887493/1‑145 (MQ=25)
                                                                                                     cccctggtaaaaaaaTGCCTCCGGGGGGGTTGCCTGCAAGTTATTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAAt                                                   <  2:468096/135‑1 (MQ=255)
                                                                                                     cGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAAAGAAAAAAAATTATAACCCCCCTCTTGTGTGGGGAACGCTGCTGGTGACaaaatgg                                               >  1:1242600/1‑144 (MQ=21)
                                                                                                       ccttgtaaaaaatTGCCTGCGGGCGGGCTGCCGGCAAGTTTTTTTCCCCCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAAtgtg                                                <  2:626541/138‑1 (MQ=255)
                                                                                                        cATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGAAAAAAATTTTATCCCCTCCTGTTGTGTGGGGAACGCTGCTGGTGACAAtgtg                                                >  1:1613105/1‑145 (MQ=37)
                                                                                                             gaaaaaaaGCCTTCCGGGGCGGCTGTCTCCCGTTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCGTGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAg                                      <  2:2076678/149‑1 (MQ=255)
                                                                                                                aaaaTGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGAAAAAACCTTTAAGCCCCCCTGTTGTTTTGGGATCCCTGCTGGTGAAAAATTGGGAAGGGGAggggg                                  >  1:810743/1‑147 (MQ=255)
                                                                                                                     tcccTTCGCGGGGGTTGTCTCCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTgcggcg                               <  1:1498079/148‑1 (MQ=255)
                                                                                                                       ccTCCGCGGCGCCTTTCGCCAAGTTTTTTTCCCCCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAAt                            <  2:434068/151‑1 (MQ=255)
                                                                                                                        cTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGAAAAAAGTTTTAACCCCCCCAGTGTTGTGGGGATAGCTGCTTGGGAAAAAGTGGGGAGGGAAGAgcgg                                 >  1:2049355/1‑145 (MQ=255)
                                                                                                                            cgcgGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGAAAAAAATTTTAAACCCCCCTTTTTTGTGGGGATACCTGCTGGGTACAAAGTGGGAAggggaaggggggggatc                           >  1:307131/1‑134 (MQ=255)
                                                                                                                               ggcggcTGTCCCCCGGTTTTTTTCCCCCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAAt                            <  2:1993011/144‑1 (MQ=255)
                                                                                                                                gcggcTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACTGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCa                    <  2:89198/151‑1 (MQ=255)
                                                                                                                                     tttcgccgagtttttttcCCGCCAAAACAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCa                    <  2:1370889/133‑1 (MQ=255)
                                                                                                                                        ctcccAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGTGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGaa                 <  2:506444/143‑1 (MQ=255)
                                                                                                                                           ccAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACcgc             <  1:1550277/147‑1 (MQ=255)
                                                                                                                                           ccAGTTTTTTTCCCCCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCt          <  1:932212/150‑1 (MQ=255)
                                                                                                                                           ccAGTTATTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCa                    <  1:1800706/140‑1 (MQ=255)
                                                                                                                                           ccAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGCAGGGAAAAAAACTTAAACCCCCCCCGTTCTGTGGCGATCGCTCCGGGGGACCAAGTGGGAAGGGGAGGGCGGCGGAACAGCTccccaaacccgc           >  2:814227/1‑139 (MQ=255)
                                                                                                                                                  tttttCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGAAAAAACTTTTAACCCCCCCTTTTTTGTGGGGATCGCTGCTGGGAAAAAAGTGGAAAAGGAAGGGGGGGGGAACAGCCCCCCAACCGCGCCAATTTAg   >  1:135170/1‑150 (MQ=255)
                                                                                                                                                  ttttcCCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAcc               <  2:1308814/138‑1 (MQ=255)
                                                                                                                                                  ttttcCCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGa    <  1:395976/149‑1 (MQ=255)
                                                                                                                                                   tttcCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGt  <  1:1516533/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                    
CGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGT  >  NC_000913/357733‑358026

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