Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913357,8730CT42.5% 105.2 / 41.0 80S273N (AGC→AAC) cynRtranscriptional activator of cyn operon; autorepressor
Reads supporting (aligned to +/- strand):  ref base C (33/13);  new base T (0/34);  total (33/47)
Fisher's exact test for biased strand distribution p-value = 4.64e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑‑GGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGA  >  NC_000913/357725‑358010
                                                                                                                                                       |                                                                                                                                             
aCTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTa                                                                                                                                             <  1:955844/151‑1 (MQ=255)
 cTGACCTCGTAGTCTTGATACGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTa                                                                                                                                             <  4:584492/150‑1 (MQ=255)
 cTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTa                                                                                                                                             <  4:1188186/150‑1 (MQ=255)
 cTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTAt                                                                                                                                            >  4:602084/1‑151 (MQ=255)
  tGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTAtt                                                                                                                                           >  1:1204611/1‑151 (MQ=255)
   gAACTCGTAGGCCTGCTAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGAGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGACATGTGCAAAAATGCCTTCGCGGCGGCTGTTTGCCAGCTAttt                                                                                                                                          >  1:1252609/1‑151 (MQ=255)
   gAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTAttt                                                                                                                                          <  4:621991/151‑1 (MQ=255)
   gAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTAt                                                                                                                                            <  1:319183/149‑1 (MQ=255)
     aCTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTAttttt                                                                                                                                        <  4:1066164/151‑1 (MQ=255)
     aCTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTAtttt                                                                                                                                         <  4:331863/150‑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)
                                                                                                                                                       |                                                                                                                                             
ACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCT‑‑GGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGA  >  NC_000913/357725‑358010

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