Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 624,590 T→A 5.4% intergenic (‑80/‑295) 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,5900TA 5.4% 709.1 / 22.6 204intergenic (‑80/‑295)fepB/entCferrienterobactin ABC transporter periplasmic binding protein/isochorismate synthase 1
Reads supporting (aligned to +/- strand):  ref base T (95/95);  new base A (4/7);  total (101/103)
Fisher's exact test for biased strand distribution p-value = 5.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.78e-01

Only 100 of 208 total aligned reads displayed.

TGCGGCTATTCCTGAAAGCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTT  >  NC_000913/624442‑624738
                                                                                                                                                    |                                                                                                                                                    
 gCGGCTATTCCTGAAAGCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTg                                                                                                                                                    <  2:168418/150‑1 (MQ=255)
  cgtcTATTCCTGAAAGCAACAGTCCTGAGAATAGAAGGGCGGTGCGGTAGAGCGTGGCGAGTCTCCCAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTgg                                                                                                                                                   <  1:524819/147‑1 (MQ=255)
   ggCTATTCCTGAAAGCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTggg                                                                                                                                                  <  1:407462/150‑1 (MQ=255)
   ggCTATTCCTGAAAGCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTggg                                                                                                                                                  <  1:4224858/150‑1 (MQ=255)
      tATTCCTGAAAGCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGgag                                                                                                                                               <  1:3528023/150‑1 (MQ=255)
      tATTCCTGAAAGCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGgag                                                                                                                                               <  1:2744513/150‑1 (MQ=255)
       aTTCCTGAAAGCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGgtga                                                                                                                                              <  2:2970967/150‑1 (MQ=255)
       aTTCCTGAAAGCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGgaga                                                                                                                                              <  1:2911410/150‑1 (MQ=255)
        ttCCTGAAAGCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGagag                                                                                                                                             <  2:3276653/150‑1 (MQ=255)
        ttCCTGAAAGCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGagag                                                                                                                                             <  2:3372016/150‑1 (MQ=255)
           cTGAAAGCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGt                                                                                                                                          >  1:977559/1‑150 (MQ=255)
           cTGAAAGCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGt                                                                                                                                          >  2:881952/1‑150 (MQ=255)
            tGAAAGCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCTCCCCAAAGGGGCGAgggga                                                                                                                                         <  2:3336355/150‑3 (MQ=255)
            tGAAAGCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCTCCCCAAAGGGGCGAgggga                                                                                                                                         <  2:693187/150‑3 (MQ=255)
            tGAAAGCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGtt                                                                                                                                         >  1:1123595/1‑150 (MQ=255)
             gAAAGCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTa                                                                                                                                        >  2:2019347/1‑150 (MQ=255)
               aaGCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAgg                                                                                                                                      <  1:394245/150‑1 (MQ=255)
               aaGCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAgg                                                                                                                                      <  1:1295899/150‑1 (MQ=255)
               aaGCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAgg                                                                                                                                      <  2:3689267/150‑1 (MQ=255)
                 gCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGt                                                                                                                                    >  2:1311420/1‑150 (MQ=255)
                   aaaaGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGa                                                                                                                                  >  2:3627780/1‑150 (MQ=255)
                     aaGTCCTGTTAATAGAAGGGCGGTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATCAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAgg                                                                                                                                <  2:4113673/150‑1 (MQ=255)
                           tGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAAca                                                                                                                          >  2:410732/1‑150 (MQ=255)
                             ttAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAAcagc                                                                                                                        >  1:1077960/1‑150 (MQ=255)
                             ttAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAAcagc                                                                                                                        >  2:2450063/1‑150 (MQ=255)
                               aaTAGATGGGCGTTGCGGTAGAGGGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGTTTATGGGGTTAATGACAATTTCGAAAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGTTGTGGGTTTGGGTTTGGGGAAcagcca                                                                                                                      <  2:3059509/150‑1 (MQ=255)
                               aaTAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAAcagcca                                                                                                                      >  2:3470116/1‑150 (MQ=255)
                               aaTAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAAcagcca                                                                                                                      >  1:26253/1‑150 (MQ=255)
                                atatAAGGGCGTTGCGGTTGAGCGGGGCTTGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTACGGGCGTAATGACAAATTCGACAACGCGCACATGCCGTCCCCTCGCCCCTTTGGGGTGAGGTTTTGGGTGAGGGGAAcagccag                                                                                                                     <  2:1384896/150‑1 (MQ=38)
                                   gacggggGTTGCGGTAGTGCGGGGCGTGTCTCACAAATCAGCTTCCTTTTATTAATATGGTTATGGGCTTAATTACAAATTCTAAAACGGGCACAATCCGTCCCCTCGCCCCTTTGGGGTGTGGTTTTGGGTGTGGGGAACAGCCAGCAc                                                                                                                  <  2:211532/147‑1 (MQ=255)
                                            tGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGaa                                                                                                         >  1:1792025/1‑150 (MQ=255)
                                                   gagCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAAc                                                                                                  <  2:2306490/150‑1 (MQ=255)
                                                     gCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAAccc                                                                                                >  1:31751/1‑150 (MQ=255)
                                                        gggCGAGTCTCACAAATCAGCTTCCTTTTATTAATAATGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCTTCCCCTCGCCCCTTTGGGGAGAGGGTTTGGGTGTGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCa                                                                                             <  2:528417/150‑1 (MQ=38)
                                                         ggCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTATGGGCGTAATGAAAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCAc                                                                                            <  2:1363232/150‑1 (MQ=255)
                                                           cGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAACTTCGACAAAGCGCACCATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCAccc                                                                                          <  2:2776033/150‑1 (MQ=255)
                                                                ctcACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGTTGAGGGTTAGGGTGTGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGcc                                                                                     <  2:493720/150‑1 (MQ=255)
                                                                ctcACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGTAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGcc                                                                                     <  2:2653269/150‑1 (MQ=255)
                                                                 tcACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGccc                                                                                    >  1:2866827/1‑150 (MQ=255)
                                                                 tcACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGccc                                                                                    >  1:849883/1‑150 (MQ=255)
                                                                 tcACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGccc                                                                                    >  1:917974/1‑150 (MQ=255)
                                                                  cacaAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCt                                                                                   >  1:3132420/1‑150 (MQ=255)
                                                                  cacaAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCt                                                                                   >  1:175829/1‑150 (MQ=255)
                                                                  cacaAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCCCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCt                                                                                   >  2:3540548/1‑150 (MQ=255)
                                                                    aaaaTCCGCTTCCTTTTTTTAATATGTTTATGGGCTTAATCAAAATTTCTAAAAAGCGAAAAATCCGCCCCCCCGCCCCTTGGGGGTTTGGTTTTGGGTGTGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCa                                                                                 <  1:1017774/149‑1 (MQ=255)
                                                                      aaTCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCAcc                                                                               >  1:3872352/1‑150 (MQ=255)
                                                                            cTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGac                                                                         >  2:664044/1‑149 (MQ=255)
                                                                            cTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGaa                                                                         >  1:3285587/1‑150 (MQ=255)
                                                                              tCCTTTTATTAATATGTTTAAGGGCTTAATGAAAAATTCGAAAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGATAGTTTTTGGGTGTGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAgg                                                                       <  2:3177441/150‑1 (MQ=255)
                                                                              tCCTTTTATTAATAAGGTTATGGGCTTATTGAAAAATTCGACAACGCGCACAAGCCGTCCCCTCGCCCCTTTGGGTAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAgg                                                                       <  2:3702957/150‑1 (MQ=25)
                                                                              tCCTGTTTTTAATATGTTTAAGGGCTTAATCACAAATTCGACAACGCGCACAATCCGGCCCCTCGCCCGTTTGGGGATTGGTTTTGGGTGCGGGCAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAgg                                                                       <  2:2014166/150‑1 (MQ=255)
                                                                              tCCTGTTTTAAATATGTTTATGGGCGTACCCCAAAATCCGACACAGCCCACAACCCGCCCCCCCCCCCGTTTGGGTTGGGGGTTGGGGTGCGGGCAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAgg                                                                       <  1:3331129/150‑1 (MQ=255)
                                                                                 tGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGgag                                                                    >  2:1289689/1‑150 (MQ=255)
                                                                                 tGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGgag                                                                    >  2:1240744/1‑150 (MQ=255)
                                                                                 tGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGgag                                                                    >  1:1677291/1‑150 (MQ=255)
                                                                                  gTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGgaga                                                                   >  2:2782111/1‑150 (MQ=255)
                                                                                         atacGGTTACGGCCTTAATCACAAATCCGAAAAAGCGCAAAATCCGTCCCCGCGCCCCTTTGGGGAGTGGTTTTGGGTGTGGGCAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCa                                                            <  1:2328801/150‑1 (MQ=255)
                                                                                         ataaGGTTAGGGGCGTAATCAAAAATTCTAAAACCCGCACAATCCTTCCCCCCGCCCCTTTGGGTTGCGGGTTTGGGTGTGGGCAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCa                                                            <  1:1779507/150‑1 (MQ=255)
                                                                                          taaGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAg                                                           <  2:641265/150‑1 (MQ=255)
                                                                                               ttAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACgg                                                      >  1:1933889/1‑150 (MQ=255)
                                                                                                  aGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGACCAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGgcgc                                                   >  2:4070736/1‑150 (MQ=255)
                                                                                                  aGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGgcgc                                                   >  1:1144261/1‑150 (MQ=255)
                                                                                                   gggCTTAAGGACAAATGCGACAAAGCGCACAATCCGTCCCCGCGCCCCTTTGGGGTGTGGTTTTGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCa                                                  <  2:2057302/150‑1 (MQ=255)
                                                                                                   gggCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCa                                                  <  2:2147767/150‑1 (MQ=255)
                                                                                                   gggCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCa                                                  <  2:662088/150‑1 (MQ=255)
                                                                                                    ggCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAg                                                 >  2:4327328/1‑150 (MQ=255)
                                                                                                    ggCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAg                                                 <  2:348342/150‑1 (MQ=255)
                                                                                                           tGAAAAATTCGACAAAGCGCACAATCCGTCCCCGCGCCCGTTTGGGGAGTGGGTTCGGGTGCGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATca                                          <  1:3103038/150‑1 (MQ=255)
                                                                                                             aCAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATcaca                                        >  1:4110056/1‑150 (MQ=255)
                                                                                                              cAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACAt                                       >  1:2370354/1‑150 (MQ=255)
                                                                                                               acaTTCGCCACACCGCACAATCCGTCCCCTGGCCCCTTTGGGGCGTGGGTTCGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACAtt                                      <  2:941972/148‑1 (MQ=255)
                                                                                                               aaaTTCGACAAAGCGCACAATCCGGCCCCCCGCCCCTTGGGGGTGTGGGTTTGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACAtt                                      <  1:1408975/150‑1 (MQ=255)
                                                                                                               aaaTTCGACAAAGCGCACAATCCGACCCCTCGCCCGTTTGGGTTGTGGGTTTGGGTGTGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACAtt                                      <  1:3014816/150‑1 (MQ=255)
                                                                                                                 aTTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTgc                                    >  1:4330408/1‑150 (MQ=255)
                                                                                                                  ttCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTgcg                                   >  1:4212314/1‑150 (MQ=255)
                                                                                                                   tCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTgcgc                                  <  2:3798884/150‑1 (MQ=255)
                                                                                                                     gACAACGCGCACAAGCCGTCCCCCCGCCCCTTTGGGGAGAGGGTTTGGGTGCGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCtt                                <  2:902742/150‑1 (MQ=255)
                                                                                                                      aCAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACCTTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTa                               >  2:2856560/1‑150 (MQ=255)
                                                                                                                      aCAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTa                               >  1:3303033/1‑150 (MQ=255)
                                                                                                                      aCAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTa                               >  2:2980488/1‑150 (MQ=255)
                                                                                                                        aaaGCGCACACTCCGTCCCCTCGCCCCTTGGGGGTGTGGGTTTGGTTGTGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATg                             <  1:2221692/150‑1 (MQ=255)
                                                                                                                        aaaGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTTGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATg                             <  1:668460/150‑1 (MQ=255)
                                                                                                                         caGCGCACACTCGGGCCCCTCGCCCCGTGGGGGAGTGGGTTAGGGTGCGGGGCACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGc                            <  2:2593653/149‑1 (MQ=255)
                                                                                                                         acgcgcACAACCCGTCCCCTCGCCCCTTGGGGGAGTGGGTGAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGc                            <  2:3848947/148‑1 (MQ=255)
                                                                                                                         acgcgcAAAAGCCGCCCCATCGCCCCTTTGGGGTGAGGGTGAGGTTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCATGACATCACATTGCGCTTATGc                            <  2:1317065/148‑1 (MQ=255)
                                                                                                                         aaGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGc                            <  2:4253561/150‑1 (MQ=255)
                                                                                                                            cgcACAATCCGTCCCCTCGCCCCTTTGGGGTGTGGGTTTGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGaa                         <  1:3107053/150‑1 (MQ=255)
                                                                                                                             gcACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAAt                        >  1:563099/1‑150 (MQ=255)
                                                                                                                              cacaATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATc                       >  1:927268/1‑150 (MQ=255)
                                                                                                                              cacaATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATc                       >  1:1580339/1‑150 (MQ=255)
                                                                                                                              cacaATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATc                       <  2:634866/150‑1 (MQ=255)
                                                                                                                              cacaATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATc                       >  1:1162230/1‑150 (MQ=255)
                                                                                                                                      gtcgcgtCGCCCTTTGGGGGTGTGGGTTTGGTTTAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGTGAGAGGGGTCAGAACGGCGCAGTACATCACATTGCGCTTATGCGAATCCATCaata               <  1:2375994/144‑1 (MQ=255)
                                                                                                                                      gTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCaata               >  2:2395370/1‑150 (MQ=255)
                                                                                                                                       tCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCaataa              >  1:3902317/1‑150 (MQ=255)
                                                                                                                                         cccTCGCCCCTTTGGGGAGTGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATg            <  2:3132825/150‑1 (MQ=255)
                                                                                                                                            tCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCtt         >  1:4305462/1‑150 (MQ=255)
                                                                                                                                            tCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCtt         >  1:3698295/1‑150 (MQ=255)
                                                                                                                                                  cTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCAtt   >  1:4283639/1‑150 (MQ=255)
                                                                                                                                                   tttGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCAttt  >  1:1389342/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TGCGGCTATTCCTGAAAGCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTT  >  NC_000913/624442‑624738

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: