Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,560,9960GT24.0% 647.9 / 61.2 339intergenic (+98/‑372)yffL/yffMCPZ‑55 prophage; uncharacterized protein/CPZ‑55 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (122/135);  new base T (30/51);  total (152/187)
Fisher's exact test for biased strand distribution p-value = 1.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 342 total aligned reads displayed.

AAAGGGAGTGCAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGC  >  NC_000913/2560846‑2561146
                                                                                                                                                       |                                                                                                                                                      
aaaGGGAGTGCAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCg                                                                                                                                                        <  1:1716800/151‑1 (MQ=255)
aaaGGGAGTGCAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCg                                                                                                                                                        >  4:1418715/1‑151 (MQ=255)
  aGGGAGTGCAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGTGTCAAAAGTTGCCGtt                                                                                                                                                      >  1:1338825/1‑151 (MQ=255)
     gactGCAACGAATACCAGCTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGTGAAATCTTCCCGCCTTTTCTTTGCCATAAATCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTACTGGTGGGTCAAAACTTGCCGTTAAg                                                                                                                                                   <  5:409926/148‑1 (MQ=255)
      aGTGCAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGACGTTACAGTTGGGTAAAATGATGCCTTTAAgg                                                                                                                                                  >  3:946689/1‑151 (MQ=255)
           aaCGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGt                                                                                                                                             <  8:819754/151‑1 (MQ=255)
            aCGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTc                                                                                                                                            <  2:1724961/151‑1 (MQ=255)
             cGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGATTTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTcc                                                                                                                                           >  6:230076/1‑150 (MQ=255)
             cGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCa                                                                                                                                           <  8:1685789/151‑1 (MQ=255)
               aaTACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAATGTTGCCTTTAAGGTGGGTCaa                                                                                                                                          >  3:1248700/1‑150 (MQ=255)
               aaTACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCaa                                                                                                                                          >  3:1550882/1‑150 (MQ=255)
               aaTACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGATAGTGTTGGGTAAAAAGTTGACGTTAAGCTGGGTCaa                                                                                                                                          >  1:927115/1‑150 (MQ=255)
                aTACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTGTCTGCTGGGTCAAAAGTTGCCGTGAAGGGGGGTGaatt                                                                                                                                        >  1:403233/1‑149 (MQ=255)
                  cccATTTATTATTCAAGCCGTAAAACATGAAATTACTGAATCAGTCGGGCAATGTTCCCGCGTTTTCTTTGCCTTAATTCCATGAATTCTGGTGGGGCAAAAGTTGCGGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAg                                                                                                                                       <  7:887850/149‑1 (MQ=255)
                   ccAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGtt                                                                                                                                     <  3:1297278/151‑1 (MQ=255)
                      gTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGcc                                                                                                                                  <  3:427934/151‑1 (MQ=255)
                         attattCAAGCCGTACAACATGAAATAACGGAATCAGGCTGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGtt                                                                                                                               <  7:1689632/151‑1 (MQ=255)
                         attattCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGtt                                                                                                                               <  4:1070652/151‑1 (MQ=255)
                         attattCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGtt                                                                                                                               <  4:987264/151‑1 (MQ=255)
                         attattCAAGCCGTACAACATGAAATAACGGAATCAGGAGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGtt                                                                                                                               <  2:1510222/151‑1 (MQ=255)
                             ttCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTGAAgg                                                                                                                           >  8:1669001/1‑151 (MQ=255)
                              tCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAgg                                                                                                                           <  3:754907/150‑1 (MQ=255)
                               cAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGTTGGGTAAAAATTTGACTTTAAGGTGGGTCAAAATTTGAAGTAAAGGTg                                                                                                                         >  5:78321/1‑151 (MQ=255)
                                  gCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTTACTTTAAGGTGGGt                                                                                                                      >  6:1289051/1‑151 (MQ=255)
                                   ccGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTc                                                                                                                     <  8:1101502/151‑1 (MQ=255)
                                        caacaTGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGCGCaaaat                                                                                                                >  7:983749/1‑150 (MQ=255)
                                          acaTGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGtt                                                                                                              <  3:555256/151‑1 (MQ=255)
                                           caTGAAATAACGGAATCAGGCGGGCAATGTTCCAGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTg                                                                                                             <  8:1714542/151‑1 (MQ=255)
                                            aTGAAATAACGGAATCAGGCGGGCAGTGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAATGGTGTGTAAAAATTTGCATTTAAGGAGAGTCAAAATTTGAAGTAAAGGTgtcggaaattt                                                                                                               >  5:1465447/1‑138 (MQ=255)
                                            aTGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAATGTTGCCGTTAAGGTGGGTGTAATGTTGc                                                                                                            >  2:262985/1‑151 (MQ=255)
                                              gAAATAAGGGAATCAGTCGGGCAATTTTCCCGCCTTTTCTTTTCCTTAATTCCATGAATTCAGGTGGCTCAAAAGAAGCCCTTAGTGGTGGGTAAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGcc                                                                                                           <  4:283757/150‑1 (MQ=37)
                                               aaaTCACGGAATCCGGAGCGCACTCTTCCCGCCTTTTCTATGCCTTATTTCCATGAATTCAGGGGGGTCAAAAGTGGGCCTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCg                                                                                                          <  5:1301715/150‑1 (MQ=25)
                                               aaaTAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGACAAAAGTTCCCGTTAGTGTAGTGTAAAATGTTGCATTCTAGGTGTGTCAAAATTTTACGATAACGTCTGAGAAAAGTTGCCg                                                                                                          >  1:1572731/1‑150 (MQ=255)
                                                     cGGAATCCGTCCGGCATTTTTCCAGCCTTTTCTTTGCCTTAAATCCATGAATTCTGGTGGGTCAAAAGCTGCCTTTAGTGGTGGGTAAAAAGTTGCCTTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGa                                                                                                   <  6:1597602/151‑1 (MQ=25)
                                                      ggAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCTAAAGTTGCCGAGTa                                                                                                      >  4:1729850/1‑147 (MQ=255)
                                                       gAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTAAAAATTTGCCTTTAAGGAGGGTCAAAAGTTGAAGAGAACGTGTGTGAAAATTTGACTAGTAAGCAt                                                                                                 >  1:1020644/1‑151 (MQ=255)
                                                        aagccggCGGGCAATCTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGGGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAAtt                                                                                                <  4:429602/146‑1 (MQ=255)
                                                             gccGGGCACTCTTGCAGCCCTTTCTATCCCTTAATTTAATGAATTCCTGGTGCTTTAAAGTTGCCGTTTTTGCTGGGTAAAATGTGGCCTTTTATGTGGGTGAAAAGATGCACTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTctact                                                                                           <  2:679144/149‑1 (MQ=255)
                                                                gggCACTGATCCCGGCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGGCCTTAGTGGTGGGTCAAAAGTTGCCTTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTAtt                                                                                        <  2:798473/151‑1 (MQ=37)
                                                                gggCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGACGTTAAGGTGGGTCAAAATTTGACGACTAAGAATTCTACTAtt                                                                                        >  6:520192/1‑151 (MQ=255)
                                                                gggCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGACGTTAAGGTGTGTCAAAATTTGAAGTTAAGGTGGGGCAAAAGTTGACTAGTAAGAATTCTACTAtt                                                                                        >  2:217689/1‑151 (MQ=38)
                                                                gggCAATGTTCCCGCCTTTTCTTTGCCTTAATTACATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAACGTGGCTCAAAATTTGGCTAGTAAGAATTCTAGTAtt                                                                                        >  2:1627362/1‑151 (MQ=37)
                                                                  gcattattCCAGCCTTTTCTTTTCCTTAATTCAATGAATTCATTTTGGTAAAAATTTGGCTTTACTGGTGGGTAAAAACTTGCCTTTAAGGTGGGTTAAAAGTTGCCTTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTaa                                                                                      <  1:1157745/145‑1 (MQ=21)
                                                                    actGTTCCGGTCTTTTCTTTGCCTTTAATTCATGAATTCAGGTTGGTCAATAGTTTCGGATAGTGGTTGGTTAAATGTTGCCGTTTAGGTGGGTAAAAAGTTTCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGa                                                                                    <  4:789043/149‑1 (MQ=25)
                                                                    aaTGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTTCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGa                                                                                    >  1:1472777/1‑150 (MQ=255)
                                                                    aaTGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGATAGTGGTGGGTAAAAATTTGCCGTTAAGGAGGGTCAAAATTTGAAGTTAAGGTGTCTCAAAAATTGAAGAGTAACAATTATACAATTa                                                                                       >  7:1407716/1‑148 (MQ=255)
                                                                       gatgcagCCTTTTCTTTGCCTTAATTCCATGAATTCAGGGGGGTCAAAAGATGCCGTTAGTGGTGGGTCAAAACTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAAt                                                                                  <  4:624368/144‑1 (MQ=255)
                                                                         tCCCGCCTTTTCTTTCCCTTAAATCCATGAATTCAGGTGGGTCAAAAGTTGCCGATACTGCTGGGTAAAAAGTTGCCGATAACGAGTGACAAAAGTTGCAGATAAGAAGGCTCAAAAGTTGACAAGTAAGAATTCTACAATaaatatatct                                                                               >  6:119163/1‑144 (MQ=255)
                                                                         tCCAGCCCTTTCTTATCCATAATTACATGAAGGCCGCCCGGTCTAAAGTAGCCCTTACTGCGGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCCTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTc                                                                                <  6:639971/150‑1 (MQ=18)
                                                                          gccGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAATGTTGCCGTTAGTGGTGGGTCAAATGTTGCCTTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTa                                                                              <  8:442275/150‑1 (MQ=255)
                                                                           ccGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCt                                                                               >  3:1334558/1‑149 (MQ=255)
                                                                           ccGCCTTTTCTATGCCTTCATTCAATGAATTCAGGTGGGTCATAAGTTGCCGTTATTGGGGGGTCAAAACTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTa                                                                              <  6:764681/150‑1 (MQ=37)
                                                                            ggccctagctatGCCATTATTTCATCAATTCCTGGGGGTCAAAATCTGGCCTTACTGGTGGGTAAAAAGTTGCCTTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAgg                                                                            <  3:1750044/140‑1 (MQ=21)
                                                                            cGCCTTTGCTTTGCCTTAATTGCATGAATTCAGGCGGGGAAAAAGTCGCCCTTAGTGGTGGGTCAAAACTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAgg                                                                            <  4:1227727/151‑1 (MQ=37)
                                                                             gCCTTTTCTTTGCCATAATTCCATGAATTCATGTGGGTGAAAAGTTGCCGTTAGTTGTGGGTCAAAAGTAGCCGTTAAGTTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGc                                                                           <  2:1107318/151‑1 (MQ=38)
                                                                                 tttCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAATGTTGCCGTTAAGGTGGGGCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGAc                                                                        <  2:423513/150‑1 (MQ=255)
                                                                                 tttCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAATGGTGGGTAAAAATTTGCCTTTAAGGAGGGTCAAAATTTTAAGTTAAGGTGTGTGAAAATTTGACTAGAAAGAATTATAAAATAAAGATTAATAGGCTACg                                                                       >  7:203196/1‑151 (MQ=16)
                                                                                 tttCTATGCCTTAAATACATGAATTCAGCGGGGTCAAAAGTTGCCCTTCGGGGTGGGTCAAAAGGTGCCGTTAAGGTGGGTCAAAACCTTCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGAc                                                                        <  4:694062/150‑1 (MQ=37)
                                                                                   gctataccatAAATCCATGAATTCCGGGGGGTCAAAAGAAGGCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACgg                                                                      <  8:642395/141‑1 (MQ=25)
                                                                                      tttacTTAATTAAATGAATTCAGTTGGGTATTAAAGTTGCCTTTATTTGGGGTTATTTACTTTCCTTTTATTTGGGTAAAAAGTTGCATTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGa                                                                    <  1:810236/146‑1 (MQ=255)
                                                                                      ttGCCTTAAATCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCTTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGa                                                                    <  8:1596698/149‑1 (MQ=255)
                                                                                        tccttttattcttGAATTCTGGTTGGTTAAAAGTTGGGTTTAGGGGTGGGTAAATTGTTGCCTTTTTGTTGGTTTAAATGTTGCTGTTATGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGtt                                                                 <  6:884924/138‑1 (MQ=255)
                                                                                        gCCTTAATTCCATGAATTCAGGGGGGTCAAAAGTCGCCCTTAGTGGTGGGTCAAAAGGTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGttt                                                                <  3:300500/151‑1 (MQ=255)
                                                                                           ttAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGACGTTAAGGTGGGTCAAAAGTTGGAGTTAAGGTGGCTCAAATGTTGACTAGTAAGAATTCTACTATTAAGATTTATACGCGACgtcatattttt                                                             >  3:1495260/1‑141 (MQ=25)
                                                                                             aaTTCCATGAATTCAGGGGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAg                                                           <  2:1046237/151‑1 (MQ=255)
                                                                                              aTTCCATGAATTCTTGTTGGTTTAAAGTTGGCGTTATTGGTGGGTCAAATGTTGCCTTTTAGGTGGGTCAAAAGTTGCCTTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGc                                                          <  5:601218/151‑1 (MQ=37)
                                                                                                  cATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAATTTGCCGTTAAGGTGGGGCAAAATTTGCCGTGAAGGTGGGGCAAAAGTTGCCGAGTAACAATTCTACAATTAAAAATACTAGGCGCCGTGAGATTTTTAgaatac                                                      >  4:127741/1‑146 (MQ=25)
                                                                                                   aTGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGGGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAATGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGAGCCGGGAGTTTGTTAGCCTccc                                                     >  3:1035300/1‑151 (MQ=255)
                                                                                                      aaTTCAGGTGGGTCAAAATTTGCCGTTAGTGGTGGGTCAAAAGTAGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCcggt                                                  <  4:153875/151‑1 (MQ=255)
                                                                                                      aaTTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTATGAATTCTAGGCGACGGGAGTTTTTTAGCCTCAcggg                                                  >  3:926936/1‑150 (MQ=255)
                                                                                                       aTTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCcggtg                                                 >  2:1328582/1‑150 (MQ=255)
                                                                                                       aTTCAGGTGGGTAAAAAGTTGGCGGTAGTGGTGGGGCAAAACTTGCCTTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGTTTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCcggtc                                                 <  3:659661/151‑1 (MQ=38)
                                                                                                         tCAGGTGGGTAAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCcggtcg                                                <  8:1267462/150‑1 (MQ=255)
                                                                                                          cacgTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGAAAGAATTCTACAAATAAGAATTCTAGGCGACGGGAGTTTTTTAGAGTCCcggtcggg                                              >  7:51103/4‑150 (MQ=37)
                                                                                                          cAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGGAGTTAAGGTGGGTGAAATTTTGACGAGTAAGAATTGTACTATTAAGAATTCTTGGCGACGGGAGTTTTGTAGAGTAGcggtcggt                                              >  3:377659/1‑151 (MQ=25)
                                                                                                           cagttttgtaaaaaGTTGCCGTTATTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCcggtcggt                                              <  7:746938/141‑1 (MQ=255)
                                                                                                                 gTCAAAAGTTCCCGTTACGGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGGTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATagag                                       <  3:1102366/151‑1 (MQ=255)
                                                                                                                      aaaGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGGGGGGCAAAATTTGAGGACTAAGATTTCTACTATTAAGAATTCTATGCAAATGGAGTTTTTTAGAATAACGGGCGGTCATAGATAGTc                                   >  7:1277997/1‑151 (MQ=21)
                                                                                                                        tttaggcctttttggtggGTATATATTTGCCGTTTAGTTGGGTCAAAAGTTTCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCg                                  <  1:212573/138‑1 (MQ=38)
                                                                                                                           agCCGTTATTGGTGGGTAAAAAGTTGCCTTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCt                               <  7:1546997/149‑1 (MQ=255)
                                                                                                                             ccGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCt                               >  6:431760/1‑148 (MQ=255)
                                                                                                                              gcttATTGGGGGGTAAAAAGTTGCCTTTAATGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTagag                           <  1:645721/149‑1 (MQ=38)
                                                                                                                                 tAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGtt                         >  7:314913/1‑150 (MQ=255)
                                                                                                                                 tAGTGGTGGGTCAAAAGGTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGtt                         <  2:546108/150‑1 (MQ=255)
                                                                                                                                 tAGTGGTGCTCAAAAGGAGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTa                        <  3:1588709/150‑1 (MQ=255)
                                                                                                                                  aGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTaa                       <  2:88288/151‑1 (MQ=255)
                                                                                                                                  aGTGGGGGCTCAAAAGCTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTCAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTaa                       <  8:1406912/151‑1 (MQ=255)
                                                                                                                                   gggggggctCAAAAGGTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTaaa                      <  3:1328973/143‑1 (MQ=255)
                                                                                                                                    tggtggGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTaaa                      >  4:201726/1‑150 (MQ=255)
                                                                                                                                     ggtggGTAAAATCTTGCCTTTTATGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACa                    <  5:1213910/151‑1 (MQ=38)
                                                                                                                                          gTCAACAGCTGCCGTTCAGGGGGCTCAAAAGTTGCCGTTACGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACagaaga               <  1:585650/151‑1 (MQ=255)
                                                                                                                                           tCAAATGTTGCCGTTTATGTGGGTAAAAAGTTTCCGTTAAGGTGGGTCAAAAGTTGCCGAGTATGAATTCTACTATTAAGAATTCTTGTCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACagaaga               <  7:1105813/150‑1 (MQ=38)
                                                                                                                                            aaaaaGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAat             <  8:1505130/150‑1 (MQ=255)
                                                                                                                                             aaaaGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAatt            >  4:903916/1‑151 (MQ=255)
                                                                                                                                               aaGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAattat          >  1:724272/1‑151 (MQ=255)
                                                                                                                                                    gCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTcgc     >  5:892537/1‑151 (MQ=255)
                                                                                                                                                    gCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTcgc     >  8:275522/1‑151 (MQ=255)
                                                                                                                                                     ccGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTcgcc    <  2:310030/151‑1 (MQ=255)
                                                                                                                                                      cGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTcgccg   >  6:1026775/1‑151 (MQ=255)
                                                                                                                                                      cGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTcgccg   >  5:216019/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: