Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,663,553 A→C 6.1% intergenic (+26/+337) mdtF → / ← gadW multidrug efflux pump RND permease MdtF/DNA‑binding transcriptional dual regulator GadW

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,663,5530AC 6.1% 841.8 / 35.2 259intergenic (+26/+337)mdtF/gadWmultidrug efflux pump RND permease MdtF/DNA‑binding transcriptional dual regulator GadW
Reads supporting (aligned to +/- strand):  ref base A (115/128);  new base C (4/12);  total (119/140)
Fisher's exact test for biased strand distribution p-value = 1.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.81e-01

Only 100 of 262 total aligned reads displayed.

GCAAAACGCGGTAGGTACCGGCGTGATGGGCGGGATGTTTGCCGCAACAGTGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAAACCACTCCAGCCATTCGTCATCAATTTGAACAACACAATACTGACCCACA  >  NC_000913/3663404‑3663702
                                                                                                                                                     |                                                                                                                                                     
gcAAAACGCGGTAGGTACCGGCGTGATGGGCGGGATGTTTGCCGCAACAGTGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACa                                                                                                                                                       <  1:859967/150‑1 (MQ=255)
 cAAAACGCGGTAGGTACCGGCGTGATGGGCGGGATGTTTGCCGCAACAGTGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACAt                                                                                                                                                      >  1:455410/1‑150 (MQ=255)
  aaaaCGCGGTAGGTACCGGCGTGATGGGCGGGATGTTTGCCGCAACAGTGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATg                                                                                                                                                     >  1:1938245/1‑150 (MQ=255)
    aaCGCGGTAGGTACCGGCGTGATGGGCGGGATGTTTGCCGCAACAGTGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGaa                                                                                                                                                   <  2:1096087/150‑1 (MQ=255)
    aaCGCGGTAGGTACCGGCGTGATGGGCGGGATGTTTGCCGCAACAGTGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGaa                                                                                                                                                   <  1:60198/150‑1 (MQ=255)
    aaCGCGGTAGGTACCGGCGTGATGGGCGGGATGTTTGCCGCAACAGTGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGaa                                                                                                                                                   <  1:1474017/150‑1 (MQ=255)
      cgcgGTAGGTACCGGCGTGATGGGCGGGATGTTTGCCGCAACAGTGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAgg                                                                                                                                                 >  2:1404302/1‑150 (MQ=255)
       gcgGTAGGTACCGGCGTGATGGGCGGGATGTTTGCCGCAACAGTGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGt                                                                                                                                                <  1:278392/150‑1 (MQ=255)
         ggtaggtaCCGGCGTGATGGGCGGGATGTTTGCCGCAACAGTGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTc                                                                                                                                              <  1:1688310/150‑1 (MQ=255)
           taggtaCCGGCGTGATGGGCGGGATGTTTGCCGCAACAGTGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAg                                                                                                                                            <  1:137748/150‑1 (MQ=255)
             ggtaCCGGCGTGATGGGCGGGATGTTTGCCGCAACAGTGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGcc                                                                                                                                          <  2:1430478/150‑1 (MQ=255)
             ggtaCCGGCGTGATGGGCGGGATGTTTGCCGCAACAGTGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGcc                                                                                                                                          <  2:1722789/150‑1 (MQ=255)
               taCCGGCGTGATGGACGGGATGTTTGCCGCAACAGTGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCtt                                                                                                                                        <  2:661645/150‑1 (MQ=255)
                aCCGGCGTGATGGGCGGGATGTTTGCCGCAACAGTGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGCTGAACCTTCATGTTCAACCTTa                                                                                                                                       <  2:763897/150‑1 (MQ=255)
                 ccGGCGTGATGGGCGGGATGTTTGCCGCAACAGTGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGCTGAACCTTCATGTTCAACCTTAc                                                                                                                                      <  2:1862776/150‑1 (MQ=255)
                  cGGCGTGATGGGCGGGATGTTTGCCGCAACAGTGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGATAGTAAGGTTGAACATGAAGGTTCAGCCTTAct                                                                                                                                     <  1:1482225/150‑1 (MQ=255)
                  cGGCGTGATGGGCGGGATGTTTGCCGCAACAGTGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTAct                                                                                                                                     <  2:1354821/150‑1 (MQ=255)
                  cGGCGTGATGGGCGGGATGTTTGCCGCAACAGTGCTGGCAATTTACTTCGTTCCAGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGCTGAACCTTCATGTTCAACCTTAct                                                                                                                                     <  2:574622/150‑1 (MQ=255)
                      gTGATGGGCGGGATGTTTGCCGCAACAGTGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCttt                                                                                                                                 <  2:195999/150‑1 (MQ=255)
                       tGATGGGCGGGATGTTTGCCGCAACAGTGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGCTGAACCTTCATGTTCAACCTTActctcat                                                                                                                                <  1:1002267/150‑4 (MQ=38)
                        gATGGGCGGGATGTTTGCCGCAACAGTGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTcc                                                                                                                               <  1:356550/150‑1 (MQ=255)
                          tGGGCGGGATGTTTGCCGCAACAGTGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTg                                                                                                                             >  2:244034/1‑150 (MQ=255)
                             gCGGGATGTTTGCCGCAACAGTGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTa                                                                                                                          >  2:527564/1‑150 (MQ=255)
                               gggATGTTTGCCGCAACAGTGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAAc                                                                                                                        <  1:1066481/150‑1 (MQ=255)
                               gggATGTTTGCCGCAACAGTGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAAc                                                                                                                        <  1:1697919/150‑1 (MQ=255)
                               gggATGTTTGCCGCAACAGTGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGCTGAACCTTCATGTTCAACCTTAccctcat                                                                                                                                <  2:20671/142‑4 (MQ=38)
                               gggATGTTTGCCGCAACAGTGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGCTGAACCTTCATGTTCAACCTTAccctcat                                                                                                                                >  1:20671/1‑139 (MQ=38)
                                ggATGTTTGCCGCAACAGTGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAAcc                                                                                                                       >  1:1850503/1‑150 (MQ=255)
                                ggATGTTTGCCGCAACAGTGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAAcc                                                                                                                       <  1:1368280/150‑1 (MQ=255)
                                      ttGCCGCAACAGTGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTcac                                                                                                                 <  1:291178/150‑1 (MQ=255)
                                      ttGCCGCAACAGTGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTcac                                                                                                                 >  2:308301/1‑150 (MQ=255)
                                      ttGCCGCAACAGTGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTcac                                                                                                                 >  2:192540/1‑150 (MQ=255)
                                       tGCCGCAACAGTGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTcacc                                                                                                                >  1:175648/1‑150 (MQ=255)
                                       tGCCGCAACAGTGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGCTTCAGCCTTACTCTTTCCTGCTAACCATTcacc                                                                                                                <  2:1343585/150‑1 (MQ=255)
                                        gCCGCAACAGTGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACCTCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTACCCATTcccca                                                                                                               >  2:370635/1‑150 (MQ=255)
                                          cGCAACAGTGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTCGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACt                                                                                                             >  1:800464/1‑150 (MQ=255)
                                           gCAACAGTGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTc                                                                                                            <  1:713442/150‑1 (MQ=255)
                                            cAACAGTGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTca                                                                                                           <  1:883767/150‑1 (MQ=255)
                                            cAACAGTGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTca                                                                                                           >  1:101992/1‑150 (MQ=255)
                                                 gTGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCaacaac                                                                                                      >  2:1104537/1‑150 (MQ=255)
                                                   gCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCa                                                                                                    >  2:1039210/1‑150 (MQ=255)
                                                    cTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAg                                                                                                   >  2:816200/1‑150 (MQ=255)
                                                    cTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAg                                                                                                   >  2:1774057/1‑150 (MQ=255)
                                                     tGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGc                                                                                                  >  1:1777224/1‑150 (MQ=255)
                                                          aTTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAAcc                                                                                             >  1:668971/1‑150 (MQ=255)
                                                             tACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTc                                                                                          <  2:251599/150‑1 (MQ=255)
                                                              aCTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCa                                                                                         <  1:817734/150‑1 (MQ=255)
                                                               cTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAg                                                                                        <  2:854044/150‑1 (MQ=255)
                                                                ttCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAgg                                                                                       <  2:1640443/150‑1 (MQ=255)
                                                                ttCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAgg                                                                                       <  1:728555/150‑1 (MQ=255)
                                                                   gTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGTTTCAGCCTTACTCTTTCCTTCTAACCATTCACAACTTATCAATAATCTAAGCGTGTGGAAt                                                                                    <  2:1938245/150‑1 (MQ=255)
                                                                   gTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAAt                                                                                    <  1:937147/150‑1 (MQ=255)
                                                                    ttCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTca                                                                                                           <  2:1671142/126‑1 (MQ=255)
                                                                       cGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATAccc                                                                                <  2:21442/150‑2 (MQ=255)
                                                                       cGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCa                                                                                >  2:460912/1‑150 (MQ=255)
                                                                           cTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCt                                                                            >  1:1440961/1‑150 (MQ=255)
                                                                            tttttCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCtt                                                                           >  1:1614313/1‑150 (MQ=255)
                                                                                tCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTTCTAACCATTCAGCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTa                                                                       >  2:945487/1‑150 (MQ=255)
                                                                                 cGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTTAGGAATACCAGGCTTAGTac                                                                      >  2:843454/1‑150 (MQ=255)
                                                                                   ttGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTacac                                                                    <  1:1187879/150‑1 (MQ=255)
                                                                                   ttGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTacac                                                                    <  2:1373151/150‑1 (MQ=255)
                                                                                    tGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTAcaca                                                                   <  1:996265/150‑1 (MQ=255)
                                                                                       aGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGACTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCa                                                                >  2:985222/1‑150 (MQ=255)
                                                                                        gTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCAc                                                               <  2:1792529/150‑1 (MQ=255)
                                                                                        gTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCAc                                                               >  2:1871158/1‑150 (MQ=255)
                                                                                         tGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACt                                                              <  2:379952/150‑1 (MQ=255)
                                                                                          ggAACATCTCTTTGCCCGCTTTAAAAACGCGTAACGTGTAAATGAGAGTACGGTTGACCATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCCCTCAACCACCCGCTCACCGTCAGGAATACCAGGCGTAGTACCCAGCACTg                                                             >  1:1932200/1‑150 (MQ=255)
                                                                                          ggAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTg                                                             >  1:11625/1‑150 (MQ=255)
                                                                                          ggAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTg                                                             <  2:1440961/150‑1 (MQ=255)
                                                                                           gAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGa                                                            <  2:11625/150‑1 (MQ=255)
                                                                                                tctcTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTa                                                       >  1:1236064/1‑150 (MQ=255)
                                                                                                tctcTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTa                                                       <  1:1330483/150‑1 (MQ=255)
                                                                                                tctcTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTa                                                       <  1:188197/150‑1 (MQ=255)
                                                                                                 ctctTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAg                                                      >  1:499780/1‑150 (MQ=255)
                                                                                                  tctTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGa                                                     >  1:466641/1‑150 (MQ=255)
                                                                                                  tctTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGa                                                     >  1:3461/1‑150 (MQ=255)
                                                                                                    tttGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGaaa                                                   <  1:1800193/150‑1 (MQ=255)
                                                                                                    tttGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGaaa                                                   <  2:792111/150‑1 (MQ=255)
                                                                                                     ttGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAAAc                                                  <  1:848561/150‑1 (MQ=255)
                                                                                                       gCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGATGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAAACCa                                                >  1:1116134/1‑150 (MQ=255)
                                                                                                       gCCCGCCTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAAACCa                                                >  2:499147/1‑150 (MQ=255)
                                                                                                           gCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAAACCACTcc                                            <  1:1404302/150‑1 (MQ=255)
                                                                                                           gCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAAACCACTcc                                            <  1:143124/150‑1 (MQ=255)
                                                                                                               tAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAAt                                                                                    <  2:1295930/106‑1 (MQ=255)
                                                                                                                aaaaaaGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAAACCACTCCAGCCa                                       <  2:722748/150‑1 (MQ=255)
                                                                                                                   aaaGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAAACCACTCCAGCCATTc                                    >  1:1280643/1‑150 (MQ=255)
                                                                                                                      gCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAAACCACTCCAGCCATTCGtc                                 >  1:34007/1‑150 (MQ=255)
                                                                                                                      gCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAAACCACTCCAGCCATTCGtc                                 <  2:967819/150‑1 (MQ=255)
                                                                                                                         tAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAAACCACTCCAGCCATTCGtcatc                              <  1:194759/150‑1 (MQ=255)
                                                                                                                                tAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAAACCACTCCAGCCATTCGTCATCAATTTGa                       <  2:1520337/150‑1 (MQ=255)
                                                                                                                                tAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAAACCACTCCAGCCATTCGTCATCAATTTGa                       <  2:1614313/150‑1 (MQ=255)
                                                                                                                                 aaaTGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAAACCACTCCAGCCATTCGTCATCAATTTGaa                      >  2:802868/1‑150 (MQ=255)
                                                                                                                                  aaTGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAAACCACTCCAGCCATTCGTCATCAATTTGaac                     <  1:1719354/150‑1 (MQ=255)
                                                                                                                                      agagTAAGGCTGAACCTTCATGTTCAACCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAAACCACTCCAGCCATTCGTCATCAATTTGAacaaca                 >  1:1052544/1‑150 (MQ=39)
                                                                                                                                             ggTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACAACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAAACCACTCCAGCCATTCGTAATCAATTTGAACAACACAATACt          >  2:180783/1‑150 (MQ=255)
                                                                                                                                              gTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAAAACCAGGCTTAGTACACAGCACTGAAAGTAGAAACCACTCCAGCCATTCGTCATCAATTTGAACAACACAATACTg         >  2:400375/1‑150 (MQ=255)
                                                                                                                                                 gAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAAACCACTCCAGCCATTCGTCATCAATTTGAACAACACAATACTGAcc      <  1:237853/150‑1 (MQ=255)
                                                                                                                                                    cATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAAACCACTCCAGCCATTCGTCATCAATTTGAACAACACAATACTGACCcac   <  2:702700/150‑1 (MQ=255)
                                                                                                                                                    cATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAAACCACTCCAGCCATTCGTCATCAATTTGAACAACACAATACTGACCcac   <  1:579479/150‑1 (MQ=255)
                                                                                                                                                     aTGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAAACCACTCCAGCCATTCGTCATCAATTTGAACAACACAATACTGACCcaca  <  2:1231963/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCAAAACGCGGTAGGTACCGGCGTGATGGGCGGGATGTTTGCCGCAACAGTGCTGGCAATTTACTTCGTTCCGGTCTTTTTCGTTGTAGTGGAACATCTCTTTGCCCGCTTTAAAAAAGCGTAACGTGTAAATGAGAGTAAGGTTGAACATGAAGGTTCAGCCTTACTCTTTCCTGCTAACCATTCACCACTCAACAACCAGCTAACCGTCAGGAATACCAGGCTTAGTACACAGCACTGAAAGTAGAAACCACTCCAGCCATTCGTCATCAATTTGAACAACACAATACTGACCCACA  >  NC_000913/3663404‑3663702

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