Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,560,9960GT27.2% 584.0 / 68.4 332intergenic (+98/‑372)yffL/yffMCPZ‑55 prophage; uncharacterized protein/CPZ‑55 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (108/130);  new base T (32/57);  total (142/190)
Fisher's exact test for biased strand distribution p-value = 1.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 339 total aligned reads displayed.

AAAGGGAGTGCAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGC  >  NC_000913/2560846‑2561146
                                                                                                                                                      |                                                                                                                                                      
    ggAGTGCAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTaa                                                                                                                                                    <  4:242914/151‑1 (MQ=255)
    ggAGTGCAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTaa                                                                                                                                                    <  6:112618/151‑1 (MQ=255)
    ggAGTGCAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTACTGGTGGGTAAAATTTTGCCTTTa                                                                                                                                                     >  3:1500769/1‑150 (MQ=255)
       gggCAACGAATACCAGGTATTATTCAAGCCGTACAACATGAAATAACGGACTCAGGCGGGCACTGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGCTGCCGTTAAGGt                                                                                                                                                 <  4:445347/149‑1 (MQ=255)
        tGCAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAATGATGTGTAAAAATTTGCATTTAAggat                                                                                                                                                >  7:1059920/1‑149 (MQ=255)
             cGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTAAAAATTTGCCGTTTAGGAGGGTc                                                                                                                                            >  3:1161108/1‑150 (MQ=255)
              gAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGGGGGTCa                                                                                                                                           >  3:710745/1‑150 (MQ=255)
               aaTACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCATAAGTTGCCTTGAAGGTGGGTCaa                                                                                                                                          >  2:526596/1‑150 (MQ=255)
               aaTACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCaaa                                                                                                                                         <  3:1302457/151‑1 (MQ=255)
                aTACCAGTTATTATTCAAGCCGTACAACATGAAATACCGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGGCAAAAGTTGCCGTTAGTGGTGGGGCAAAAGTTGCCGTTAAGGTGGGTCaaaa                                                                                                                                        <  8:1505581/151‑1 (MQ=255)
                aTACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCaaaa                                                                                                                                        <  6:610298/151‑1 (MQ=255)
                  aCCAGTTATTATTCAAGCAGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGt                                                                                                                                      <  4:1331701/151‑1 (MQ=255)
                   ccAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAATTTGCCTGTAAGGTTGTTCAAAAGtt                                                                                                                                     >  3:581608/1‑151 (MQ=255)
                   ccAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGtt                                                                                                                                     <  7:1177877/151‑1 (MQ=255)
                    cAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGCTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTg                                                                                                                                    <  7:418137/151‑1 (MQ=255)
                     aGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGTTGGGTAAAAAGTTGCCGTTAAGGAGGGTCAAAATTTg                                                                                                                                    >  1:1339402/1‑150 (MQ=255)
                     aGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCACGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAATGTTGTCTAAAAATTTGACATAAAGGAGAGTCAAAATTTg                                                                                                                                    >  4:205669/1‑150 (MQ=255)
                        tattatTCAAGCCGTACAAGATGAAATACAGGAATCAGGAGGGCAATGTTCCCGCATTTTCTTTGACTTAATTGCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTTCCTTTAAGGTGGGTCAAAAGTTGCCg                                                                                                                                 <  1:856516/150‑1 (MQ=255)
                        tattatTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCg                                                                                                                                 <  1:1023212/150‑1 (MQ=255)
                        tattatTCAAGCCGTACAACATGAAATAACGGAATCAGCCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCg                                                                                                                                 <  1:694674/150‑1 (MQ=255)
                         attattCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGtt                                                                                                                               <  1:155428/151‑1 (MQ=255)
                               cAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGGTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGTCTCAAAAGTTGCCGTGAAGGTg                                                                                                                         >  2:942475/1‑151 (MQ=255)
                                aaGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAATTTGCCTTTAAGGAGGGTAAAAATTTTACGTTTAGGTgg                                                                                                                        >  1:1446774/1‑151 (MQ=255)
                                aaGCCGTACAACATGAAATAACGGAATCAGGAGTGCAATGTTGCATCCTTTTCTTTGCCTTAATTGCATGAATTCATGTGGGTGAAAAGTTGGCGTTATTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTgg                                                                                                                        <  3:346893/151‑1 (MQ=37)
                                   ccGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGAAGTTAAGGTGGGtg                                                                                                                     >  2:333156/1‑150 (MQ=255)
                                   ccGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCCTTAAGGTGGGTCAAAAGTTGCCGTTAAgggggggc                                                                                                                     >  6:253949/1‑149 (MQ=255)
                                   ccGTACAACATGAAATAAAGTAATCAGTCGGGCAATTTTCCCGCCTTTTCTTTGCCATAAATCCATGAATTCAGGTGGGTCAAAAGCTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTc                                                                                                                     <  8:1004706/151‑1 (MQ=37)
                                    cGTACAACATGAAATAACGGAATCAGGCGGGCACTGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCa                                                                                                                    <  4:115061/151‑1 (MQ=255)
                                    cGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCa                                                                                                                    <  5:468093/151‑1 (MQ=255)
                                     gTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAATGTTGCCGTGAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGGCa                                                                                                                    >  2:417483/1‑150 (MQ=255)
                                     gTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAATGTTGCCGTGAAGGTGGGGCTATAGTTTCCGGTTCGGTGGGTct                                                                                                                    >  2:743488/1‑149 (MQ=255)
                                       acaacaTGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTAAATATTTGCCGTTAAGGTgggggcaat                                                                                                                 >  5:558561/1‑145 (MQ=255)
                                        gaacaTGAAATAATGGAATCCGGCGGGCAATGTTGCCGCCTTTTCTTTGCCTTTATTCCATGAATTCAGGTGGGTCAAAAGTTGCGGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCaaaa                                                                                                                 <  8:891769/149‑1 (MQ=255)
                                         aacaTGAAATAACGGAATCAGGCGGGCAATGTTCCCGGCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTAAAAAGTTGGCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTTCCGTTAAGGTGGGTCAAAAGt                                                                                                               <  3:1036788/151‑1 (MQ=255)
                                         aacaTGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGATATTGTTTGGTAAAAAGTTGAATTTAAGGAGGGTAAAAATTTAAATTTAAAGTGTGTCaaaatt                                                                                                               >  7:100696/1‑149 (MQ=255)
                                          acaTGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTGAAGGTGGGTCAAATGTTGGCGTTAAGGTGGGTCACAAGtt                                                                                                              >  5:301315/1‑151 (MQ=255)
                                           caTGAAATAACTGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTg                                                                                                             <  6:1503062/151‑1 (MQ=255)
                                                 aTAACGGAATCAGCCGGGCAATTTTCCAGCCTTTTCTTTGCCATAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCCTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGt                                                                                                       <  1:516700/151‑1 (MQ=255)
                                                   aaCGGACTCCGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGGGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTa                                                                                                      <  3:201882/150‑1 (MQ=255)
                                                    atggAATCAGTGGGGCATTTTTGCAGGCTTTTGTTTTCCTTAATTTAATGAATTCTTTTGGGTTAAATGTAGCGGTTATTGGTGGGTAAATTGTTGCCTTTTAGGTGGGTCAAAAGTTTCCGTTAAGGTGGGTCAAAAGTTGCCGAGTaa                                                                                                     <  7:209809/148‑1 (MQ=12)
                                                    aagtaaTCAGACGGGCAATTTTCCAGCCTTTTGTTTGCCTTAAATCCATGAATTCAGGTGGGTCAAAAGTAGCCCTTAGTGGTGGGTCAAAAGGGGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAg                                                                                                    <  6:40964/147‑1 (MQ=37)
                                                      ggAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGGCAAAAGTTGCGGAGTa                                                                                                      >  5:1533222/1‑147 (MQ=255)
                                                      ggAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTGAAGGTGGGTCAAATGTTGCCGAGTAAgt                                                                                                   >  5:202999/1‑149 (MQ=255)
                                                      ggAATCACTCTGGCAATGTTCCCTCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTAAAAAGTTGCCCTTAGTGTTGGGTAAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGGCGAGTAAGaa                                                                                                  <  5:1266652/151‑1 (MQ=37)
                                                        aaTCCGTCGCGCAATTATCCCGCCTTTTCTTTTCCTTAATTCCATGAATTCAGGTGGGTCAAAAGATGCCGTTAGTGGTGGGTCAAAAGTTGCCTTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAAtt                                                                                                <  4:1498218/151‑1 (MQ=37)
                                                        aaTCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGGGGGTCAAAAGTTGCCGAGTAAGAAtt                                                                                                >  3:1414015/1‑151 (MQ=255)
                                                         aTCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGGGGGTCAAAAGTTGGCGTGAAGGTGGGGGACAATTTGCCGCGTAAggagtt                                                                                               >  6:1062049/1‑146 (MQ=37)
                                                         aTCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCAGTTAAGGAGGGTCAAAATTTGCTTTTAACGTGTCTCAAAATTTGACGAGAAAGAAtt                                                                                                >  1:199659/1‑150 (MQ=255)
                                                           ccggCGGGCAATGTTCCCGCCGTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCCTTAGGGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTcta                                                                                             <  5:1175928/149‑1 (MQ=255)
                                                           cAGGCGGGCAATGTTCCCGCCTTTGGTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTcta                                                                                             <  5:1403762/151‑1 (MQ=255)
                                                           aagtcGGGCATTGTTCCCGCCTTTTGTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGGCTTTTGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTct                                                                                              <  3:539781/146‑1 (MQ=38)
                                                             ggCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATCCATGAATTCAGGTGGGTCAAAAGTTTGCTTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTct                                                                                              <  5:1028494/147‑1 (MQ=255)
                                                               ggggCAATGTTCCAGCCTTTTCTTTTCCTTAATTTCATGAATTCATGTGTGTCAAAAGTTGCCCTTAGTGGTGGGTCAAAACTTGCCTTTAAGGTGGGTCAAAAGGTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTAt                                                                                         <  2:217905/150‑1 (MQ=37)
                                                                    aaTGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAATAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTa                                                                                       >  3:478984/1‑148 (MQ=255)
                                                                         tCCCTCCTTTTCTTTTCCTTAATTCCATGAATTCAGGTTCGTCAAAAGATGCCTTTACTGGTGGGTCAAATGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGaa                                                                                   <  2:862959/147‑1 (MQ=37)
                                                                          cccGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTTAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTa                                                                              <  6:933280/151‑1 (MQ=255)
                                                                          cccGCCTTTAGTTTTCCTTTATTTCATGAATTCAGGTGGGATAAAAGTAGGCGTTACTGGTGGGTAAAATGTTGCCGTTAAGGTGGGTCAAATGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCt                                                                               <  8:398713/150‑1 (MQ=25)
                                                                          cccGCCTTATGTTTTCCTTAAATCCATGAACTCCGGGGGCTTAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTTAGGTGGGTCAAAAGTTTCCTTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCt                                                                               <  1:109648/150‑1 (MQ=25)
                                                                          cccGCCTTATCTTTGCCTTTATTCCATGAATTCAGGGCGGTCAAAAGGTGCCGTTAGTGGTGGCTCACAAGTTGCCGTTAAGGGGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCt                                                                               <  6:32233/150‑1 (MQ=37)
                                                                             gCCTTATCTTTGCCTTAAATCCATGAATTCAGGTGGGTCAAAAGTTGCCTTTACTGGGGGGTAAAAAGTTGCCTTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGc                                                                           <  2:129873/151‑1 (MQ=38)
                                                                              gcTTTTGTTTGGCTTATTTCCATGTATTCTGGTTGGTTAATAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCg                                                                          <  1:579017/150‑1 (MQ=25)
                                                                                ttttCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGGCGTTAAGGTGGGGCAAAAGTAGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGAc                                                                        >  6:146715/1‑151 (MQ=255)
                                                                                ttttCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCAGAGTAAGAATTCTACAATTAACAATTCTACGCGAc                                                                        >  1:957349/1‑151 (MQ=255)
                                                                                 ttggattGCCTTAAATCCATGAATTCCGGCGGGTCAAAAGTTGCCGTTAGTGCGGGGTCAAAAGTTGCCGTTACGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGAc                                                                        <  2:61393/145‑1 (MQ=37)
                                                                                 ttgcTTTCCCTTTATTCCATGAATTCCGGGGGGACAAAAGTTGCCCTTAGTGGTGGGTCAAAACTTGCCGTTAAGGTGGGTCAAATGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGAc                                                                        <  7:154954/147‑1 (MQ=37)
                                                                                  ttgtttGCCATAAATCCATGAATTCAGGGGGGTAAAAAGGTGCCCTTAGTGGTGGGTCAAAAGATGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACg                                                                       <  3:94852/147‑1 (MQ=37)
                                                                                  ttCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCGAAAGTTGCCGAGTAAGAATTCTGCTATTAAGAATTCTAGGCGACgg                                                                      >  5:245738/1‑151 (MQ=255)
                                                                                      aggCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGGGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAg                                                                   <  3:147079/148‑1 (MQ=255)
                                                                                         ccTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAATTTGACGTTAAGGTGGCTCAAAAGTTGACTAGTAAGAATTATAATATTAAGAATTCTAGGCGACGGGATtttt                                                               >  4:714042/1‑151 (MQ=37)
                                                                                             aaTTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAATTTGCCGAGTAAGAATTGTACTATTAAGATTTCTAGGCGACGGCAGtttttt                                                             >  7:914224/1‑149 (MQ=255)
                                                                                               tgcgatataTTCATGTTGGTAAAACGTTTCGCTTAATTGTTGGTTAATATTTTCCCTTAAGGTGGGTCAAAAGTTGCAGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGc                                                          <  1:1209527/142‑1 (MQ=21)
                                                                                                       aTTCCGTTGCGTCAAAATTTGCCCTTCCTGGGGGGTAAAAAGTTGCCTTTAAGGTGGGTAAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCcggtc                                                 <  7:1433343/151‑1 (MQ=37)
                                                                                                          cAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGCATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCcggtcggt                                              >  5:1093709/1‑151 (MQ=255)
                                                                                                          cAGCTGGCTAAAAAGTTGCCTTTAGTGGGGGGTCAACAGTTGCCTTTAAGGTGGGTCAAATGTTTCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCcggtcgg                                               <  7:59436/150‑1 (MQ=38)
                                                                                                                ggtttaaaTTTGGCTTTACTTGTGGGTAATTAGTTGCCTTTACGGTGGGTAAAATGTTGCCTTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATaga                                        <  1:256987/146‑1 (MQ=255)
                                                                                                                  tCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATag                                         >  2:675721/1‑148 (MQ=255)
                                                                                                                    taaatttcgggTTATTGGTGGGTAATTTGTTGCCGTTTATTTGGGTAAATAGTTGCCTTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATagagag                                     <  8:166717/140‑1 (MQ=255)
                                                                                                                      aagatccccttAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCg                                  <  5:1589904/142‑1 (MQ=255)
                                                                                                                            cccttccgggtggCTAAAATGGTGCCGTTTCGGTGGGTCAAAAGCTGCCCTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTaga                            <  2:1587568/143‑1 (MQ=255)
                                                                                                                             cGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTag                             >  7:698120/1‑149 (MQ=255)
                                                                                                                             cGTTAGTGGGGGGTCAAAAGTTGCCGTTACGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTagag                           <  6:222380/151‑1 (MQ=255)
                                                                                                                             cGTTAGTGGGGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTagag                           <  5:780538/151‑1 (MQ=255)
                                                                                                                              gTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGt                          >  6:1340538/1‑151 (MQ=255)
                                                                                                                               ttAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGtt                         >  5:435033/1‑151 (MQ=255)
                                                                                                                                 actggtggGTAATAAGTTGCCTTTTAGGTGGGTCAATTGTTTCATTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTaa                       <  4:883852/149‑1 (MQ=255)
                                                                                                                                    ggtggGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAAc                     >  7:1421933/1‑150 (MQ=255)
                                                                                                                                     gttggTCTAATGTTGCCGTTACGTTGGGTCAAAAGTTTCAGTTCAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACag                   <  7:322962/148‑1 (MQ=255)
                                                                                                                                        gttaaattgTTGCCTTTTAGTTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACagaag                <  8:1082029/143‑1 (MQ=39)
                                                                                                                                          tCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCGTAGAGTTAAACagaaga               >  1:349755/1‑150 (MQ=255)
                                                                                                                                            aaaaGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAatt            >  5:768368/1‑151 (MQ=255)
                                                                                                                                             ataGTTGCCTTTTAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATATAGAGTCGCCTAGAGTTAAACAGAAGAatta           <  8:717705/149‑1 (MQ=255)
                                                                                                                                             aaaGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAatt            >  2:825533/1‑150 (MQ=255)
                                                                                                                                               aGTTGCCCTTAAGTTGCGTCAAAAGGTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCAGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAattatt         <  2:1028584/151‑1 (MQ=255)
                                                                                                                                                gatGCCTTTAAGGTGGGTCAAAAGTTGCCGTTCAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTAttt        <  5:850328/149‑1 (MQ=255)
                                                                                                                                                  tGCCTTTTAGTTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTc       <  6:883710/150‑1 (MQ=255)
                                                                                                                                                  tGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTcg      >  3:293283/1‑151 (MQ=255)
                                                                                                                                                    ccGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTcgcc    >  8:1399826/1‑151 (MQ=255)
                                                                                                                                                     cGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCGCGGTCGGTCATAGATAGTAGCCTAGAGTTAAACAGAAGAATTATTTcgccg   >  8:1592301/1‑151 (MQ=255)
                                                                                                                                                      gTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTcgccgc  >  7:740762/1‑151 (MQ=255)
                                                                                                                                                      gTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTcgccgc  >  6:722123/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 < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: