Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,673,545 C→A 5.3% intergenic (‑44/‑368) mdtJ ← / → tqsA multidrug/spermidine efflux pump membrane subunit MdtJ/autoinducer 2 exporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,673,5450CA 5.3% 277.1 / 3.5 83intergenic (‑44/‑368)mdtJ/tqsAmultidrug/spermidine efflux pump membrane subunit MdtJ/autoinducer 2 exporter
Reads supporting (aligned to +/- strand):  ref base C (40/38);  new base A (3/2);  total (43/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.56e-01

AAATAAAGCCGCCATTTCCCTCACTGACGCTCGCCCATTTCATTGACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACTAAAGCGCGTTGCGCCAGTGCTGCTCACCTACGAGCAAGATAGTGGATGAGGTACGAAGATGGGTAAATGTAGAAAATAACGT  >  NC_000913/1673396‑1673691
                                                                                                                                                     |                                                                                                                                                  
aaaTAAAGCCGCCATTTCCCTCACTGACGCTCGCCCATTTCATTGACAGCGTACCGGTAATTTCTGTCGCAAGAGCCAGCCCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTc                                                                                                                                                    >  1:1000797/1‑150 (MQ=255)
aaaTAAAGCCGCCATTTCCCTCACTGACGCTCGCCCATTTCATTGACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTc                                                                                                                                                    >  1:25377/1‑150 (MQ=255)
    aaaGCCGCCATTTCCCTCACTGACGCTCGCCCATTTCATTGACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCg                                                                                                                                                <  1:336102/150‑1 (MQ=255)
      aGCCGCCATTTCCCTCACTGACGCTCGCCCATTTCATTGACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGaa                                                                                                                                              <  1:187417/150‑1 (MQ=255)
      aGCCGCCATTTCCCTCACTGACGCTCGCCCATTTCATTGACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGaa                                                                                                                                              >  2:1331591/1‑150 (MQ=255)
        ccgccATTTCCCTCACTGACGCTCGCCCATTTCATTGACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAAc                                                                                                                                            <  2:463592/150‑1 (MQ=255)
          gccATTTCCCTCACTGACGCTCGCCCATTTCATTGACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTGAATTAAGCGaaaatt                                                                                                                                          <  1:1062948/150‑4 (MQ=255)
          gccATTTCCCTCACTGACGCTCGCCCATTTCATTGACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTGAATTAAGCGaaaatt                                                                                                                                          <  2:1245828/150‑4 (MQ=255)
             aTTTCCCTCACTGACGCTCGCCCATTTCATTGACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCaaag                                                                                                                                       <  1:464237/150‑1 (MQ=255)
             aTTTCCCTCACTGACGCTCGCCCATTTCATTGACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCaaag                                                                                                                                       <  1:851255/150‑1 (MQ=255)
             aTTTCCCTCACTGACGCTCGCCCATTTCATTGACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCaaag                                                                                                                                       <  2:125198/150‑1 (MQ=255)
                 cccTCACTGACGCTCGCCCATTTCATTGACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCaaagtaaa                                                                                                                                   >  2:830255/1‑150 (MQ=255)
                       cTGACGCTCGCCCATTTCATTGACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTc                                                                                                                             <  1:980456/150‑1 (MQ=255)
                            gCTCGCCCATTTCATTGACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGcc                                                                                                                        <  2:1352997/150‑1 (MQ=255)
                                 cccATTTCATTGACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCcgcg                                                                                                                   >  1:20567/1‑150 (MQ=255)
                                       tCATTGACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAgc                                                                                                             >  2:1222475/1‑150 (MQ=255)
                                           tGACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCaa                                                                                                         >  1:542406/1‑150 (MQ=255)
                                              cAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCt                                                                                                      >  2:1258681/1‑150 (MQ=255)
                                                gCGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCttt                                                                                                    >  2:574763/1‑150 (MQ=255)
                                                 cGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTg                                                                                                   >  2:621934/1‑150 (MQ=255)
                                                   tACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGc                                                                                                 <  1:1301325/150‑1 (MQ=255)
                                                   tACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGc                                                                                                 >  2:876892/1‑150 (MQ=255)
                                                    aCCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCa                                                                                                <  1:1331591/150‑1 (MQ=255)
                                                      cGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCaga                                                                                              >  2:925199/1‑150 (MQ=255)
                                                      cGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCATTCTCATGAAAGAGAATTATTTTAATTTTCGCTTAATTAAGATAAACCAAAGTAAAGTGTTAAAGCCACGCGAATAGCGCAAGCTTTGGCaga                                                                                              >  2:644643/1‑150 (MQ=255)
                                                       ggTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCagaa                                                                                             >  2:455325/1‑150 (MQ=255)
                                                          aaTTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAgaagaa                                                                                          >  1:724438/1‑150 (MQ=255)
                                                            tttCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAg                                                                                        <  2:1241143/150‑1 (MQ=255)
                                                             ttCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCTAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAgg                                                                                       >  1:131814/1‑150 (MQ=255)
                                                             ttCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCTAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAgg                                                                                       <  2:131814/150‑1 (MQ=255)
                                                               cTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGAc                                                                                     >  1:485972/1‑150 (MQ=255)
                                                               cTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGAc                                                                                     >  2:1346660/1‑150 (MQ=255)
                                                                    gCAATAGCCAGACCTAATAAAATCCAATAAATATACTTTGTGCTTCTCCTGCAAGATAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACTAAAg                                                                                <  2:419002/150‑1 (MQ=255)
                                                                      aaTAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACTAAAgcg                                                                              <  2:672029/150‑1 (MQ=255)
                                                                       aTAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACTAAAgcgc                                                                             <  2:1021927/150‑1 (MQ=255)
                                                                       aTAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACTAAAgcgc                                                                             >  2:801887/1‑150 (MQ=255)
                                                                           ccAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACTAAAGCGCGTTg                                                                         >  1:1371405/1‑150 (MQ=255)
                                                                               aCCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACTAAAGCGCGTTGCGcc                                                                     <  1:1411483/150‑1 (MQ=255)
                                                                               aCCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACTAAAGCGCGTTGCGcc                                                                     <  1:621934/150‑1 (MQ=255)
                                                                               aCCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACTAAAGCGCGTTGCGcc                                                                     <  1:59299/150‑1 (MQ=255)
                                                                               aCCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACTAAAGCGCGTTGCGcc                                                                     <  2:140895/150‑1 (MQ=255)
                                                                                ccTAATAAATTCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACTAAAGCGCGTTGCGCCa                                                                    >  1:2293/1‑150 (MQ=255)
                                                                                ccTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACTAAAGCGCGTTGCGCCa                                                                    >  1:1342891/1‑150 (MQ=255)
                                                                                ccTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACTAAAGCGCGTTGCGCCa                                                                    >  2:313113/1‑150 (MQ=255)
                                                                                   aataaAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACTAAAGCGCGTTGCGCCAGtg                                                                 <  1:925199/150‑1 (MQ=255)
                                                                                   aataaAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACTAAAGCGCGTTGCGCCAGtg                                                                 >  1:1424465/1‑150 (MQ=255)
                                                                                   aataaAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACTAAAGCGCGTTGCGCCAGtg                                                                 <  1:1165333/150‑1 (MQ=255)
                                                                                    ataaAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACTAAAGCGCGTTGCGCCAGtgc                                                                >  2:766352/1‑150 (MQ=255)
                                                                                          tCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACTAAAGCGCGTTGCGCCAGTGCTGCTCa                                                          <  1:830255/150‑1 (MQ=255)
                                                                                          tCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACTAAAGCGCGTTGCGCCAGTGCTGCTCa                                                          <  1:574763/150‑1 (MQ=255)
                                                                                          accAATAAATATACATTGTCCTTCTCCTGAAAGAGAATTATTTTAATTTTCGCTTAATTCACCGAAACCAAAGTAAAGTGGTCAAGCACCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACTAAAGCGCGTAGCGCCAGTCCTGCTCa                                                          >  2:1031094/2‑150 (MQ=255)
                                                                                             aataaataTACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACTAAAGCGCGTTGCGCCAGTGCTGCTCACCt                                                       >  2:1067496/1‑150 (MQ=255)
                                                                                                aaataTACATTGTCCTTCTCCTGCAAGAGATTTATTTTAATTTTAGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATATCGCAAGCTTTGGCAGAAGAAAGGACTAATGCGCGTTGGGCCAGTGCTGCTCACCTACg                                                    <  2:135193/150‑1 (MQ=255)
                                                                                                aaataTACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACTAAAGCGCGTTGCGCCAGTGCTGCTCACCTACg                                                    <  1:801887/150‑1 (MQ=255)
                                                                                                  atatACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACTAAAGCGCGTTGCGCCAGTGCTGCTCACCTACGAg                                                  <  1:832276/150‑1 (MQ=255)
                                                                                                       cATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACTAAAGCGCGTTGCGCCAGTGCTGCTCACCTACGAGCAAGa                                             <  2:25377/150‑1 (MQ=255)
                                                                                                        aTTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACTAAAGCGCGTTGCGCCAGTGCTGCTCACCTACGAGCAAGAt                                            <  2:724438/150‑1 (MQ=255)
                                                                                                           gTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACTAAAGCGCGTTGCGCCAGTGCTGCTCACCTACGAGCAAGATAGt                                         >  1:157235/1‑150 (MQ=255)
                                                                                                            tCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACTAAAGCGCGTTGCGCCAGTGCTGCTCACCTACGAGCAAGATAGTg                                        <  2:542406/150‑1 (MQ=255)
                                                                                                             ccTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACTAAAGCGCGTTGCGCCAGTGCTGCTCACCTACGAGCAAGATAGTgg                                       >  2:1052188/1‑150 (MQ=255)
                                                                                                                tctcCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACTAAAGCGCGTTGCGCCAGTGCTGCTCACCTACGAGCAAGATAGTGGATg                                    >  2:149160/1‑150 (MQ=255)
                                                                                                                   ccTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACTAAAGCGCGTTGCGCCAGTGCTGCTCACCTACGAGCAAGATAGTGGATGAgg                                 >  1:83384/1‑150 (MQ=255)
                                                                                                                   ccTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACTAAAGCGCGTTGCGCCAGTGCTGCTCACCTACGAGCAAGATAGTGGATGAgg                                 >  2:846688/1‑150 (MQ=255)
                                                                                                                    cTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACTAAAGCGCGTTGCGCCAGTGCTGCTCACCTACGAGCAAGATAGTGGATGAGGt                                >  1:782313/1‑150 (MQ=255)
                                                                                                                     tGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACTAAAGCGCGTTGCGCCAGTGCTGCTCACCTACGAGCAAGATAGTGGATGAGGTa                               >  2:1293117/1‑150 (MQ=255)
                                                                                                                     tGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACTAAAGCGCGTTGCGCCAGTGCTGCTCACCTACGAGCAAGATAGTGGATGAGGTa                               <  1:756603/150‑1 (MQ=255)
                                                                                                                     tGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACTAAAGCGCGTTGCGCCAGTGCTGCTCACCTACGAGCAAGATAGTGGATGAGGTa                               >  2:871453/1‑150 (MQ=255)
                                                                                                                     tGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACTAAAGCGCGTTGCGCCAGTGCTGCTCACCTACGAGCAAGATAGTGGATGAGGTa                               >  2:1136571/1‑150 (MQ=255)
                                                                                                                      gCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACTAAAGCGCGTTGCGCCAGTGCTGCTCACCTACGAGCAAGATAGTGGATGAGGTAc                              <  1:805418/150‑1 (MQ=255)
                                                                                                                       cAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACTAAAGCGCGTTGCGCCAGTGCTGCTCACCTACGAGCAAGATAGTGGATGAGGTACg                             <  1:1264961/150‑1 (MQ=255)
                                                                                                                        aaGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACTAAAGCGCGTTGCGCCAGTGCTGCTCACCTACGAGCAAGATAGTGGATGAGGTACGa                            <  1:461840/150‑1 (MQ=255)
                                                                                                                            gaATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACTAAAGCGCGTTGCGCCAGTGCTGCTCACCTACGAGCAAGATAGTGGATGAGGTACGAAGAt                        >  2:1018758/1‑150 (MQ=255)
                                                                                                                                 attTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACTAAAGCGCGTTGCGCCAGTGCTGCTCACCTACGAGCAAGATAGTGGATGAGGTACGAAGATGGGTa                   >  1:832095/1‑150 (MQ=255)
                                                                                                                                  ttttAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACTAAAGCGCGTTGCGCCAGTGCTGCTCACCTACGAGCAAGATAGTGGATGAGGTACGAAGATGGGTaa                  >  1:1355029/1‑150 (MQ=255)
                                                                                                                                     ttggtttCGCTGAATTAAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACTAAAGCGCGTTGCGCCAGTGCTGCTCACCTACGAGCAAGATAGTGGATGAGGTACGAAGATGGGTAAATg               >  1:995110/5‑150 (MQ=255)
                                                                                                                                      aaTTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACTAAAGCGCGTTGCGCCAGTGCTGCTCACCTACGAGCAAGATAGTGGATGAGGTACGAAGATGGGTAAATGt              >  2:483622/1‑150 (MQ=255)
                                                                                                                                       ggtttCGCTGAATTAAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACTAAAGCGCGTTGCGCCAGTGCTGCTCACCTACGAGCAAGATAGTGGATGAGGTACGAAGATGGGTAAATGTa             >  1:168997/3‑150 (MQ=255)
                                                                                                                                       aTTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCTAGCTTTGGCAGAAGAAAGGACTAAAGCGCGTTGCGCCAGTGCTGCTCACCTACGAGCAAGATAGTGGATGAGGTACGAAGATGGGTAAATGTa             <  1:1400335/150‑1 (MQ=255)
                                                                                                                                           tCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACTAAAGCGCGTTGCGCCAGTGCTGCTCACCTACGAGCAAGATAGTGGATGAGGTACGAAGATGGGTAAATGTAGaaa         >  2:1284618/1‑150 (MQ=255)
                                                                                                                                            cGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACTAAAGCGCGTTGCGCCAGTGCTGCTCACCTACGAGCAAGATAGTGGATGAGGTACGAAGATGGGTAAATGTAGaaaa        >  1:564406/1‑150 (MQ=255)
                                                                                                                                             gCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACTAAAGCGCGTTGCGCCAGTGCTGCTCACCTACGAGCAAGATAGTGGATGAGGTACGAAGATGGGTAAATGTAGAAAAt       >  2:879862/1‑150 (MQ=255)
                                                                                                                                              cTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACTAAAGCGCGTTGCGCCAGTGCTGCTCACCTACGAGCAAGATAGTGGATGAGGTACGAAGATGGGTAAATGTAGAAAATa      <  1:1222475/150‑1 (MQ=255)
                                                                                                                                               ttAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACTAAAGCGCGTTGCGCCAGTGCTGCTCACCTACGAGCAAGATAGTGGATGAGGTACGAAGATGGGTAAATGTAGAAAATaa     <  1:313113/150‑1 (MQ=255)
                                                                                                                                               ttAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACTAAAGCGCGTTGCGCCAGTGCTGCTCACCTACGAGCAAGATAGTGGATGAGGTACGAAGATGGGTAAATGTAGAAAATaa     <  1:1052188/150‑1 (MQ=255)
                                                                                                                                                  aTTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACTAAAGCGCGTTGCGCCAGTGCTGCTCACCTACGAGCAAGATAGTGGATGAGGTACGAAGATGGGTAAATGTAGAAAATAACGt  <  1:849056/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
AAATAAAGCCGCCATTTCCCTCACTGACGCTCGCCCATTTCATTGACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACTAAAGCGCGTTGCGCCAGTGCTGCTCACCTACGAGCAAGATAGTGGATGAGGTACGAAGATGGGTAAATGTAGAAAATAACGT  >  NC_000913/1673396‑1673691

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: