Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,560,9960GT24.6% 727.0 / 75.9 379intergenic (+98/‑372)yffL/yffMCPZ‑55 prophage; uncharacterized protein/CPZ‑55 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (148/137);  new base T (31/62);  total (180/199)
Fisher's exact test for biased strand distribution p-value = 1.89e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 385 total aligned reads displayed.

AAAGGGAGTGCAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGC  >  NC_000913/2560846‑2561146
                                                                                                                                                      |                                                                                                                                                      
 aaGGGAGTGCAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGt                                                                                                                                                       <  5:1577793/151‑1 (MQ=255)
        tGCAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGTGTCAAAAGTTGCCGTTAAgg                                                                                                                                                  >  2:840838/1‑149 (MQ=255)
         gCAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTgg                                                                                                                                               <  1:1808925/151‑1 (MQ=255)
         gCAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTgg                                                                                                                                               <  2:288975/151‑1 (MQ=255)
          cAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTggg                                                                                                                                              >  7:1292162/1‑151 (MQ=255)
           aaCGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGt                                                                                                                                             <  3:1610098/151‑1 (MQ=255)
             cGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTc                                                                                                                                            >  6:419702/1‑150 (MQ=255)
             cGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCa                                                                                                                                           <  8:1200889/151‑1 (MQ=255)
              gAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCaa                                                                                                                                          <  7:970925/151‑1 (MQ=255)
              gAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCaa                                                                                                                                          <  1:811235/151‑1 (MQ=255)
                aTACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTGAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCaat                                                                                                                                         >  3:1090692/1‑149 (MQ=255)
                 tACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGAGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAATTTGCCTTTAAGGAGTGTCaaa                                                                                                                                         >  4:820667/1‑149 (MQ=255)
                   ccAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGtt                                                                                                                                     <  7:1035443/151‑1 (MQ=255)
                         attattCAAGCCGTACAACATGAAATAACGGAATCAGCCGGGCAATGTTCACGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGCGTCAAAAGTTGCCGTTAGTGCTGTCTCAAAAGATGACGTTAAGGAGGGACAAAAGTTgaacat                                                                                                                               >  6:97602/1‑146 (MQ=255)
                           tattCAATCCGTAGAAACTGAAATTACGGAATCAGGCGGGGATTGTTCCGGCCTTTTGTTTGCGTTAATTCGATGAATTCAGGTGGGTCAAAAGTTGGCGTTTGTGGTGGGTGAAAAGTTCGCGTTAAGGTGGGTCAAAAGTTGCCGTTa                                                                                                                              <  1:1015656/150‑1 (MQ=25)
                            attCAAGCCGTACAACATGAAATAACGGAATCAGTCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAg                                                                                                                            <  8:804635/151‑1 (MQ=255)
                             ttCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAATGGTGTGTCAAAAGTTGCCGTTAAGGTGGCTCAAAAGTTGAAGTTAAgg                                                                                                                           >  6:232743/1‑151 (MQ=255)
                                aaGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGgggg                                                                                                                        <  1:327764/151‑1 (MQ=255)
                                  gCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGt                                                                                                                      <  3:1309237/151‑1 (MQ=255)
                                   ccGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTGTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTc                                                                                                                     <  7:380675/151‑1 (MQ=255)
                                    cGTACAACATGAAATAACGGAATCAGTCGTGCATTGTTACCGCCTTTTCTTTGCCTTAATTACATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTAAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCa                                                                                                                    <  3:957309/151‑1 (MQ=38)
                                    cGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCa                                                                                                                    <  6:319312/151‑1 (MQ=255)
                                    cGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGATAGTGGTGGGTCAAAAGTTGCCTTTAAGGTGTGTCAAAAGATgaagtgtagggggggc                                                                                                                     >  4:1081852/1‑135 (MQ=255)
                                       acaacaTGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCaaa                                                                                                                  <  1:1374144/150‑1 (MQ=255)
                                       acaacaTGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGGAAAAATTTGCCGTTAAGCTTGGTCAAAATTTGACGtgtagctggctgaaa                                                                                                                  >  1:1636001/1‑136 (MQ=255)
                                        caacaTGAAATAACGGAATCCGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAg                                                                                                                <  5:821218/151‑1 (MQ=255)
                                        caacaTGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTATTGGTGGGTAAAAATTTGCAGTTAAGGTGGGTCGAAAGTTTAAGTTAAGGTGGGTCaaaat                                                                                                                >  2:740350/1‑150 (MQ=255)
                                         cacaTGAAATAACGGAATCAGGCGGCCAATTTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGt                                                                                                               <  1:16789/150‑1 (MQ=255)
                                                 aTAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTGAAAAGTTGCCGTGAAGGTGGGGCTAAAGTTGCCTAgg                                                                                                       >  7:270643/1‑150 (MQ=255)
                                                  tAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGGGGGTCAAAAGTTGCCGGTAAGGTGGGTCTAAAGTTGCCGAGt                                                                                                       >  5:769608/1‑150 (MQ=255)
                                                   aaCGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTGAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGACGAGTa                                                                                                      >  5:1024273/1‑150 (MQ=255)
                                                     cGGAATCAGGCGGGCAATTTTCCCGCCTTTTGTTTGCCTTAATTACATGAATTCAGGGGGGTCAAAAGTAGCCGTTAGGGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGa                                                                                                   <  5:781684/151‑1 (MQ=38)
                                                           cAGGCGTGCATTGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTct                                                                                              <  6:1819578/150‑1 (MQ=255)
                                                            cgtagggCATTTTTCCCGGCCTTTGTATGCCATTATTTAATGAATTCTGGTGGGTAAAAATTTGGCCTTACTGGTGGGTCAATAGTTGCCTTTTAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTctac                                                                                            <  6:1487906/147‑1 (MQ=16)
                                                            aGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAATTTGCCGTGAAGGTGGGTCAAAAGTTGCCGTTAAGGTGTCGCAAAAGTTGCCTAGTAAGCAGTctac                                                                                            >  2:1336450/1‑151 (MQ=38)
                                                            aGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCAGTTAAGGTGGGGCAAAATTTGACGAGTAAGAATTcga                                                                                             >  4:1749050/1‑148 (MQ=255)
                                                            aGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGACGTTAAGGTGGGTCAAAATTTGCCTAGTAAGAATTcta                                                                                             >  3:47102/1‑150 (MQ=255)
                                                              ccGGGCACTCTTCCAGCCGTATGTTTTCCTTAAATCAATGAATTCAGGTTGGTAAAAATTTGCCGTTAGTGGTGGGTCAAAAGTTGCCTTTAAGGTGGGTTAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTctacta                                                                                          <  2:377076/150‑1 (MQ=21)
                                                               cGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGCTCAAAAGTTGCCGTTACTGTTGGGTAAAACTTTGCCGTTAAGGAGGCTAAAAATTTGAAGTTTACGGGGGTCAAAAGTTGACGAGTAAGAATTcttcaat                                                                                         >  4:693050/1‑146 (MQ=255)
                                                                  gcagtcttCCAGCCTTTTCTTTGCCTTAATTTACTGCAGTGTGTGGGGTCAAAAGTAGCCGTTCTTGGTGGGTCTAATGATGCAGTTAAGTTGGGTCAAAAGTTGCAGTTCAGGTGTGTCAAAAGTTTCCGAGTAAGAATTCTACTATTaa                                                                                      <  6:703290/145‑1 (MQ=255)
                                                                  gatttcttcctgtattttcttttcCTTTAATCAATGAATTCATGTGTGTAAAAAGTTGCCTTTAGTGTTGGGTCAAAAGTTGCCTTTATGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTa                                                                                       <  2:599972/136‑1 (MQ=21)
                                                                  gCAATGTTCCCGCCTTTTGTTTGCCTTAATTCCATGAATTCAGGGGGGTCAAAAGTAGCCGTTAGTGGTGGGTAAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTaa                                                                                      <  1:341176/151‑1 (MQ=255)
                                                                  gCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGTGTGAAAAGTTGCCGAGAAAGAATTCTACTATTa                                                                                       >  3:269888/1‑150 (MQ=255)
                                                                  gCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTAAAAATTTGCCGTTAAGGTGGGTCAAAATTTGACGTTAAGGTGGGGCAAAATTTGCAGAGAAAGAATGCTACAAtat                                                                                       >  3:740362/1‑148 (MQ=37)
                                                                      tGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAATGTAGCCGTTAAGCTGCGTCAATAGTTGGCGTTAAGGTGGGTCAAAAATTGCCGAGTAAGAATTCTAATATTaatatt                                                                                  >  3:606400/1‑147 (MQ=38)
                                                                        ttCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAATGGTGGGTCAAAAGTTGCCTTATAGGAGGCTCAATAGTTTAAGTGAAGGTGGGTCAAATTTTGACTAGTAACGAATCTACAATTAAAAAtt                                                                                 >  6:835487/1‑150 (MQ=19)
                                                                        ttCCAGCCTTACCTTTGCCTTAATTACATGAATTCAGGTGGGTAAAAAGTTGCCGTTAGTGGTGGCTAAAAAGTTGCCTTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAAtt                                                                                 <  1:446397/150‑1 (MQ=38)
                                                                         tCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTATTGTTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAATTTGCATTAAAGGTGGGGCAAATTTTGAAGAGTAACAATTGTAATATTATGAATTCt                                                                               >  2:1558438/1‑151 (MQ=21)
                                                                          cccGCCTTTTCTTTCCCTTAATTCAATGAATTCCTGTGGCTCAAAAGTAGCCCTTCGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTa                                                                              <  3:381463/151‑1 (MQ=37)
                                                                           aaggcTTTTGTTTTCCTTTATTTAATGAATTCATTTTGGTTAAAATTTGGCTTTATTGGTGGGTATATTGTTGCCTTTTATTTGGGTCAAAAGTTGCCTTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTa                                                                              <  4:445416/146‑1 (MQ=16)
                                                                            agCCTTTTCTTTGCCTTAATTACATGAATTCAGGTGTGTCAAAAGTTGCCTTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAgg                                                                            <  6:1017831/150‑1 (MQ=255)
                                                                                ttttCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTGAGTGGTGGGTAAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGGCAAAAGTTGCCGAGTAAGAATTCTACGATTATCATTTATAGGCgcc                                                                        >  7:1091787/1‑149 (MQ=37)
                                                                                     tttGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGCGTCAAAAGTTGCAGTTAAGGTGGGGCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCCACGGGat                                                                   >  7:610609/1‑150 (MQ=255)
                                                                                       tGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAATGTTGCCGTTAAGGTGGGTCAATAGTTGCCGTGAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACGATTAAGAATTCTAGGCGACGGGAGtt                                                                 >  6:1644120/1‑151 (MQ=255)
                                                                                          cataaataaaTGAATTCCGGGGGGTAAAAAGCAGGCTTTAGTGGTGGGTAAATACCTGCCTTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGttttt                                                              <  5:1688794/142‑1 (MQ=25)
                                                                                          cTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGATAGTGGTGTGTCAAAATTTGCCTTGAAGGAGAGTCCAAATTTGCAGATAACCGGGGGGAAAATTTGAATAGAAAGAATTGTAAGATAAAAATTTCTAggacacgtgatatttt                                                              >  8:1496110/1‑137 (MQ=255)
                                                                                            tAAATCCAGGAATTCAGGGGGGCCAAAAGCTGCCGTTAGGGTTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGGTGCCCTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGtttttt                                                             <  3:83257/150‑1 (MQ=37)
                                                                                                tcaaTGAATTCAGGTGGGTCAAAAGCTGCCGTTAGTGGTGGGTCAAAAGCTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGcc                                                         <  3:1403626/147‑1 (MQ=255)
                                                                                                 ccATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCAGTTAAGGAGGGTCAAAATTTTAAGTTAAGGTGTCGGTAATTTTGCCGAGGAACCCGTCTACTATTAAGATTTCAAGCCGCGGCCAGTATTTAAGCCTc                                                       >  5:1866195/1‑151 (MQ=12)
                                                                                                  cATGAATTCAGGTGGGTCAAAAGCTGGCGTTAGGGGTGGGTCAAAAGGTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTcc                                                      <  1:748478/151‑1 (MQ=255)
                                                                                                      aaTTCAGGTGGGTAAAAATAAGCCGTTACGGGTGGGTAAATAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCcggt                                                  <  1:1096177/151‑1 (MQ=37)
                                                                                                       aTTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCTTTAAGGTGTGTCAAAATTTGGCGTTAAGGTGGGTCTAATTTTGCATAGTAAGAATTCTACAATTAAGAATTGTAGGCGACGTGAGTTTTTTAGcgttacgg                                                   >  4:344552/1‑142 (MQ=255)
                                                                                                       aTTCAGGAGGGTTAAAATCTGCCGTTACTGGTGGGTAAATTGTTGCCTTTTAGGTGGGTCAAAAGTTGCCTTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCcggtc                                                 <  1:207296/151‑1 (MQ=37)
                                                                                                         tCAGTTGGGGCAAAAGTTGCCGTTACTGGGGGGTCAAAAGCTGCCTTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCcggtcg                                                <  1:1479130/150‑1 (MQ=38)
                                                                                                              ggggTGAAAATTAGCGGTTATGGGTGGGTATAATGTTGCCTTTTATTTGCGTCAAATGTTGCCTTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATa                                          <  3:1241870/150‑1 (MQ=255)
                                                                                                                ggTCAAAATTAGCCTTTAGTGGTGGGTCAAAAGTTGCCTTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATag                                         <  8:649492/150‑1 (MQ=255)
                                                                                                                 gTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATagag                                       <  5:300156/151‑1 (MQ=255)
                                                                                                                   cAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTGAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATagag                                       >  8:1687091/1‑149 (MQ=255)
                                                                                                                   cAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGCCGGGAGTTTTTTAGCCTCCCGGTCGGTCATagag                                       >  1:1263589/1‑149 (MQ=255)
                                                                                                                   cAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATagag                                       >  6:582326/1‑149 (MQ=255)
                                                                                                                   cAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGACTAACGGTCGGTCATagag                                       >  6:78460/1‑149 (MQ=255)
                                                                                                                    aaaaGTTGCCCTTAGGGGTGGGTCAAAAGGTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGt                                    <  1:1022005/151‑1 (MQ=255)
                                                                                                                     aaaGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTc                                   >  4:1668204/1‑151 (MQ=255)
                                                                                                                      aagtatgctttATTGGTGGGTAAAATGTTGCCTTTAAGTTGGGTAAAAAGTTGCCTTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCg                                  <  4:1460421/142‑1 (MQ=38)
                                                                                                                      aaGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTGAAAATTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACTGAAGTTTTTAAGCCTACCGGTCGGTCATAGAGAGTCg                                  >  8:377429/1‑151 (MQ=38)
                                                                                                                          tGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGGTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCt                               >  6:1190285/1‑150 (MQ=255)
                                                                                                                          cgCCCTTAGTGGTGGCTCAAAAGGTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTa                              <  5:600648/150‑1 (MQ=255)
                                                                                                                           ggcGTTAGGGGTGGGTCAAATCTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTag                             <  5:491731/149‑1 (MQ=255)
                                                                                                                            gctttAGGGGGGGGTCAAATGTTGCCGTTAATGTGGGTCAAAAGTTGCCGTTCAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTaga                            <  8:134838/148‑1 (MQ=255)
                                                                                                                             cGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGGCGGGAGTTTTTTAGCCTCGCGGTCGCTCATAGAGAGTCTCGTag                             >  4:1281380/1‑149 (MQ=255)
                                                                                                                              gTTACTGGTGGGTCAAAAGTTGCCTTTAAGGTGGGTCAAAAGGTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGt                          <  5:1372343/151‑1 (MQ=255)
                                                                                                                                  gtggtgGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTaa                       >  6:134052/1‑150 (MQ=255)
                                                                                                                                    gttgggtattttgTTGCCTTTTAGTTGGGTTAAATGTTTGATTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAAc                     <  2:1194624/138‑1 (MQ=255)
                                                                                                                                    ggtgggtaatttgTTGCCTTTAAGTTGGGTCAAAAGTTGCAGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAAc                     <  8:1063823/138‑1 (MQ=39)
                                                                                                                                    ggtggGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAAc                     >  4:1544873/1‑150 (MQ=255)
                                                                                                                                     gtggGTCAATAGTTGGCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACa                    <  3:1778563/150‑1 (MQ=255)
                                                                                                                                      ggggTCAGATGCTGCCGTTACTGTGCGTTAAAACCTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACaga                  <  7:1488518/150‑1 (MQ=255)
                                                                                                                                          tCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCGGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACagaaga               >  8:2581/1‑150 (MQ=255)
                                                                                                                                            aaaaGTTGCCGTTAAGGTGGGTAAAAAGTTGCCGTTAATGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAatt            <  6:1445519/151‑1 (MQ=255)
                                                                                                                                             ataGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAatta           <  4:268403/149‑1 (MQ=255)
                                                                                                                                              aaGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAattat          >  5:1216674/1‑151 (MQ=255)
                                                                                                                                               aGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAattatt         >  7:1843006/1‑151 (MQ=255)
                                                                                                                                                 ttGCCGTTAAGGTGGGTCAAAAGTTTGCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTAttt        <  2:667882/150‑1 (MQ=255)
                                                                                                                                                  tGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCCGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTcg      >  7:1637446/1‑151 (MQ=255)
                                                                                                                                                   gCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTcgc     >  2:1478865/1‑151 (MQ=255)
                                                                                                                                                    ccGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTcgcc    >  3:1021093/1‑151 (MQ=255)
                                                                                                                                                     cGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTcgccg   <  8:58241/151‑1 (MQ=255)
                                                                                                                                                     cGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTAGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCCCTCATAGAGAGTCGCCTAGAGTGAAACAGAACAATTATTTcgcc  >  7:661530/1‑150 (MQ=255)
                                                                                                                                                      gTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTcgccgc  >  2:1043535/1‑151 (MQ=255)
                                                                                                                                                      gTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTcgccgc  >  8:1774500/1‑151 (MQ=255)
                                                                                                                                                      |                                                                                                                                                      
AAAGGGAGTGCAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGC  >  NC_000913/2560846‑2561146

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: