Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,560,9960GT26.3% 587.2 / 70.0 343intergenic (+98/‑372)yffL/yffMCPZ‑55 prophage; uncharacterized protein/CPZ‑55 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (136/116);  new base T (41/49);  total (177/166)
Fisher's exact test for biased strand distribution p-value = 1.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 347 total aligned reads displayed.

AAGGGAGTGCAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGC  >  NC_000913/2560847‑2561146
                                                                                                                                                      |                                                                                                                                                      
 aGGGAGTGCAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGtt                                                                                                                                                      <  4:631105/151‑1 (MQ=255)
      gTGCAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGGTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGt                                                                                                                                                 <  7:139671/151‑1 (MQ=255)
        gCAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTg                                                                                                                                                <  7:278557/150‑1 (MQ=255)
          aaCGAATACCAGTTATTATTCAAGCCGTACAACATGAAATTACTGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTATTTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGt                                                                                                                                             <  5:1002122/151‑1 (MQ=255)
          aaCGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCGTGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGTTGGGTAAAAAGTTGCCGTTAAGGAGGGt                                                                                                                                             >  2:830524/1‑151 (MQ=255)
          aaCGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAACCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTgggg                                                                                                                                             >  3:1356144/1‑150 (MQ=255)
             gAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTAAAAAGTTGCCGTTAAGTTGGGTCa                                                                                                                                           >  5:1353425/1‑150 (MQ=255)
              aaTACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCaaa                                                                                                                                         <  8:113827/151‑1 (MQ=255)
               aTACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGATAAGGAGGGTCaaa                                                                                                                                         >  6:958849/1‑150 (MQ=255)
                 aCCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTTAGGTGGGTCAAAAGt                                                                                                                                      >  1:1325996/1‑151 (MQ=255)
                 aCCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGt                                                                                                                                      <  6:1340839/151‑1 (MQ=255)
                  acAGTTAATATTCAACCCCTACATCAACAAATAACGGATTAAGGATGGCTTATTTCCCTCATTTTATTTGCCTTAATTGCATCATTTCAGGTTGTTCAAAAATTTGCGTTACTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGtt                                                                                                                                     <  4:142430/150‑1 (MQ=12)
                    aGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTATTGGTGGGTAAAAAGTAGACTTTAAGGAGCGTCAAAAGTTGc                                                                                                                                   >  2:1271215/1‑151 (MQ=255)
                    aGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGACGTTAAGGTGGGTCAAATGTTGc                                                                                                                                   >  1:1627497/1‑151 (MQ=255)
                     gTTATTATTCAAGCCGTACAAGATGTAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGcc                                                                                                                                  <  6:85366/151‑1 (MQ=255)
                     gTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAATTTGCCGTTAAGGAGTGTAAAAATTTga                                                                                                                                   >  2:777223/1‑149 (MQ=255)
                      ttattaTTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAATTTGc                                                                                                                                   >  6:1257174/1‑149 (MQ=255)
                       tattatTCAATCCGTACAACATTAAATAACGGAATCATTCTGGCAATGTTCCCGCGTTTTCTTTGCCTTAATTCCATGAATTCAGGTTGGTCAAAAGTTGCCTTTAGTGGTGGGTCAAAAGTTTCCGTTAAGGTGGGTCAAAAGTTGCCg                                                                                                                                 <  2:1547925/150‑1 (MQ=255)
                       tattatTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTAAAAAGTTGCCTTTAAGGAGGGTCAATAGtttcagt                                                                                                                                >  5:342082/1‑146 (MQ=255)
                         ttattCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGACTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCTTTAAGGGGGGTCAAAATTAGACGtt                                                                                                                               >  3:32987/1‑150 (MQ=255)
                          tattCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTa                                                                                                                              <  8:165310/150‑1 (MQ=255)
                          tattCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTa                                                                                                                              <  8:463475/150‑1 (MQ=255)
                          tattCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGCTGGGTCAAAAGTTGCGGTTAAGGTGGCTCAAATTTTGCAGTTa                                                                                                                              >  1:268447/1‑150 (MQ=255)
                           attCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTaa                                                                                                                             <  3:420992/150‑1 (MQ=255)
                           attCAAGCAGTACAACATGAAATAACGGAATCAGTCGGGCATTCTTTCCGCCTTTTCTTTTCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGGCGTTTGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAg                                                                                                                            <  7:1049938/151‑1 (MQ=255)
                             tCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGATAGTGGTGGGTCAAATTTTGCCGGTAAGGTGGGGGAAAAGTTGGCGTTAAgcg                                                                                                                          >  8:906755/1‑149 (MQ=255)
                              cAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGt                                                                                                                          <  4:1579046/150‑1 (MQ=255)
                                 gCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCTTTAAGGTGTGTCAAAAGTTGCAGTTAAGGTggcg                                                                                                                      >  3:355418/1‑149 (MQ=255)
                                  ccGTACAACATGAAATTACGGAATCCGTCTGGCATTTTTCCCGCATTTTCTTTGCCTTTATTCCATGAATTCAGGTGGGTCAAAAGTTGGCGTTATTGGTGGGTAAATTGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTc                                                                                                                     <  7:409162/151‑1 (MQ=25)
                                   cGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGATGCCGTAAAGGTGGCTCAAAAGTTGCAGTTTACGTgtggc                                                                                                                     >  3:21767/1‑146 (MQ=255)
                                    gTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCaa                                                                                                                   <  6:1703190/151‑1 (MQ=255)
                                       caacaTGAAATACCGTAATCTGGCGTGCATTGTTCCCGGCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAATGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAg                                                                                                                <  2:222851/151‑1 (MQ=37)
                                           aTGAAATAAAGGAATCAGTCGGGCATTGTTGCCGCCTTTTCTTTTCCATAAATCAATGAATTCAGGGGGGTCAAAAGTAGGCGTTAGTGGTGGGTAAAAAGGTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGc                                                                                                            <  6:1531011/151‑1 (MQ=25)
                                             gAAATAACCGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGGCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGcc                                                                                                           <  6:1066542/150‑1 (MQ=255)
                                              aaaTAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGGCGTTAAGGTGGGTCAAAAGTTGCCg                                                                                                          >  1:1054961/1‑150 (MQ=255)
                                                  taCGTAATCAGTAGGGCATTGATGCAGCCTTATCTTTTCCTTAAATTCATGAATTCAGGGGGGTCAAAAGTTGGCTTTAGGGGTGGGTAAAAAGTTGCATTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTa                                                                                                      <  5:601097/149‑1 (MQ=21)
                                                     ggAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTGAAGGTGGGGGAAAAGTTGGAGAGTAAgc                                                                                                   >  8:1439667/1‑149 (MQ=38)
                                                     ggAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGGCAAAAGTTGCCGTTAAGCGGGGTCAAAAGTTGCGGAGTAAgc                                                                                                   >  8:1239641/1‑149 (MQ=255)
                                                     ggAACCTGGCTGGCATTGTTACTGGCTTTTGTATCCCATAATTAAATGAATTCTGGGGGGGAAAATGCAGGCGTTATGGGTGGGTCAAAAGTTGCCTTTAAGTTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGa                                                                                                   <  3:1660195/150‑1 (MQ=12)
                                                       aaTCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTGAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGGCGTTTAGGTGGGTCAAAAGTTGCCGTGTAAGGAtt                                                                                                >  6:1013072/1‑151 (MQ=255)
                                                          cctgctgccactGTTTCAGCCTTTGGTTTGCCTTAAATCCATGAATTCCTGTGGGTCAAAAGCTGCCTTTTCTGGTGGGTCAATAGATGCCGTTAAGGTGGGTCAAAAGTTGCCGTTCAGGTGGGTCAAAAGTTGCCGAGTAAGAATTcta                                                                                             <  5:1502673/140‑1 (MQ=21)
                                                           aGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTAAAAAGTTGCCTTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGGGAAAAGTTGACGAGTAACAATTcta                                                                                             >  3:1017787/1‑150 (MQ=38)
                                                                  cAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTACTGGTGGCTAAAAATTTGACTTAAAGGAGGGTCAAAAGTTGAAGTGAAGGGGTGTGAAAATTTGCCGAGGAAGAAGTGTACTATTaat                                                                                     >  3:1631313/1‑150 (MQ=255)
                                                                   acttatatcGCCTTTTGATTTCCTTTAATTAATGAATTCTGTGGGGTTTAAATTAGGCTTTACTGGTGGGTAAAATGTTGCCTTTTAGGTGGGTCAAAAGTTTCATTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGa                                                                                    <  5:238663/143‑1 (MQ=16)
                                                                        tCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTATTGGTGGGTCAAAAGTTGCATTTAAGGTGGGTCAACAGTTGAAGTTAAGGTGGGGCAAAATTTGCCGAGTAACTATTATAATATTATGAAttgt                                                                               >  8:1513023/1‑149 (MQ=37)
                                                                        tCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAATTTGCCGTTAAGGTGTGTCAAAAGTTGACGTTAAGGTGGGTCAAAAGTTGACGAGTAAGAATTCTACAATTAAAAAttat                                                                               >  7:728457/1‑149 (MQ=38)
                                                                         cccGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCt                                                                               <  1:1099895/150‑1 (MQ=255)
                                                                         cccGCCTTTGAATTGCATTAATTCCTGAATTTCTGGTGGGTAAAATGTGCCGTTTAGGGGGGGTAAAAACGTGCCGTTTAGTTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTa                                                                              <  1:91194/150‑1 (MQ=19)
                                                                          gcGCCTTATCTATTCCATTATTCATTGAATTCATTGGGGTATAAATTTCCCGTTAGTTGTGGGTAAATAGGTGCCGTTAATTTGGGTCAAAAGTTGCCCTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAg                                                                             <  8:708781/150‑1 (MQ=16)
                                                                            gCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGTGTCAAAATTTGCCGAGTAAGGATGGTACAAGTAAGATTTATAGGc                                                                           >  7:1744775/1‑151 (MQ=25)
                                                                             ccatttctttcCCATAATAACATGCATTCAGGTGGGTAATAAGTTCCCCTTAGTGTTGGGTCAAAAGTTGCCCTTAAGCTGGGGCAAAAGTTGCCGTTAAGGAGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCg                                                                          <  4:338890/148‑1 (MQ=255)
                                                                             ccTTTTCTTTTCCTTAATTCATTGATTTCTGTGGTGTTAAAATTTGCCGTTTTTTGTGGGTAATTTGTTGCCGTTTAGTTGGGTCAAATGTTTCCGTTAAGGTGGGTCAAAAGTTGGCGAGTAAGAATTCTACTATTAAGAATTCTAGGCg                                                                          <  5:1186774/151‑1 (MQ=255)
                                                                              cTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAATGTTGGCTCAAAATTTGCCTTAAAGGTGGGTCAAAAGTTGACGAAAACCTGGCTCTAAAATAGCCTAGTAAGCATTCTAAAATAAACAATACTAAGcc                                                                          >  2:9776/1‑149 (MQ=16)
                                                                               ttttCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAATGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCGGCGTAAGAATTCTACTATTAAGAATTCTAGGCGAc                                                                        >  3:573822/1‑151 (MQ=255)
                                                                               ttttCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGAGGTGGGTCAAATGTTGCCGTTAAGGTGGGTCAAAAGTTGCGGTTAAGGTGGGTCAAATGTTGCCGAGTAAGACTTCTACTATTAACAATTCTACGCGAc                                                                        >  1:169636/1‑151 (MQ=37)
                                                                                ttttttTGCCTTAATTACATGAATTCAGTTGGGTCTATTTTTGGCTTTAGTGTTGTGTAAATAGTTGCCTTTTATGTGGGTTAAATGTTTCATTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGAc                                                                        <  2:1044148/150‑1 (MQ=255)
                                                                                tttCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTTAGGTGGGTGAAAAGTTGCCGAGTAAGAATTCTACTATTACTAATTCTTGGCGACg                                                                       >  1:1558495/1‑151 (MQ=255)
                                                                                  tCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGTTGGGTCAAAAGTTGCCGTTAAGGAGGGTCAAAAGTTGAAGATAAGGTGTGTCAAAATTTGACGAGAAACAATTCTAAAATTAAGTATTATAGGCGaag                                                                       >  2:377035/1‑147 (MQ=25)
                                                                                     ttGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAACTTGCCTTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGa                                                                    <  8:286634/149‑1 (MQ=255)
                                                                                       gCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGACGTTAAGGTGGGTCAAAAGTTGCCGTGAAGGTGGGTCAAAAGTTGCCGAGTAACAATTCTACTATTAAAAATAGTAGGCGACGGGAGttt                                                                >  6:531702/1‑151 (MQ=37)
                                                                                         cTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAATGGTGGGTCAAAAGTTGCCTTTAAGGTGGGTCAAAAGTTGACGATAAGGTGGCGCAAAATTTGAAGAGTAACAATAATAAAATAAATAATACTAGTCGaagagacttatt                                                              >  7:983875/1‑140 (MQ=12)
                                                                                          ttAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTAAAAATTTGCATTTAAGGAGGGTCAAAATTTTAAGTAAACGGGGCTCAAATTTTGACTAGTAAGAATTCTACAATTAATATTAATACGCGACGTGAGttttat                                                             >  4:284017/1‑149 (MQ=255)
                                                                                           ttaattaaTGAATTCTGTTTGGTTAAAATTTGGCTTTATTGGTGGGTAAATTGCTGCCTTTTATTTGGGTTAAAAGTTGCCTTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGttttt                                                              <  4:1368171/142‑1 (MQ=255)
                                                                                           tAAATCCATCAATTCATGGTGCTCAAAAGTAGGCGTTACGGGTGGGTCAAACGGCGCCTTTTCGGTGGGTCAAAAGGTGCCCTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGtttttt                                                             <  3:798669/150‑1 (MQ=25)
                                                                                            aaTTCCATGAATTCAGGTTGGTAAAAAGTTGGCTTTATTGGTGGGTCAAATGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTa                                                            <  7:209046/150‑1 (MQ=37)
                                                                                            aaTTCCATGAATTCAGGTGGCTCAAAAGCTGCCGTTAGTGGTGGGTAAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTa                                                            <  8:792490/150‑1 (MQ=255)
                                                                                             aTTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTGAAGGTGGGTCAAAAGTTGACGCGTAAGAATTCTACAATTATGAATTCTAGGCGACGGGAGTTTTTTAGc                                                          >  1:1759276/1‑151 (MQ=255)
                                                                                               tCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGcgt                                                        >  6:1494852/1‑149 (MQ=255)
                                                                                               tCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTAAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGcc                                                         <  1:809338/150‑1 (MQ=255)
                                                                                                 cATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCTAAAGTTGCCGAGTAAGCATTCTAATATTAAGAATTCTAGGCGAAGGGAGTTTTTTAGCCttc                                                      >  3:1088393/1‑149 (MQ=255)
                                                                                                     aaTTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCAGTTAAGGTGGGTGTAAATTTGACGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGTGAGTTTTTTAGACTACcggg                                                  >  8:722660/1‑150 (MQ=37)
                                                                                                      aTTCAGGTGGGTCAAAAGTTGCCGTTAGTTGTGGGTCATAAGGTGCCTTTAAGGTGGGTCAAAAGTTGCAGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCcggtc                                                 <  3:992930/151‑1 (MQ=255)
                                                                                                       attccgggtggTAAATAGTTGGCGTTATTGGTGGGTAAATAGTTGCCTTTTAGGTGGGTAAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCcggtc                                                 <  4:186363/142‑1 (MQ=37)
                                                                                                           aGGTGGGTAAAAAGGTGGCGTTAGTGGTGGGTCAAAACTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCggtcggtc                                             <  6:644159/151‑1 (MQ=255)
                                                                                                             gTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAATTTGACGTGAAGGGGGGTCTAATGTTGCCGAGTAAGAATTCTACTATTAAGATTTATAGGCGACGGGAGTTTTTTAGTCTCCCGGTGGGTCa                                            >  1:1189661/1‑150 (MQ=37)
                                                                                                              tGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCAt                                           >  1:1376319/1‑150 (MQ=255)
                                                                                                                      aaGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCg                                  >  4:487918/1‑151 (MQ=255)
                                                                                                                       agaagCCGTTAGTGGTGGGTCAAAACTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGc                                 <  8:13018/147‑1 (MQ=255)
                                                                                                                       aGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGc                                 >  3:96497/1‑151 (MQ=255)
                                                                                                                         ttGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCt                               >  6:1035812/1‑151 (MQ=255)
                                                                                                                         ttGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGCTCAAAAGATCCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCACCCGGTCGGTCATAGAGAGTCGCCt                               >  4:753154/1‑151 (MQ=255)
                                                                                                                           ggcGTTAGTGGTGGGTAAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTa                              <  3:1424754/148‑1 (MQ=255)
                                                                                                                            ggtttATTGGTGGGTATATTGTTGCCTTTTATGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTaga                            <  5:1775412/148‑1 (MQ=38)
                                                                                                                            ccGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCt                               >  4:1046477/1‑148 (MQ=255)
                                                                                                                            ccGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGACt                               >  7:1084984/1‑148 (MQ=255)
                                                                                                                            ccGTTAGTGGGGGGTCAAATGGGGCCTTTAAGGTGGGTCAAAAGTTCCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTaga                            <  3:669132/151‑1 (MQ=39)
                                                                                                                             cGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGACTCCCGGTCGGTCATAGAGAGTCGCCTag                             >  4:236149/1‑149 (MQ=255)
                                                                                                                               ttATTGGTGGGTAAAATGTTGCCTTTAAGGTGGGTCAAAAGTTTCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGt                          <  3:1216695/150‑1 (MQ=255)
                                                                                                                                 aGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTa                        >  2:551393/1‑150 (MQ=255)
                                                                                                                                  ggggtggCTCAAATCTTGCCTTTAAGTTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTaaa                      <  3:1788843/149‑1 (MQ=38)
                                                                                                                                         gTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACaga                  >  6:227697/1‑148 (MQ=255)
                                                                                                                                            aaaaGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAatt            >  7:1590323/1‑151 (MQ=255)
                                                                                                                                              aactggCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAattat          <  1:708981/146‑1 (MQ=255)
                                                                                                                                              aaGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAattat          <  1:511080/151‑1 (MQ=255)
                                                                                                                                                gTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTAttt        >  7:1431167/1‑151 (MQ=255)
                                                                                                                                                  tGCCTTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTc       <  5:738692/150‑1 (MQ=255)
                                                                                                                                                  tGCAGTTAAGGTGGGTAAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTc       <  5:1791060/150‑1 (MQ=255)
                                                                                                                                                   gCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCggtcggtc                                             <  8:64524/111‑1 (MQ=255)
                                                                                                                                                    ccGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTcgcc    >  4:1617257/1‑151 (MQ=255)
                                                                                                                                                      gTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTcgccgc  >  8:320674/1‑151 (MQ=255)
                                                                                                                                                      |                                                                                                                                                      
AAGGGAGTGCAACGAATACCAGTTATTATTCAAGCCGTACAACATGAAATAACGGAATCAGGCGGGCAATGTTCCCGCCTTTTCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGC  >  NC_000913/2560847‑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: