Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 624,610 A→T 12.1% intergenic (‑100/‑275) fepB ← / → entC ferrienterobactin ABC transporter periplasmic binding protein/isochorismate synthase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913624,6100AT12.1% 833.9 / 13.4 329intergenic (‑100/‑275)fepB/entCferrienterobactin ABC transporter periplasmic binding protein/isochorismate synthase 1
Reads supporting (aligned to +/- strand):  ref base A (173/83);  new base T (2/46);  total (177/152)
Fisher's exact test for biased strand distribution p-value = 1.48e-17
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 344 total aligned reads displayed.

AAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGAT  >  NC_000913/624461‑624759
                                                                                                                                                     |                                                                                                                                                     
aaaaGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGa                                                                                                                                                       >  1:4755921/1‑150 (MQ=255)
aaaaGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGa                                                                                                                                                       <  1:1820544/150‑1 (MQ=255)
aaaaGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGa                                                                                                                                                       <  1:3950062/150‑1 (MQ=255)
aaaaGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGa                                                                                                                                                       <  1:317007/150‑1 (MQ=255)
  aaGTCCTTTTAATAGAAGGGCGTTGCTGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAgg                                                                                                                                                     <  1:2714605/150‑1 (MQ=255)
   aGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAggg                                                                                                                                                    <  2:4452357/150‑1 (MQ=255)
   aGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAggg                                                                                                                                                    <  2:4272370/150‑1 (MQ=255)
       cTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAAc                                                                                                                                                >  1:1546389/1‑150 (MQ=255)
       cTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAAc                                                                                                                                                >  2:5239723/1‑150 (MQ=255)
        tGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAAca                                                                                                                                               >  1:3757192/1‑150 (MQ=255)
         gTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAAcag                                                                                                                                              <  2:3951691/150‑1 (MQ=255)
           tAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCGGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAAcagcc                                                                                                                                            >  2:4129967/1‑150 (MQ=255)
              tAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAAcagccagc                                                                                                                                         >  1:2082977/1‑150 (MQ=255)
                    ggCGTTGCGGTCGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTATGGGCGTAAGGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGt                                                                                                                                   <  2:6024386/150‑1 (MQ=255)
                    ggCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGt                                                                                                                                   >  2:2886391/1‑150 (MQ=255)
                       gTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGGGGGGAACAGCAAGCACAGGTGCg                                                                                                                                >  2:1455842/1‑150 (MQ=255)
                         tGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATACGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGaa                                                                                                                              <  2:6097668/150‑1 (MQ=255)
                          gCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAAc                                                                                                                             >  2:3339943/1‑150 (MQ=255)
                           cGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACa                                                                                                                            >  1:3067698/1‑150 (MQ=255)
                            ggTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACAt                                                                                                                           >  1:2652898/1‑150 (MQ=255)
                                  gCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAAcca                                                                                                                     >  1:1875408/1‑149 (MQ=255)
                                   gggggCGAGTCTAAAAAATCTGCTTCCTGTTATTATTAAGGTTAAGGGCTTAATCACAAATTCTACAAAGCGCACAATCCGTCCCCTCGCCCATTTGGGTAGCGGGTTTGGGTGCGGGGAACAGCCAGCACTGGTGCGAACATTAACCCt                                                                                                                    <  2:5231861/149‑1 (MQ=255)
                                     gggCGAGTCTCAAAAATCAGCTTCCTTTTATTAATAAGGTTAAGGGCGTATTGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGTGTGGGTTTGGGTGCGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCa                                                                                                                  <  2:2737472/150‑1 (MQ=38)
                                        cGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCAccc                                                                                                               >  2:4748872/1‑150 (MQ=255)
                                             caaaaaaaTCCGCTCCCTGTAATAAATACGTTTATGGGCTTAATCAAAAATCCCACAACGCGCAAAAACCGTCCCCCCGCCCTTTTGGGTTGTGGGTTTGGGTGGGGGCAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGcc                                                                                                          <  1:2207930/145‑1 (MQ=255)
                                              tcACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGccc                                                                                                         >  1:174769/1‑150 (MQ=255)
                                                acaAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTc                                                                                                       >  1:2381217/1‑150 (MQ=255)
                                                  aaaTCAGCTTCCTTTTATTATTACGTTTATGGGCTTAATTACAAATTCGAAAAAGCGCACAAGCCGTCCCCTCGCCCTTTGGGGGTGTGGTTTTGGGTGTGGGTAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCAc                                                                                                     <  1:6675379/150‑1 (MQ=255)
                                                  aaaTCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCAc                                                                                                     >  1:1677872/1‑150 (MQ=255)
                                                  aaaTCAGATTCGTTTTATTAATATGGTTAAGGGCGTAATGACAAATTCGACAACGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCAc                                                                                                     <  2:2608264/150‑1 (MQ=37)
                                                   aaTCAGCTTCCTGTTATTAATAAGGTTACGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGTGAGGTTTTGGGTGTGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCAcc                                                                                                    <  2:2677394/150‑1 (MQ=255)
                                                   aaTCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCAcc                                                                                                    >  1:1321177/1‑150 (MQ=255)
                                                        gCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGa                                                                                               >  2:5217142/1‑150 (MQ=255)
                                                         cTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGac                                                                                              >  1:5582455/1‑149 (MQ=255)
                                                         cTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGaa                                                                                              <  2:4220663/150‑1 (MQ=255)
                                                         cTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGaa                                                                                              <  2:1637358/150‑1 (MQ=255)
                                                          ttCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAg                                                                                             >  1:5262280/1‑150 (MQ=255)
                                                           tCCTTTTATTAATAAGGTTATGGGCTTAATTACAAATTCGACAAAGCGCACAATCCTTCCCCTCGCCCCTTTGGGGTGTGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAgg                                                                                            <  2:5396949/150‑1 (MQ=25)
                                                            ccTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAggg                                                                                           <  2:3805451/150‑1 (MQ=255)
                                                             cttttattaATAAGGTTACGGGCTTAATGACAAATTCGACAACGCGCAAAATCCGTCCCCTCCCCCCTTTGGGGAGTGGGTTTGGGTGCGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGga                                                                                          <  1:4941691/147‑1 (MQ=255)
                                                              tGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGgag                                                                                         >  2:6391766/1‑150 (MQ=255)
                                                                  attaATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAggg                                                                                     >  2:2125235/1‑150 (MQ=255)
                                                                       taaGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAg                                                                                >  1:5300109/1‑150 (MQ=255)
                                                                       taaGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAg                                                                                >  2:2569495/1‑150 (MQ=255)
                                                                         aGGTTATGGGCGTATTGACAAATTCGAAAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGTGTGGTTTCGGGTGTGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGaa                                                                              <  2:791859/150‑1 (MQ=25)
                                                                           gTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACg                                                                            >  2:1927705/1‑150 (MQ=255)
                                                                               cgggCGTATTGACAAATCCGCCACCGCGCACCCTCCGTCCCCTCGCCCCTTTGGGGAGCGGGTTCGGGTGCGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGgcgc                                                                        <  2:6455682/149‑1 (MQ=21)
                                                                               aGGGCTTAATGAAAAATTCGACAAAGCGCACAATCCGCCCCCTCGCCCCTTGGGGTAGTGGTTTTGGGTGTGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGgcgc                                                                        <  1:3562290/150‑1 (MQ=255)
                                                                                gggCTTAATGAAAAATTCTAAAACGCGCACAATCCGCCCCCGCGCCCGTTGGGGTTGTGTTTTTGGGTGTGGGTAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCa                                                                       <  1:3221493/150‑1 (MQ=255)
                                                                                gggCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGTTGTGGGTTTGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCa                                                                       <  1:6756153/150‑1 (MQ=255)
                                                                                  gCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAgg                                                                     >  1:1179002/1‑150 (MQ=255)
                                                                                  gCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAgg                                                                     >  1:1061258/1‑150 (MQ=255)
                                                                                  gCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAgg                                                                     >  1:4302226/1‑150 (MQ=255)
                                                                                   cttAATGACAAATTCCACAAAGCGCAAAATCCGTCCCCGCGCCCCTTTGGGGAGTGGGTTTGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGa                                                                    <  1:4827652/148‑1 (MQ=255)
                                                                                    gTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGGGAGGGGGCAGAACGGCGCAGGAc                                                                   >  1:4647082/1‑150 (MQ=255)
                                                                                     tAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGCTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACa                                                                  >  2:2216418/1‑150 (MQ=255)
                                                                                            aaaTTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACAtt                                                           >  2:6280512/1‑150 (MQ=255)
                                                                                             aaTCCGCCAAAGCCCCCAACCCGTCCCCTCGCCCCTTGGGGGAGCGGTTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTg                                                          <  1:5982428/150‑1 (MQ=255)
                                                                                              aTTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTgc                                                         <  2:4548206/150‑1 (MQ=255)
                                                                                               ttCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTgcg                                                        >  1:3012000/1‑150 (MQ=255)
                                                                                                  gACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCtt                                                     >  2:2067243/1‑150 (MQ=255)
                                                                                                   aCAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTa                                                    >  2:6759860/1‑150 (MQ=255)
                                                                                                    cAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTAt                                                   >  1:2954954/1‑150 (MQ=255)
                                                                                                     aaaGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATg                                                  >  2:2438991/1‑150 (MQ=255)
                                                                                                     aaaGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATg                                                  >  1:3689568/1‑150 (MQ=255)
                                                                                                     aaaGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATg                                                  >  1:4042498/1‑150 (MQ=255)
                                                                                                            acaATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATcc                                           >  2:3120828/1‑150 (MQ=255)
                                                                                                              aaTCCGTCCCCTCGCCCCTTTGGGGTGTGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCAt                                         <  2:5307595/150‑1 (MQ=255)
                                                                                                              aaTCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCAt                                         <  2:4908871/150‑1 (MQ=255)
                                                                                                               aTCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCg                                                                                                                                >  1:3523112/1‑62 (MQ=255)
                                                                                                               aTCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCg                                                                                                                                <  2:3523112/62‑1 (MQ=255)
                                                                                                                tCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCa                                       <  2:1102740/150‑1 (MQ=255)
                                                                                                                 cctaacccTCTCCCCAAAGGGGCGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCGGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCaa                                      >  1:5134831/6‑150 (MQ=37)
                                                                                                                  cGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCaat                                     >  2:5611272/1‑150 (MQ=255)
                                                                                                                   taacccTCTCCCCAAAGGGGCGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATca                                                               >  1:2742013/4‑123 (MQ=25)
                                                                                                                   gTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCaata                                    >  2:5842810/1‑150 (MQ=255)
                                                                                                                     ccccTCGCCCCTTTGGGGAGAGGTTTAGGGTGTGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGAGCAGGACATCACATTGCGCTTATGCGAATCCATCaataat                                  <  1:6615668/150‑1 (MQ=255)
                                                                                                                      cccTCTCCCCAAAGGGGCGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATg                                 >  1:736415/1‑150 (MQ=37)
                                                                                                                        cTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCt                               >  1:5204131/1‑150 (MQ=255)
                                                                                                                            ccccaaaggggCGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTtctc                           >  2:3606847/8‑150 (MQ=255)
                                                                                                                               cTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCAtt                        >  1:5760290/1‑150 (MQ=255)
                                                                                                                                tttGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCAttt                       >  1:3629593/1‑150 (MQ=255)
                                                                                                                                tttGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCAttt                       >  1:4859006/1‑150 (MQ=255)
                                                                                                                                   ggggAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCa                    >  2:294265/1‑150 (MQ=255)
                                                                                                                                     ggtgtgggTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCAtt                  <  2:1546482/145‑1 (MQ=255)
                                                                                                                                     ggAGTGGGTTTGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCAtt                  <  2:1864626/150‑1 (MQ=255)
                                                                                                                                     ggAGTGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCAtt                  <  2:2986960/150‑1 (MQ=255)
                                                                                                                                     ggAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCAtt                  >  2:1275517/1‑150 (MQ=255)
                                                                                                                                       agagGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGt                >  2:6698676/1‑150 (MQ=255)
                                                                                                                                         cgggTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGATCATTAACCCTCACCCCAGCCCTCTCCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGTATCCATCAATAATGCTTCTCATTTTCATTGTaa              <  1:5432812/149‑1 (MQ=255)
                                                                                                                                          gggTTAGGGTGTGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAAc             <  2:5891330/150‑1 (MQ=255)
                                                                                                                                           ggTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAAcc            >  2:87847/1‑150 (MQ=255)
                                                                                                                                               tgggTGCGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACaa        <  2:2357034/149‑1 (MQ=255)
                                                                                                                                                gggTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAAc       <  1:6498006/150‑1 (MQ=255)
                                                                                                                                                gggTGAGGGCAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAAc       <  1:6336418/150‑1 (MQ=255)
                                                                                                                                                    gAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGa   <  2:5262284/150‑1 (MQ=255)
                                                                                                                                                     aGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGAt  <  1:6071128/150‑1 (MQ=255)
                                                                                                                                                     aGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGAt  <  1:982162/150‑1 (MQ=255)
                                                                                                                                                     aGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGAt  <  1:255663/150‑1 (MQ=255)
                                                                                                                                                     aGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGAt  <  1:1927705/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGAT  >  NC_000913/624461‑624759

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: