Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,588,928 T→G 5.6% I445S (ATT→AGT)  acrD → aminoglycoside/multidrug efflux system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,588,9280TG 5.6% 1976.4 / 47.3 588I445S (ATT→AGT) acrDaminoglycoside/multidrug efflux system
Reads supporting (aligned to +/- strand):  ref base T (311/236);  new base G (2/39);  total (313/275)
Fisher's exact test for biased strand distribution p-value = 9.55e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 600 total aligned reads displayed.

TTCGCGATGGTGCTGGCGATCGGTCTGCTGGTGGATGACGCCATCGTGGTGGTGGAAAACGTCGAACGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTCGGCGGCACCACCGGTGCCATCTATCGCCAGTTCTCTATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCA  >  NC_000913/2588780‑2589077
                                                                                                                                                    |                                                                                                                                                     
ttCGCGATGGTGCTGGCGATCGGTCTGCTGGTGGATGACGCCATCGTGGTGGTGGAAAACGTCGAACGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGGGGGCACTGGTCGGGAtt                                                                                                                                                      >  1:3413213/1‑150 (MQ=255)
gtCGCGATGGTGCGGGCGATCGGTCTGCTGGTGGATGACGCCATCGTGGTGGTGGAAAACGTCGAACGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGGGGGCACTGGTCGGGAtt                                                                                                                                                      <  2:1814388/149‑1 (MQ=255)
 tCGCGATGGTGCTGGCGATCGGTCTGCTGGTGGATGACGCCATCGTGGTGGTGGAAAAAGTCGAACGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGGGGGCACTGGTCGGGATTg                                                                                                                                                     >  2:9211714/1‑150 (MQ=255)
    cgATGGTGCTGGCGATCGGTCTGCTGGTGGATGACGCCATCGTGGTGGTGGAAAACGTCGAACGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGGGGGCACTGGTCGGGATTGCGa                                                                                                                                                  >  2:9338966/1‑150 (MQ=255)
    cgATGGTGCTGGCGATCGGTCTGCTGGTGGATGACGCCATCGTGGTGGTGGAAAACGTCGAACGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGGGGGCACTGGTCGGGATTGCGa                                                                                                                                                  <  2:6672141/150‑1 (MQ=255)
    cgAGGGTGCTGGCGATCGGTCTGCTGGTGGATGACGCCATCGTGGTGGTGGAAAACGTCGAACGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGGGGGCACTGGTCGGGATTGCGa                                                                                                                                                  <  2:454810/150‑1 (MQ=255)
       tGGTGCTGGCGATCGGTCTGCTGGTGGATGACGCCATCGTGGTGGTGGAAAACGTCGAACGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGGGGGCACTGGTCGGGATTGCGATgg                                                                                                                                               >  2:8055838/1‑150 (MQ=255)
           gCTGGCGATCGGTCTGCTGGTGGATGACGCCATCGTGGTGGTGGAAAACGTCGAACGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGGGGGCACTGGTCGGGATTGCGATGGTTCt                                                                                                                                           >  1:9020124/1‑150 (MQ=255)
            cTGGCGATCGGTCTGCTGGTGGATGACGCCATCGTGGTGGTGGAAAACGTCGAACGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGGGGGCACTGGTCGGGATTGCGATGGTTCtt                                                                                                                                          >  1:1331900/1‑150 (MQ=255)
             tGGCGATCGGTCTGCTGGTGGATGACGCCATCGTGGTGGTGGAAAACGTCGAACGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGGGGGCACTGGTCGGGATTGCGATGGTTCttt                                                                                                                                         <  1:5657229/150‑1 (MQ=255)
             tGGCGATCGGTCTGCTGGTGGATGACGCCATCGTGGTGGTGGAAAACGTCGAACGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGGGGGCACTGGTCGGGATTGCGATGGTTCttt                                                                                                                                         <  1:3981377/150‑1 (MQ=255)
              ggCGATCGGTCTGCTGGTGGATGACGCCATCGTGGTGGTGGAAAACGTCGAACGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGGGGGCACTGGTCGGGATTGCGATGGTTCTTTc                                                                                                                                        <  2:2160339/150‑1 (MQ=255)
                 gATCGGTCTGCTGGTGGATGACGCCATCGTGGTGGTGGAAAACGTCGAACGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGGGGGCACTGGTCGGGATTGCGATGGTTCTTTcggc                                                                                                                                     >  1:9400925/1‑150 (MQ=255)
                  aTCGGTCTGCTGGTGGATGACGCCATCGTGGTGGTGGAAAACGTCGAACGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGGGGGCACTGGTCGGGATTGCGATGGTTCTTTcggcg                                                                                                                                    >  2:9723690/1‑150 (MQ=255)
                         tgctgGTGGATGACGCCATCGGGGTGGTGGAAAACGTCGAACGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTg                                                                                                                             >  2:2637250/1‑150 (MQ=255)
                           ctgGTGGATGACGCCATCGTGGTGGTGGAAAACGTCGAACGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGGGGGCACTGGTCGGGATTGCGATGGtt                                                                                                                                             <  2:570108/132‑1 (MQ=255)
                             ggtggATGACGCCATCGTGGTGGTGGAAAACGTCGAACGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTAcc                                                                                                                         <  1:2565446/150‑1 (MQ=255)
                                ggATGACGCCATCGTGGTGGTGGAAAACGTCGAACGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAAt                                                                                                                      <  1:7917016/150‑1 (MQ=255)
                                 gATGACGCCATCGTGGTGGTGGAAAACGTCGAACGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATg                                                                                                                     >  1:1053083/1‑150 (MQ=255)
                                   tGACGCCATCGTGGTGGTGGAAAACGTCGAACGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGc                                                                                                                   <  2:5625435/150‑1 (MQ=255)
                                   tGACGCCATCGTGGTGGTGGAAAACGTCGAACGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGc                                                                                                                   <  1:5760581/150‑1 (MQ=255)
                                    gACGCCATCGTGGTGGTGGAAAACGTCGAACGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGcc                                                                                                                  >  2:7014937/1‑150 (MQ=255)
                                      cGCCATCGTGGTGGTGGAAAACGTCGAACGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCctt                                                                                                                >  1:8253064/1‑150 (MQ=255)
                                      cGCCATCGTGGTGGTGGAAAACGTCGAACGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCctt                                                                                                                >  2:2402555/1‑150 (MQ=255)
                                       gCCATCGTGGTGGTGGAAAACGTCGAACGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCcttc                                                                                                               >  1:149317/1‑150 (MQ=255)
                                       gCCATCGTGGTGGTGGAAAACGTCGAACGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCcttc                                                                                                               >  1:8659149/1‑150 (MQ=255)
                                              tggtggtggAAAACGTCGAACGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTcggcg                                                                                                        <  1:7282018/150‑1 (MQ=255)
                                               ggtggtggAAAACGTCGAACGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGGGGGCACTGGTCGGGAGTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTcggcgg                                                                                                       <  1:4148148/150‑1 (MQ=255)
                                               ggtggtggAAAACGTCGAACGTATGATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAAGCGATGGGGCAGATCCAGGGGGCACTGGTCGGGAGTGCGATGGTTCTTTCGGCGGTATTTGCACCAATGGCCTTCTTcggcgg                                                                                                       <  1:3233390/150‑1 (MQ=255)
                                                   gtggAAAACGTCGAACGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTCGGCGGcacc                                                                                                   >  1:7680005/1‑150 (MQ=255)
                                                     ggAAAACGTCGAACGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTCGGCGGcaccac                                                                                                 <  2:7145630/150‑1 (MQ=255)
                                                      gAAAACGTCGAACGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCATGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTCGGCGGCaccacc                                                                                                >  1:8792511/1‑150 (MQ=255)
                                                      gAAAACGTCGAACGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTCGGCGGCaccacc                                                                                                >  2:2989027/1‑150 (MQ=255)
                                                        aaaCGTCGAACGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTCGGCGGCACCACCgg                                                                                              >  2:1992197/1‑150 (MQ=255)
                                                        aaaCGTCGAACGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTCGGCGGCACCACCgg                                                                                              <  2:7591247/150‑1 (MQ=255)
                                                         aaCGTCGAACGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTCGGCGGCACCACCGGt                                                                                             >  2:5266303/1‑150 (MQ=255)
                                                         aaCGTCGAACGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTCGGCGGCACCACCGGt                                                                                             <  2:4991947/150‑1 (MQ=255)
                                                         aaCGTCGAACGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGGGGGCACTGGTCGGGATGGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTCGGCGGCACCACCGGt                                                                                             <  1:3788805/150‑1 (MQ=255)
                                                            gTCGAACGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTCGGCGGCACCACCGGTGcc                                                                                          >  1:801336/1‑150 (MQ=255)
                                                             tCGAACGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTCGGCGGCACCACCGGTGCCa                                                                                         >  1:541707/1‑150 (MQ=255)
                                                               gAACGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTCGGCGGCACCACCGGTGCCATc                                                                                       >  1:6298381/1‑150 (MQ=255)
                                                               gAACGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTCGGCGGCACCACCGGTGCCATc                                                                                       >  1:9899762/1‑150 (MQ=255)
                                                               gAACGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTCGGCGGCACCACCGGTGCCATc                                                                                       >  1:2340888/1‑150 (MQ=255)
                                                               gAACGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTCGGCGGCACCACCGGTGCCATc                                                                                       >  1:7453573/1‑150 (MQ=255)
                                                                aaCGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTCGGCGGCACCACCGGTGCCATCt                                                                                      <  2:3429691/150‑1 (MQ=255)
                                                                aaCGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTACTTGTACCAATGGCCTTCTTCGGCGGCACCACCGGTGCCATCt                                                                                      >  1:6402563/1‑150 (MQ=255)
                                                                 aCGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTCGGCGGCACCACCGGTGCCATCTa                                                                                     <  1:3615681/150‑1 (MQ=255)
                                                                          gagtgagGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGGGGGCACTGGTCGGGATGGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTCGGCGGCACCACCGGTGCCATCTATCGCCAGtt                                                                            <  1:2857737/150‑1 (MQ=255)
                                                                            gtgagGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTCGGCGGCACCACCGGTGCCATCTATCGCCAGTtct                                                                          >  1:9216435/1‑150 (MQ=255)
                                                                             tgagGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTCGGCGGCACCACCGGTGCCATCTATCGCCAGTtctc                                                                         >  1:1809006/1‑150 (MQ=255)
                                                                                 gaaggaCTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTCGGCGGCACCACCGGTGCCATCTATCGCCAGTTCTCTAtt                                                                     >  1:6642813/1‑150 (MQ=255)
                                                                                      actcactcCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTCGGCGGCACCACCGGTGCCATCTATCGCCAGTTCTCTATTACCAt                                                                >  1:8459579/1‑150 (MQ=255)
                                                                                       ctcactcCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTCGGCGGCACCACCGGTGCCATCTATCGCCAGTTCTCTATTACCAtt                                                               >  1:6923152/1‑150 (MQ=255)
                                                                                       ctcactcCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTCGGCGGCACCACCGGTGCCATCTATCGCCAGTTCTCTATTACCAtt                                                               >  1:2665725/1‑150 (MQ=255)
                                                                                        tcactcCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTCGGCGGCACCACCGGTGCCATCTATCGCCAGTTCTCTATTACCAttg                                                              >  1:5557032/1‑150 (MQ=255)
                                                                                          actcCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTCGGCGGCACCACCGGTGCCATCTATCGCCAGTTCTCTATTACCAttgtt                                                            >  2:3906317/1‑150 (MQ=255)
                                                                                            tcctcGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTCGGCGGCACCACCGGTGCCATCTATCGCCAGTTCTCTATTACCATTGTTgc                                                          <  2:6259658/150‑1 (MQ=255)
                                                                                             cctcGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGGGGGCACTGGTCGGGAGGGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTCGGCGGCACCACCGGTGCCATCTATCGCCAGTTCTCTATTACCATTGTTgcg                                                         <  1:305650/150‑1 (MQ=255)
                                                                                              ctcGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTCGGCGGCACCACCGGTGCCATCTATCGCCAGTTCTCTATTACCATTGTTgcgg                                                        >  1:5949201/1‑150 (MQ=255)
                                                                                               tcGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTCGGCGGCACCACCGGTGCCATCTATCGCCAGTTCTCTATTACCATTGTTgcggc                                                       <  1:9637094/150‑1 (MQ=255)
                                                                                               tcGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTCGGCGGCACCACCGGTGCCATCTATCGCCAGTTCTCTATTACCATTGTTgcggc                                                       <  1:7016790/150‑1 (MQ=255)
                                                                                                cgcgAAGCCACACGTAAATCGATGGGGCAGATCCAGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTCGGCGGCACCACCGGTGCCATCTATCGCCAGTTCTCTATTACCATTGTTgcggcg                                                      >  1:1680759/1‑150 (MQ=255)
                                                                                                 gcgAAGCCACACGTAAATCGATGGGGCAGATCCAGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTCGGCGGCACCACCGGTGCCATCTATCGCCAGTTCTCTCTTACCATTGTTGCGGCGa                                                     >  1:2214991/1‑150 (MQ=255)
                                                                                                 gcgAAGCCACACGTAAATCGATGGGGCAGATCCAGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTCGGCGGCACCACCGGTGCCATCTATCGCCAGTTCTCTATTACCATTGTTGCGGCGa                                                     <  2:4303336/150‑1 (MQ=255)
                                                                                                  cgAAGCCACACGTAAATCGATGGGGCAGATCCAGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTCGGCGGCACCACCGGTGCCATCTATCGCCAGTTCTCTATTACCATTGTTGCGGCGAt                                                    <  2:5934339/150‑1 (MQ=255)
                                                                                                  cgAAGCCACACGTAAATCGATGGGGCAGATCCAGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTCGGCGGCACCACCGGTGCCATCTATCGCCAGTTCTCTATTACCATTGTTGCGGCGAt                                                    <  2:10186634/150‑1 (MQ=255)
                                                                                                   gAAGCCACACGTAAATCGATGGGGCAGATCCAGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTCGGCGGCACCACCGGTGCCATCTATCGCCAGTTCTCTATTACCATTGTTGCGGCGATg                                                   <  1:1129963/150‑1 (MQ=255)
                                                                                                     aGCCACACGTAAATCGATGGGGCAGATCCAGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTCGGCGGCACCACCGGTGCCATCTATCGCCAGTTCTCTATTACCATTGTTGCGGCGATGGt                                                 <  2:1419991/150‑1 (MQ=255)
                                                                                                       ccACACGTAAATCGATGGGGCAGATCCAGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTCGGCGGCACCACCGGTGCCATCTATCGCCAGTTCTCTATTACCATTGTTGCGGCGATGGTGc                                               >  1:5050433/1‑150 (MQ=255)
                                                                                                       ccACACGTAAATCGATGGGGCAGATCCAGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTCGGCGGCACCACCGGTGCCATCTATCGCCAGTTCTCTATTACCATTGTTGCGGCGATGGTGc                                               >  2:4789832/1‑150 (MQ=255)
                                                                                                        cacaCGTAAATCGATGGGGCAGATCCAGGGGGCACTGGTCGGGATTGCGATGGTTCCTTCGGCGGTATTTGTACCAATGGCCTTCTTCGGCGGCACCACCGGTGCCATCTATCGCCAGTTCTCTATTACCATTGTTGCGGCGATGGTGCt                                              <  2:4666304/150‑1 (MQ=255)
                                                                                                          cacGTAAATCGATGGGGCAGATCCAGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTCGGCGGCACCACCGGTGCCATCTATCGCCAGTTCTCTATTACCATTGTTGCGGCGATGGTGCTGt                                            <  2:8781292/150‑1 (MQ=255)
                                                                                                          cacGTAAATCGATGGGGCAGATCCAGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTCGGCGGCACCACCGGTGCCATCTATCGCCAGTTCTCTATTACCATTGTTGCGGCGATGGTGCTGt                                            <  2:1288590/150‑1 (MQ=255)
                                                                                                          cacGTAAATCGATGGGGCAGATCCAGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTCGGCGGCACCACCGGTGCCATCTATCGCCAGTTCTCTATTACCATTGTTGCGGCGATGGTGCTGt                                            <  2:536275/150‑1 (MQ=255)
                                                                                                           acGTAAATCGATGGGGCAGATCCAGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTCGGCGGCACCACCGGTGCCATCTATCGCCAGTTCTCTATTACCATTGTTGCGGCGATGGTGCTGTc                                           >  1:4649965/1‑150 (MQ=255)
                                                                                                            cGTAAATCGATGGGGCAGATCCAGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTCGGCGGCACCACCGGTGCCATCTATCGCCAGTTCTCTATTACCATTGTTGCGGCGATGGTGCTGTCa                                          <  2:8854986/150‑1 (MQ=255)
                                                                                                            cGTAAATCGATGGGGCAGATCCAGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTCGGCGGCACCACCGGTGCCATCTATCGCCAGTTCTCTATTACCATTGTTGCGGCGATGGTGCTGTCa                                          <  2:3760373/150‑1 (MQ=255)
                                                                                                                aaTCGATGGGGCAGATCCAGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTCGGCGGCACCACCGGTGCCATCTATCGCCAGTTCTCTATTACCATTGTTGCGGCGATGGTGCTGTCAGTAc                                      <  2:7407579/150‑1 (MQ=255)
                                                                                                                 aTCGATGGGGCAGATCCAGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTCGGCGGCACCACCGGTGCCATCTATCGCCAGTTCTCTATTACCATTGTTGCGGCGATGGTGCTGTCAGTACt                                     <  2:9587361/150‑1 (MQ=255)
                                                                                                                    gATGGGGCAGATCCAGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTCGGCGGCACCACCGGTGCCATCTATCGCCAGTTCTCTATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGt                                  >  2:6839926/1‑150 (MQ=255)
                                                                                                                     aTGGGGCAGATCCAGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCACTGGCCTTCTTCGGCGGCACCACCGGTGCCATCTATCGCCAGTTCTCTATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGtc                                 >  2:3153290/1‑149 (MQ=255)
                                                                                                                     aTGGGGCAGATCCAGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTCGGCGGCACCACCGGTGCCATCTATCGCCAGTTCTCTATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTa                                 >  2:1100058/1‑150 (MQ=255)
                                                                                                                       ggggCAGATCCAGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTCGGCGGCACCACCGGTGCCATCTATCGCCAGTTCTCTATTACCACTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGc                               >  1:10154472/1‑150 (MQ=255)
                                                                                                                         ggCAGATCCAGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTCGGCGGCACCACCGGTGCCATCTATCGCCAGTTCTCTATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCga                             <  2:7742386/150‑1 (MQ=255)
                                                                                                                         ggCAGATCCAGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTCGGCGGCACCACCGGTGCCATCTATCGCCAGTTCTCTATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCga                             <  1:2466676/150‑1 (MQ=255)
                                                                                                                          gCAGATCCAGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTCGGCGGCACCACCGGTGCCATCTATCGCCAGTTCTCTATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCgat                            >  2:3613925/1‑150 (MQ=255)
                                                                                                                          gCAGATCCAGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTCGGCGGCACCACCGGTGCCATCTATCGCCAGTTCTCTATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCgat                            >  1:6544264/1‑150 (MQ=255)
                                                                                                                            aGATCCAGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTCGGCGGCACCACCGGTGCCATCTATCGCCAGTTCTCTATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCgatga                          >  2:8227154/1‑150 (MQ=255)
                                                                                                                              aTCCAGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTCGGCGGCACCACCGGTGCCATCTATCGCCAGTTCTCTATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATc                        <  1:4846909/150‑1 (MQ=255)
                                                                                                                              aTCCAGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTCGGCGGCACCACCGGTGCCATCTATCGCCAGTTCTCTATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATc                        <  2:9039851/150‑1 (MQ=255)
                                                                                                                                     gggCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTCGGCGGCACCACCGGTGCCATCTATCGCCAGTTCTCTATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTc                 <  1:4449565/150‑1 (MQ=255)
                                                                                                                                      ggCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTCGGCGGCACCACCGGTGCCATCTATCGCCAGTTCTCTATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTcc                >  2:1694503/1‑150 (MQ=255)
                                                                                                                                            ggTCGGGAGGGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTCGGCGGCACCACCGGTGCCATCTATCGCCAGTTCTCTATTTCCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGctct          <  1:6081720/150‑1 (MQ=255)
                                                                                                                                                gggATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTCGGCGGCACCACCGGTGCCATCTATCGCCAGTTCTCTATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCtgtgt      <  2:1910548/150‑1 (MQ=255)
                                                                                                                                                 ggATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTCGGCGGCACCACCGGTGCCATCTATCGCCAGTTCTCTATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAg                                >  2:8527110/1‑123 (MQ=255)
                                                                                                                                                  gATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTCGGCGGCACCACCGGTGCCATCTATCGCCAGTTCTCTATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGc    <  2:5461302/150‑1 (MQ=255)
                                                                                                                                                  gATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTCGGCGGCACCACCGGTGCCATCTATCGCCAGTTCTCTATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGc    <  1:189831/150‑1 (MQ=255)
                                                                                                                                                  gATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTCGGCGGCACCACCGGTGCCATCTATCGCCAGTTCTCTATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGc    <  1:5305084/150‑1 (MQ=255)
                                                                                                                                                    ttGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTCGGCGGCACCACCGGTGCCATCTATCGCCAGTTCTCTATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCca  <  2:6485708/150‑1 (MQ=255)
                                                                                                                                                    ttGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTCGGCGGCACCACCGGTGCCATCTATCGCCAGTTCTCTATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCca  >  1:10025315/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TTCGCGATGGTGCTGGCGATCGGTCTGCTGGTGGATGACGCCATCGTGGTGGTGGAAAACGTCGAACGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTCGGCGGCACCACCGGTGCCATCTATCGCCAGTTCTCTATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCA  >  NC_000913/2588780‑2589077

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