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

TCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGA  >  NC_000913/357732‑358024
                                                                                                                                                         |                                                                                                                                                  
tCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTcc                                                                                                                                                  >  2:1083342/1‑151 (MQ=255)
 cGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTcc                                                                                                                                                  <  3:1044853/150‑1 (MQ=255)
   tAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTccgcc                                                                                                                                               >  3:330078/1‑151 (MQ=255)
   tAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTccgcc                                                                                                                                               >  2:1072839/1‑151 (MQ=255)
    aGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTccgccg                                                                                                                                              >  2:264657/1‑151 (MQ=255)
    aGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTccgccg                                                                                                                                              >  2:577806/1‑151 (MQ=255)
       ccTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGca                                                                                                                                            >  2:991986/1‑150 (MQ=255)
         tGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGcaac                                                                                                                                          >  3:490141/1‑150 (MQ=255)
             aaGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAc                                                                                                                                     <  3:1425133/151‑1 (MQ=255)
              aGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAcc                                                                                                                                    <  2:896948/151‑1 (MQ=255)
              aGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAcc                                                                                                                                    <  2:1446198/151‑1 (MQ=255)
               gCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAcc                                                                                                                                    <  2:1232859/150‑1 (MQ=255)
                 gTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAccgcc                                                                                                                                 <  3:24959/151‑1 (MQ=255)
                     ccaataaGGCATTCCCAGCCGCAGCCTTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACACCCGCGCATTTATCCAACGCCTTGTGAAAAAATCCCTCCGGGGCGGCTGTCTCCCAGCTTTTTTTCCGCCGCAACAAAAccgccg                                                                                                                                <  3:976774/142‑1 (MQ=18)
                         tCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTtctctc                                                                                                                         >  2:1422117/1‑151 (MQ=255)
                           aGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTcc                                                                                                                        >  3:303283/1‑150 (MQ=255)
                           aGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTcc                                                                                                                        >  2:1591011/1‑150 (MQ=255)
                             gCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCagt                                                                                                                     >  2:115315/1‑151 (MQ=255)
                             gCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCagt                                                                                                                     >  3:678947/1‑151 (MQ=255)
                               aaTTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCagt                                                                                                                     >  3:1626660/1‑149 (MQ=255)
                                aTTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCagtagt                                                                                                                  >  2:619217/1‑151 (MQ=255)
                                aTTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCagtagt                                                                                                                  >  2:752390/1‑151 (MQ=255)
                                 ttCCAGCCGCAGCCCTGTGTCAGCGGCTCCCGTGATTCATTCCCGCCAACACCCGCGCTTTTACCCAACGCCATGTGAAAAAATGCCTTCGCGGCGGCTGTCTCCCATTTATTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCagtagt                                                                                                                  <  2:394864/150‑1 (MQ=18)
                                 ttCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTg                                                                                                                 >  3:1270600/1‑151 (MQ=255)
                                 ttCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTg                                                                                                                 >  3:536020/1‑151 (MQ=255)
                                      gCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCggg                                                                                                            >  3:636234/1‑151 (MQ=255)
                                       ccGCAGCCTTGTTCCAGCGCCTCCCGTGTTCCTTTCCCCCAAACACCCGCGCTTTTACCCACCGCAGTGTGAAAAAATGCCTTCGGGGGGGCTGTCTGCCGGTTTTTTTCCCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCgggg                                                                                                           <  3:70444/151‑1 (MQ=255)
                                       ccGCAGCCTTGTGTCGGGGCCTCCCGTTATTCTTCCCCCCCAACACCCCCTCTTTTCCCCACCGCCTTGCAAAAAAATGCCTTCGGGGGGTCTGTCTGCCTGCTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCgggg                                                                                                           <  3:165028/151‑1 (MQ=255)
                                       ccGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggg                                                                                                           >  2:94718/1‑151 (MQ=255)
                                        cGCAGCCTTGTGCCGGCGGCTCCGTTATTTCTTTCCCCCAACCACCCGCGCTTTTCCCAACCCCCTTTTAAAAAAATGCTTCCGCGGCGGTTGTCGCCCACTTTTTTTCCCGCCCCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGc                                                                                                          <  2:475221/151‑1 (MQ=255)
                                        cGCAAACCTGTGGCAGCGCCTCCGTTATTTCTTCCCCGCAAACAACCGCCCTTTTTCCCACCGCCTTGTCCAAAAGTGCTTCGGGGGCGGCTGTGTGCCGGTTTTTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGc                                                                                                          <  3:71883/151‑1 (MQ=255)
                                         gCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggg                                                                                                           >  2:784155/1‑149 (MQ=255)
                                         gCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggg                                                                                                           >  2:995855/1‑149 (MQ=255)
                                            gACCTGTGTCGGCGCCTCCCGTGTTTCATTCCCGCCAACACCCGCGCATTTACCCACCGCCTTGTGAAAAAATGCCTTCGCGGGGGTTGTCGCCCAGTTTTTTTCCCCCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCaa                                                                                                        <  2:903562/149‑1 (MQ=255)
                                                tgtgTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggg                                                                                                           >  2:396177/1‑142 (MQ=255)
                                                tgtgTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGaa                                                                                                  >  2:322155/1‑150 (MQ=255)
                                                tgggccgggggCTCCCGTGTTTCTTTCCCCCCAACACCCGCGCTTTTACCCACCCCCATGTGAAAAATTCCCTCCGCGGCGGTTGTCTGCAAGTTTTTTTCCCGCCACAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGaa                                                                                                  <  3:148192/141‑1 (MQ=255)
                                                     cAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGCAGAGaaaaaaa                                                                                             >  2:669123/1‑146 (MQ=255)
                                                       gCGGCTACCGTGTTTCATTCCCCCCAACCCCCGCGCTTTTCCCCACCGCCATGGAAAAAAATGCCTGCGGGGGGGCTGCCCGCCATTTTTTTTCCCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGc                                                                                             <  3:1105348/148‑1 (MQ=255)
                                                         ggCCACCGTAATTTTTTCCCGCCAACACCCGCGCTTTTATCAACCGACTTGTGAAAAAATGCCTTCGCGGGGGTTGTCTGCCATTTTTTTTGCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCtt                                                                                          <  2:1397977/148‑1 (MQ=255)
                                                          gCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGCAGGGaaaaaaaattt                                                                                         >  2:1147582/1‑141 (MQ=25)
                                                              ccgtgttttCTTTTCCGCCAACACCCGCGCATTTGCCAACCGCCATATACAAAATTGCCTCCGCGGCGACTCCCTCCAACCTTTTTTCCCGCCACAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTaa                                                                                       <  2:1067470/140‑1 (MQ=255)
                                                              ccGTTATTCTTTCCCGCCAACACCCCCCCATTTATCCACCCCATTGGGCAAAAACGCCTCCGCGGCGGTTGTCTCCCAGCTTTTTTCCCGCCGCAATAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCt                                                                                           <  2:1634742/142‑1 (MQ=255)
                                                               cGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGAGgaaaaagatatagca                                                                                     >  3:991055/1‑136 (MQ=25)
                                                                     ttCATTCCCGCCAACACCCGCGCTTTTACCCACCGCCATGTGAAAAAATGCCTCCGCGGCGGCTGTCTGCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCAtg                                                                             <  2:884911/150‑1 (MQ=255)
                                                                               ccAACACCCGCGTTTTTACCAACCGCCATGTAAAAAATTGCCTGCGGGGGGGCTTTCTGCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCaa                                                                  <  2:1595156/151‑1 (MQ=255)
                                                                                       cccGCGCTTTTACCCACCCCAATGTGAAAAAATGCCTCCGCGGGGGCTTTCTCCCAGTTTTTTTCCCGCCGAAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATc                                                                <  3:1335253/146‑1 (MQ=255)
                                                                                          gcgcTTTTACCAACCGCAATGTGAAAAAAGCCTTCGGCGGCGGCTGTCTCCCAGTTTTTTTCCCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCg                                                               <  2:377633/145‑1 (MQ=255)
                                                                                           cgtaTTTACCCACCCCCTTGTCAAAAATGGCCTTCGGGGCGGCTTTCTGCCAGTTTTTTTCCCCCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTgg                                                        <  3:1023012/148‑1 (MQ=255)
                                                                                                    ccacccccATGTAAAAAAATGCCTCCGGGGCGGCTGTCGCCCAGTTTTTTTCCCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAAtg                                                <  3:1157216/144‑1 (MQ=255)
                                                                                                     cAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGAGAAAAAAATTTTAACCCCCCCTGTTTGGTGGGGAACGCTGCTgggtaaaa                                                  >  1:10347/1‑141 (MQ=25)
                                                                                                        cGCCATGTGAAAAAATGCCTCCGCGGCGGCTGTCTCCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGaa                                           <  2:52696/151‑1 (MQ=255)
                                                                                                         gCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGCAGGGAAAAAACTTTTAAACCCGCCTATTTTGTGGGGAACGCTGCTGGttgcaa                                                  >  2:730861/1‑138 (MQ=21)
                                                                                                               tgCAAAAATCCCTTGGCGGCGGCTGTCTGCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGaa                                     <  2:155345/150‑1 (MQ=255)
                                                                                                               gtttaataaggCCCTCGCGGCGGCGGTCTGCCAGCTATTTTTCCGCCATAAGAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAAGCGCTGCTGGTAACAATGTGGAAAGGGaa                                     <  2:1379316/141‑1 (MQ=255)
                                                                                                                  aaaatTCCCTCCGCGGCGGCTGTCTGCCAGTTTTTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTg                                  <  2:922209/150‑1 (MQ=255)
                                                                                                                  aaaagTGCCTTCGCGGCGTCTGCCGGCCGGTTTTTTTCCCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTg                                  <  3:1534811/150‑1 (MQ=255)
                                                                                                                  aaaaaTCCCTTCGGGGGGGTTGTCTGCCAGTTTTTTTCCCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAggg                                       <  2:660253/145‑1 (MQ=255)
                                                                                                                    aattGCCTCCGCGGGGGTTGTCTCCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTgcgg                               <  2:506257/148‑1 (MQ=255)
                                                                                                                    aactGCCTCGGCGGCGGCTTCCTGCAAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAg                                    <  3:1343315/143‑1 (MQ=255)
                                                                                                                    aaaaccccTCCGGGGCGCCTGTTTGCCATCTTTTTTCCCGCCCAAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTgcg                                <  2:1071955/145‑1 (MQ=255)
                                                                                                                     aattcccTTCGGGGCGGTTTTCTCCAAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTgcgg                               <  3:1253449/146‑1 (MQ=255)
                                                                                                                      aTGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAAAGAAAAAAATTTTAACCCCCCCTGTTTTGTTGGGAACGCCGCTGGGTACAAAGTGGGAAGGGGAgggggggg                             >  2:813301/1‑144 (MQ=255)
                                                                                                                      aTGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCGAGAGGAAAAAATTTTAACCCCCCCATTTTTGTGGGGAACCCTGCTGGGAAAAAAATGGGAAggggaaggggg                               >  2:1181885/1‑141 (MQ=255)
                                                                                                                          ccTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGAGAAAAAACATTAAATCCGCCATGTTTTGTGGGAAATCCTGCCGGTGAAAAAATGGGAAgggggggggggggg                            >  2:554405/1‑138 (MQ=255)
                                                                                                                                              cctttttttttcCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCgcgc         <  3:847256/142‑1 (MQ=255)
                                                                                                                                              ccAGTTTTTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTa       <  2:927245/151‑1 (MQ=255)
                                                                                                                                              ccAGTTTTTTTCCCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAcc             <  3:1442934/145‑1 (MQ=255)
                                                                                                                                              ccAGTTTTTTTCCCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACcgcg          <  2:1199181/148‑1 (MQ=255)
                                                                                                                                              ccAGTTTTTTTCCCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCGGCTCCAGAACCGCGCTa       <  2:166166/151‑1 (MQ=255)
                                                                                                                                              ccAGTTTTTTTCCCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACcg            <  2:490176/146‑1 (MQ=255)
                                                                                                                                              caattttttttcCCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGaa               <  3:1144159/136‑1 (MQ=255)
                                                                                                                                              caagtttttttcCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACcgc           <  3:881320/140‑1 (MQ=255)
                                                                                                                                                     ttttcCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGAGCGGCGAATCAGCTCCAGAACCGCGCTAATTGa  <  2:457546/149‑1 (MQ=255)
                                                                                                                                                     ttttcCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACcgcg          <  3:371708/141‑1 (MQ=255)
                                                                                                                                                         |                                                                                                                                                  
TCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGA  >  NC_000913/357732‑358024

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