Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,589,122 G→A 15.3% G510S (GGC→AGC)  acrD → aminoglycoside/multidrug efflux system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,589,1220GA15.3% 312.0 / 37.9 130G510S (GGC→AGC) acrDaminoglycoside/multidrug efflux system
Reads supporting (aligned to +/- strand):  ref base G (68/41);  new base A (14/6);  total (82/48)
Fisher's exact test for biased strand distribution p-value = 6.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.50e-01

Only 100 of 141 total aligned reads displayed.

GCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGA  >  NC_000913/2589023‑2589221
                                                                                                   |                                                                                                   
gcgATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAg                                                                                                     <  1:3915161/100‑1 (MQ=255)
 cgATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGGGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAgg                                                                                                    <  2:3222267/100‑1 (MQ=255)
 cgATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGGGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAg                                                                                                     <  1:3915089/99‑1 (MQ=255)
    tGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCtt                                                                                                 >  1:2158490/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)
         cTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTATGACACAATACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCtttttt                                                                                             >  2:2019552/1‑99 (MQ=255)
          tGTCAGTACTGGTAGCGATGATCCTCTCTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTg                                                                                            <  1:2490792/99‑1 (MQ=255)
          tGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGc                                                                                           <  2:94471/100‑1 (MQ=255)
          tGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGc                                                                                           <  1:3012951/100‑1 (MQ=255)
           gTCAGTACTGGTAGCGATGTTCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAAGCCTTTTTGcc                                                                                          <  1:2321340/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                                                                                        <  2:945844/100‑3 (MQ=255)
             cAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATAGGCAAAAAAGCCTTTTTGccca                                                                                        <  1:1864509/100‑3 (MQ=255)
              aGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTgg                                                                                       >  1:1778747/1‑100 (MQ=255)
              aGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTg                                                                                        <  2:383448/99‑1 (MQ=255)
                tACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGtt                                                                                     <  1:3455386/100‑1 (MQ=255)
                tACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGt                                                                                      >  2:3714150/1‑99 (MQ=255)
                 aCTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGttt                                                                                    <  2:3558350/100‑1 (MQ=255)
                  cTGGTCGCGATGATCCTAACTCCGGCTCTGTGTACCAGACTCCTTAAGCCACTGAAAAAAGGTGAGCCTCATGGGCAAAAAGGCTTTTTTGCCTGGTTTa                                                                                   >  2:2087767/1‑100 (MQ=255)
                    ggTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAAGTGAACATCATAAGCAAAAAGGCTTTTTTGCCTGGTTTAAc                                                                                 >  1:187733/1‑100 (MQ=255)
                    ggTAGCGATGATCCGCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTGTTTTGCCTGGTTTAAc                                                                                 >  2:2279277/1‑100 (MQ=255)
                      tAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCa                                                                               >  2:530657/1‑100 (MQ=255)
                         cGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGAt                                                                            >  2:3854550/1‑100 (MQ=255)
                         cGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGAt                                                                            >  1:2918720/1‑100 (MQ=255)
                         cGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCAACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGAt                                                                            >  2:216446/1‑100 (MQ=255)
                          gatgatCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGAt                                                                            >  1:1851557/1‑99 (MQ=255)
                          gatgatACTTACTCCGGATCTGTGTGCCACACTAATTAAGCCACTGAAAAAAGGTGAGCATCAGAGGCAAAAAGGCTTGGTTGCCTAGGTTAACCACAt                                                                            >  2:1688010/1‑99 (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)
                                   cACTCCGGCTCTGTGTGCCACACTACGTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGACTTTTTTGCCTGGTTTAACCAGATGTTTAACCGc                                                                  >  2:3702898/1‑100 (MQ=255)
                                     cTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCaa                                                                >  1:2449603/1‑100 (MQ=255)
                                       ccGGCTCTGTGTGCCACACTACTTAAGCCACAGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACg                                                              >  2:413278/1‑100 (MQ=255)
                                         ggCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGcc                                                            <  1:640548/100‑1 (MQ=255)
                                          gCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCg                                                           >  2:2662360/1‑100 (MQ=255)
                                          gCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCg                                                           <  2:3511594/100‑1 (MQ=255)
                                             ctGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACACCGaa                                                         >  1:653342/1‑99 (MQ=255)
                                                tgtgCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCt                                                     >  2:559876/1‑100 (MQ=255)
                                                 gtgCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTa                                                    <  2:50937/100‑1 (MQ=255)
                                                  tggcacaCTACTTAAGGCACTGAAAAAATGTTATCATCATTTGCAAAAAGTCTGGTTTGCCTGGTTTAACCATATGTTTAACCGCAACGCCTAACGCTAc                                                   <  1:2739478/97‑1 (MQ=255)
                                                   gCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACg                                                  >  1:3266428/1‑100 (MQ=255)
                                                     cacaCTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAAACAGATGTTTAACCGCAACGCAGAACGATACGaa                                                >  1:2563132/1‑100 (MQ=255)
                                                       cacTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGaaaa                                              >  2:2394604/1‑100 (MQ=255)
                                                            ctTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAgggg                                         >  2:3254097/1‑100 (MQ=255)
                                                            ctTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAgggg                                         >  2:251151/1‑100 (MQ=255)
                                                            ctTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAgggg                                         >  2:876175/1‑100 (MQ=255)
                                                            ctTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAgggg                                         <  2:2052475/100‑1 (MQ=255)
                                                            ctTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAggg                                          >  1:2591766/1‑99 (MQ=255)
                                                             ttAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGt                                        <  2:3888725/100‑1 (MQ=255)
                                                             ttAAGCCACTGAAAAAAGGTGAGAATCATGGACAACAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCATCTCCGAACGATACGAAAAAGGGGt                                        <  2:389357/100‑1 (MQ=255)
                                                              tAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTg                                       >  2:10270/1‑100 (MQ=255)
                                                              tAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTg                                       >  1:2596769/1‑100 (MQ=255)
                                                               aaGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTgg                                      <  2:2527315/100‑1 (MQ=255)
                                                                 gCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCg                                    >  1:1029073/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)
                                                                          aaaaGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATtctc                           <  2:1198797/100‑1 (MQ=255)
                                                                           aaaGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTcc                          <  2:784680/100‑1 (MQ=255)
                                                                           aaaGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCAAAAAGATAAAAAAAAGGGGTGGCGAAAATTCTcc                          <  2:1481787/100‑1 (MQ=255)
                                                                           aaaGGTGAGCATCACGAGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGTAAAAAAAAAAAAAAAAGGGGTGGCGAAAATTCTcc                          <  2:1480667/100‑1 (MQ=255)
                                                                            caGGTGAGCATCATGGACAAAAAGGCTTTTTTGCCTGGTTTAACCAGCGGTTTAACCGCAACGCCGAACGCTACGAAAAATGTGTGGCGACAATtct                            >  2:3660698/2‑97 (MQ=255)
                                                                                tGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGt                     >  2:584415/1‑100 (MQ=255)
                                                                                 gAGCATCATGGGCAAAAAGGTTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAAAAAAAAAAGAAAAGAAAAAGGGGTGGCGAAAATTCTCCACCGt                     <  1:817758/99‑1 (MQ=255)
                                                                                    taaaccaggCAAAAAAGCCTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGcc                 >  1:1599683/8‑100 (MQ=38)
                                                                                     atcatGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCt                >  1:2960774/1‑100 (MQ=255)
                                                                                      tcatGGGCAAAAAGGCTTTTTTGCCTGGTTTAATCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTg               <  1:963356/100‑1 (MQ=255)
                                                                                      aaccaggCAAAAAAGCCTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTg               >  2:2886158/6‑100 (MQ=255)
                                                                                       accaggCAAAAAAGCCTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTgc              >  2:833999/5‑100 (MQ=255)
                                                                                       accaggCAAAAAAGCCTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTgc              >  2:3233611/5‑100 (MQ=255)
                                                                                        ccaggCAAAAAAGCCTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGAGGTGGCGAAAATTCTCCACCGTAGCCTgcg             >  1:2861821/4‑100 (MQ=255)
                                                                                         caggCAAAAAAGCCTTTTTGTTTTTTTTTAACAGATGTTTAACCGCAACGTCGAACGCTATGAAAAAGGTGTTGCGAAAATTCTTCACCATAGCCt                >  2:2273319/3‑96 (MQ=255)
                                                                                         caggCAAAAAAGCCTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTgcgc            >  1:308805/3‑100 (MQ=255)
                                                                                         caggCAAAAAAGCCTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGCAGAGAAAATGCTCAACCGTAGCCTgcgc            >  2:2484411/3‑100 (MQ=255)
                                                                                         caggCAAAAAAGCCTTTTTGCCTGCTTTAACCAGATGTTTAACCGCAACGCCGAACGCTCCGAAAAAGGGCGGGCGAAAATTCTCCACCGTAGCCTgc              >  2:1506267/3‑98 (MQ=255)
                                                                                          gggCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCt           >  1:3562125/1‑100 (MQ=255)
                                                                                          gggCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCt           <  1:3203512/100‑1 (MQ=255)
                                                                                          gggCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCt           <  2:3443979/100‑1 (MQ=255)
                                                                                          aggCAAAAAAGCCTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCt           >  2:2198489/2‑100 (MQ=255)
                                                                                          aggCAAAAAAGCCTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCt           >  1:1178419/2‑100 (MQ=255)
                                                                                          aggCAAAAAAGCCTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCt           >  1:2291042/2‑100 (MQ=255)
                                                                                           ggCAAAAAGGCTTTTTTGCCTTGTTTAACCAAATGTTTAACCGCAATGCCGTATGCTACGAAAAAGTGGTGGCGAAAATTCTCTACTGTAGTCttcact           >  2:634226/1‑94 (MQ=255)
                                                                                           ggCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTATCCACCGTAGCCTGCGCTg          >  1:2311236/1‑100 (MQ=255)
                                                                                               aaaaGGCTTTTTTGCCTGGTTTAACCGGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGAt       <  2:3048210/99‑1 (MQ=255)
                                                                                                  aGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATtgtg   <  1:10270/100‑1 (MQ=255)
                                                                                                   ggCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGa  <  1:965300/100‑1 (MQ=255)
                                                                                                   ggCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCAACCGTAGCCTGCGCTGGATTGTGa  >  1:1910890/1‑100 (MQ=255)
                                                                                                   ggCTTTTTTGCCTGGTTTAAACAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCCCCGTAGCCTGCGCTGGATtgtg   <  2:2534863/99‑1 (MQ=255)
                                                                                                   |                                                                                                   
GCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGA  >  NC_000913/2589023‑2589221

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: