Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,589,125 T→C 16.6% F511L (TTT→CTT)  acrD → aminoglycoside/multidrug efflux system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,589,1250TC16.6% 270.9 / 47.3 126F511L (TTT→CTT) acrDaminoglycoside/multidrug efflux system
Reads supporting (aligned to +/- strand):  ref base T (63/42);  new base C (14/7);  total (77/49)
Fisher's exact test for biased strand distribution p-value = 6.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.15e-01

Only 100 of 136 total aligned reads displayed.

TGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGATTT  >  NC_000913/2589027‑2589224
                                                                                                  |                                                                                                   
 ggTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCttt                                                                                                   >  2:3675185/1‑100 (MQ=255)
 ggTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTACCACTCTACTTAAGCCGTTGAAAAAAGATGAGCATCATGGGCAAAAAGGCttt                                                                                                   >  2:3053135/1‑100 (MQ=255)
    gCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCtttttt                                                                                                <  1:743377/100‑1 (MQ=255)
    gCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGATGAGCATCATGGGCAAAAAGGCttttt                                                                                                 >  1:3443609/1‑99 (MQ=255)
    gCTGTCAGTACCGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCtttttt                                                                                                >  2:5698/1‑100 (MQ=255)
     cTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTATGACACAATACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCtttttt                                                                                                >  2:2019552/1‑99 (MQ=255)
      tGTCAGTACTGGTAGCGATGATCCTCTCTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTg                                                                                               <  1:2490792/99‑1 (MQ=255)
      tGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGc                                                                                              <  1:3012951/100‑1 (MQ=255)
      tGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGc                                                                                              <  2:94471/100‑1 (MQ=255)
       gTCAGTACTGGTAGCGATGTTCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAAGCCTTTTTGcc                                                                                             <  1:2321340/100‑1 (MQ=255)
       gTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAAGCCTTTTTGcc                                                                                             <  1:1375850/100‑1 (MQ=255)
       gTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAAGCCTTTTTGcc                                                                                             <  2:262511/100‑1 (MQ=255)
        tCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCt                                                                                            >  1:1853305/1‑100 (MQ=255)
        tCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAAGCCTTTTTGccc                                                                                            <  2:3540289/100‑2 (MQ=255)
         cAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTg                                                                                           >  1:2464231/1‑100 (MQ=255)
         cAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAAGCCTTTTTGccca                                                                                           >  1:945855/1‑98 (MQ=255)
         cAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAAGCCTTTTTGccca                                                                                           <  2:945844/100‑3 (MQ=255)
          aGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTgg                                                                                          >  1:1778747/1‑100 (MQ=255)
            tACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGt                                                                                         >  2:3714150/1‑99 (MQ=255)
             aCTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGttt                                                                                       <  2:3558350/100‑1 (MQ=255)
                ggTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAAGTGAACATCATAAGCAAAAAGGCTTTTTTGCCTGGTTTAAc                                                                                    >  1:187733/1‑100 (MQ=255)
                ggTAGCGATGATCCGCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTGTTTTGCCTGGTTTAAc                                                                                    >  2:2279277/1‑100 (MQ=255)
                 gTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAAcc                                                                                   >  2:3883240/1‑100 (MQ=255)
                 gTAGCGATGATCCTCACTCCGGCTATGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAAcc                                                                                   <  1:508001/100‑1 (MQ=255)
                  tAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCa                                                                                  >  2:530657/1‑100 (MQ=255)
                  tAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAACCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCa                                                                                  >  2:539212/1‑100 (MQ=255)
                     cGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGAt                                                                               >  2:1516105/1‑100 (MQ=255)
                     cGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGAt                                                                               >  2:3854550/1‑100 (MQ=255)
                     cGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCAACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGAt                                                                               >  2:216446/1‑100 (MQ=255)
                      gatgatCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGAt                                                                               >  1:1851557/1‑99 (MQ=255)
                        tgatTCTCACTCCGTCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGtt                                                                            <  1:2534079/100‑1 (MQ=255)
                        tgatCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGtt                                                                            >  2:1305744/1‑100 (MQ=255)
                         gatCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGttt                                                                           >  1:2529261/1‑100 (MQ=255)
                          atCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTa                                                                          >  2:3650160/1‑100 (MQ=255)
                          atCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGACAAAAAGGCTTTTTTGCCTGGTTGAAACAGATGTTTa                                                                          >  1:2275024/1‑100 (MQ=255)
                             cTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAAcc                                                                       >  1:1366787/1‑100 (MQ=255)
                              tCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCg                                                                      <  2:3821727/100‑1 (MQ=255)
                               cACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGACAAAAAGGCTTTTTTGCCTGGTTTAACAAGATGTTTAACCGc                                                                     >  1:100731/1‑100 (MQ=255)
                                 cTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCaa                                                                   >  1:2449603/1‑100 (MQ=255)
                                   ccGGCTCTGTGTGCCACACTACTTAAGCCACAGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACg                                                                 >  2:413278/1‑100 (MQ=255)
                                      gCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCg                                                              <  2:3511594/100‑1 (MQ=255)
                                         ctGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACACCGaa                                                            >  1:653342/1‑99 (MQ=255)
                                              tggcacaCTACTTAAGGCACTGAAAAAATGTTATCATCATTTGCAAAAAGTCTGGTTTGCCTGGTTTAACCATATGTTTAACCGCAACGCCTAACGCTAc                                                      <  1:2739478/97‑1 (MQ=255)
                                               gCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACg                                                     >  1:3266428/1‑100 (MQ=255)
                                                ccACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGa                                                    >  2:3817554/1‑100 (MQ=255)
                                                  acacTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGaaa                                                  <  1:3018592/100‑1 (MQ=255)
                                                     ctactTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTTTTTTTTTTTCTCAGATGTATAACCGCAACGCCGAACGCTACGaaaa                                                 >  2:212085/1‑98 (MQ=255)
                                                        ctTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGCAAAAgggg                                            >  2:1036239/1‑100 (MQ=255)
                                                        ctTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAgggg                                            >  2:251151/1‑100 (MQ=255)
                                                        ctTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAgggg                                            <  2:2052475/100‑1 (MQ=255)
                                                        ctTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAgggg                                            >  2:876175/1‑100 (MQ=255)
                                                        ctTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAgggg                                            >  2:3254097/1‑100 (MQ=255)
                                                        ctTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAggg                                             >  1:2591766/1‑99 (MQ=255)
                                                         ttAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGt                                           <  2:3888725/100‑1 (MQ=255)
                                                          tAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTg                                          >  1:2596769/1‑100 (MQ=255)
                                                          tAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTg                                          >  2:10270/1‑100 (MQ=255)
                                                           aaGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTgg                                         <  2:2527315/100‑1 (MQ=255)
                                                             gCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCg                                       >  1:1029073/1‑100 (MQ=255)
                                                             gCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCg                                       >  1:227735/1‑100 (MQ=255)
                                                             gCCACTGAAAAAAGGTGAGCATAATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAAGGATACGAAAAAGGGGTGGCg                                       >  2:2916280/1‑100 (MQ=255)
                                                               cACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGaa                                     >  1:2894000/1‑100 (MQ=255)
                                                                 cTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGaaaa                                   <  2:2811704/100‑1 (MQ=255)
                                                                 cTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGaaaa                                   >  1:494624/1‑100 (MQ=255)
                                                                  tGAAAAAAGGTGATTATCATGGTGTAAATGGTTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGaaaa                                   <  1:3543988/99‑1 (MQ=255)
                                                                  tGAAAAAAGGTGAGCATAATAAGAAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAAt                                  <  2:1103533/100‑1 (MQ=255)
                                                                   gAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAGGGGGTGGCGAAAAtt                                 >  2:543906/1‑100 (MQ=255)
                                                                   gAAAAAAGGTGAGCATCATGGGAAAAAAGGCTTTTTTGCATAGTTTAACCAAATGTTTAACAGAAACGACAAACACTAAGAAAAAAAGGTGGCGaaaa                                   >  2:3599381/1‑98 (MQ=255)
                                                                    aaaaaaGGTGAGCATCATCGGGAAAAACGCCTTTTTGCCTGGTTTAACCATATGTTTATCCGCATTGAAAAAAAAAAAAAAAAAGGGGTGGCGAAAATtc                                <  2:165712/100‑1 (MQ=255)
                                                                     aaaaaGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATtct                               <  1:1829258/100‑1 (MQ=255)
                                                                       aaaGGTGAGCATCACGAGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGTAAAAAAAAAAAAAAAAGGGGTGGCGAAAATTCTcc                             <  2:1480667/100‑1 (MQ=255)
                                                                        caGGTGAGCATCATGGACAAAAAGGCTTTTTTGCCTGGTTTAACCAGCGGTTTAACCGCAACGCCGAACGCTACGAAAAATGTGTGGCGACAATtct                               >  2:3660698/2‑97 (MQ=255)
                                                                          ggTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCAcc                          >  2:191724/1‑100 (MQ=255)
                                                                            tGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGt                        >  2:584415/1‑100 (MQ=255)
                                                                               gCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGc                     >  1:352144/1‑100 (MQ=255)
                                                                                 atcatGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCt                   >  1:2960774/1‑100 (MQ=255)
                                                                                  tcatGGGCAAAAAGGCTTTTTTGCCTGGTTTAATCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTg                  <  1:963356/100‑1 (MQ=255)
                                                                                  aaccaggCAAAAAAGCCTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTg                  >  2:2886158/6‑100 (MQ=255)
                                                                                   accaggCAAAAAAGCCTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTgc                 >  2:3233611/5‑100 (MQ=255)
                                                                                   accaggCAAAAAAGCCTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTgc                 >  2:833999/5‑100 (MQ=255)
                                                                                    ccaggCAAAAAAGCCTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGAGGTGGCGAAAATTCTCCACCGTAGCCTgcg                >  1:2861821/4‑100 (MQ=255)
                                                                                     tGGGCAAAAAGGCTTTTTTGACTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTgcgc               <  1:679112/100‑1 (MQ=255)
                                                                                     caggCAAAAAAGCCTTTTTGTTTTTTTTTAACAGATGTTTAACCGCAACGTCGAACGCTATGAAAAAGGTGTTGCGAAAATTCTTCACCATAGCCt                   >  2:2273319/3‑96 (MQ=255)
                                                                                     caggCAAAAAAGCCTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGCAGAGAAAATGCTCAACCGTAGCCTgcgc               >  2:2484411/3‑100 (MQ=255)
                                                                                     caggCAAAAAAGCCTTTTTGCCTGCTTTAACCAGATGTTTAACCGCAACGCCGAACGCTCCGAAAAAGGGCGGGCGAAAATTCTCCACCGTAGCCTgc                 >  2:1506267/3‑98 (MQ=255)
                                                                                      gggCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCt              <  2:3443979/100‑1 (MQ=255)
                                                                                      gggCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCt              >  1:3562125/1‑100 (MQ=255)
                                                                                      gggCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCt              <  1:3203512/100‑1 (MQ=255)
                                                                                      aggCAAAAAAGCCTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCt              >  2:2198489/2‑100 (MQ=255)
                                                                                      aggCAAAAAAGCCTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCt              >  1:2291042/2‑100 (MQ=255)
                                                                                      aggCAAAAAAGCCTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCt              >  1:1178419/2‑100 (MQ=255)
                                                                                       ggCAAAAAGGCTTTTTTGCCTTGTTTAACCAAATGTTTAACCGCAATGCCGTATGCTACGAAAAAGTGGTGGCGAAAATTCTCTACTGTAGTCttcact              >  2:634226/1‑94 (MQ=255)
                                                                                         cAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGa           <  2:3927633/100‑1 (MQ=255)
                                                                                              aGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATtgtg      <  2:2311214/100‑1 (MQ=255)
                                                                                              aGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATtgtg      <  1:10270/100‑1 (MQ=255)
                                                                                              aGGCTTTTTTGCCTGGTTTAACAAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATtgtg      >  2:824357/1‑100 (MQ=255)
                                                                                               ggCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCAACCGTAGCCTGCGCTGGATTGTGa     >  1:1910890/1‑100 (MQ=255)
                                                                                               ggCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGATGGCGAAAATTCTCAACCGTAACCTGCGCTGGATTGTGa     >  1:1861741/1‑100 (MQ=255)
                                                                                               ggCTTTTTTGCCTGGTTTAAACAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCCCCGTAGCCTGCGCTGGATtgtg      <  2:2534863/99‑1 (MQ=255)
                                                                                                  ttttttGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGAttt  <  1:836908/100‑1 (MQ=255)
                                                                                                  ttttttGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGAttn  <  1:1733048/100‑2 (MQ=255)
                                                                                                  |                                                                                                   
TGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGATTT  >  NC_000913/2589027‑2589224

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