Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,673,545 C→A 6.0% 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 6.0% 335.1 / 7.6 100intergenic (‑44/‑368)mdtJ/tqsAmultidrug/spermidine efflux pump membrane subunit MdtJ/autoinducer 2 exporter
Reads supporting (aligned to +/- strand):  ref base C (39/55);  new base A (4/2);  total (43/57)
Fisher's exact test for biased strand distribution p-value = 3.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.39e-01

Only 100 of 105 total aligned reads displayed.

AAATAAAGCCGCCATTTCCCTCACTGACGCTCGCCCATTTCATTGACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACTAAAGCGCGTTGCGCCAGTGCTGCTCACCTACGAGCAAGATAGTGGATGAGGTACGAAGATGGGTAAATGTAGAAAATAACGTCCT  >  NC_000913/1673396‑1673694
                                                                                                                                                     |                                                                                                                                                     
aaaTAAAGCCGCCATTTCCCTCACTGACGCTCGCCCATTTCATTGACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTc                                                                                                                                                       >  1:549547/1‑150 (MQ=255)
aaaTAAAGCCGCCATTTCCCTCACTGACGCTCGCCCATTTCATTGACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTc                                                                                                                                                       <  2:1388809/150‑1 (MQ=255)
aaaTAAAGCCGCCATTTCCCTCACTGACGCTCGCCCATTTCATTGACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTc                                                                                                                                                       <  1:1407790/150‑1 (MQ=255)
   tAAAGCCGCCATTTCCCTCACTGACGCTCGCCCATTTCATTGACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGc                                                                                                                                                    <  1:601533/150‑1 (MQ=255)
   tAAAGCCGCCATTTCCCTCACTGACGCTCGCCCATTTCATTGACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGc                                                                                                                                                    <  1:546133/150‑1 (MQ=255)
    aaaGCCGCCATTTCCCTCACTGACGCTCGCCCATTTCATTGACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCg                                                                                                                                                   <  2:67179/150‑1 (MQ=255)
      aGCCGCCATTTCCCTCACTGACGCTCGCCCATTTCATTGACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGaa                                                                                                                                                 >  2:1549392/1‑150 (MQ=255)
        ccgccATTTCCCTCACTGACGCTCGCCCATTTCATTGACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAAc                                                                                                                                               <  2:221262/150‑1 (MQ=255)
         cgccATTTCCCTCACTGACGCTCGCCCATTTCATTGACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTGAATTAAGCGaaaat                                                                                                                                              <  1:429253/150‑3 (MQ=255)
          gccATTTCCCTCACTGACGCTCGCCCATTTCATTGACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCa                                                                                                                                             >  1:1481216/1‑150 (MQ=255)
             aTTTCCCTCACTGACGCTCGCCCATTTCATTGACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCaaag                                                                                                                                          >  1:1197194/1‑150 (MQ=255)
             aTTTCCCTCACTGACGCTCGCCCATTTCATTGACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCaaag                                                                                                                                          <  1:32208/150‑1 (MQ=255)
             aTTTCCCTCACTGACGATCGCCCATTTCATTGACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCaaag                                                                                                                                          <  1:657342/150‑1 (MQ=255)
               ttCCCTCACTGACGCTCGCCCATTTCATTGACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCaaagta                                                                                                                                        <  2:935996/150‑1 (MQ=255)
                tCCCTCACTGACGCTCGCCCATTTCATTGACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCaaagtaa                                                                                                                                       >  1:838389/1‑150 (MQ=255)
                 cccTCACTGACGCTCGCCCATTTCATTGACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCaaagtaaa                                                                                                                                      <  2:340127/150‑1 (MQ=255)
                 cccTCACTGACGCTCGCCCATTTCATTGACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCaaagtaaa                                                                                                                                      >  2:833261/1‑150 (MQ=255)
                  ccTCACTGACGCTCGCCCATTTCATTGACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCaaagtaaag                                                                                                                                     >  1:273111/1‑150 (MQ=255)
                    tCACTGACGCTCGCCCATTTCATTGACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTg                                                                                                                                   <  1:1170424/150‑1 (MQ=255)
                     cACTGACGCTCGCCCATTTCATTGACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTgg                                                                                                                                  <  2:582722/150‑1 (MQ=255)
                     cACTGACGCTCGCCCATTTCATTGACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTgg                                                                                                                                  <  2:1567416/150‑1 (MQ=255)
                           cGCTCGCCCATTTCATTGACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCACGc                                                                                                                            >  2:604735/1‑150 (MQ=255)
                              tCGCCCATTTCATTGACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGcccc                                                                                                                         >  2:529373/1‑150 (MQ=255)
                               cGCCCATTTCATTGACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCcg                                                                                                                        <  2:1481216/150‑1 (MQ=255)
                                 cccATTTCATTGACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCcgcg                                                                                                                      >  1:1373576/1‑150 (MQ=255)
                                 cccATTTCATTGACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCcgcg                                                                                                                      >  2:785033/1‑150 (MQ=255)
                                 cccATTTCATTGACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCcgcg                                                                                                                      <  2:1065554/150‑1 (MQ=255)
                                  ccATTTCATTGACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCgcgc                                                                                                                     <  2:802897/150‑1 (MQ=255)
                                    aTTTCATTGACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCAt                                                                                                                   <  2:1305672/150‑1 (MQ=255)
                                     tttCATTGACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATa                                                                                                                  >  1:118678/1‑150 (MQ=255)
                                       tCATTGACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAgc                                                                                                                >  1:1424738/1‑150 (MQ=255)
                                             aCAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGc                                                                                                          >  2:720844/1‑150 (MQ=255)
                                              cAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCt                                                                                                         <  1:147396/150‑1 (MQ=255)
                                                 cGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGTGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTg                                                                                                      >  1:368936/1‑150 (MQ=255)
                                                   tACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGc                                                                                                    >  2:107537/1‑150 (MQ=255)
                                                   tACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGc                                                                                                    <  2:149477/150‑1 (MQ=255)
                                                   tACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGc                                                                                                    <  2:549547/150‑1 (MQ=255)
                                                   tACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGc                                                                                                    <  2:970572/150‑1 (MQ=255)
                                                   tACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGc                                                                                                    <  1:1283147/150‑1 (MQ=255)
                                                   tACAGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAACCCCCGCCCATAGCGCAAGCTTTGGc                                                                                                    >  2:508937/1‑150 (MQ=255)
                                                    aCCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCa                                                                                                   >  1:678115/1‑150 (MQ=255)
                                                    aCCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCa                                                                                                   <  1:1457736/150‑1 (MQ=255)
                                                       ggTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCagaa                                                                                                >  1:1320461/1‑150 (MQ=255)
                                                        gTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCagaag                                                                                               <  2:310245/150‑1 (MQ=255)
                                                            tttCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAg                                                                                           <  2:343221/150‑1 (MQ=255)
                                                            tttCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAg                                                                                           <  2:149425/150‑1 (MQ=255)
                                                                 gTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACTa                                                                                      >  2:69552/1‑150 (MQ=255)
                                                                    gCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACTAAAg                                                                                   >  2:780658/1‑150 (MQ=255)
                                                                    gCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACTAAAg                                                                                   <  1:1549392/150‑1 (MQ=255)
                                                                     cAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATATTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACTAAAgc                                                                                  >  1:542814/1‑150 (MQ=255)
                                                                      aaTAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACTAAAgcg                                                                                 >  1:257301/1‑150 (MQ=255)
                                                                         aGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACTAAAGCGCGt                                                                              <  2:434566/150‑1 (MQ=255)
                                                                            cAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACTAAAGCGCGTTgc                                                                           >  1:100981/1‑150 (MQ=255)
                                                                            cAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACTAAAGCGCGTTgc                                                                           >  2:490493/1‑150 (MQ=255)
                                                                               aCCTAATAAAATCCAATAAATATACATTTTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACTAAAGCGCGTTGCGcc                                                                        <  2:851329/150‑1 (MQ=255)
                                                                               aCCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGTTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCTCCGCGCATAGCGCAATCTTTGGCAGAAGAAAGGACTAAAGCGCGTTGCGcc                                                                        <  2:762277/150‑1 (MQ=255)
                                                                               aCCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACTAAAGCGCGTTGCGcc                                                                        <  2:1417717/150‑1 (MQ=255)
                                                                               aCCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACTAAAGCGCGTTGCGcc                                                                        <  2:273111/150‑1 (MQ=255)
                                                                                ccTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACTAAAGCGCGTTGCGCCa                                                                       >  2:31057/1‑150 (MQ=255)
                                                                                  taataaAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACTAAAGCGCGTTGCGCCAGt                                                                     >  1:1579035/1‑150 (MQ=255)
                                                                                   aataaAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCACGCTTTGGCAGAAGAAAGGACTAACGCGCGTTGCGCCAGtg                                                                    >  2:1036988/1‑150 (MQ=255)
                                                                                    ataaAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACTAAAGCGCGTTGCGCCAGtgc                                                                   <  2:894923/150‑1 (MQ=255)
                                                                                     taaAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGGATAGCGAAAGCTTTGGCAGATGAAAGGACTAAAGCGCGTTGCGGAAGtgct                                                                  <  2:1373576/150‑1 (MQ=255)
                                                                                          tCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACTAAAGCGCGTTGCGCCAGTGCTGCTCa                                                             >  1:969568/1‑150 (MQ=255)
                                                                                           ccAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACTAAAGCGCGTTGCGCCAGTGCTGCTCAc                                                            >  2:473209/1‑150 (MQ=255)
                                                                                             aataaataTACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGACGAAAGGACTAAAGCGCGTTGCGCCAGTGCTGCTCACCt                                                          <  1:833261/150‑1 (MQ=255)
                                                                                              ataaataTACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACTAAAGCGCGTTGCGCCAGTGCTGCTCACCTa                                                         <  2:1424738/150‑1 (MQ=255)
                                                                                                 aataTACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACTAAAGCGCGTTGCGCCAGTGCTGCTCACCTACGa                                                      <  1:31057/150‑1 (MQ=255)
                                                                                                  atatACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACTAAAGCGCGTTGCGCCAGTGCTGCTCACCTACGAg                                                     <  1:369027/150‑1 (MQ=255)
                                                                                                   tataCATTGTCCTTCTCCTGCAAGAGAATTATTTTATTTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACTAAAGCGCGTTGCGCCAGTGCTGCTCACCTACGAGc                                                    <  2:1561802/150‑1 (MQ=255)
                                                                                                   tataCATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACTAAAGCGCGTTGCGCCAGTGCTGCTCACCTACGAGc                                                    >  2:793825/1‑150 (MQ=255)
                                                                                                       cATTGTCCTTCTCCTGCAAGTGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACTAAAGCGCGTTGCGCCAGTGCTGCTCACCTACGAGCAAGa                                                <  2:368936/150‑1 (MQ=255)
                                                                                                       cATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACTAAAGCGCGTTGCGCCAGTGCTGCTCACCTACGAGCAAGa                                                <  2:1288777/150‑1 (MQ=255)
                                                                                                          tGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACTAAAGCGCGTTGCGCCAGTGCTGCTCACCTACGAGCAAGATAg                                             <  1:1211089/150‑1 (MQ=255)
                                                                                                           gTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACTAAAGCGCGTTGCGCCAGTGCTGCTCACCTACGAGCAAGATAGt                                            >  2:269072/1‑150 (MQ=255)
                                                                                                            tCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACTAAAGCGCGTTGCGCCAGTGCTGCTCACCTACGAGCAAGATAGTg                                           <  1:720844/150‑1 (MQ=255)
                                                                                                             ccTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACTAAAGCGCGTTGCGCCAGTGCTGCTCACCTACGAGCAAGATAGTgg                                          >  2:326863/1‑150 (MQ=255)
                                                                                                                tctcCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACTAAAGCGCGTTGCGCCAGTGCTGCTCACCTACGAGCAAGATAGTGGATg                                       <  2:118678/150‑1 (MQ=255)
                                                                                                                    cTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACTAAAGCGCGTTGCGCCAGTGCTGCTCACCTACGAGCAAGATAGTGGATGAGGt                                   <  1:529373/150‑1 (MQ=255)
                                                                                                                    cTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACTAAAGCGCGTTGCGCCAGTGCTGCTCACCTACGAGCAAGATAGTGGATGAGGt                                   >  1:1076047/1‑150 (MQ=255)
                                                                                                                      gCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACTAAAGCGCGTTGCGCCAGTGCTGCTCACCTACGAGCAAGATAGTGGATGAGGTAc                                 >  1:146365/1‑150 (MQ=255)
                                                                                                                       cAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACAAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACTAAAGCGCGTTGCGCCAGAGCTGCTCACCTACGAGCAAGATAGTGGATGAGGTACg                                >  2:933974/1‑150 (MQ=255)
                                                                                                                             aaTTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACTAAAGCGCGTTGCGCCAGTGCTGCTCACCTACGAGCAAGATAGTGGATGAGGTACGAAGATg                          <  1:107537/150‑1 (MQ=255)
                                                                                                                                 attTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACTAAAGCGCGTTGCGCCAGTGCTGCTCACCTACGAGCAAGATAGTGGATGAGGTACGAAGATGGGTa                      <  1:1156768/150‑1 (MQ=255)
                                                                                                                                  ttttAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACTAAAGCGCGTTGCGCCAGTGCTGCTCACCTACGAGCAAGATAGTGGATGAGGTACGAAGATGGGTaa                     >  2:160364/1‑150 (MQ=255)
                                                                                                                                       ggtttCGCTGAATTAAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACTAAAGCGCGTTGCGCCAGTGCTGCTCACCTACGAGCAAGATAGTGGATGAGGTACGAAGATGGGTAAATGTa                >  1:1620481/3‑150 (MQ=255)
                                                                                                                                       ggtttCGCTGAATTAAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACTAAAGCGCGTTGCGCCAGTGCTGCTCACCTACGAGCAAGATAGTGGATGAGGTACGAAGATGGGTAAATGTa                >  1:1373013/3‑150 (MQ=255)
                                                                                                                                       aTTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACTAAAGCGCGTTGCGCCAGTGCTGCTCACCTACGAGCAAGATAGTGGATGAGGTACGAAGATGGGTAAATGTa                <  1:780658/150‑1 (MQ=255)
                                                                                                                                        ttttCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACTAAAGCGCGTTGCGCCAGTGCTGCTCACCTACGAGCAAGATAGTGGATGAGGTACGAAGATGGGTAAATGTAg               >  2:931976/1‑150 (MQ=255)
                                                                                                                                        gtttCGCTGAATTAAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACTAAAGCGCGTTGCGCCAGTGCTGCTCACCTACGAGCAAGATAGTGGATGAg                                     >  1:1192399/2‑128 (MQ=255)
                                                                                                                                        gtttCGCTGAATTAAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACTAAAGCGCGTTGCGCCAGTGCTGCTCACCTACGAGCAAGATAGTGGATGAg                                     <  2:1192399/127‑1 (MQ=255)
                                                                                                                                        gtttCGCTGAATTAAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACTAAAGCGCGTTGCGCCAGTGCTGCTCACCTACGAGCAAGATAGTGGATGAGGTACGAAGATGGGTAAATGTAg               >  2:662402/2‑150 (MQ=255)
                                                                                                                                              cTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACTAAAGCGCGTTGCGCCAGTGCTGCTCACCTACGAGCAAGATAGTGGATGAGGTACGAAGATGGGTAAATGTAGAAAATa         >  1:1216446/1‑150 (MQ=255)
                                                                                                                                                 aaTTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACTAAAGCGCGTTGCGCCAGTGCTGCTCACCTACGAGCAAGATAGTGGATGAGGTACGAAGATGGGTAAATGTAGAAAATAACg      <  2:633384/150‑1 (MQ=255)
                                                                                                                                                    tCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGGCTAAAGCGCGTTGCGCCAGTGCTGCTCACCTACGAGCAAGATAGTGGATGAGGTACGAAGATGGGTAAATGTAGAAAATAACGTcc   <  2:157482/150‑1 (MQ=255)
                                                                                                                                                    tCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACTAAAGCGCGTTGCGCCAGTGCTGCTCACCTACGAGCAAGATAGTGGATGAGGTACGAAGATGGGTAAATGTAGAAAATAACGTcc   <  2:1579035/150‑1 (MQ=255)
                                                                                                                                                    tCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACTAAAGCGCGTTGCGCCAGTGCTGCTCACCTACGAGCAAGATAGTGGATGAGGTACGAAGATGGGTAAATGTAGAAAATAACGTcc   <  2:1197194/150‑1 (MQ=255)
                                                                                                                                                    tCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACTAAAGCGCGTTGCGCCAGTGCTGCTCACCTACGAGCAAGATAGTGGATGAGGTACGAAGATGGGTAAATGTAGAAAATAACGTcc   <  1:473209/150‑1 (MQ=255)
                                                                                                                                                     cAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACTAAAGCGCGTTGCGCCAGTGCTGCTCACCTACGAGCAAGATAGTGGATGAGGTACGAAGATGGGTAAATGTAGAAAATAACGTCCt  >  2:1510745/1‑150 (MQ=255)
                                                                                                                                                     cAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACTAAAGCGCGTTGCGCCAGTGCTGCTCACCTACGAGCAAGATAGTGGATGAGGTACGAAGATGGGTAAATGTAGAAAATAACGTCCt  >  2:1505072/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AAATAAAGCCGCCATTTCCCTCACTGACGCTCGCCCATTTCATTGACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACTAAAGCGCGTTGCGCCAGTGCTGCTCACCTACGAGCAAGATAGTGGATGAGGTACGAAGATGGGTAAATGTAGAAAATAACGTCCT  >  NC_000913/1673396‑1673694

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: