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

ACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑‑GGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTT  >  NC_000913/357730‑358027
                                                                                                                                                         |                                                                                                                                                       
aCTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTAttttt                                                                                                                                                         <  4:1066164/151‑1 (MQ=255)
  tCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTcc                                                                                                                                                       >  2:189856/1‑151 (MQ=255)
   cGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTccg                                                                                                                                                      <  4:169584/151‑1 (MQ=255)
    gTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTccgc                                                                                                                                                     <  4:141328/151‑1 (MQ=255)
     tAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTccgcc                                                                                                                                                    >  4:1328856/1‑151 (MQ=255)
      aGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTccgccg                                                                                                                                                   <  1:1292118/151‑1 (MQ=255)
          cTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGcaac                                                                                                                                               >  1:1453369/1‑151 (MQ=255)
            gATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGcaaca                                                                                                                                              <  1:1055590/150‑1 (MQ=255)
               aaGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAc                                                                                                                                          <  4:311546/151‑1 (MQ=255)
                aGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAATTGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAcc                                                                                                                                         >  1:960258/1‑151 (MQ=255)
                  cGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAccg                                                                                                                                        <  4:842140/150‑1 (MQ=255)
                    tAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAccgccg                                                                                                                                     >  1:80790/1‑151 (MQ=255)
                     aGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGt                                                                                                                                    >  1:841046/1‑151 (MQ=255)
                       cgcATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTtc                                                                                                                                  >  4:1024163/1‑151 (MQ=255)
                       cgcATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTtc                                                                                                                                  >  2:174553/1‑151 (MQ=255)
                                 aaTTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCagt                                                                                                                          >  4:709733/1‑149 (MQ=255)
                                   ttCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTg                                                                                                                      >  4:1091658/1‑151 (MQ=255)
                                     ccccccGCAGCCGTGTGTCGGGGCCTCCCGTGTTTCTTTCCCCCCAACACCCGCGCTTTTATCAAACGCCTTGTAAAAAATTCCCTCCGGGGCGGTTGTCGGCCGGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCagtag                                                                                                                        <  3:180308/143‑1 (MQ=255)
                                      cAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCg                                                                                                                   >  4:643699/1‑151 (MQ=255)
                                       aGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgg                                                                                                                  >  4:265081/1‑151 (MQ=255)
                                       aGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgg                                                                                                                  >  1:1248575/1‑151 (MQ=255)
                                         ccGCAGCCTTGTGTCGGCGCCTCCCGTTTTTCTTTCCCGCAAACACCCGCGTTTTTACCCACCGCCGTGTAAAAAAATGCCTCCGCGGCGGTTGTCTGCCAGTTTTTTTCCCGCCGCAAAAAACCCGCCGTTCTCTCCAGTAGTGGCggg                                                                                                                 <  4:992550/150‑1 (MQ=255)
                                         ccGCAGCCTTGTGTCAGCGGCTCCCGTGTTTCTTTCCCCCCAACACCCGCGCTTTTCCCCACCCCCTTGTGCAAAAATCCCTCCGGGGGGGCTGTCTGCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCg                                                                                                                   <  1:983614/148‑1 (MQ=255)
                                         ccGCAGCCCTGTGTCAGCGGCTACCGTGTTTCTTTCCCGCCAACACCCGCGCTTTTACCCAACGCCATGTGAAAAAATCCCTGCGCGGCGGCTGTCTCCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCgg                                                                                                                  <  1:994694/149‑1 (MQ=255)
                                         ccGCAGCCCTGTGCCGGCGCCCCCCGTGTTTCTTTCCCCCAAACACCCGCGTTTTTACCCACCCCCATGTAAAAAATTGCCTCCGGGGGGGCGGCCGCCCAGTTTTTTTCCCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCggg                                                                                                                 <  4:292756/150‑1 (MQ=255)
                                         ccGCAGACTTGTGTCGGGGCCTCCCGTGATTCTTTCCCCCCAACACCCGCCCTTTTACCCAACGCCTTGTGAAAAATTCCTTCCGGGGCGGTTGTCGGCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCgg                                                                                                                  <  1:56401/149‑1 (MQ=255)
                                          cGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCggggg                                                                                                               >  4:496844/1‑150 (MQ=255)
                                          cGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggg                                                                                                                >  1:372179/1‑150 (MQ=255)
                                            cagccgtgtgCCAGCGCCTCCCGTGATTCTTCCCCCCCAACACCCGCGCTTTTCCCCAACGCCATGTGAAAAAATGCCTCCGGGGCGTCTGTCTGCCAGTTTTTTTCCCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCaa                                                                                                             <  4:1232255/145‑1 (MQ=255)
                                              gcctgtgcccAGCGGCAACCGTGTTTCTTTCCCGCAAACACCCGCACTTTTACCAAACGCCAGGTAAAAAATGCCCGCCGGGGGGGCTGTCTGCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCaa                                                                                                             <  4:1202687/140‑1 (MQ=255)
                                              gcccTGTGCCAGGGCCTCCCGTGTTTCTTTCCCCCCAACACCCGCGCTTTTACCCACCGCCTTGTGAAAAAATCCTTCCGGGGGGGCTGTCTGCCATTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAaga                                                                                                           <  1:156131/149‑1 (MQ=255)
                                               aCCTGTGTCAGCGGCTCCCGTGATTCATTTCCCCCAACACCCGCGATTTTACCCAACCCCATGTCAAAAAATGCCTTCGCGGGGGCTGCCTCCCAGTTTTTTTCCCCCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCgg                                                                                                                  <  1:161173/143‑1 (MQ=255)
                                                ccttgtccggCGGCTCCCGTGATTCATTCCCCCCAACACCCGCGCTTTTACCCACCCCCTTGTGCAAAATTCCCTCCGCGGCGGTTGCCTGCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAag                                                                                                            <  1:569720/139‑1 (MQ=255)
                                                 cTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCggggggggag                                                                                                          >  4:15904/1‑143 (MQ=255)
                                                 cTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggg                                                                                                                >  1:878203/1‑143 (MQ=255)
                                                  tgtgTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCggggg                                                                                                               >  3:31396/1‑142 (MQ=255)
                                                    ttccAGCGGCTCCGTTGATTATTTCCCGCCAACACCCGCGCATTTACCCACCCCCTTGTGAAAAATCGCCTCGGGGGGGGCTGTCTGCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCa                                                                                                              <  1:882375/139‑1 (MQ=255)
                                                    tgccAGCGCCTCCCGTAATTCTTTCCCCCCAACACCCGCCCTTTTACCCACCCCAATGTGAAAAAATGCCTCCGCGGCGGCTGTCTCCCACTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAAt                                                                                                      <  1:924682/147‑1 (MQ=255)
                                                     gtCGGCGGCTCCCGTTATTCTTTCCCGCCAACCCCCGCGCATTTACCCACCGCCATGTGCAAAAAGCCCTCCGGGGGGGCTTTCTGCCGGTTATTTTCCCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAg                                                                                                    <  4:221312/151‑1 (MQ=255)
                                                       cAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggggcagcgcaaa                                                                                                      >  4:297361/1‑137 (MQ=255)
                                                         gCGGCTCCCGTGTTTCTTTCCCGCCACCACCCGCGCTTTTACCCACCGCCTTGTAAAAAAATGCTTCGGGGGCGGCTGTCTCCGAGTTTTTTTCCCCCCCCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAg                                                                                                    <  1:618516/147‑1 (MQ=255)
                                                         gCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAAAGAAAAAACttt                                                                                               >  1:1373059/1‑151 (MQ=255)
                                                         gCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCGAGAGgaaaaagatt                                                                                                >  4:745333/1‑144 (MQ=255)
                                                           ggCTCCCGTGATTCATTCCCGCCAACACCCGCGCTTTTACCCACCGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAg                                                                                                    <  1:52958/145‑1 (MQ=25)
                                                            gCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGAAAAAAATTTaa                                                                                             >  1:903117/1‑150 (MQ=255)
                                                               aCCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAAAGAAAAAGCCTTTAAAcc                                                                                          >  1:77570/1‑151 (MQ=255)
                                                                ccGTGTTTCTTTCCCCCCAACACCCGCGCTTTTACCCACCCCCATTTGAAAAATTCCCTCCGCGGCGGCTGTCTCCCAGTTTTTTTCCCGCCGAAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCg                                                                                        <  1:262317/151‑1 (MQ=255)
                                                                ccGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGCGGAGaaaaaaattttaacccc                                                                                         >  4:268662/1‑137 (MQ=25)
                                                                 ccgtaTTCCTTTCCCGCCACCAACCGCCATTTTACCCACCGCCATGTCAAAAAATGCCTCCGCGGGGGCTGTCTGCCAGTTATTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGccc                                                                                         <  4:313013/146‑1 (MQ=255)
                                                                    taattcttcccccccAACACCCGCGCTTTTACCCACCCCCTTGTAAAAAAAGCCCTCCGCGGGGGCTGCCTGCCAGTTTTTTTCCCGCGAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCt                                                                                                 <  4:723792/126‑1 (MQ=255)
                                                                      ttCTTTCCCCCCAACACCCGCGCTTTTACCCACCCCCTTGTGAAAAAATGCTTCCGGGGCGGCTGTCTCCCAGTTTTTTTCCCGCCGAAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCAtg                                                                                   <  1:324729/150‑1 (MQ=255)
                                                                           ttccccccAACACCCGCGCTTTTACCCACCCCCTTGTGAAAAAATCCCTCCGGGGCGGTTGTCGCCCAGTTTTTTTCCCCCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCAtgt                                                                                  <  1:663391/141‑1 (MQ=255)
                                                                             ccccccaccaCCCGCGCTTTCACCCAACCCAATTTGAAAAAATGCCGCCGCGGCGGCTGTCTGCCAGTTTTTTTCCCGCCCAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCAtgttgt                                                                               <  4:931769/144‑1 (MQ=255)
                                                                                        ccgcccttttccccACCGCCATTTCAAAAAATCCTTCCGCGGGGGCTTTCTCCCATTTTTTTTCCCCCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgctg                                                                 <  1:1440791/138‑1 (MQ=255)
                                                                                         cgctcttttACCCAACGCAATGTGAAAAAGTCCCTTGGCGGCGGCGTTCGGCCATTTTTTTTCCCGCCCCAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgctgct                                                               <  4:1030254/145‑1 (MQ=255)
                                                                                           cgcTTTTACCCACCCCCATGTAAAAAAATCCCTCCGGGGGGGCTTTCTGCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGctgct‑‑g                                                            <  1:1295545/150‑1 (MQ=255)
                                                                                               tttATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGCAGAGAAAAAGCATTTAACCCCCCCTGTTGTGTGGCAATCgc                                                                   >  1:399145/1‑141 (MQ=37)
                                                                                               tttACCCACCCCCATGTAAAAAATTCCCCCCGGGGGGGCTGTCTCCCATTTTTTTTCCCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑‑GGt                                                          <  4:146337/148‑1 (MQ=255)
                                                                                                        cGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGAAAAAACATTAAACCCCCCATCTTTTGTGGGAAACGCtggt‑‑gggtaaaaattg                                                 >  1:1154509/1‑138 (MQ=25)
                                                                                                             tgtgCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCGAGGGAAAAAACATTTAACCCCTCCTTTTGTGTGGGGAACGCTGCT‑‑GGGGACAAAATGGGAAggg                                          >  1:568071/1‑150 (MQ=255)
                                                                                                                  aaaaaTGCCTTCGGGGCGGTTGTCCGCCACTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑‑GGTAACAATGTGGAAAGGGAAg                                       <  1:462100/148‑1 (MQ=255)
                                                                                                                  aaaaaTCCCTTGGCGGCGGCTGTCGCCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑‑GGTAACAATGTGGAAAGGGaa                                        <  4:1336957/147‑1 (MQ=255)
                                                                                                                   aaaaTGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGAGAAATATATTTTACACCCCCTTGTTTTGTGGGCATCCCCGCTGGGGTAAAAATGGGAAAGGGagtggg                                    >  2:19597/1‑146 (MQ=255)
                                                                                                                    aaaTGCCTCCGGGGGGGCTGTCGCCCATTTTTTTTCCCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑‑GGTAACAATGTGGAAAGGGAAGTgcgg                                  <  4:429547/151‑1 (MQ=255)
                                                                                                                    aaaTCCCTTCGCGGCGGCTGTCTCCCAGTTTTTTTCCCCCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑‑GGTAACAATGTGGAAAGGGAAGTgcgg                                  <  1:1350146/151‑1 (MQ=255)
                                                                                                                     aaTGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAAAAAATTTTAACCCCTCCTGTTGTGTGGGAATCGCTGCT‑‑GGTGAAAAAGTGGAAAGGGGAGGgcggg                                 >  1:842602/1‑150 (MQ=255)
                                                                                                                      aTGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGGGAAAAAACATTAAACCCCTCCTGTTTTGTTGGGAACCGTGCT‑‑GGGTAAAATTTGGGAAGGGGAgggggggg                                >  3:91003/1‑144 (MQ=255)
                                                                                                                          ccTTGGGGGCGGCTTTCTGCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑‑GGTAACAATGTGGAAAGGGAAGTGCGGCGa                               <  1:888597/149‑1 (MQ=255)
                                                                                                                                  ggcggcTGTCTGCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑‑GGTAACAATGTGGAAAGGGAAGTgcggcg                                <  2:93248/141‑1 (MQ=255)
                                                                                                                                  ggcgcttttCTGCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑‑GGTAACAATGTGGAAAGGGAAGTGCGGCGaa                              <  4:1252503/135‑1 (MQ=255)
                                                                                                                                          tctcCCAGTTTTTTTCCCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑‑GGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGa                   <  1:645096/146‑1 (MQ=255)
                                                                                                                                            tgcccttttttttcCCGCCCAAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑‑GGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCa                     <  4:807373/133‑1 (MQ=255)
                                                                                                                                             cccattttttttcCCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑‑GGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACcgcg             <  3:103474/141‑1 (MQ=255)
                                                                                                                                              caagtttttttcCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑‑GGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGaa                  <  1:1193358/136‑1 (MQ=255)
                                                                                                                                                   ttttttcCCCCCGAAAAAAAGCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑‑GGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTg      <  4:28885/148‑1 (MQ=255)
                                                                                                                                                        ttCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGAGCCAGAGAAAAAGATTTTAAAACCTCAACTTGTGTGGGAATCCCTGCT‑‑GGGAACAAAGTGGGAAGGGGAGGGCGGCGGAACAGCTCCCGAACCCCGCCCATTggggt  >  2:135547/1‑145 (MQ=255)
                                                                                                                                                         |                                                                                                                                                       
ACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑‑GGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTT  >  NC_000913/357730‑358027

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: