Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,589,125 T→C 15.1% 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,1250TC15.1% 257.2 / 37.4 113F511L (TTT→CTT) acrDaminoglycoside/multidrug efflux system
Reads supporting (aligned to +/- strand):  ref base T (54/41);  new base C (9/8);  total (63/50)
Fisher's exact test for biased strand distribution p-value = 7.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.35e-01

Only 100 of 120 total aligned reads displayed.

TGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGATT  >  NC_000913/2589027‑2589223
                                                                                                  |                                                                                                  
tGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCGACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCtt                                                                                                   >  1:1626638/1‑100 (MQ=255)
tGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCtt                                                                                                   <  2:538509/100‑1 (MQ=255)
 ggTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAGGTGAGCATCATGGGCAAAAAGGCtttt                                                                                                 >  2:555855/1‑100 (MQ=255)
    gCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCtttttt                                                                                               >  2:2257557/1‑100 (MQ=255)
     cTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTg                                                                                              >  1:1997352/1‑100 (MQ=255)
      tGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAGGTGAGCATCATGGGCAAAAAAGCCTTTTTGcc                                                                                            <  2:394025/100‑1 (MQ=255)
       gTCAGTACTGGTAGCGATGATCCTCACTCCGGTTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAAGCCTTTTTGcc                                                                                            <  1:965861/100‑1 (MQ=255)
       gTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGcc                                                                                            >  1:941136/1‑100 (MQ=255)
       gTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAAGCCTTTTTGcc                                                                                            <  1:738361/100‑1 (MQ=255)
       gTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAAGCCTTTTTGcc                                                                                            <  2:592141/100‑1 (MQ=255)
       gTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAAAAAAAAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCt                                                                                           <  1:119621/100‑1 (MQ=255)
        tCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCt                                                                                           <  2:2167364/100‑1 (MQ=255)
        tCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCt                                                                                           >  1:995216/1‑100 (MQ=255)
         cAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTg                                                                                          >  1:1872769/1‑100 (MQ=255)
         cAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTg                                                                                          >  2:2382875/1‑100 (MQ=255)
         cAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCACTGAAAAAAGGTGAGCATCATGGGCAAAAAAGCCTTTTTGcccat                                                                                         <  1:1870254/100‑4 (MQ=255)
         cAGTACTGGTAGAGATGATACTCACTACGGATCTGTGTGCCACACTACAAAAAAAAAAAAAAAAAGGTGAGCATCATGGGCAAAAAAGCCTTTTTGcc                                                                                            <  1:1545741/98‑1 (MQ=255)
          aGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGACAAAAGGTGAGCATCATCGGCAACAAGGCTTTTTTGCCTgg                                                                                         >  1:686506/1‑100 (MQ=255)
          aGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTgg                                                                                         <  2:1172311/100‑1 (MQ=255)
            tACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGtt                                                                                       <  2:617737/100‑1 (MQ=255)
             aCTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGtt                                                                                       <  1:233121/99‑1 (MQ=255)
              cTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTa                                                                                     >  1:2592976/1‑100 (MQ=255)
                ggTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAAc                                                                                   >  2:1928922/1‑100 (MQ=255)
                ggTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGAGCAAAAAGGCTTTTTTGCCTGGTTTAAc                                                                                   >  2:2006207/1‑100 (MQ=255)
                 gTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAAcc                                                                                  <  2:2436660/100‑1 (MQ=255)
                  tAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCa                                                                                 <  2:1266754/100‑1 (MQ=255)
                  tAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCa                                                                                 <  1:595608/100‑1 (MQ=255)
                   aGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGa                                                                               <  2:1167828/100‑1 (MQ=255)
                      gatgatCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATg                                                                             >  1:850031/1‑100 (MQ=255)
                       atgatCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGtt                                                                           <  1:2266895/100‑1 (MQ=255)
                         gatCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGttt                                                                          >  2:178289/1‑100 (MQ=255)
                         gatCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGttt                                                                          <  2:1996260/100‑1 (MQ=255)
                         gatCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGttt                                                                          <  2:293520/100‑1 (MQ=255)
                          atCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTa                                                                         >  2:718151/1‑100 (MQ=255)
                            ccTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAAc                                                                       <  1:178288/100‑1 (MQ=255)
                             cTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAAcc                                                                      <  2:1602626/100‑1 (MQ=255)
                              tCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCg                                                                     <  2:2126809/100‑1 (MQ=255)
                              tCACTCCGGCTCTGTGTGCCACACTACTTAAGACAAAGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAAcc                                                                      <  1:577975/99‑1 (MQ=255)
                                 cTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTTTTTTGTTTAACCAGATGTTTAACCGCaa                                                                  >  1:1350537/1‑100 (MQ=255)
                                 cTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAAc                                                                 >  1:1313708/1‑100 (MQ=255)
                                    cGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGc                                                               >  2:2155232/1‑100 (MQ=255)
                                     ggCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGcc                                                              <  2:2318688/100‑1 (MQ=255)
                                      gCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCg                                                             <  2:1421776/100‑1 (MQ=255)
                                         ctGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAAc                                                          <  2:1608643/100‑1 (MQ=255)
                                          tgtgtgCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACg                                                         <  2:2448232/100‑1 (MQ=255)
                                           gtgtgCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGc                                                        >  2:654417/1‑100 (MQ=255)
                                           gtgtgCCACACTACTTAAGCCACTGAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCt                                                       >  1:2583270/1‑100 (MQ=255)
                                              tgCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTAc                                                     <  2:1254336/100‑1 (MQ=255)
                                               gCCACACTACTTAAGCCACTGAACAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCATGGATTAACCAGATGTTTAACCGCAACGCATAACGCTAc                                                     >  1:985862/1‑99 (MQ=255)
                                                ccACACTACTTAAGCCACTGAATAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGa                                                   <  1:2323701/100‑1 (MQ=255)
                                                 cacaCTACTTAAGCCACTGAGAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTTCTTTTTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGaa                                                  >  1:918835/1‑100 (MQ=255)
                                                      tactTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAgg                                             <  2:2441771/100‑1 (MQ=255)
                                                       actTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTTTTTTTTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAggg                                            >  1:1612611/1‑100 (MQ=255)
                                                       actTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAggg                                            >  1:1035861/1‑100 (MQ=255)
                                                          tAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTg                                         >  1:2530776/1‑100 (MQ=255)
                                                            aGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGc                                       <  2:1997363/100‑1 (MQ=255)
                                                             gCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCg                                      <  2:138158/100‑1 (MQ=255)
                                                                aCTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCgca                                    >  2:548274/1‑97 (MQ=255)
                                                                 cTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTTTTTTTTTTTTTCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGaaa                                   >  1:880557/1‑99 (MQ=255)
                                                                 cTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGaaaa                                  >  2:2239785/1‑100 (MQ=255)
                                                                  tGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCTGAACGCTACGAAAAAGGGGTGGCGAAAAt                                 >  1:1184068/1‑100 (MQ=255)
                                                                  tGAAAAAAGGTGAGCATCATGAGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAAt                                 <  1:2006196/100‑1 (MQ=255)
                                                                   gAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAAtt                                >  1:1722521/1‑100 (MQ=255)
                                                                   gAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAAtt                                >  1:358881/1‑100 (MQ=255)
                                                                    gaaaaaGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATtc                               >  2:119133/2‑100 (MQ=255)
                                                                    aaaaaaGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATtc                               >  2:2216906/1‑100 (MQ=255)
                                                                     aaaaaGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATtct                              <  2:1052181/100‑1 (MQ=255)
                                                                       aaaGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGAGaattttcgcc                            <  2:377321/100‑9 (MQ=255)
                                                                            tGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCAccct                       >  1:2401258/1‑98 (MQ=255)
                                                                            tGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGt                       >  1:530655/1‑100 (MQ=255)
                                                                               gCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGc                    >  2:2474651/1‑100 (MQ=255)
                                                                               gCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGc                    >  2:2473730/1‑100 (MQ=255)
                                                                                catcatGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGcc                   <  2:586016/100‑1 (MQ=255)
                                                                                 atcatGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCt                  >  2:1686189/1‑100 (MQ=255)
                                                                                 atcatGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCt                  >  1:2241371/1‑100 (MQ=255)
                                                                                 atcatGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCt                  >  1:1958116/1‑100 (MQ=255)
                                                                                 atcatGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCt                  >  2:1267031/1‑100 (MQ=255)
                                                                                  tcatGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTg                 >  1:1651653/1‑100 (MQ=255)
                                                                                  aaccaggCAAAAAAGCCTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTg                 >  1:2571271/6‑100 (MQ=255)
                                                                                   accaggCAAAAAAGCCTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTgc                >  1:2150246/5‑100 (MQ=255)
                                                                                    ccaggCAAAAAAGCCTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTgcg               >  2:1186571/4‑100 (MQ=255)
                                                                                    ccaggCAAAAAAGCCTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTgcg               >  2:412852/4‑100 (MQ=255)
                                                                                    ccaggCAAAAAAGCCTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTgcgc              >  1:430232/4‑100 (MQ=255)
                                                                                    atGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTgcg               >  2:2587565/1‑100 (MQ=255)
                                                                                     tGGGCAAAAAGGCTTTTTTTTTTTTTTTTTCCAGATGTTTAACCGCAACGCTTAACGTTATGACAATGGGGGGGGGAAACTTCTTCATCGTAGCCTgc                >  1:2433915/1‑98 (MQ=255)
                                                                                     tGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTgcgc              >  1:152103/1‑100 (MQ=255)
                                                                                     caggCAAAAAAGCCTTTTTTTCTGGTTTAATTAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTgcgc              >  1:1020207/3‑100 (MQ=255)
                                                                                     caggCAAAAAAGCCTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTgcgc              >  2:421391/3‑100 (MQ=255)
                                                                                     caggCAAAAAAGCCTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTgcgc              >  2:792133/3‑100 (MQ=255)
                                                                                       ggCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTg            >  2:1596625/1‑100 (MQ=255)
                                                                                        gCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTgg           >  1:85822/1‑100 (MQ=255)
                                                                                         cAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGa          >  1:2562730/1‑100 (MQ=255)
                                                                                         cAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGa          >  1:287168/1‑100 (MQ=255)
                                                                                         cAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGa          >  2:779884/1‑100 (MQ=255)
                                                                                            aaaGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATtg       >  1:1068280/1‑100 (MQ=255)
                                                                                            aaaGGCTTTTGTGCATGGTTTCGCCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAATGGGTGGCGAAAATTCGCCACCGTAGCCTGCGCTGGATtg       <  1:2447868/100‑1 (MQ=255)
                                                                                             aaGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATtgt      <  2:354439/100‑1 (MQ=255)
                                                                                               ggCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGa    <  1:2491633/100‑1 (MQ=255)
                                                                                                gCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGAt   <  2:2331115/100‑1 (MQ=255)
                                                                                                 cTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGAt   <  1:1309199/99‑1 (MQ=255)
                                                                                                  |                                                                                                  
TGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGATT  >  NC_000913/2589027‑2589223

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