Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 624,598 A→T 16.7% intergenic (‑88/‑287) 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,5980AT16.7% 830.6 / 40.7 333intergenic (‑88/‑287)fepB/entCferrienterobactin ABC transporter periplasmic binding protein/isochorismate synthase 1
Reads supporting (aligned to +/- strand):  ref base A (166/90);  new base T (2/60);  total (168/165)
Fisher's exact test for biased strand distribution p-value = 8.59e-21
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 351 total aligned reads displayed.

ATTCCTGAAAGCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAA  >  NC_000913/624449‑624747
                                                                                                                                                     |                                                                                                                                                     
aTTCCTGAAAGCAAAAGTCCTGTTAATAGTAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTGCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCTCACAATCCGTCCCCTCGCCCCTTTGGGgaga                                                                                                                                                       <  1:1627154/150‑1 (MQ=255)
aTTCCTGAAAGCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGgaga                                                                                                                                                       <  2:2460460/150‑1 (MQ=255)
 ttCCTGAAAGCAAAAGTCCTGTTAATAGACGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGagag                                                                                                                                                      <  2:936476/150‑1 (MQ=255)
  tCCTGAAAGCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAgg                                                                                                                                                     <  1:5148802/150‑1 (MQ=255)
  tCCTGAAAGCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAgg                                                                                                                                                     <  1:1682902/150‑1 (MQ=255)
    cTGAAAGCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGt                                                                                                                                                   >  2:5551798/1‑150 (MQ=255)
     tGACAGCAACAGGCCGGTCAAGCGCAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGtt                                                                                                                                                  <  1:5195642/150‑1 (MQ=255)
     tGAAAGCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGtt                                                                                                                                                  >  2:199923/1‑150 (MQ=255)
     tGAAAGCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGtt                                                                                                                                                  >  2:2300148/1‑150 (MQ=255)
        aaGCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAgg                                                                                                                                               <  1:3138507/150‑1 (MQ=255)
        aaGCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAgg                                                                                                                                               <  1:2786375/150‑1 (MQ=255)
        aaGCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAgg                                                                                                                                               <  2:4380914/150‑1 (MQ=255)
         aGCAAAATTACTGTTAAAAGAAGGGCGTTGCGGTAGAGAGTCGCGAGTCTAACAAATCAGCTACCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCACCTGTGGGGAGAGGGTTAggg                                                                                                                                              >  2:5382228/1‑150 (MQ=255)
               aGTCCTGTTAAGAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAggg                                                                                                                                        <  2:2777841/150‑1 (MQ=255)
                  ccTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGaa                                                                                                                                     <  2:1540126/150‑1 (MQ=255)
                     gTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACTAATCAGCTTCGTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAAcag                                                                                                                                  >  2:3469510/1‑150 (MQ=255)
                       tAATAGACGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTTGGGTGAGGGGAAcagcc                                                                                                                                <  2:5628229/150‑1 (MQ=255)
                         aTAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAAcagccag                                                                                                                              >  2:1317101/1‑150 (MQ=255)
                            gAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCAc                                                                                                                           >  1:3296423/1‑150 (MQ=255)
                              aGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAATGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGTGGGGAACAGCCAGCACTg                                                                                                                         <  2:6006673/150‑1 (MQ=255)
                               gggCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAATGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTgg                                                                                                                        >  1:3510945/1‑150 (MQ=255)
                                ggCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGt                                                                                                                       >  2:4711368/1‑150 (MQ=255)
                                 gCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTg                                                                                                                      <  2:2983721/150‑1 (MQ=255)
                                         gTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACAtt                                                                                                              >  2:2448232/1‑150 (MQ=255)
                                         gTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACAtt                                                                                                              >  2:1206343/1‑150 (MQ=255)
                                         gTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACAtt                                                                                                              >  1:1065791/1‑150 (MQ=255)
                                            gagCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAACGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGTGGGGAACAGCCAGCACTGGTGCGAACATTAAc                                                                                                           <  2:5324851/150‑1 (MQ=255)
                                                 gggCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGACCATTAACCCTCa                                                                                                      >  1:1311514/1‑150 (MQ=255)
                                                    ctcgtaccacaAATCAGCTGCCTGTTATTAATAAGGTTATGGGGGTAATCACAATTCCGCAAAAGCGCAAAATCCGCCCCCAGGCCCCTTTGGGTAGCGGTTTGGGGTGCGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCAccc                                                                                                   <  2:2641815/143‑1 (MQ=255)
                                                    cGAGTCTCACAACTCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGTAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCAccc                                                                                                   <  2:3685541/150‑1 (MQ=255)
                                                          tcACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGccc                                                                                             >  1:2866061/1‑150 (MQ=255)
                                                          tcACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGccc                                                                                             >  1:1432052/1‑150 (MQ=255)
                                                          tcACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGccc                                                                                             >  2:4582531/1‑150 (MQ=255)
                                                           cacaAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCt                                                                                            >  2:1709524/1‑150 (MQ=255)
                                                            acaAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTc                                                                                           >  2:5243269/1‑150 (MQ=255)
                                                            acaAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTc                                                                                           >  1:75911/1‑150 (MQ=255)
                                                             caAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCa                                                                                          >  1:5788537/1‑150 (MQ=255)
                                                                  cAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTg                                                                                     >  1:5506734/1‑150 (MQ=255)
                                                                    gCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGa                                                                                   <  2:1405498/150‑1 (MQ=255)
                                                                     cTTCCTTTTATTAATAAGGTTAAGGGCTTAATGAAAAATTCGAAAAAGCGCACAATCCGGCCCCTCGCCCATTTGGGTAGTGGGTTTGGGTTAGGGTAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGaa                                                                                  <  2:2245395/150‑1 (MQ=255)
                                                                       tCCTTTTTTTATTATGGTTATGGGGTTATTGAAAATTTCTACAAAGCGCACAATCGGTCCCCTCGCCCCTTTGGGTAGTGGGTTAGGGTGTGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAgg                                                                                <  2:2143858/150‑1 (MQ=37)
                                                                        ccTGTTATTAATACGGTTATGGGCGTAATGAAAAATTCCACAAAGCGCACAATCCGGCCCCCCGCCCCTTTGGGGTGAGGGTTTGGGTGAGGGGCACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAggg                                                                               <  1:5969826/150‑1 (MQ=255)
                                                                        ccTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAggg                                                                               >  2:3444057/1‑150 (MQ=255)
                                                                         cTGTTATTAATAAGGTTAAGGGCGTAATGACAAATCCGCCAAAGCGCACAATCCGTCCCCCCGCCCCTTTGGGGCGAGGGTTAGGGTTAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGga                                                                              <  2:339423/150‑1 (MQ=255)
                                                                          tGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGgag                                                                             >  1:3204956/1‑150 (MQ=255)
                                                                           ttttttaATAAGGTTACGGGCGTAATGAAAAATTCTACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGTAGAGGGTTTGGGTTAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGgaga                                                                            <  2:3587574/146‑1 (MQ=255)
                                                                            ttattaATATGGTTATGGGTTTATTCAAAACTTCCAAAAAGCGCACAATCCGCCCCCGCGCCCTTTTGGGTTGTGGTTTCGGGTGCGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGagag                                                                           <  1:1186633/150‑1 (MQ=255)
                                                                            ttattaATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGtgag                                                                           >  1:1021042/1‑150 (MQ=255)
                                                                              attaATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAggg                                                                         >  2:3103098/1‑150 (MQ=255)
                                                                               ttaATATGGTTACGGGCGTAATGACAAATTCGACAACGCGAACAAGCGGTCCCCGCGCCCGTTTGGGTAGTGGTTTTGGGGGCGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAgggg                                                                        <  2:521060/150‑1 (MQ=255)
                                                                                taatacGGTTACGGGCTTAATGACAAATACGACAAAGCGCACAATCCGACCCCACGCCCCTTTGGGGAGTGGGTTAGGGTGAGGGCAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAggggg                                                                       <  1:3615204/150‑1 (MQ=255)
                                                                                aaataaGGTAAAGGCCGTAATCACAAATACCACAAACCCCACAAACCGTCCCCACGCCCCTTTGGGGAGAGTGTAAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACAATAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAggggg                                                                       <  2:806274/149‑1 (MQ=255)
                                                                                   tatggTTATGGGCTTACTGAAAAATTCGACAACGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTTGGGTGCGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAg                                                                    <  1:3817305/147‑1 (MQ=25)
                                                                                       gTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACg                                                                >  1:916935/1‑150 (MQ=255)
                                                                                        ttAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGACGGGAGAGGGGGCAGAACgg                                                               >  1:2546640/1‑150 (MQ=255)
                                                                                        ttAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACgg                                                               >  1:5328722/1‑150 (MQ=255)
                                                                                        tactgggCGTAATCACAACTTCGACAACGCGCACAATCCGCCCCCGCGCCCCTTGGGGGTGTGGGTTTGGGTGCGGGCAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACgg                                                               <  1:4627189/146‑1 (MQ=255)
                                                                                          acgggCTTATTGACAAATTCGACAACGCGCACAATCCGACCCCTCGCCCCTTTGGGTTGTGGTTTTGGGTGTGGGTAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGgcg                                                             <  2:1097718/148‑1 (MQ=255)
                                                                                          acgggCGTAATGACAAATCCGACAACGCGCACAATCCGTCCCCCCGCCCCTTTGGGGAGCGGGTTAGGGTGCGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGgcg                                                             <  1:3460786/148‑1 (MQ=255)
                                                                                          aaGGGCGTAATTAAAAATTCGAAAACGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGTGGTTTTGGGTGTGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGgcg                                                             <  1:4988638/150‑1 (MQ=25)
                                                                                             ggCGTATTGCCACATTCGACACCGCGCCCACGCCGGCCCCTCGCCCGTTGGGGGCGCGGGTTTGGGGGCGGGGAACCGCCAGCACTGGTGCGAACCGGAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAg                                                          <  2:1799947/150‑1 (MQ=255)
                                                                                             ggCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAg                                                          >  1:5669328/1‑150 (MQ=255)
                                                                                             ggCGTAATGACAAATACGACAACGCGCACACTCCGTCCCCTCGCCCCTTTGGGTCGTGGGTTTGGGTGTGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAg                                                          <  2:2113090/150‑1 (MQ=255)
                                                                                              gCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAgg                                                         >  2:4038215/1‑150 (MQ=255)
                                                                                               cttAATTACAATTTCGACCAAGCGCAAAACCCGCCCCCTCGCCCCTTTGGGGTGTGGGTTTGGGTGTGGGCAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGa                                                        <  1:87323/148‑1 (MQ=255)
                                                                                                gTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGAc                                                       >  2:4843287/1‑150 (MQ=255)
                                                                                                gTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGAc                                                       >  2:4201319/1‑150 (MQ=255)
                                                                                                  aaTGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACAt                                                     >  1:4188752/1‑150 (MQ=255)
                                                                                                    tGACAAATGCGCCACCGCGCACAATCCGGCCCCGCGCCCCTTGGGGGAGCGGGTTCGGGTGAGGGGAACAGCCAGCACTGGTGCGACCATGAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACCTca                                                   <  2:4239791/150‑1 (MQ=255)
                                                                                                      aCAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATcaca                                                 >  2:4535953/1‑150 (MQ=255)
                                                                                                          agtCGCCAACGCGCACACTCCGTCCCCTCGCCCCTTGGGGGTGTGGGTTGGGGTGTGGGGAACCGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTgc                                             <  1:2088521/148‑1 (MQ=255)
                                                                                                           ttCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTTGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTgcg                                            <  1:5383610/150‑1 (MQ=255)
                                                                                                           ttCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTgcg                                            >  1:4142145/1‑150 (MQ=255)
                                                                                                              gACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCtt                                         >  2:5165894/1‑150 (MQ=255)
                                                                                                               aCAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTa                                        <  2:5432242/150‑1 (MQ=255)
                                                                                                                      gcACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAAt                                 >  1:2505525/1‑150 (MQ=255)
                                                                                                                        acaATCCGTCCCCTCGCCCCTTTGGGGAGTGGGTTTGGGTGTGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATcc                               <  1:1266146/150‑1 (MQ=255)
                                                                                                                          aaTCCGCCCCCCCGCCCGTTGGGGGCGGGGGTTTGGGGGGGGGCAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCAt                             <  2:2553105/150‑1 (MQ=255)
                                                                                                                            tCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCa                           >  1:656678/1‑150 (MQ=255)
                                                                                                                             ccGTCCCCTCGCCCCTTTGGGGAGTGGGTTCGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCaa                          <  1:5254662/150‑1 (MQ=255)
                                                                                                                               gTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCaata                        >  1:2154549/1‑150 (MQ=255)
                                                                                                                                     tCGCCCCTTTGGGTTGTGGTTTTGGGTGAGGGCAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCtt                  <  1:5599695/150‑1 (MQ=255)
                                                                                                                                     tCGCCCCTTTGGGTAGAGGGTTAGGGTGTGGGTAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCtt                  <  2:67692/150‑1 (MQ=255)
                                                                                                                                        ccccTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTtctc               >  2:4930140/1‑150 (MQ=255)
                                                                                                                                        ccccTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTtctc               >  2:5075516/1‑150 (MQ=255)
                                                                                                                                         cccTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCa              <  1:4524831/150‑1 (MQ=255)
                                                                                                                                          ccTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCAt             <  2:2866061/150‑1 (MQ=255)
                                                                                                                                            tttGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCAttt           >  2:4077168/1‑150 (MQ=255)
                                                                                                                                               ggggTGTGGGTTTGGGTGTGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCa        <  2:1582247/150‑1 (MQ=255)
                                                                                                                                               ggggAGAGGGTTTGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCa        <  2:3270711/150‑1 (MQ=255)
                                                                                                                                               ggggAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCa        <  2:4465163/150‑1 (MQ=255)
                                                                                                                                                gggAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCAt       >  2:2282646/1‑150 (MQ=255)
                                                                                                                                                gggAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCAt       >  2:5115507/1‑150 (MQ=255)
                                                                                                                                                 ggtgagGGTTAGGGTGCGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCAtt      <  1:2968821/147‑1 (MQ=255)
                                                                                                                                                 ggAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCAtt      >  2:5133807/1‑150 (MQ=255)
                                                                                                                                                 ggAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCAtt      >  2:380096/1‑150 (MQ=255)
                                                                                                                                                 ggAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCAtt      <  2:75911/150‑1 (MQ=255)
                                                                                                                                                  gagaGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTg     >  2:3856426/1‑150 (MQ=255)
                                                                                                                                                     agGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTaa  >  1:774645/1‑150 (MQ=255)
                                                                                                                                                     agGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTaa  <  2:1919622/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ATTCCTGAAAGCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAA  >  NC_000913/624449‑624747

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: