Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,987,103 C→T 5.3% intergenic (+27/‑433) yqeG → / → yqeH putative transporter/putative LuxR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,987,1030CT 5.3% 911.6 / 25.7 264intergenic (+27/‑433)yqeG/yqeHputative transporter/putative LuxR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (115/135);  new base T (4/10);  total (119/145)
Fisher's exact test for biased strand distribution p-value = 2.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.83e-01

Only 100 of 265 total aligned reads displayed.

TTTTCATCATGCCTACGCTGTCAACGTATCTCATCCCGGCGCTTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATCTTAAATAAGACTATTTAAGATGCATAACTT  >  NC_000913/2986955‑2987252
                                                                                                                                                    |                                                                                                                                                     
  ttCATCATGCCTACGCTGTCAACGTATCTCATCCCGGCGCTTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTAc                                                                                                                                                    >  2:4585527/1‑150 (MQ=255)
     atcatGCCTACGCTGTCAACGTATCTCATCCCGGCGCTTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAAt                                                                                                                                                 <  1:2778946/150‑1 (MQ=255)
      tcatGCCTACGCTGTCAACGTATCTCATCCCGGCGCTTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAAtt                                                                                                                                                >  1:4267051/1‑150 (MQ=255)
      tcatGCCTACGCCGTCAACGTATCTCATCCCGGCGCTTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAAtt                                                                                                                                                >  2:2579157/1‑150 (MQ=255)
        atGCCTACGCTGTCAACGTATCTCATCCCGGCGCTTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAAtttt                                                                                                                                              >  1:301566/1‑150 (MQ=255)
          gCCTACGCTGTCAACGTATCTCATCCCGGCGCTTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGAATGGTTCAAAATTGTAGCAATTTTGa                                                                                                                                            <  2:4364174/150‑1 (MQ=255)
          gCCTACGCTGTCAACGTATCTCATCCCGGCGCTTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGAATGGTTCAAAATTGTAGCAATTTTGa                                                                                                                                            <  2:1356768/150‑1 (MQ=255)
           ccTACGCTGTCAACGTATCTCATCCCGGCGCTTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGaa                                                                                                                                           <  2:4644194/150‑1 (MQ=255)
           ccTACGCTGTCAACGTATCTCATCCCGGCGCTTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGaa                                                                                                                                           <  1:1820544/150‑1 (MQ=255)
            cTACGCTGTCAACGTATCTCATCCCGGCGCTTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAAc                                                                                                                                          >  1:2968573/1‑150 (MQ=255)
             tACGCTGTCAACGTATCTCATCCCGGCGCTTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAAcc                                                                                                                                         >  2:2542647/1‑150 (MQ=255)
                 cTGTCAACGTATCTCATCCCGGCGCTTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAAccattc                                                                                                                                     <  1:1961876/150‑1 (MQ=255)
                   gTCAACGTATCTCATCCCGGCGCTTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAAccattcca                                                                                                                                   <  2:4685801/150‑1 (MQ=255)
                   gTCAACGTATCTCATCCCGGCGCTTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAAccattcca                                                                                                                                   <  1:2282939/150‑1 (MQ=255)
                   gTCAACGTATCTCATCCCGGCGCTTAAACCCTGGCGTTCCATCGGAAATCTGATGACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAAccattcca                                                                                                                                   <  2:4130028/150‑1 (MQ=255)
                    acAACTTATCTCATCCCGGCGCTTAAACCCTGGCGTTCCATCGGAAATCTTAATACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAAccattccat                                                                                                                                  <  1:4416769/149‑1 (MQ=255)
                     cAACGTATCTCATCCCGGCGCTTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAAccattccatt                                                                                                                                 <  2:1340377/150‑1 (MQ=255)
                     cAACGTATCTCATCCCGGCGCTTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAAccattccatt                                                                                                                                 <  1:3457421/150‑1 (MQ=255)
                        cGTATCTCATCCCGGCGCTTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTTCTTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATTTTTCAAAATTGCTACAATTTTGAACCATTCCAttttt                                                                                                                              <  1:2582765/150‑1 (MQ=255)
                           aTCTCATCCCGGCGCTTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCat                                                                                                                           <  2:4496814/150‑1 (MQ=255)
                           aTCTCATCCCGGCGCTTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCat                                                                                                                           <  2:1273503/150‑1 (MQ=255)
                            tctcATCCCGGCGCTTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCata                                                                                                                          >  2:1891093/1‑150 (MQ=255)
                                aTCCCGGCGCTTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAatat                                                                                                                      >  2:393320/1‑150 (MQ=255)
                                   ccGGCGCTTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCa                                                                                                                   <  2:674314/150‑1 (MQ=255)
                                   ccGGCGCTTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCa                                                                                                                   <  2:3022525/150‑1 (MQ=255)
                                   ccGGCGCTTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCa                                                                                                                   <  2:4000120/150‑1 (MQ=255)
                                      gcgcTTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCAttt                                                                                                                <  2:4640294/150‑1 (MQ=255)
                                      gcgcTTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCAttt                                                                                                                <  1:3106132/150‑1 (MQ=255)
                                         ctttatcccTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCTGTAATGTTCTTTAGCTAATATCAGTTGTTGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGtt                                                                                                             <  1:2371946/144‑1 (MQ=255)
                                           tAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGttat                                                                                                           <  1:2542647/150‑1 (MQ=255)
                                               cccTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTattatt                                                                                                       >  2:2216367/1‑150 (MQ=255)
                                                ccTGGCTGTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATTGTTCAAAATTGCTACAATTTTGTACCATTCCATTTTTCATAATATTCATTTGTTATTATTc                                                                                                      <  1:1541685/150‑1 (MQ=255)
                                                ccTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTATGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTc                                                                                                      >  2:857262/1‑150 (MQ=255)
                                                ccTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATACCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTc                                                                                                      >  1:2503631/1‑150 (MQ=255)
                                                 cTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTcc                                                                                                     <  1:1859245/150‑1 (MQ=255)
                                                  tGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCa                                                                                                    >  2:3492488/1‑150 (MQ=255)
                                                  tGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCa                                                                                                    >  2:3371103/1‑150 (MQ=255)
                                                  tGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCa                                                                                                    >  2:3063982/1‑150 (MQ=255)
                                                    gCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCAcc                                                                                                  <  2:2293656/150‑1 (MQ=255)
                                                        tCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTAt                                                                                              <  1:4666018/150‑1 (MQ=255)
                                                        tCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTAt                                                                                              <  1:518031/150‑1 (MQ=255)
                                                          cATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCa                                                                                            <  1:136048/150‑1 (MQ=255)
                                                          cATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCa                                                                                            >  2:1482482/1‑150 (MQ=255)
                                                             cGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACtt                                                                                         <  2:2712802/150‑1 (MQ=255)
                                                             cGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACtt                                                                                         >  1:4129227/1‑150 (MQ=255)
                                                             cGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACtt                                                                                         <  2:1756908/150‑1 (MQ=255)
                                                               gAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATCTCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTtct                                                                                       >  2:1795151/1‑150 (MQ=255)
                                                               gAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTtct                                                                                       <  2:206229/150‑1 (MQ=255)
                                                                  aTCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTtctctc                                                                                    >  2:1436265/1‑150 (MQ=255)
                                                                   tCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCg                                                                                   >  1:2239418/1‑150 (MQ=255)
                                                                     tGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTc                                                                                 <  1:861767/150‑1 (MQ=255)
                                                                      gATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCa                                                                                <  1:903212/150‑1 (MQ=255)
                                                                          aCGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATc                                                                            <  2:2428889/150‑1 (MQ=255)
                                                                            gCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCt                                                                          >  2:3859442/1‑150 (MQ=255)
                                                                            gCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCt                                                                          >  1:2881574/1‑150 (MQ=255)
                                                                            gCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCt                                                                          >  2:3753782/1‑150 (MQ=255)
                                                                            gCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCt                                                                          >  2:589967/1‑150 (MQ=255)
                                                                               gATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGt                                                                       <  2:1963035/150‑1 (MQ=255)
                                                                                  cGTGGGTATCCTGTGCGTATCGGTTATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATTATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGtatt                                                                    <  2:1655888/150‑1 (MQ=255)
                                                                                      ggTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCa                                                                <  1:744353/150‑1 (MQ=255)
                                                                                       gTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCag                                                               >  1:4048967/1‑150 (MQ=255)
                                                                                       gTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCag                                                               >  1:4049015/1‑150 (MQ=255)
                                                                                         aTCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCagtt                                                             >  1:3989458/1‑150 (MQ=255)
                                                                                         aTCCTGTGCGTATCGGAAATGTTCTTTAGCTAATATCAGTTGTGGATGCTACAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCagtt                                                             <  2:3469036/150‑1 (MQ=255)
                                                                                          tCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCagtta                                                            <  2:1397106/150‑1 (MQ=255)
                                                                                            cTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTATCTCGTCATATCCTAGTATTATCagttagt                                                          >  1:4250973/1‑150 (MQ=255)
                                                                                             tgtgCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGCTGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCagttagtt                                                         <  2:3845554/150‑1 (MQ=255)
                                                                                              gtgCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAgttagtta                                                        <  2:4408078/150‑1 (MQ=255)
                                                                                                gCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTaaa                                                      <  2:2909547/150‑1 (MQ=255)
                                                                                                gCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTaaa                                                      <  1:1572442/150‑1 (MQ=255)
                                                                                                gCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTaaa                                                      <  1:145656/150‑1 (MQ=255)
                                                                                                 cGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAac                                                     <  2:4207866/150‑1 (MQ=255)
                                                                                                  gTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAaca                                                    <  1:3591534/150‑1 (MQ=255)
                                                                                                     tCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAacacac                                                 <  2:961679/150‑1 (MQ=255)
                                                                                                     tCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAacacac                                                 <  1:2147666/150‑1 (MQ=255)
                                                                                                     tCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAacacac                                                 <  1:3217903/150‑1 (MQ=255)
                                                                                                       ggTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATCTTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACaa                                               >  2:111046/1‑150 (MQ=255)
                                                                                                       ggTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACaa                                               >  2:708126/1‑150 (MQ=255)
                                                                                                       ggTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTACAAAATTTCTACAATTTTGAACCATTCCATTTTTAATAATAATCATTTGTTATTATTCAACCCTAAAACTAATCTCGTCATATCCTAGTATTATAAGTTAGTTAAACACACaa                                               >  1:717284/1‑150 (MQ=255)
                                                                                                        gTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACaaa                                              >  2:2636063/1‑150 (MQ=255)
                                                                                                            tGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAgaga                                          <  2:4609602/150‑1 (MQ=255)
                                                                                                                 tttAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGa                                     <  1:996512/150‑1 (MQ=255)
                                                                                                                  ttAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGaa                                    <  1:3342801/150‑1 (MQ=255)
                                                                                                                      cTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATc                                >  2:849736/1‑150 (MQ=255)
                                                                                                                       tAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATCt                               <  1:3735075/150‑1 (MQ=255)
                                                                                                                           atCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATCTTaaa                           <  1:1436265/150‑1 (MQ=255)
                                                                                                                                ttGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATCTTAAATAAGa                      <  2:214683/150‑1 (MQ=255)
                                                                                                                                aatggaaTGGTTCAAAATTGTAGCAATTTTGAACCATCCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATCTTAAATAAGa                      >  1:3153461/7‑150 (MQ=39)
                                                                                                                                      aTGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATCTTAAATAAGACTAttt                <  2:360278/150‑1 (MQ=255)
                                                                                                                                      aTGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATCTTAAATAAGACTAttt                >  2:3324514/1‑150 (MQ=255)
                                                                                                                                       tGGTTCAAAATTGTAGCAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATCTTAAATAAGACTATTTa               >  1:1069895/1‑150 (MQ=255)
                                                                                                                                        ggTTCAAAATTGCTACAAATTTGAACCATTCCATTTTTCAAAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTAATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATATGAACATCTTAAATAAGACTATTTaa              >  1:3867164/1‑150 (MQ=255)
                                                                                                                                         gTTCAACATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATCTTAAATAAGACTATTTAAg             <  1:2205872/150‑1 (MQ=255)
                                                                                                                                         gTTCAAAATTGTAGCAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATCTTAAATAAGACTATTTAAg             >  2:2403695/1‑150 (MQ=255)
                                                                                                                                          ttCAAAATTGCTACACTTATGAACCAGTCCATTTTTCCTCATATTCAGTTGTTATTATTCCACCCTATCACGTCTCTCGTCATCTCCTAGTATTATCAGTTAGTTACACACACAAAGAGATCTGAACATCTTAAATAGTCTTATTTAAGa            <  2:4045014/150‑1 (MQ=255)
                                                                                                                                          ttCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATCTTAAATAAGACTATTTAAGa            >  2:4766606/1‑150 (MQ=255)
                                                                                                                                          ttCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATCTTAAATAAGACTATTTAAGa            >  2:1121633/1‑150 (MQ=255)
                                                                                                                                          ttCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATCTTAAATAAGACTATTTAAGa            <  1:1011898/150‑1 (MQ=255)
                                                                                                                                           tCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATCTTAAATAAGACTATTTAAGAt           >  1:4707316/1‑150 (MQ=255)
                                                                                                                                                  tGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATCTTAAATAGTCTTATTTAAGATgttcaga    <  2:899463/150‑7 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TTTTCATCATGCCTACGCTGTCAACGTATCTCATCCCGGCGCTTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATCTTAAATAAGACTATTTAAGATGCATAACTT  >  NC_000913/2986955‑2987252

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: