Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,301,071 C→T A670T (GCA→ACA)  mdtO ← membrane translocase (MDR) of MdtNOP efflux pump, PET family

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,301,0710CT97.0% 1299.3 / 19.1 403A670T (GCA→ACA) mdtOmembrane translocase (MDR) of MdtNOP efflux pump, PET family
Reads supporting (aligned to +/- strand):  ref base C (6/6);  new base T (203/181);  total (211/192)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.92e-01

Only 100 of 403 total aligned reads displayed.

CTGTTGATGAGGTGCAGAATCCTTACGTACCAGGGCACAGCCGGAAATCAGCGTCGTGCTGCCGAGAATGCTGCACAACAGCAGACGTGAAAGTTGACGATTGATCATTGCGTGGCTCCTTGTGCCTGTTCCGTGGCGGGCGTTAATGCCGGGGCTGTCCAGTCTGGCAAGCGGGCGAAAAGCTGGCAGACGTGTTGCTCCTGTTTTAATAAGGTGGGCAGGATGGCCTGAGAGGCGGGTGTCTCTTCAAGCGTTATTTGAGGAGAACGGCTGAGTGCGGTTGCCAGACCGGCAGC  >  NC_000913/4300923‑4301218
                                                                                                                                                    |                                                                                                                                                   
  gTTGATGAGGTGCAGAATCCTTACGTACCAGGGCACAGCCGGAAATCAGCGTCGTGCTGCCGAGAATGCTGCACAACAGCAGACGTGAAAGTTGACGATTGATCATTGCGTGGCTCCTTGTGCCTGTTCCGTGGCGGGCGTTAATGTCggg                                                                                                                                                 <  8:1471551/151‑1 (MQ=255)
  gTTGATGAGGTGCAGAATCCTTACGTACCAGGGCACAGCCGGAAATCAGCGTCGTGCTGCCGAGAATGCTGCACAACAGCAGACGTGAAAGTTGACGATTGATCATTGCGTGGCTCCTTGTGCCTGTTCCGTGGCGGGCGTTAATGTCggg                                                                                                                                                 <  3:1198526/151‑1 (MQ=255)
           gggCAGAATCCTTACGTACCAGGGCACAGCCGGAAATCAGCGTCGTGCTGCCGAGAATGCTGCACAACAGCAGACGTGAAAGTTGACGATTGATCATTGCGTGGCTCCTTGTGCCTGTTCCGTGGCGGGCGTTAATGTCGGGGCTGTCCAg                                                                                                                                        <  2:419708/149‑1 (MQ=255)
             gCAGAATCCTTACGTACCAGGGCACAGCCGGAAATCAGCGTCGTGCTGCCGAGAATGCTGCACAACAGCAGACGTGAAAGTTGACGATTGATCATTGCGTGGCTCCTTGTGCCTGTTCCGTGGCGGGCGTTAATGTCGGGGCTGTCCAGt                                                                                                                                       <  5:77180/150‑1 (MQ=255)
             gCAGAATCCTTACGTACCAGGGCACAGCCGGAAATCAGCGTCGTGCTGCCGAGAATGCTGCACAACAGCAGACGTGAAAGTTGACGATTGATCATTGCGTGGCTCCTTGTGCCTGTTCCGTGGCGGGCGTTAATGTCGGGGCTGTCCAGTc                                                                                                                                      <  1:1276922/151‑1 (MQ=255)
              cAGAATCCTTACGTACCAGGGCACAGCCGGAAATCAGCGTCGTGCTGCCGAGAATGCTGCACAACAGCAGACGTGAAAGTTGACGATTGATCATTGCGTGGCTCCTTGTGCCTGTTCCGTGGCGGGCGTTAATGTCGGGGCTGTCCAGTCt                                                                                                                                     <  4:42608/151‑1 (MQ=255)
               aGAATCCTTACGTACCAGGGCACAGCCGGAAATCAGCGTCGTGCTGCCGAGAATGCTGCACAACAGCAGACGTGAAAGTTGACGATTGATCATTGCGTGGCTCCTTGTGCCTGTTCCGTGGCGGGCGTTAATGTCGGGGCTGTCCAGTCTg                                                                                                                                    >  4:73972/1‑151 (MQ=255)
               aGAATCCTTACGTACCAGGGCACAGCCGGAAATCAGCGTCGTGCTGCCGAGAATGCTGCACAACAGCAGACGTGAAAGTTGACGATTGATCATTGCGTGGCTCCTTGTGCCTGTTCCGTGGCGGGCGTTAATGTCGGGGCTGTCCAGTCTg                                                                                                                                    >  3:268492/1‑151 (MQ=255)
                 aaTCCTTACGTACCAGGGCACAGCCGGAAATCAGCGTCGTGCTGCCGAGAATGCTGCACAACAGCACACGGGAAAGTTGACCATTGAACATTGCGTGGCGCCTTGTGCCGGTTCCTTGCCGGGCGTTAATGACGGGGCTGACCAGTCTgcc                                                                                                                                  >  1:1209037/1‑149 (MQ=255)
                    ccTTACGTACCAGGGCACAGCCGGAAATCAGCGTCGTGCTGCCGAGAATGCTGCACAACAGCAGACGTGAAAGTTGACGATTGATCATTGCGTGGCTCCTTGTGCCTGTTCCGTGGCGGGCGTTAATGTCGGGGCTGTCCAGTCTGGCaa                                                                                                                                <  3:90113/150‑1 (MQ=255)
                    ccTTACGTACCAGGGCACAGCCGGAAATCAGCGTCGTGCTGCCGAGAATGCTGCACAACAGCAGACGTGAAAGTTGACGATTGATCATTCGGTGGCTCCTTGTGCCTGTTCCGTGGCGGGCGTTAATGTCGGGGCTGTCCAGTCTGGCa                                                                                                                                 >  8:738082/1‑149 (MQ=255)
                      ttACGTACCAGGGCACAGCCGGAAATCAGCGTCGTGCTGCCGAGAATGCTGCACAACAGCAGACGTGAAAGTTGACGATTGATCATTGCGTGGCTCCTTGTGCCTGTTCCGTGGCGGGCGTTAATGTCGGGGCTGTCCAGTCTGGCAAGCg                                                                                                                             >  8:942174/1‑151 (MQ=255)
                       tACGTACCAGGGCACAGCCGGAAATCAGCGTCGTGCTGCCGAGAATGCTGCACAACAGCAGACGTGAAAGTTGACGATTGATCATTGCGTGGCTCCTTGTGCCTGTTCCGTGGCGGGCGTTAATGTCGGGGCTGTCCAGTCTGGCAAGc                                                                                                                              <  4:807167/149‑1 (MQ=255)
                       tACGTACCAGGGCACAGCCGGAAATCAGCGTCGTGCTGCCGAGAATGCTGCACAACAGCAGACGTGAAAGTTGACGATTGATCATTGCGTGGCTCCTTGTGCCTGTTCCGTGGCGGGCGTTAATGTCGGGGCTGTCCAGTCTGGCAAGCg                                                                                                                             <  7:795697/150‑1 (MQ=255)
                         cGTACCAGGGCACAGCCGGAAATCAGCGTCGTGCTGCCGAGAATGCTGCACAACAGCAGACGTGAAAGTTGACGATTGATCATTGCGTGGCTCCTTGTGCCTGTTCCGTGGCGGGCGTTAATGTCGGGGCTGTCCAGTCTGGCAAGCGGGc                                                                                                                          <  1:1342313/151‑1 (MQ=255)
                          gTACCAGGGCACAGCCGGAAATCAGCGTCGTGCTGCCGAGAATGCTGCACAACAGCAGACGTGAAAGTTGACGATTGATCATTGCGTGGCTCCTTGTGCCTGTTCCGTGGCGGGCGTTAATGTCGGGGCTGTCCAGTCTGGCAAGCGGGCg                                                                                                                         <  1:1031504/151‑1 (MQ=255)
                           tACCAGGGCACAGCCGGAAATCAGCGTCGTGCTGCCGAGAATGCTGCACAACAGCAGACGTGAAAGTTGACGATTGATCATTGCGTGGCTCCTTGTGCCTGTTCCGTGGCGGGCGTTAATGTCGGGGCTGTCCAGTCTGGCAAGCGGGCg                                                                                                                         >  1:657227/1‑150 (MQ=255)
                            aCCAGGGCACAGCCGGAAATCAGCGTCGTGCTGCCGAGAATGCTGCACAACAGCAGACGTGAAAGTTGACGATTGATCATTGCGTGGCTCCTTGTGCCTGTTCCGTGGCGGGCGTTAATGTCGGGGCTGTCCAGTCTGGCAAGCGGGCGaa                                                                                                                       <  5:1543692/151‑1 (MQ=255)
                            aCCAGGGCACAGCCGGAAATCAGCGTCGTGCTGCCGAGAATGCTGCACAACAGCAGACGTGAAAGTTGACGATTGATCATTGCGTGGCTCCTTGTGCCTGTTCCGTGGCGGGCGTTAATGTCGGGGCTGTCCAGTCTGGCAAGCGGGCGa                                                                                                                        <  2:1675028/150‑1 (MQ=255)
                             ccAGGGCACAGCCGGAAATCAGCGTCGTGCTGCCGAGAATGCTGCACACCAGCAGACGTGAAAGTTGACGATTGATCATTGCGTGGCTCCTTGTGCCTGTTCCGTGGCGGGCGTTAATGTCGGGGCTGTCCAGTCTGGCAAGCGGGCGaa                                                                                                                       >  8:825206/1‑150 (MQ=255)
                             ccAGGGCACAGCCGGAAATCAGCGTCGTGCTGCCGAGAATGCTGCACAACAGCAGACGTGAAAGTTGACGATTGATCATTGCGTGGCTCCTTGTGCCTGGTCCGTGGCGGGCGTTAATGTCGGGGCTGTCCAGTCTGGCAAGCGGGCGaa                                                                                                                       >  3:999325/1‑150 (MQ=255)
                               aGGGCAGAGCCGGAAATCAGCGTCGTGCTGCCGAGAATGCTGCACAACAGCAGACGTGAAAGTTGACGATTGATCATTGCGTGGCTCCTTGTGCCTGTTCCGTGGCGGGCGTTAATGTCGGGGCTGTCCAGTCTGGCAAGCGGGCGAAAAg                                                                                                                    <  6:1366774/151‑1 (MQ=255)
                               aGGGCACAGCCGGAAATCAGCGTCGTGCTGCCGAGAATGCTGCACAACAGCAGACGTGAAAGTTGACGATTGATCATTGCGTGGCTCCTTGTGCCTGTTCCGTGGCGGGCGTTAATGTCGGGGCTGTCCAGTCTGGCAAGCGGGCGAAAAg                                                                                                                    <  5:913765/151‑1 (MQ=255)
                                gggCAGAGCCGGAAATGAGCGTCGTGCTGCCGTGAATGCTGCACAACAGCAGACGTGAAAGTTGACGATTGATCATTGCGTGGCTCCTTGTGCCTGTTCCGTGGCGGGCGTTAATGTCGGGGCTGTCCAGTCTGGCAAGGGGGCGAAAAGc                                                                                                                   <  3:1460443/151‑1 (MQ=255)
                                 ggCACAGCCGGAAATCAGGGTCGTGCTGCCGAGAATGCTGCACAACAGCAGACGTGAAAGTTGACGATTGATCATTGCGGGGCTCCTTGTGCCTGTTCCGTGGCGGGCGTTAATGCCGGGGCTGTCCAGTCTGGCAAGCGGGCGAAAAGCt                                                                                                                  <  3:538462/151‑1 (MQ=255)
                                 ggCACAGCCGGAAATCAGCGTCGTGCCGCCGAGAATGCTGCACAACAGCAGACGTGAAAGTTGACGATTGATCATTGCGTGGCTCCTTGTGCCTGTTCCGTGGCGGGCGTTAATGTCGGGGCTGTCCAGTCTGGCAAGCGGGCGAAAAGCt                                                                                                                  <  7:1344095/151‑1 (MQ=255)
                                  gCACAGCCGGAAATCAGCGTCGTGCTGCCGAGAATGCTGCACAACAGCAGACGTGAAAGTTGACGATTGATCATTGCGTGGCTCCTTGTGCCTGTTCCGTGGCGGGCGTTAATGTCGGGGCTGTCCAGTCTGGCAAGCGGGCGAAAAGCTg                                                                                                                 <  4:544423/151‑1 (MQ=255)
                                    acaGCCGGAAATCAGCGTGGTGCTGCCGATAATGGTGGACAACAGGATACGTGAAAGTTGTCGATTGATCATTGCGTTGGTGCTTGTGGCTGTTCCGTGGCGGGCGTTAATGTGGGGGGTGTCCAGTCTGGCAAGCGGGCGAAAAGCTGGc                                                                                                               <  8:1278968/151‑1 (MQ=255)
                                       gCCGGAAATCAGCGTCGTGCTGCCGAGAATGCTGCACAACAGCAGACGTGAAAGTTGACGATTGATCATTGCGTGGCTCCTTGTGCCTGTTCCGTGGCGGGCGTTAATGTCGGGGCTGTCCAGTCTGGCAAGCGGGCGAAAAGCTGGCAGa                                                                                                            <  5:1625692/151‑1 (MQ=255)
                                         cGGAAATCAGCGTCGTGCTGCCGAGAATGCTGCACAACAGCAGACGTGAAAGTTGACGATTGATCATTGCGTGGCTCCTTGTGCCTGTTCCGTGGCGGGCGTTAATGTCGGGGCTGTCCAGTCTGGCAAGCGGGCGAAAAGCTGGCAGACg                                                                                                          >  5:309328/1‑151 (MQ=255)
                                          ggAAATCAGCGTCGTGCTGCCGAGAATGCTGCACAACAGCAGACGTGAAAGTTGACGATTGATCATTGCGTGGCTCCTTGTGCCTGTTCCGTGGCGGGCGTTAATGTCGGGGCTGTCCAGTCTGGCAAGCGGGCGAAAAGCTGGCAGACgt                                                                                                         >  3:776860/1‑151 (MQ=255)
                                           gAAATCAGCGTCGTGCTGCCGAGAATGCTGCACAACAGCAGACGTGAAAGTTGACGATTGATCATTGCGTGGCTCCTTGTGCCTGTTCCGTGGCGGGCGTTAATGTCGGGGCTGTCCAGTCTGGCAAGCGGGCGAAAAGCTGGCAGACgtg                                                                                                        <  2:99925/151‑1 (MQ=255)
                                           gAAATCAGCGTCGTGCTGCCGAGAATGCTGCACAACAGCAGACGTGAAAGTTGACGATTGATCATTGCGTGGCTCCTTGTGCCTGTTCCGTGGCGGGCGTTAATGTCGGGGCTGTCCAGTCTGGCAAGCGGGCGAAAAGCTGGCAGACgtg                                                                                                        <  8:321077/151‑1 (MQ=255)
                                             aaTCAGCGTCGTGCTGCCGAGAATGCTGCACAACAGCAGACGTGAAAGTTGACGATTGATCATTGCGTGGCTCCTTGTGCCTGTTCCGTGGCGGGCGTTAATGTCGGGGCTGTCCAGTCTGGCAAGCGGGCGAAAAGCTGGCAGACGTGtt                                                                                                      <  3:457532/151‑1 (MQ=255)
                                                  gCGTCGTGCTGCCGAGAATGCTGCACAACAGCAGACGTGAAAGTTGACGATTGATCATTGCGTGGCTCCTTGTGCCTGTTCCGTGGCGGGCGTTAATGTCGGGGCTGTCCAGTCTGGCAAGCGGGCGAAAAGCTGGCAGACGTGTTGCTcc                                                                                                 >  3:108284/1‑151 (MQ=255)
                                                   cgtcgtGCTGCCGAGAATGCTGCACAACAGCAGACGTGAAAGTTGACGATTGATCATTGCGTGGCTCCTTGTGCCTGTTCCGTGGCGGGCGTTAATGTCGGGGCTGTCCAGTCTGGCAAGCGGGCGAAAAGCTGGCAGACGTGTTGCTCCt                                                                                                <  8:445138/151‑1 (MQ=255)
                                                       gtGCTGCCGAGAATGCTGCACAACAGCAGACGTGAAAGTTGACGATTGATCATTGCGTGGCTCCTTGTGCCTGTTCCGTGGCGGGCGTTAATGTCGGGGCTGTCCAGTCTGGCAAGCGGGCGAAAAGCTGGCAGACGTGTTGCTCCTGttt                                                                                            <  8:1599770/151‑1 (MQ=255)
                                                       gtGCTGCCGAGAATGCTGCACAACAGCAGACGTGAAAGTTGACGATTGATCATTGCGTGGCTCCTTGTGCCTGTTCCGTGGCGGGCGTTAATGTCGGGGCTGTCCAGTCTGGCAAGCGGGCGAAAAGCTGGCAGACGTGTTGCTCCTGttt                                                                                            <  4:908882/151‑1 (MQ=255)
                                                        ggctgcCTAGAATGCGGCACAACAGCAGACGTGAAAGGTGACGATTGATCATTGCGTGGCTCCTTGTGCCTGTTCCGTGGCGGGCGTTAATGTCGGGGCTGTCCAGTCTGGCAAGCGGGCGAAAAGCTGGCAGACGTGTTGCTCCTGtttt                                                                                           <  2:60744/150‑1 (MQ=255)
                                                         gctgcCGAGAATGCTGCACAACAGCAGACGTGAAAGTTGACGATTGATCATTGCGTGGCTCCTTGTGCCTGTTCCGTGGCGGGCGTTAATGTCGGGGCTGTCCAGTCTGGCAAGCGGGCGAAAAGCTGGCAGACGTGTTGCTCCTGtttt                                                                                           >  1:873479/1‑150 (MQ=255)
                                                         gctgcCGAGAATGCTGCACAACAGCAGACGTGAAAGTTGACGATTGATCATTGCGTGGCTCCTTGTGCCTGTTCCGTGGCGGGCGTTAATGTCGGGGCTGTCCAGTCTGGCAAGCGGGCGAAAAGCTGGCAGACGTGTTGCTCCTGtttt                                                                                           >  2:850966/1‑150 (MQ=255)
                                                           ggcCGAGAATGCTGCACAACAGCAGACGTGAAAGTTGACGATTGATCATTGCGTGGCTCCTTGTGCCTGTTCCGTGGCGGGCGTTAATGTCGGGGCTGTCCAGTCTGGCAAGCGGGCGAAAAGCTGGCAGACGTGTTGCTCCTGTTTtaat                                                                                        <  4:1603131/150‑1 (MQ=255)
                                                            gcCGAGAATGCTGCACAACAGCAGACGTGAAAGTTGACGATTGATCATTGCGTGGCTGCTTGTGCCTGTTCCGTGGCGGGCGTTAATGTCGGGGCAGTCCAGTCTGGCAAGCGGGCGAAAAGCTGGCAGACGTGTTGCGCCTGTTTtaat                                                                                        >  6:141294/1‑150 (MQ=255)
                                                            gcCGAGAATGCTGCACAACAGCAGACGTGAAAGTTGACGATTGATCATTGCGTGGCTCCTTGTGCCTGTTCCGTGGCGGGCGTTAATGCCGGGGCTGTCCAGTCTGGCAAGCGGGCGAAAAGCTGGCAGACGTGTTGCTCCTGTTTtaata                                                                                       <  7:1495484/151‑1 (MQ=255)
                                                             ccGAGAATGCTGCACAACAGCAGACGTGAAAGTTGACGATTGATCATTGCGTGGCTCCTTGTGCCTGTTCCGTGGCGGGCGTTAATGTCGGGGCTGTCCAGTCTGGCAAGCGGGCGAAAAGCTGGCAGACGTGTTGCTCCTGTTTtaataa                                                                                      <  1:505794/151‑1 (MQ=255)
                                                             ccGAGAATGCTGCACAACAGCAGACGTGAAAGTTGACGATTGATCATTGCGTGGCTCCTTGTGCCTGTTCCGTGGCGGGCGTTAATGTCGGGGCTGTCCAGTCTGGCAAGCGGGCGAAAAGCTGGCAGACGTGTTGCTCCTGTTTtaata                                                                                       >  8:1149434/1‑150 (MQ=255)
                                                              cGAGAATGCTGCACAACAGCAGACGTGAAAGTTGACGATTGATCATTGCGTGGCTCCTTGTGCCTGTTCCGTGGCGGGCGTTAATGTCGGGGCTGTCCAGTCTGGCAAGCGGGCGAAAAGCTGGCAGACGTGTTGCTCCTGTTTTAATAAg                                                                                     <  2:1147499/151‑1 (MQ=255)
                                                                agaATGCTGCACAACAGCAGACGTGAAAGTTGACGATTGATCATTGCGTGGCTCCTTGTGCCTGTTCCGTGGCGGGCGTTAATGTCGGGGCTGTCCAGTCTGGCAAGCGGGCGAAAAGCTGGCAGACGTGTTGCTCCTGTTTTAATAAGGt                                                                                   >  3:161334/1‑151 (MQ=255)
                                                                agaATGCTGCACAACAGCAGACGTGAAAGTTGACGATTGATCATTGCGTGGCTCCTTGTGCCTGTTCCGTGGCGGGCGTTAATGTCGGGGCTGTCCAGTCTGGCAAGCGGGCGAAAAGCTGGCAGACGTGTTGCTCCTGTTTTAATAAGGt                                                                                   >  1:1562220/1‑151 (MQ=255)
                                                                  aaTGCTGCACAACAGCAGACGTGAAAGTTGACGATTGATCATTGCGTGGCTCCTTGTGCCTGTTCCGTGGCGGGCGTTAATGTCGGGGCTGTCCAGTCTGGCAAGCGGGCGAAAAGCTGGCAGACGTGTTGCTCCTGTTTTAATAAGGTgg                                                                                 >  7:33491/1‑151 (MQ=255)
                                                                   aTGCTGCACAACAGCAGACGTGAAAGTTGACGATTGATCATTGCGTGGCTCCTTGTGCCTGTTCCGTGGCGGGCGTTAATGTCGGGGCTGTCCAGTCTGGCAAGCGGGCGAAAAGCTGGCAGACGTGTTGCTCCTGTTTTAATAAGGTggg                                                                                >  6:456223/1‑151 (MQ=255)
                                                                     gctgcACAACAGCAGACGTGAAAGTTGACGATTGATCATTGCGTGGCTCCTTGTGCCTGTTCCGTGGCGGGCGTTAATGGCGGGGCTGGCCAGTCTGGCAAGCGGGCGAAAAGCTGGCAGACGTGTTGCTCCTGTTTTAATAAGGTGGGCa                                                                              <  7:973519/151‑1 (MQ=255)
                                                                     gctgcACAACAGCAGACGTGAAAGTTGACGATTGATCATTGCGTGGCTCCTTGTGCCTCTTCCGTGGCGGGCGTTAATGTCGGGGCTGTCCAGTCTGGCAAGCGGGCGAAAAGCTGGCAGACGTCTTGCTCCTGTTTTAATAAGGTGGGc                                                                               >  6:386561/1‑150 (MQ=255)
                                                                         cacaACAGCAGACGTGAAAGTTGACGATTGATCATTGCGTGGCTCCTTGTGCCTGTTCCGTGGCGGGCGTTAATGTCGGGGCTGTCCAGTCTGGCAAGCGGGCGAAAAGCTGGCAGACGTGTTGCTCCTGTTTTAATAAGGTGGGCAGGAt                                                                          <  5:642691/151‑1 (MQ=255)
                                                                          acaacaGCAGACGTGAAAGTTGACGATTGATCATTGCGTGGCGCCTTGTGCCTGTTCCGTGGCGGGCGTTAATGTCGGGGCTGTCCAGTCTGGCAAGCGGGCCAAAAGCTGGCAGACGTGTTGCTCCTGTTTTAATAAGGTGGGCAGGATg                                                                         <  6:1287971/151‑1 (MQ=255)
                                                                             acaGCAGACGTGAAAGTTGACGATTGATCATTGCGTGGCTCCTTGTGCCTGTTCCGTGGCGGGCGTTAATGTCGGGGCTGTCCAGTCTGGCAAGCGGGCGAAAAGCTGGCAGACGTGTTGCTCCTGTTTTAATAAGGTGGGCAGGATGGcc                                                                      >  1:765497/1‑151 (MQ=255)
                                                                                   gACGTGAAAGTTGACGATTGATCATTGCGTGGCTCCTTGTGCCTGTTCCGTGGCGGGCGTTAATGTCGGGGCTGTCCAGTCTGGCAAGCGGGCGAAAAGCTGGCAGACGTGTTGCTCCTGTTTTAATAAGGTGGGCAGGATGGCCTgag                                                                  >  7:230659/1‑149 (MQ=255)
                                                                                    aCGTGAAAGTTGACGATTGATCATTGCGTGGCTCCTTGTGCCTGTTCCGTGGCGGGCGTTAATGTCGGGGCTGTCCAGTCTTGCAAGCGGGCGAAAATCGGGCAGACGTGTTGCTCCTGTTTTAATAAGGTGGGCAGGATGGCCTGagag                                                                <  6:1232198/150‑1 (MQ=255)
                                                                                    aCGTGAAAGTTGACGATTGATCATTGCGTGGCTCCTTGTGCCTGTTCCGTGGCGGGCGTTAATGTCGGGGCTGTCCAGTCTGGCAAGCGGGCGAAAAGCTGGCAGACGTGTTGCTCCTGTTTTAATAAGGTGGGCAGGATGGCCTGAGAgg                                                               <  2:26092/151‑1 (MQ=255)
                                                                                     cGTGAAAGTTGACGATTGATCATTGCGTGGCTCCTTGTGCCTGTTCCGTGGCGGGCGTTAATGTCGGGGCTGTCCAGTCTGGCAAGCGGGCGAAAAGCTGGCAGACGTGTTGCTCCTGTTTTAATAAGGTGGGCAGGATGGCCTGAGCggg                                                              >  2:1497199/1‑150 (MQ=255)
                                                                                     cGTGAAAGTTGACGATTGATCATTGCGTGGCTCCTTGTGCCTGTTCCGTGGCGGGCGTTAATGTCGGGGCTGTCCAGTCTGGCAAGCGGGCGAAAAGCTGGCAGACGTGTTGCTCCTGTTTTAATAAGGTGGGCAGGATGGCCTGAGAGGc                                                              >  7:936756/1‑151 (MQ=255)
                                                                                       tGAAAGTTGACGATTGATCATTGAGTGGCTCATTGTGCCTGTTCCGTGGCGGGCGTTAATGTCGGGGCTGTCCAGTCTGGCAAGCGGCCGAAAAGCTGGCAGACGTGTTGCTCCTGTTTTAATAAGGTGGGCAGGATGGCCTGAGAGGCgg                                                            >  5:211558/1‑151 (MQ=255)
                                                                                         aaaGTTGACGATTGAACATTGCGTGGCTCCTTGTGCCTGTTCCGTGGCGGGCGATAATGTCGCGGCTCTCCAGTCTGGCAAGCGGGCGAAAAGCTGGCACACGTGTAGCTCCTGTAATAATAAGGTGCGCACGATGGCCTGAGAGGCGGgt                                                          >  7:398060/1‑151 (MQ=255)
                                                                                          aaGTTGACGATTGATCATTGCGTGGCTCCTTGTGCCTGTTCCGTGGCGGGCGTTAATGTCGGGGCTGTCCAGTCTGGCAAGCGGGCGAAAAGCTGGCAGACGTGTTGCTCCTGTTTTAATAAGGTGGGCAGGATGGCCTGAGAGGCGGgtg                                                         >  8:460752/1‑151 (MQ=255)
                                                                                               gACGATTGATCATTGCGTGGCTCCTTGTGCCTGTTCCGTGGCGGGCGTTAATGTCGGGGCTGTCCAGTCTGGCAAGCGGGCGAAAAGCTGGCAGACGTGTTGCTCCTGTTTTAATAAGGTGGGCAGGATGGCCTGAGAGGCGGGTGTctct                                                    >  6:1671023/1‑151 (MQ=255)
                                                                                                aCGATTGATCATTGCGTGGCTCCTTGTGCCTGTTCCGTGGCGGGCGTTAATGTCGGGGCTGTCCAGTCTGGCAAGCGGGCGAAAAGCTGGCAGACGTGTTGCTCCTGTTTTAATAAGGTGGGCAGGATGGCCTGAGAGGCGGGTGTCTCtt                                                   <  7:1544200/151‑1 (MQ=255)
                                                                                                aCGATTGATCATTGCGTGGCTCCTTGTGCCTGTTCCGTGGCGGGCGTTAATGTCGGGGCTGTCCAGTCTGGCAAGCGGGCGAAAAGCTGGCAGACGTGTTGCTCCTGTTTTAATAAGGTGGGCAGGATGGCCTGAGAGGCGGGTGTCTCtt                                                   >  4:962513/1‑151 (MQ=255)
                                                                                                 ggATTGATCATTGCGTGGCTCCTTGTGCCTGTTCCGTGGCGGGCGTTAATGTCGGGGCTGTCCAGTCTGGCAAGCGGGCGAAAAGCTGGCAGACGTGTTGCTCCTGTTTTAATAAGGTGGGCAGGATGGCCTGAGAGGCGGGTGTCTCTTc                                                  <  6:1482304/150‑1 (MQ=255)
                                                                                                 cGATTGATCATTGCGTGGCTCCTTGTGCCTGTTCCGTGGCGGGCGTTAATGTCGGGGCTGTCCAGTCTGGCAAGCGGGCGAAAAGCTGGCAGACGTGTTGCTCCTGTTTTAATAAGGTGGGCAGGATGGCCTGAGAGGCGGGTGTCTCTTc                                                  >  3:186763/1‑151 (MQ=255)
                                                                                                 cGATTGATCATTGCGTGGCTCCTTGTGCCTGTTCCGTGGCGGGCGTTAATGTCGGGGCTGTCCAGTCTGGCAAGCGGGCGAAAAGCTGGCAGACGTGTTGCTCCTGTTTTAATAAGGTGGGCAGGATGGCCTGAGAGGCGGGTGTCTCTTc                                                  >  5:1299437/1‑151 (MQ=255)
                                                                                                  gATTGATCATTGCGTGGCTCCTTGTGCCTGTTCCGTGGCGGGCGTTAATGCCGGGGCTGTCCAGTCTGGCAAGCGGGCGAAAAGCTGGCAGACGTGTTGCTCCTGTTTTAATAAGGTGGGCAGGATGGCCTGAGAGGCGGGTGTCTCTTc                                                  >  4:247657/1‑150 (MQ=255)
                                                                                                     tGATCATTGGGTGGCTCCTTGTGCCTGTTCCGTGGCGGGCGTTAATGTCGGGGCTGTCCAGTCTGGCAAGCGGGCGAAAAGCTGGCAGAGGTGTTGCTCCTGTTTTAATAAGGTGGGCAGGATGGCCTGAGAGGCGGGTGTCTCTTCAAg                                               <  6:827318/150‑1 (MQ=255)
                                                                                                     tGATCATTGCGTGGCTCCTTGTGCCTGTTCCGTGGCGGGCGTTAATGTCGGGGCTGTCCAGTCTGGCAAGCGGGCGAAAAGCTGGCAGACGTGTTGCTCCTGTTTTAATAAGGTGGGCAGGATGGCCTGAGAGGCGGGTGTCTCTTCAAGc                                              >  3:1380744/1‑151 (MQ=255)
                                                                                                      gATCATTGCGTGGCTCCTTGTGCCTGTTCCGTGGCGGGCGTTAATGTCGGGGCTGTCCAGTCTGGCAAGCGGGCGAAAAGCTGGCAGACGTGTTGCTCCTGTTTTAATAAGGTGGGCAGGATGGCCTGAGAGGCGGGTGTCTCTTCAAGCg                                             >  8:197491/1‑151 (MQ=255)
                                                                                                         cATTGCGTGGCTCCTTGTGCCTGTTCCGTGGCGGGCGTTAATGTCGGGGCTGTCCAGTCTGGCAAGCGGGCTAAAAGCTGGCAGACGTGTTGCTCCTGTTTTAATAAGGTGGGCAGGATGGCCTGAGAGGCGGGTGTCTCTTCAAGCGtt                                           <  4:1207303/150‑1 (MQ=255)
                                                                                                         cATTGCGTGGCTCCTTGTGCCTGTTCCGTGGCGGGCGTTAATGTCGGGGCTGTCCAGTCTGGCAAGCGGGCGAAAAGCTGGCAGACGTGTTGCTCCTGTTTTAATAAGGTGGTCAGGATGGCCTGAGAGGCGGGTGTCTCTTCAAGCGtt                                           >  1:370489/1‑150 (MQ=255)
                                                                                                           ttGCGTGGCTCCTTGTGCCTGTTCCGTGGCGGGCGTTAATGTCGGGGCTGTCCAGTCTGGCAAGCGGGCGAAAAGCTGGCAGACGTGTTGCTCCTGTTTTAATAAGGTGGGCAGGATGGCCTGAGAGGCGGGTGTCTCTTCAAGCGTTAtt                                        >  7:104655/1‑151 (MQ=255)
                                                                                                           ttGCGTGGCTCCTTGTGCCTGTTCCGTGGCGGGCGTTAATGTCGGGGCTGTCCAGTCTGGCAAGCGGGCGAAAAGCTGGCAGACGTGTTGCTCCTGTTTTAATAAGGTGGGCAGGATGGCCTGAGAGGCGGGTGTCTCTTCAAGCGTTAtt                                        >  7:1403897/1‑151 (MQ=255)
                                                                                                            tGCGTGGCTCCTTGTGCCTGTTCCGTGGCGGGCGTTAATGTCGGGGCTGTCCAGTCTGGCAAGCGGGCGAAAAGCTGGCAGACGTGTTGCTCCTGTTTTAATAAGGTGGGCAGGATGGCCTGAGAGGCGGGTGTCTCTTCAAGCGTTAttt                                       >  4:303411/1‑151 (MQ=255)
                                                                                                            tGCGTGGCTCCTTGTGCCTGTTCCGTGGCGGGCGTTAATGTCGGGGCTGTCCAGTCTGGCAAGCGGGCGAAAAGCTGGCAGACGTGTTGCTCCTGTTTTAATAAGGTGGGCAGGATGGCCTGAGAGGCGGGTGTCTCTTCAAGCGTTAttt                                       >  7:580009/1‑151 (MQ=255)
                                                                                                            tGCGTGGCTCCTTGTGCCTGTTCCGTGGCGGGCGTTAATGTCGGGGCTGTCCAGTCTGGCAAGCGGGCGAAAAGCTGGCAGACGTGTTGCTCCTGTTTTAATAAGGTGGGCAGGATGGCCTGAGAGGCGGGTGTCTCTTCAAGCGTTAtt                                        <  6:898241/150‑1 (MQ=255)
                                                                                                             gCGTGGCTCCTTGTGCCTGTTCCGTGGCGGGCGTTAATGTCGGGGCTGTCCAGTCTGGCAAGCGGGCGAAAAGCTGGCAGACGTGTTGCTCCTGTTTTAATAAGGTGGGCAGGATGGCCTGAGAGGCGGGTGTCTGTTCAAGCGTTATTTg                                      >  7:125658/1‑151 (MQ=255)
                                                                                                               gTGGGTCCGTGTGCCTGATCCGTGGCGGGCGTTAATGTCGGGGCTGGCCAGTCTCGGAAGCGGGCGAAAAGCTGGCAGACGTGTTGCTCCTGTTTTAATAAGGTGGGCAGGATGGCCTGAGAGGCGGGTGTCTCTTCAAGCGTTATTTgag                                    <  6:1092007/151‑1 (MQ=255)
                                                                                                                  gCTCCTTGTGCCTGTTCCGTGGCGGGCGTTAATGCCGGGGCTGTCCAGTCTGGCAAGCGGGCGAAAAGCTGGCAGACGTGTTGCTCCTGTTTTAATAAGGTGGGCAGGATGGCCTGAGAGGCGGGTGTCTCTTCAAGCGTTATTTgagg                                   >  4:614409/1‑149 (MQ=255)
                                                                                                                    tCCTTGTGCCTGTTCCGTGGCGGGCGTTAATGTCGGGGCTGTCCAGTCTGGCAAGCGGGCGAAAAGCTGGCAGACGTGTTGCTCCTGTTTTAATAAGGTGGGCAGGATGGCCTGAGAGGCGGGTGTCTCTTCAAGCGTTATTTGAGgaga                                >  5:816782/1‑150 (MQ=255)
                                                                                                                    tCCTTGTGCCTGTTCCGTGGCGGGCGTTAATGTCGGGGCTGTCCAGTCTGGCAAGCGGGCGAAAAGCTGGCAGACGTGTTGCTCCTGTTTTAATAAGGTGGGCAGGATGGCCTGAGAGGCGGGTGTCTCTTCAAGCGTTATTTGAGgaga                                >  7:174649/1‑150 (MQ=255)
                                                                                                                    gccTTGTGCCTGTTCCGTGGCGGGCGTTAATGTCGGGGCTGTCCAGTCTGGCAAGCGGGCGAAAAGCTGGCAGACGTGTTGCTCCTGTTTTAATAAGGTGGGCAGGATGGCCTGAGAGGCGGGTGTCTCTTCAAGCGTTATTTGAGGAGaa                               <  1:139278/150‑1 (MQ=255)
                                                                                                                     ccTTGTGCCTGTTCCGTGGCGGGCGTTAATGTCGGGGCTGTCCAGTCTGGCAAGCGGGCGAAAAGCTGGCAGACGTGTTGCTCCTGTTTTAATAAGGTGGGCAGGATGGCCTGAGAGGCGGGTGTCTCTTCAAGCGTTATTTGAGGAGAAc                              <  5:1156556/151‑1 (MQ=255)
                                                                                                                     ccTTGTGCCTGTTCCGTGGCGGGCGTTAATGTCGGGGCTGTCCAGTCTGCCAAGCGGGCGAAAAGCTGGCAGACGTGTTGCTCCTGTTTTAATAAGGTGGGCAGGATGGCCTGAGAGGCGGGTGTCTCTTCAAGCGTTATTTGAGGAGAAc                              <  6:143727/151‑1 (MQ=255)
                                                                                                                     ccTTGTGCCTGGTCCGTGGCGGGCGTTAATGTCGGGGCTGTCCAGTCTGGCAAGCGGGCGAAAAGCTGGCAGACGTGTTGCTCCTGTTTTAATAAGGTGGGCAGGATGGCCTGAGAGGCGGGTGTCTCTTCAAGCGTTATTTCAGGAGAAc                              <  6:218965/151‑1 (MQ=255)
                                                                                                                           gCCTGTTCCGTGGCGGGCGTTAATGTCGGGGCTGTCCAGTCTGGCAAGCGGGCGAAAAGCTGGCAGACGTGTTGCTCCTGTTTTAATAAGGTGGGCAGGATGGCCTGAGAGGCGGGTGTCTCTTCAAGCGTTATTTGAGGAGAACGGCTg                         <  1:1150724/150‑1 (MQ=255)
                                                                                                                            ccTGTTCCGTGGCGGGCGTTAATGTCGGGGCTGTCCAGTCTGGCAAGCGGGCGAAAAGCTGGCAGACGTGTTGCTCCTGTTTTAATAAGGTGGGCAGGATGGCCTGAGAGGCGGGTGTCTCTTCAAGCGTTATTTGAGGAGACCGGCTg                         >  2:315716/1‑149 (MQ=255)
                                                                                                                               gTTCCGTGGCGGGCGTTAATGTCGGGGCTGTCCAGTCTGGCAAGCGGGCGAAAAGCTGGCAGACGTGTTGCTCCTGTTTTAATAAGGTGGGCAGGATGGCCTGAGAGGCGGGTGTCTCTTCAAGCGTTATTTGAGGAGAACGGCTGAGTGc                    >  5:480029/1‑151 (MQ=255)
                                                                                                                                   cGTGGCGGGCGTTAATGTCGGGGCTGTCCAGTCTGGCAAGCGGGCGAAAAGCTGGCAGACGTGTTGCTCCTGTTTTAATAAGGTGGGCAGGATGGCCTGAGAGGCGGGTGTCTCTTCAAGCGTTATTTGAGGAGAACGGCTGAGTGCGGtt                >  3:348680/1‑151 (MQ=255)
                                                                                                                                     tGGCGGGCGTTAATGTCGGGGCTGTCCAGTCTGGCAAGCGGGCGAAAAGCTGGCAGACGTGTTGCTCCTGTTTTAATAAGGTGGGCAGGATGGCCTGAGAGGCGGGTGTCTCTTCAAGCGTTATTTGAGGAGAACGGCTGAGTGCGGTTGc              >  2:1393081/1‑151 (MQ=255)
                                                                                                                                      ggcgggcgTTAATGTCGGGGCTGTCCAGTCTGGCAAGCGGGCGAAAAGCTGGCAGACGTGTTGCTCCTGTTTTAATAAGGTGGGCAGGATGGCCTGAGAGGCGGGTGTCTCTTCAAGCGTTATTTGAGGAGAACGGCTGAGTGCGGTTGcc             >  2:687191/1‑151 (MQ=255)
                                                                                                                                          ggcgTTAATGTCGGGGCTGTCCAGTCTGGCAAGCGGGCGAAAAGCTGGCAGACGTGTTGCTCCTGTTTTAATAAGGTGGGCAGGATGGCCTGAGAGGCGGGTGTCTCTTCAAGCGTTATTTGAGGAGAACGGCTGAGTGCGGTTGGCAGAc         >  7:615428/1‑151 (MQ=255)
                                                                                                                                          ggcgTTAATGTCGGGGCTGTCCAGTCTGGCAAGCGGGCGAAAAGCTGGCAGACGTGTTGCTCCTGTTTTAATAAGGTGGGCAGGATGGCCTGAGAGGCGGGTGTCTCTTCAAGCGTTATTTGAGGAGAACGGCTGAGTGCGGTTGCCAGAc         <  2:547891/151‑1 (MQ=255)
                                                                                                                                            cgTTAATGTCGGGGCTGTACAGTCTGGCAAGCGGGCGAAAAGCTGGCAGACGTGTTGCTCCTGTTTTAATAAGGTGGGCAGGATGGCCTGAGAGGCGGGTGTCTCTTCAAGCGTTATTTGAGGAGAACGGCTGAGTGCGGTTGCCAGAccc       >  4:1428030/1‑150 (MQ=255)
                                                                                                                                             gTTAATGTCGGGGCTGTCCAGTCTGGCAAGCGGGCGAAAAGCTGGCAGACGTGTTGCTCCTGTTTTAATAAGGTGGGCAGGATGGCCTGAGAGGCGGGTGTCTCTTCAAGCGTTATTTGAGGAGAACGGCTGAGTGCGGTTGCCAGACCgg      <  5:1433044/151‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CTGTTGATGAGGTGCAGAATCCTTACGTACCAGGGCACAGCCGGAAATCAGCGTCGTGCTGCCGAGAATGCTGCACAACAGCAGACGTGAAAGTTGACGATTGATCATTGCGTGGCTCCTTGTGCCTGTTCCGTGGCGGGCGTTAATGCCGGGGCTGTCCAGTCTGGCAAGCGGGCGAAAAGCTGGCAGACGTGTTGCTCCTGTTTTAATAAGGTGGGCAGGATGGCCTGAGAGGCGGGTGTCTCTTCAAGCGTTATTTGAGGAGAACGGCTGAGTGCGGTTGCCAGACCGGCAGC  >  NC_000913/4300923‑4301218

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