Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,891,768 G→A A52T (GCA→ACA)  mdtL → multidrug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,891,7680GA100.0% 560.1 / NA 143A52T (GCA→ACA) mdtLmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (72/71);  total (72/71)

Only 100 of 143 total aligned reads displayed.

CCGCTTTTTGATTTGTAGTTTTGCCCTGGTTTTACTTTATCCCGCCGGGATTGATATGTACCTCGTTGGTTTACCGCGCATCGCCGCCGATCTCAATGCCAGCGAAGCGCAGTTGCATATTGCGTTCTCCGTATATCTGGCGGGGATGGCAGCTGCGATGTTATTTGCCGGTAAAGTGGCCGATCGTTCAGGGAGAAAGCCGGTCGCCATACCCGGCGCGGCGCTATTTATTATTGCCTCGGTGTTCTGTTCACTGGCTGAAACCAGCACGTTATTTCTTGCAGGCCGATTTCTAC  >  NC_000913/3891620‑3891915
                                                                                                                                                    |                                                                                                                                                   
   cTTTTTGATTTGTAGTTTTGCCCTGGTTTTACTTTATCCCGCCGGGATTGATATGTACCTCGTTGGTTTACCGCGCATCGCCGCCGATCTCAATGCCAGCGAAGCGCAGTTGCATATTGCGTTCTCCGTATATCTGGCGGGGATGACAGCt                                                                                                                                                >  4:417508/1‑151 (MQ=255)
     ttttGATTTGTAGTTTTGCCCTGGTTTTCCTTTAGCCCGCCGGGATTGATATGTACCTCGTTGGTTTACCGCGCATCGCCGCCGATCTCAATGCCAGCGAAGCGCAGTTGCATATTGCGTTCTCCGTATATCTGGCGGGGATGACAGCTg                                                                                                                                               <  1:1740859/150‑1 (MQ=255)
     ttttGATTTGTAGTTTTGCCCTGGTTTTACTTTATCCCGCCGGGATTGATATGTACCTCGTTGGTTTACCGCGCATCGCCGCCGATCTCAATGCCAGCGAAGCGCAGTTGCATATTGCGTTCTCCGTATATCTGGCGGGGATGACAGCt                                                                                                                                                <  4:399806/149‑1 (MQ=255)
     ttttGATTTGTAGTTTTGCCCTGGTTTTACTTTATCCCGCCGGGATTGATATGTACCTCGTTGGTTTACCGCGCATCGCCGCCGATCTCAATGCCAGCGAAGCGCAGTTGCATATTGCGTTCTCCGTATATCTGGCGGGGATGACAGCt                                                                                                                                                <  1:973990/149‑1 (MQ=255)
     ttttGATTTGTAGTTTTGCCCTGGTTTTACTTTATCCCGCCGGGATTGATATGTACCTCGTTGGTTTACCGCGCATCGCCGCCGATCTCAATGCCAGCGAAGCGCAGTTGCATATTGCGTTCTCCGTATATCTGGCGGGGATGACAGCTg                                                                                                                                               >  4:706154/1‑150 (MQ=255)
     ttttGATTTGTAGTTTTGCCCTGGTTTTACTTTATCCCGCCGGGATTGATATGTACCTCGTTGGTTTACCGCGCATCGCCGCCGATCTCAATGCCAGCGAAGCGCAGTTGCATATTGCGTTCTCCGTATATCTGGCGGGGATGACAGCTGc                                                                                                                                              >  4:714015/1‑151 (MQ=255)
     ttttGATTTGTAGTTTTGCCCTGGTTTTACTTTATCCCGCCGGGATTGATATGTACCTCGTTGGTTTACCGCGCATCGCCGCCGATCTCAATGCCAGCGAAGCGCAGTTGCATATTGCGTTCTCCGTATATCTGGCGGGGATGACAGCTGc                                                                                                                                              >  3:46008/1‑151 (MQ=255)
      tttGATTTGTAGTTTTGCCCTGGTTTTACTTTATCCCGCCGGGATTGATATGTACCTCGTTGGTTTACCGCGCATCGCCGCCGATCTCAATGCCAGCGAAGCGCAGTTGCATATTGCGTTCTCCGTATATCTGGCGGGGATGACAGCTGCg                                                                                                                                             >  2:367161/1‑151 (MQ=255)
      tttGATTTGTAGTTTTGCCCTGGTTTTACTTTATCCCGCCGGGATTGATATGTACCTCGTTGGTTTACCGCGCATCGCCGCCGATCTCAATGCCAGCGAAGCGCAGTTGCATATTGCGTTCTCCGTATATCTGGCGGGGATGACAGCTGCg                                                                                                                                             >  4:1169902/1‑151 (MQ=255)
       ttGATTTGTAGTTTTGCCCTGGTTTTACTTTATCCCGCCGGGATTGATATGTACCTCGTTGGTTTACCGCGCATCGCCGCCGATCTCAATGCCAGCGAAGCGCAGTTGCATATTGCGTTCTCCGTATATCTGGCGGGGATGACAGCTGCg                                                                                                                                             <  4:169338/150‑1 (MQ=255)
        tGATTTGTAGTTTTGCCCTGGTTTTACTTTATCCCGCCGGGATTGATATGTACCTCGTTGGTTTACCGCGCATCGCCGCCGATCTCAATGCCAGCGAAGCGCAGTTGCATATTGCGTTCTCCGTATATCTGGCGGGGATGACTGCTGCg                                                                                                                                             >  1:1883092/1‑149 (MQ=255)
        tGATTTGTAGTTTTGCCCTGGTTTTACTTTATCCCGCCGGGATTGATATGTACCTCGTTGGTTTACCGCGCATCGCCGCCGATCTCAATGCCAGCGAAGCGCAGTTGCATATTGCGTTCTCCGTATATCTGGCGGGGATGACAGCTGCGAt                                                                                                                                           >  1:1579323/1‑151 (MQ=255)
          aTTTGTAGTTTTGCCCTGGTTTTACTTTATCCCGCCGGGATTGATATGTACCTCGTTGGTTTACCGCGCATCGCCGCCGATCTCAATGCCAGCGAAGCGCAGTTGCATATTGCGTTCTCCGTATATCTGGCGGGGATGACAGCTGCGATGt                                                                                                                                         <  4:543591/151‑1 (MQ=255)
            ttGTAGTTTTGCCCTGGTTTTACTTTATCCCGCCGGGATTGATATGTACCTCGTTGGTTTACCGCGCATCGCCGCCGATCTCAATGCCAGCGAAGCGCAGTTGCATATTGCGTTCTCCGTATATCTGGCGGGGATGACAGCTGCGATGtt                                                                                                                                        <  4:489443/150‑1 (MQ=255)
             tGTAGTTTTGCCCTGGTTTTACTTTATCCCGCCGGGATTGATATGTACCTCGTTGGTTTACCGCGCATCGCCGCCGATCTCAATGCCAGCGAAGCGCAGTTGCATATTGCGTTCTCCGTATATCTGGCGGGGATGACAGCTGCGATGTTa                                                                                                                                       <  2:353291/150‑1 (MQ=255)
             tGTAGTTTTGCCCTGGTTTTACTTTATCCCGCCGGGATTGATATGTACCTCGTTGGTTTACCGCGCATCGCCGCCGATCTCAATGCCAGCGAAGCGCAGTTGCATATTGCGTTCTCCGTATATCTGGCGGGGATGACAGCTGCGATGTTAt                                                                                                                                      >  1:1575610/1‑151 (MQ=255)
              gTAGTTTTGCCCTGGTTTTACTTTATCCCGCCGGGATTGATATGTACCTCGTTGGTTTACCGCGCATCGCCGCCGATCTCAATGCCAGCGAAGCGCAGTTGCATATTGCGTTCTCCGTATATCTGGCGGGGATGACAGCTGCGATGTTAt                                                                                                                                      <  2:157542/150‑1 (MQ=255)
               tAGTTTTGCCTTGGTTTTACTTTATCCCGCCGGGATTGATATGTACCTCGTTGGTTTACCGCGCATCGCCGCCGATCTCAATGCCAGCGAAGCGCAGTTGCATATTGCGTTCTCCGTATATCTGGCGGGGATGACAGCTGCGATGTTAtt                                                                                                                                     <  1:1621579/150‑1 (MQ=255)
                aGTTTTGCCCTGGTTTTACTTTATCCCGCCGGGATTGATATGTACCTCGTTGGTTTACCGCGCATCGCCGCCGATCTCAATGCCAGCGAAGCGCAGTTGCATATTGCGTTCTCCGTATATCTGGCGGGGATGACAGCTGCGATGTTATTTg                                                                                                                                   >  1:1222195/1‑151 (MQ=255)
                  ttttGCCCTGGTTTTCCTTTACCCCGCCGGGATTGATATGTACCTCGTTGGTTTACCGCGCATCGCCGCCGATCTCAATGCCAGCGAAGCGCAGTTGCATATTGCGTTCTCCGTATATCTGGCGGGGATGACAGCTGCGATGTTATTTg                                                                                                                                   <  1:1081085/149‑1 (MQ=255)
                  ttttGCCCTGGTTTTACTTTATCCCGCCGGGATTGATATGTACCTCGTTGGTTTACCGCGCATCGCCGCCGATCTCAATGCCAGCGAAGCGCAGTTGCATATTGCGTTCTCCGTATATCTGGCGGGGATGACAGCTGCGATGTTATTTGc                                                                                                                                  <  1:1548098/150‑1 (MQ=255)
                   tttGCCCTGGTTTTCCTTTATCCCGCCGGGATTGATATGTACCTCGTTGGTTTACCGCGCATCGCCGCCGATCTCAATGCCAGCGAAGCGCAGTTGCATATTGCGTTCTCCGTATATCTGGCGGGGATGACAGCTGCGATGTTATTTGcc                                                                                                                                 <  3:211010/150‑1 (MQ=255)
                     tGCCCTGGTTTTACTTTATCCCGCCGGGATTGATATGTACCTCGTTGGTTTACCGCGCATCGCCGCCGATCTCAATGCCAGCGAAGCGCAGTTGCATATTGCGTTCTCCGTATATCTGGCGGGGATGACAGCTGCGATGTTATTTGCCgg                                                                                                                               >  1:1643491/1‑150 (MQ=255)
                        ccTGGTTTTACTTTATCCCGCCGGGATTGATATGTACCTCGTTGGTTTACCGCGCATCGCCGCCGATCTCAATGCCAGCGAAGCGCAGTTGCATATTGCGTTCTCCGTATATCTGGCGGGGATGACAGCTGCGATGTTATTTGCCGGTaaa                                                                                                                           <  2:70407/151‑1 (MQ=255)
                         cTGGTTTTACTTTATCCCGCCGGGATTGATATGTACCTCGTTGGTTTACCGCGCATCGCCGCCGATCTCAATGCCAGCGAAGCGCAGTTGCATATTGCGTTCTCCGTATATCTGGCGGGGATGACAGCTGCGATGTTATTTGCCGGTaa                                                                                                                            >  1:271017/1‑149 (MQ=255)
                         cTGGTTTTACTTTATCCCGCCGGGATTGATATGTACCTCGTTGGTTTACCGCGCATCGCCGCCGATCTCAATGCCAGCGAAGCGCAGTTGCATATTGCGTTCTCCGTATATCTGGCGGGGATGACAGCTGCGATGTTATTTGCCGGTaa                                                                                                                            >  4:1552082/1‑149 (MQ=255)
                          tGGTTTTACTTTATCCCGCCGGGATTGATATGTACCTCGTTGGTTTACCGCGCATCGCCGCCGATCTCAATGCCAGCGAAGCGCAGTTGCATATTGCGTTCTCCGTATATCTGGCGGGGATGACAGCTGCGATGTTATTTGCCGGTAAAGt                                                                                                                         >  4:920816/1‑151 (MQ=255)
                           ggTTTTACTTTATCCCGCCGGGATTGATATGTACCTCGTTGGTTTACCGCGCATCGCCGCCGATCTCAATGCCAGCGAAGCGCAGTTGCATATTGCGTTCTCCGTATATCTGGCGGGGATGACAGCTGCGATGTTATTTGCCGGTAAAGTg                                                                                                                        >  4:1627057/1‑151 (MQ=255)
                                 aCTTTATCCCGCCGGGATTGATATGTACCTCGTTGGTTTACCGCGCATCGCCGCCGATCTCAATGCCAGCGAAGCGCAGTTGCATATTGCGTTCTCCGTATATCTGGCGGGGATGACAGCTGCGATGTTATTTGCCGGTAAAGTGGCCGa                                                                                                                   <  4:1644283/150‑1 (MQ=255)
                                 aCTTTATCCCGCCGGGATTGATATGTACCTCGTTGGTTTACCGCGCATCGCCGCCGATCTCAATGCCAGCGAAGCGCAGTTGCATATTGCGTTCTCCGTATATCTGGCGGGGATGACAGCTGCGATGTTATTTGCCGGTAAAGTGGCCGAt                                                                                                                  >  4:265754/1‑151 (MQ=255)
                                   tttATCCCGCCGGGATTGATATGTACCTCGTTGGTTTACCGCGCATCGCCGCCGATCTCAATGCCAGCGAAGCGCAGTTGCATATTGCGTTCTCCGTATATCTGGCGGGGATGACAGCTGCGATGTTATTTGCCGGTAAAGTGGCCGATCg                                                                                                                >  1:1266379/1‑151 (MQ=255)
                                    ttATCCCGCCGGGATTGATATGTACCTCGTTGGTTTACCGCGCATCGCCGCCGATCTCAATGCCAGCGAAGCGCAGTTGCATATTGCGTTCTCCGTATATCTGGCGGGGATGACAGCTGCGATGTTATTTGCCGGTAAAGTGGCCGATc                                                                                                                 <  3:440375/149‑1 (MQ=255)
                                    ttATCCCGCCGGGATTGATATGTACCTCGTTGGTTTACCGCGCATCGCCGCCGATCTCAATGCCAGCGAAGCGCAGTTGCATATTGCGTTCTCCGTATATCTGGCGGGGATGACAGCTGCGATGTTATTTGCCGGTAAAGTGGCCGATCGt                                                                                                               >  4:1401291/1‑151 (MQ=255)
                                         ccgccgGGATTGATATGTACCTCGTTGGTTTACCGCGCATCGCCGCCGATCTCAATGCCAGCGAAGCGCAGTTGCATATTGCGTTCTCCGTATATCTGGCGGGGATGACAGCTGCGATGTTATTTGCCGGTAAAGTGGCCGATCGTTCAgg                                                                                                          <  1:175932/151‑1 (MQ=255)
                                         ccgccgGGATTGATATGTACCTCGTTGGTTTACCGCGCATCGCCGCCGATCTCAATGCCAGCGAAGCGCAGTTGCATATTGCGTTCTCCGTATATCTGGCGGGGATGACAGCTGCGATGTTATTTGCCGGTAAAGTGGCCGATCGTTCAgg                                                                                                          >  3:96826/1‑151 (MQ=255)
                                          cgccgGGATTGATATGTACCTCGTTGGTTTACCGCGCATCGCCGCCGATCTCAATGCCAGCGAAGCGCAGTTGCATATTGCGTTCTCCGTATATCTGGCGGGGATGACAGCTGCGATGTTATTTGCCTGTAAAGTGGCCGATCGTTCAggg                                                                                                         <  1:1410672/151‑1 (MQ=255)
                                           gccgGGATTGATATGTACCTCGTTGGTTTACCGCGCATCGCCGCCGATCTCAATGCCAGCGAAGCGCAGTTGCATATTGCGTTCTCCGTATATCTGGCGGGGATGACAGCTGCGATGTTATTTGCCGGTAAAGTGGCCGATCGTTCAggg                                                                                                         >  4:447531/1‑150 (MQ=255)
                                             cgGGATTGATATGTACCTCGTTGGTTTACCGCGCATCGCCGCCGATCTCAATGCCAGCGAAGCGCAGTTGCATATTGCGTTCTCCGTATATCTGGCGGGGATGACAGCTGCGATGTTATTTGCCGGTAAAGTGGCCGATCGTTCAggg                                                                                                         >  4:1491185/1‑148 (MQ=255)
                                              gggATTGATATGTACCTCGTTGGTTTACCGCGCATCGCCGCCGATCTCAATGCCAGCGAAGCGCAGTTGCATATTGCGTTCTCCGTATATCTGGCGGGGATGACAGCTGCGATGTTATTTGCCGGTAAAGTGGCCGATCGTTCAGGGAGaa                                                                                                     <  4:819549/151‑1 (MQ=255)
                                               ggATTGATATGTACCTCGTTGGTTTACCGCGCATCGCCGCCGATCTCAATGCCAGCGAAGCGCAGTTGCATATTGCGTTCTCCGTATATCTGGCGGGGATGACAGCTGCGATGTTATTTGCCGGTAAAGTGGCCGATCGTTCAGGGAGaaa                                                                                                    <  3:273553/151‑1 (MQ=255)
                                                gATTGATATGTACCTCGTTGGTTTACCGCGCATCGCCGCCGATCTCAATGCCAGCGAAGCGCAGTTGCATATTGCGTTCTCCGTATATCTGGCGGGGATGACAGCTGCGATGTTATTTGCCGGTAAAGTGGCCGATCGTTCAGGGAGaa                                                                                                     >  4:1353025/1‑149 (MQ=255)
                                                    gATATGTACCTCGTTGGTTTACCGCGCATCGCCGCCGATCTCAATGCCAGCGAAGCGCAGTTGCATATTGCGTTCTCCGTATATCTGGCGGGGATGACAGCTGCGATGTTATTTGCCGGTAAAGTGGCCGATCGTTCAGGGAGAAAGCCgg                                                                                               <  4:1368850/151‑1 (MQ=255)
                                                     atatCTACCTCGTTGGTTTACCGCGAATCGCCGCCGATCTCATTGCCAGCGAAGCGCAGTTGCATATTGCGTTCTCCGTATATCTGGCGGGGATGACAGCTGCGATGATATTAGCCGGTAAAGTGGCCGATCGTTCCGGGAGAATGCCTGt                                                                                              >  1:1824518/1‑151 (MQ=255)
                                                      tatGTACCTCGTTGGTTTACCGCGCATCGCCGCCGATCTCAATGCCAGCGAAGCGCAGTTGCATATTGCGTTCTCCGTATATCTGGCGGGGATGACAGCTGCGATGTTATTTGCCGGTAAAGTGGCCGATCGTTCAGGGAGAAAGCCGGt                                                                                              <  1:1047696/150‑1 (MQ=255)
                                                       atGTACCTCGTTGGTTTACCGCGCATCGCCGCCGATCTCAATGCCAGCGAAGCGCAGTTGCATATTGCGTTCTCCGTATATCTGGCGGGGATGACAGCTGCGATGTTATTTGCCGGTAAAGTGGCCGATCGTTCAGGGAGAAAGCCGGTc                                                                                             <  1:1341599/150‑1 (MQ=255)
                                                       atGTACCTCGTTGGTTTACCGCGCATCGCCGCCGATCTCAATGCCAGCGAAGCGCAGTTGCATATTGCGTTCTCCGTATATCTGGCGGGGATGACAGCTGCGATGTTATTTGCCGGTAAAGTGGCCGATCGTTCAGGGAGAAAGCCGGTCg                                                                                            >  3:232096/1‑151 (MQ=255)
                                                        tGTACCTCGTTGGTTTACCGGGCACCGCCGCCGATCTCAATGCCAGCGAAGCGCAGTTGCATATTGCGTTCGCCGTATATCTGGCGGGGATGACAGCTGCGATGTTATTTGCCGGTAAAGTGGCCGATCGTTCAGGGAGAAAGCCGGTc                                                                                             <  1:1529884/149‑1 (MQ=255)
                                                            ccTCGTTGGTTTACCGCGCATCGCCGCCGATCTCAATGCCAGCGAAGCGCAGTTGCATATTGCGTTCTCCGTATATCTGGCGGGGATGACAGCTGCGATGTTATTTGCCGGTAAAGTGGCCGATCGTTCAGGGAGAAAGCCGGTCGCCAt                                                                                        <  4:1303064/150‑1 (MQ=255)
                                                                gTTGGTTTACCGCGCATCGCCGCCGATCTCAATGCCAGCGAAGCGCAGTTGCATATTGCGTTCTCCGTATATCTGGCGGGGATGACAGCTGCGATGTTATTTGCCGGTAAAGTGGCCGATCGTTCAGGGAGAAAGCCGGTCGCCATACCcg                                                                                   >  4:188487/1‑151 (MQ=255)
                                                                gTTGGTTTACCGCGCATCGCCGCCGATCTCAATGCCAGCGAAGCGCAGTTGCATATTGCGTTCTCCGTATATCTGGCGGGGATGACAGCTGCGATGTTATTTGCCGGTAAAGTGGCCGATCGTTCAGGGAGAAAGCCGGTCGCCATACCcg                                                                                   >  1:1682537/1‑151 (MQ=255)
                                                                  tGGTTTACCGCGCATCGCCGCCGATCTCAATGCCAGCGAAGCGCAGTTGCATATTGCGTTCTCCGTATATCTGGCGGGGATGACAGCTGCGATGTTATTTGCCGGTAAAGTGGCCGATCGTTCAGGGAGAAAGCCGGTCGCCATACCcgg                                                                                  <  4:775541/150‑1 (MQ=255)
                                                                   ggTTTACCGCGCATCGCCGCCGATCTCAATGCCAGCGAAGCGCAGTTGCATATTGCGTTCTCCGTATATCTGGCGGGGATGACAGCTGCGATGTTATTTGCCGGTAAAGTGGCCGATCGTTCAGGGAGAAAGCCGGTCGCCATACCcggcg                                                                                >  1:323115/1‑151 (MQ=255)
                                                                   ggTTTACCGCGCATCGCCGCCGATCTCAATGCCAGCGAAGCGCAGTTGCATATTGCGTTCTCCGTATATCTGGCGGGGATGACAGCTGCGATGTTATTTGCCGGTAAAGTGGCCGATCGTTCAGGGAGAAAGCCGGTCGCCATACCcggc                                                                                 <  3:499291/150‑1 (MQ=255)
                                                                    gTTTACCGCGCATCGCCGCCGATCTCAATGCCAGCGAAGCGCAGTTGCATATTGCGTTCTCCGTATATCTGGCGGGGATGACAGCTGCGATGTTATTTGCCGGTAAAGTGGCCGATCGTTCAGGGAGAAAGCCGGTCGCCATACCcggcgc                                                                               >  4:723443/1‑151 (MQ=255)
                                                                      ttACCGCGCATCGCCGCCGATCTCAATGCCAGCGAAGCGCAGTTGCATATTGCGTTCTCCGTATATCTGGCGGGGATGACAGCTGCGATGTTATTTGCCGGTAAAGTGGCCGATCGTTCAGGGAGAAAGCCGGTCGCCATACCcggcgcgg                                                                             >  2:422381/1‑151 (MQ=255)
                                                                         ccGCGCATCGCCGCCGATCTCAATGCCAGCGAAGCGCAGTTGCATATTGCGTTCTCCGTATATCTGGCGGGGATGACAGCTGCGATGTTATTTGCCGGTAAAGTGGCCGATCGTTCAGGGAGAAAGCCGGTCGCCATACCcggcgcggcg                                                                           <  1:555468/150‑1 (MQ=255)
                                                                          cgcgCATCGCCGCCGATCTCAATGCCAGCGAAGCGCAGTTGCATATTGCGTTCTCCGTATATCTGGCGGGGATGACAGCTGCGATGTTATTTGCCGGTAAAGTGGCCGATCGTTCAGGGAGAAAGCCGGTCGCCATACCCGGCGCGGCGCt                                                                         >  4:782664/1‑151 (MQ=255)
                                                                           gcgcATCGCCGCCGATCTCAATGCCAGCGAAGCGCAGTTGCATATTGCGTTCTCCGTATATCTGGCGGGGATGACAGCTGCGATGTTATTTGCCGGTAAAGTGGCCGATCGTTCAGGGAGAAAGCCGGTCGCCATACCCGGCGCGGCGCta                                                                        <  1:155584/151‑1 (MQ=255)
                                                                            cgcATCGCCGCCGATCTCAATGCCAGCGAAGCGCAGTTGCATATTGCGTTCTCCGTATATCTGGCGGGGATGACAGCTGCGATGTTATTTGCCGGTAAAGTGGCCGATCGTTCAGGGAGAAAGCCGGTCGCCATACCCGGCGCGGCGCta                                                                        <  1:1730388/150‑1 (MQ=255)
                                                                               aTCGCCGCCGATCTCAATGCCAGCGAAGCGCAGTTGCATATTGCGTTCTCCGTATATCTGGCGGGGATGACAGCTGCGATGTTATTTGCCGGTAAAGTGGCCGATCGTTCAGGGAGAAAGCCGGTCGCCATACCCGGCGCGGCGCtattt                                                                     >  4:118515/1‑150 (MQ=255)
                                                                               aTCGCCGCCGATCTCAATGCCAGCGAAGCGCAGTTGCATATTGCGTTCTCCGTATATCTGGCGGGGATGACAGCTGCGATGTTATTTGCCGGTAAAGTGGCCGATCGTTCAGGGAGAAAGCCGGTCGCCATACCCGGCGCGGCGCtattt                                                                     >  2:18740/1‑150 (MQ=255)
                                                                                  gccgccGATCTCAATGCCAGCGAAGCGCAGTTGCATATTGCGTTCTCCGTATATCTGGCGGGGATGACAGCTGCGATGTTATTTGCCGGTAAAGTGGCCGATCGTTCAGGGAGAAAGCCGGTCGCCATACCCGGCGCGGCGCtatttatt                                                                  >  4:893818/1‑150 (MQ=255)
                                                                                   ccgccgATCTCAATGCCAGCGAAGCGCAGTTGCATATTGCGTTCTCCGTATATCTGGCGGGGATGACAGCTGCGATGTTATTTGCCGGTAAAGTGGCCGATCGTTCAGGGAGAAAGCCGGTCGCCATACCCGGCGCGGCGCTATTtattat                                                                >  2:245032/1‑151 (MQ=255)
                                                                                   ccgccgATCTCAATGCCAGCGAAGCGCAGTTGCATATTGCGTTCTCCGTATATCTGGCGGGGATGACAGCTGCGATGTTATTTGCCGGTAAAGTGGCCGATCGTTCAGGGAGAAAGCCGGTCGCCATACCCGGCGCGGCGCTATTtattat                                                                <  2:401441/151‑1 (MQ=255)
                                                                                    cgccgATCTCAATGCCAGCGAAGCGCAGTTGCATATTGCGTTCTCCGTATATCTGGCGGGGATGACAGCTGCGATGTTATTTGCCGGTAAAGTGGCCGATCGTTCAGGGAGAAAGCCGGTCGCCATACCCGGCGCGGCGCTATTTattatt                                                               <  4:34654/151‑1 (MQ=255)
                                                                                      ccgATCTCAATGCCAGCGAAGCGCAGTTGCATATTGCGTTCTCCGTATATCTGGCGGGGATGACAGCTGCGATGTTATTTGCCGGTAAAGTGGCCGATCGTTCAGGGAGAAAGCCGGTCGCCATACCCGGCGCGGCGCTATTTATTATTGc                                                             >  1:730388/1‑151 (MQ=255)
                                                                                       cgATCTCAATGCCAGCGAAGCGCAGTTGCATATTGCGTTCTCCGTATATCTGGCGGGGATGACAGCTGCGATGTTATTTGCCGGTAAAGTGGCCGATCGTTCAGGGAGAAAGCCGGTCGCCATACCCGGCGCGGCGCTATTTATTATTGcc                                                            >  4:1152858/1‑151 (MQ=255)
                                                                                        gATCTCAATGACAGCGAAGCGCAGTTGCATATTGCGTTCTCCGTATATCTGGCGGGGATGACAGCTGCGATGTTATTTGCCGGTAAAGTGGCCGATCGTTCAGGGAGAAAGCCGGTAGCCATACCCGGAGCGGCGCTATTTATTATTGCCt                                                           >  3:13283/1‑151 (MQ=255)
                                                                                         aTCTCAATGCCAGCGCAGCGCAGTTGCATATTGCGTTCTCCGTATATCTGGCGGGGGTGAAAGCTGCGGTGTTATTTGCCGGCAAAAAGGCCGATCGTTCAGGGGGAAAACCGGCCAGCATACCCGTCGCGGCGCGATTTATTATTGCCTc                                                          >  3:399141/1‑151 (MQ=255)
                                                                                           ctcAATGCCAGCGAAGCGCAGTTGCATATTGCGTTCTCCGTATATCTGGCGGGGATGACAGCTGCGATGTTATTTGCCGGTAAAGTGGCCGATCGTTCAGGGAGAAAGCCGGTCGCCATACCCGGCGCGGCGCTATTTATTATTGCCTCgg                                                        >  2:540812/1‑151 (MQ=255)
                                                                                             cAATGCCAGCGAAGCGCAGTTGCATATTGCGTTCTCCGTATATCTGGCGGGGATGACAGCTGCGATGTTATTTGCCGGTAAAGTGGCCGATCGTTCAGGGAGAAAGCCGGTCGCCATACCCGGCGCGGCGCTATTTATTATTGc                                                             <  4:981151/144‑1 (MQ=255)
                                                                                             cAATGCCAGCGAAGCGCAGTTGCATATTGCGTTCTCCGTATATCTGGCGGGGATGACAGCTGCGATGTTATTTGCCGGTAAAGTGGCCGATCGTTCAGGGAGAAAGCCGGTCGCCATACCCGGCGCGGCGCTATTTATTATTGCCTc                                                          <  4:620698/147‑1 (MQ=255)
                                                                                             cAATGCCAGCGAAGCGCAGTTGCATATTGCGTTCTCCGTATATCTGGCGGGGATGACAGCTGCGATGTTATTTGCCGGTAAAGTGGCCGATCGTTCAGGGAGAAAGCCGGTCGCCATACCCGGCGCGGCGCTATTTATTATTGCCTCgg                                                        <  1:845383/149‑1 (MQ=255)
                                                                                             cAATGCCAGCGAAGCGCAGTTGCATATTGCGTTCTCCGTATATCTGGCGGGGATGACAGCTGCGATGTTATTTGCCGGTAAAGTGGCCGATCGTTCAGGGAGAAAGCCGGTCGCCATACCCGGCGCGGCGCTATTTATTATTGCCTCGgt                                                       <  2:477723/150‑1 (MQ=255)
                                                                                              aaTGCCAGCGAAGCGCAGTTGCATATTGCGTTCTCCGTATATCTGGCGGGGATGACAGCTGCGATGTTATTTGCCGGTAAAGTGGCCGATCGTTCAGGGAGAAAGCCGGTCGCCATACCCGGCGCGGCGCTATTTATTATTGCCTCGgtgt                                                     <  4:1883819/151‑1 (MQ=255)
                                                                                              aaTGCCAGCGAAGCGCAGTTGCATATTGCGTTCTCCGTATATCTGGCGGGGATGACAGCTGCGATGTTATTTGCCGGTAAAGTGGCCGATCGTTCAGGGAGAAAGCCGGTCGCCATACCCGGCGCGGCGCTATTTATTATTGCCTCGgtgt                                                     >  4:1411518/1‑151 (MQ=255)
                                                                                              aaTGCCAGCGAAGCGCAGTTGCATATTGCGTTCTCCGGATATCTGTCGGGGATGACAGCTGCGGTGTTATTTGCCGGTAAAGTGGCCGATCGGTCCGGGAGAAAGCCGGTCGCCATACCCGGCGCAGCGCTATTTATTATTGCCTCGgcgt                                                     >  3:309835/1‑151 (MQ=255)
                                                                                               aTGCCAGCGAAGCGCAGTTGCATATTGCGTTCTCCGTATATCTGGCGGGGATGACAGCTGCGATGTTATTTGCCGGTAAAGTGGCCGATCGTTCAGGGAGAAAGCCGGTCGCCATACCCGGCGCGGCGCTATTTATTATTGCCTCGGtgtt                                                    <  1:425798/151‑1 (MQ=255)
                                                                                                 gCCAGCGAAGCGCAGTTGCATATTGCGTTCTCCGTATATCTGGCGGGGATGACAGCTGCGATGTTATTTGCCGGTAAAGTGGCCGATCGTTCAGGGAGAAAGCCGGTCGCCATACCCGGCGCGGCGCTATTTATTATTGCCTCGGtgttct                                                  >  1:555998/1‑151 (MQ=255)
                                                                                                     gCGAAGCGCAGTTGCATATTGCGTTCTCCGTATATCTGGCGGGGATGACAGCTGCGATGTTATTTGCCGGTAAAGTGGCCGATCGTTCAGGGAGAAAGCCGGTCGCCATACCCGGCGCGGCGCTATTTATTATTGCCTCGGtgttctgtt                                               <  1:674475/150‑1 (MQ=255)
                                                                                                        aaGCGCAGTTGCATATTGCGTTCTCCGTATATCTGGCGGGGATGACAGCTGCGATGTTATTTGCCGGTAAAGTGGCCGATCGTTCAGGGAGAAAGCCGGTCGCCATACCCGGCGCGGCGCTATTTATTATTGCCTCGGTGTTCTGTTCAc                                            >  1:1022047/1‑150 (MQ=255)
                                                                                                            gcAGTTGCATATTGCGTTCTCCGTATATCTGGCGGGGATGACAGCTGCGATGTTATTTGCCGGTAAAGTGGCCGATCGTTCAGGGAGAAAGCCGGTCGCCATACCCGGCGCGGCGCTATTTATTATTGCCTCGGTGTTCTGTTCACTGGCt                                       >  1:850456/1‑151 (MQ=255)
                                                                                                             cAGTTGCATATTGCGTTCTCCGTATATCTGGCGGGGATGACAGCTGCGATGTTATTTGCCGGTAAAGTGGCCGATCGTTCAGGGAGAAAGCCGGTCGCCATACCCGGCGCGGCGCTATTTATTATTGCCTCGGTGTTCTGTTCACTGGCTg                                      >  4:1306100/1‑151 (MQ=255)
                                                                                                             cAGTTGCATATTGCGTTCTCCGTATATCTGGCGGGGATGACAGCTGCGATGTTATTTGCCGGTAAAGTGGCCGATCGTTCAGGGAGAAAGCCGGTCGCCATACCCGGCGCGGCGCTATTTATTATTGCCTCGGTGTTCTGTTCACTGGCTg                                      >  1:707020/1‑151 (MQ=255)
                                                                                                              aGTTGCATATTGCGTTCTCCGTATATCTGGCGGGGATGACAGCTGCGATGTTATTTGCCGGTAAAGTGGCCGATCGTTCAGGGAGAAAGCCGGTCGCCATACCCGGCGCGGCGCTATTTATTATTGCCTCGGTGTTCTGTTCACTGGCTg                                      <  4:1410507/150‑1 (MQ=255)
                                                                                                              aGTTGCATATTGCGTTCTCCGTATATCTGGCGGGGATGACAGCTGCGATGTTATTTGCCGGTAAAGTGGCCGATCGTTCAGGGAGAAAGCCGGTCGCCATACCCGGCGCGGCGCTATTTATTATTGCCTCGGTGTTCTGTTCACTGGCTGa                                     <  1:232011/151‑1 (MQ=255)
                                                                                                                    atatTGCGTTCTCCGTATATCTGGCGGGGATGACAGCTGCGATGTTATTTGCCGGTAAAGTGGCCGATCGTTCAGGGAGAAAGCCGGTCGCCATACCCGGCGCGGCGCTATTTATTATTGCCTCGGTGTTCTGTTCACTGGCTGAAACCAg                               <  3:372083/151‑1 (MQ=255)
                                                                                                                      atTGCGTTCTCCGTATATCTGGCGGGGATGACAGCTGCTATGTTATTTGCCGGTAAAGTGGCCGATCGTTCAGGGAGAAAGCCGGTCGCCATACCCGGCGCGGCGCTATTTATTATTGCCTCGGTGTTCTGTTCACTGGCTGAAACCAGCa                             <  1:1440854/151‑1 (MQ=255)
                                                                                                                        tGCGTTCTCCGTATATCTGGCGGGGATGACAGCTGCGATGTTATTTGCCGGTAAAGTGGCCGATCGTTCAGGGAGAAAGCCGGTCGCCATACCCGGCGCGGCGCTATTTATTATTGCCTCGGTGTTCTGTTCACTGGCTGAAACCAGCACg                           >  2:305350/1‑151 (MQ=255)
                                                                                                                         gCGTTCTCCGTATATCTGGCGGGGATGACAGCTGCGATGTTATTTGCCGGTAAAGTGGCCGATCGTTCAGGGAGAAAGCCGGTCGCCATACCCGGCGCGGCGCTATTTATTATTGCCTCGGTGTTCTGTTCACTGGCTGAAACCAGCACGt                          >  1:936013/1‑151 (MQ=255)
                                                                                                                           gTTCTCCGTATATCTGGCGGGGATGACAGCTGCGATGTTATTTGCCGGTAAAGTGGCCGATCGTTCAGGGAGAAAGCCGGTCGCCATACCCGGCGCGGCGCTTTTTATTATTGCCTCGGAGTTCTGTTCACTGGCTGAAACCAGCACGtt                         >  1:1346691/1‑150 (MQ=255)
                                                                                                                                ccGTATATCTGGCGGGGATGACAGCTGCGATGTTATTTGCCGGTAAAGTGGCCGATCGTTCAGGGAGAAAGCCGGTCGCCATACCCGGCGCGGCGCTATTTATTATTGCCTCGGTGTTCTGTTCACTGGCTGAAACCAGCACGTTAtt                      <  3:19188/148‑1 (MQ=255)
                                                                                                                                  gTATAGCTGGCGGTGATGACAGCTGCGATGTTATTTGCCGGTAAAGTGGCCGATCGTTCAGGGAGAAATCCGGTCGCCATACCCGGCGCTGCGCTATTTATTATTGCCTCGGTGTGCTGTCCACTGGCTGAAACCAGCACGTTATTTCtag                 <  3:99739/151‑3 (MQ=255)
                                                                                                                                   tataTCTGGCGGGGATGACAGCTGCGATGTTATTTGCCGGTAAAGTGGCCGATCGTTCAGGGAGAAAGCCGGTCGCCATACCCGGCGCGGCGCTATTTATTATTGCCTCGGTGTTCTGTTCACTGGCTGAAACCAGCACGTTATTTCTTGc                >  3:513731/1‑151 (MQ=255)
                                                                                                                                      atCTGGCGGGGATGACAGCTGCGATGTTATTTGCCGGTAAAGTGGCCGATCGTTCAGGGAGAAAGCCGGTCGCCATACCCGGCGCGGCGCTATTTATTATTGCCTCGGTGTTCTGTTCACTGGCTGAAACCAGCACGTTATTTCTTGCAg              <  2:309725/150‑1 (MQ=255)
                                                                                                                                        cTGGCGGGGATGACAGCTGCGATGTTATTTGCCGGTAAAGTGGCCGATCGTTCAGGGAGAAAGCCGGTCGCCATACCCGGCGCGGCGCTATTTATTATTGCCTCGGTGTTCTGTTCACTGGCTGAAACCAGCACGTTATTTCTTGCAGGcc           >  4:1664333/1‑151 (MQ=255)
                                                                                                                                             ggggATGACAGCTGCGATGTTATTTGCCGGTAAAGTGGCCGATCGTTCAGGGAGAAAGCCGGTCGCCATACCCGGCGCGGCGCTATTTATTATTGCCTCGGTGTTCTGTTCACTGGCTGAAACCAGCACGTTATTTCTTGCAGGCCGAttt      >  4:1699882/1‑151 (MQ=255)
                                                                                                                                               ggATGACAGCTGCGATGTTATTTGCCGGTAAAGTGGCCGATCGTTCAGGGAGAAAGCCGGTCGCCATACCCGGCGCGGCGCTATTTATTATTGCCTCGGTGTTCTGTTCACTGGCTGAAACCAGCACGTTATTTCTTGCAGGCCGATTTCt    <  2:422399/151‑1 (MQ=255)
                                                                                                                                                gATGACAGCTGCGATGTTATTTGCCGGTAAAGTGGCCGATCGTTCAGGGAGAAAGCCGGTCGCCATACCCGGCGCGGCGCTATTTATTATTGCCTCGGTGTTCTGTTCACTGGCTGAAACCAGCACGTTATTTCTTGCAGGCCGATTTCTa   <  4:601554/151‑1 (MQ=255)
                                                                                                                                                 aTGACAGCTGCGATGTTATTTGCCGGTAAAGTGGCCGATCGTTCAGGGAGAAAGCCGGTCGCCATACCCGGCGCGGCGCTATTTATTATTGCCTCGGTGTTCTGTTCACTGGCTGAAACCAGCACGTTATTTCTTGCAGGCCGATTTCTAc  >  3:521322/1‑151 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CCGCTTTTTGATTTGTAGTTTTGCCCTGGTTTTACTTTATCCCGCCGGGATTGATATGTACCTCGTTGGTTTACCGCGCATCGCCGCCGATCTCAATGCCAGCGAAGCGCAGTTGCATATTGCGTTCTCCGTATATCTGGCGGGGATGGCAGCTGCGATGTTATTTGCCGGTAAAGTGGCCGATCGTTCAGGGAGAAAGCCGGTCGCCATACCCGGCGCGGCGCTATTTATTATTGCCTCGGTGTTCTGTTCACTGGCTGAAACCAGCACGTTATTTCTTGCAGGCCGATTTCTAC  >  NC_000913/3891620‑3891915

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