Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,589,125 T→C 15.8% 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.8% 334.3 / 53.6 146F511L (TTT→CTT) acrDaminoglycoside/multidrug efflux system
Reads supporting (aligned to +/- strand):  ref base T (57/66);  new base C (13/10);  total (70/76)
Fisher's exact test for biased strand distribution p-value = 4.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.94e-01

Only 100 of 154 total aligned reads displayed.

ATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGATT  >  NC_000913/2589026‑2589223
                                                                                                   |                                                                                                  
aTGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCt                                                                                                    <  2:708834/100‑1 (MQ=255)
aTGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCt                                                                                                    >  2:533725/1‑100 (MQ=255)
    tGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTTAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCttttt                                                                                                >  2:2448869/1‑100 (MQ=255)
    tGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGAGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCttttt                                                                                                >  2:2084001/1‑100 (MQ=255)
       tGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAGAGGCTTTTTTGc                                                                                             <  1:2506605/100‑1 (MQ=255)
       tGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGc                                                                                             >  2:719685/1‑100 (MQ=255)
       tGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGc                                                                                             >  1:2777689/1‑100 (MQ=255)
        gTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGcc                                                                                            <  2:2522436/100‑1 (MQ=255)
        gTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAAGCCTTTTTGcc                                                                                            <  1:1250767/100‑1 (MQ=255)
         tCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTCTTTGCCt                                                                                           <  1:1999572/100‑1 (MQ=255)
         tCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAAGCCTTTTTGccc                                                                                           <  2:2952727/100‑2 (MQ=255)
          cAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTg                                                                                          >  2:1299487/1‑100 (MQ=255)
          cAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAAGCCTTTTTGccc                                                                                           <  1:2700603/99‑2 (MQ=255)
           aGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTgn                                                                                         <  1:490929/100‑2 (MQ=255)
           aGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTgg                                                                                         >  1:1311005/1‑100 (MQ=255)
           aGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAACCCTTTTTGcccat                                                                                         <  1:2566362/100‑4 (MQ=255)
             tACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGtt                                                                                       <  1:2497533/100‑1 (MQ=255)
              aCTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATAGGCAAAAAAGCCTTTTTGcccatgat                                                                                      <  1:101965/100‑7 (MQ=38)
              aCTGGTAGCGATCATCCTCGGTACGGCTCTGTGTGCCACACTAACTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTCCCTGGtt                                                                                       <  2:453431/99‑1 (MQ=255)
                tGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATAGGCAAAAAAGCCTTTTTGcccatgatgc                                                                                    <  1:3093919/100‑9 (MQ=37)
                  gTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAAcc                                                                                  <  2:1936262/100‑1 (MQ=255)
                   tAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCa                                                                                 >  1:3549623/1‑100 (MQ=255)
                    aGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAg                                                                                >  1:1311809/1‑100 (MQ=255)
                     gCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGa                                                                               <  1:2159065/100‑1 (MQ=255)
                       gatgatCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATg                                                                             >  1:755781/1‑100 (MQ=255)
                        atgatCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGt                                                                            >  1:1570895/1‑100 (MQ=255)
                         tgatCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATg                                                                             <  2:1833005/98‑1 (MQ=255)
                           atCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTa                                                                         <  1:263475/100‑1 (MQ=255)
                            tCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTa                                                                         <  1:1452939/99‑1 (MQ=255)
                               tCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCg                                                                     <  2:1570885/100‑1 (MQ=255)
                               tCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCg                                                                     <  2:1570852/100‑1 (MQ=255)
                                cattCCGGCTCTGTCTGCCACACTACTTAAGTCACTTAAAAAACGTGAGCATCATGGGCAAAAATGCTTTTTTGCCTGGTTTAACCAGATGTTTAAc                                                                       <  2:207225/94‑1 (MQ=255)
                                cACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGc                                                                    >  1:3184716/1‑100 (MQ=255)
                                cACTCCGGCTCTGTGTACCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGc                                                                    >  2:2243846/1‑100 (MQ=255)
                                 aCTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCa                                                                   >  2:386498/1‑100 (MQ=255)
                                 aCTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCa                                                                   >  1:545585/1‑100 (MQ=255)
                                  cTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACAGCaa                                                                  >  2:3426270/1‑100 (MQ=255)
                                      ggCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGcc                                                              <  2:3262059/100‑1 (MQ=255)
                                       gCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCg                                                             >  1:1446382/1‑100 (MQ=255)
                                       gCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCg                                                             <  1:3249581/100‑1 (MQ=255)
                                          ctGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAAc                                                          <  2:1735181/100‑1 (MQ=255)
                                           tgtgtgCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGTTTTTTTTAACAAGATGTTTAACAGCAACGCCGAACg                                                         >  2:682721/1‑100 (MQ=255)
                                           tgtgtgCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACg                                                         <  2:201688/100‑1 (MQ=255)
                                           tgtgtgCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACg                                                         <  2:2283739/100‑1 (MQ=255)
                                           tgtgtgCCACACTACTTAAGCCAAAGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTAGTTTAACAAGATGTTTAAACGCAACACCGaa                                                           >  2:3009400/1‑98 (MQ=255)
                                            gtgtgCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGc                                                        >  1:2054654/1‑100 (MQ=255)
                                                gCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACg                                                    <  2:2866484/100‑1 (MQ=255)
                                                gCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACg                                                    >  2:2515404/1‑100 (MQ=255)
                                                 ccACACTACTTAAGCCACTGAAAAAAGGTTAGCATCATGGTCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACg                                                    <  2:1131419/99‑1 (MQ=255)
                                                  cacaCTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGaa                                                  >  1:2366230/1‑100 (MQ=255)
                                                    cacTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGaaaa                                                <  2:3454755/100‑1 (MQ=255)
                                                        actTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAggg                                            >  2:3667118/1‑100 (MQ=255)
                                                         ctTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAgggg                                           <  1:3428973/100‑1 (MQ=255)
                                                           tAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTg                                         >  1:1872437/1‑100 (MQ=255)
                                                           tAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTg                                         >  2:2746520/1‑100 (MQ=255)
                                                             aGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTgg                                        <  1:414264/99‑1 (MQ=255)
                                                             aGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGc                                       <  1:879050/100‑1 (MQ=255)
                                                             aGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGc                                       <  2:3326092/100‑1 (MQ=255)
                                                             aGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGc                                       <  2:3785088/100‑1 (MQ=255)
                                                                 aCTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGaaa                                   >  2:2801048/1‑100 (MQ=255)
                                                                 aCTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGaaa                                   <  1:2377124/100‑1 (MQ=255)
                                                                 aCTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGaaa                                   >  2:1766514/1‑100 (MQ=255)
                                                                  cTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTTTTTTTTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGaaaa                                  >  2:175540/1‑100 (MQ=255)
                                                                  cTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGaaaa                                  >  1:46286/1‑100 (MQ=255)
                                                                   tGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAAt                                 <  1:409105/100‑1 (MQ=255)
                                                                   tGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAAt                                 <  1:2181707/100‑1 (MQ=255)
                                                                    gAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAAtt                                >  1:776095/1‑100 (MQ=255)
                                                                     aaaaaaGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATtc                               >  2:2804152/1‑100 (MQ=255)
                                                                     aaaaaaGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACCCTACGAAAAAGGGGTGGCGAAAATtc                               <  1:2316555/100‑1 (MQ=255)
                                                                       aaaaGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATtctc                             >  2:200757/1‑100 (MQ=255)
                                                                        aaaGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAATGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTcc                            <  1:2475552/100‑1 (MQ=255)
                                                                          aGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCAc                          <  1:671410/100‑1 (MQ=255)
                                                                           ggTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCAcc                         >  1:3431622/1‑100 (MQ=255)
                                                                           ggTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCAcc                         <  1:669939/100‑1 (MQ=255)
                                                                                gCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGc                    >  2:3797717/1‑100 (MQ=255)
                                                                                gCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGc                    >  1:1004605/1‑100 (MQ=255)
                                                                                 catcatGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGcc                   <  2:2054642/100‑1 (MQ=255)
                                                                                    catGGGCAAAAAGGCTTTTTTTTTTTTTTTTACCAGATGTTTAACGCCAACAGAGAACGCTACGAAGAAGGGGTGGCGAAAATTCTCCCCCGTAGCCTg                 >  2:1698933/1‑99 (MQ=255)
                                                                                    accaggCAAAAAAGCCTTTTTGCCCGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTgc                >  2:3218775/5‑100 (MQ=255)
                                                                                     ccaggCAAAAAAGCCTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCtacg               >  1:1893800/4‑97 (MQ=255)
                                                                                     ccaggCAAAAAAGCCTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTgcg               >  1:367800/4‑100 (MQ=255)
                                                                                     ccaggCAAAAAAGCCTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTgcg               >  2:426581/4‑100 (MQ=255)
                                                                                     ccaggCAAAAAAGCCTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTgcg               >  1:3521817/4‑100 (MQ=255)
                                                                                      caggCAAAAAAGCCTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTgcgc              >  2:2827620/3‑100 (MQ=255)
                                                                                      caggCAAAAAAGCCTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTgcgc              >  1:1618862/3‑100 (MQ=255)
                                                                                      caggCAAAAAAGCCTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTgcgc              >  1:3135202/3‑100 (MQ=255)
                                                                                       gggCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCt             >  1:3097280/1‑100 (MQ=255)
                                                                                       aggCAAAAAAGCCTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCt             >  2:1400301/2‑100 (MQ=255)
                                                                                       aggCAAAAAAGCCTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCt             >  1:3282444/2‑100 (MQ=255)
                                                                                       aggCAAAAAAGCCTTTTTGCCTGGTTTAACCAGATGTTGAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCt             >  2:2348676/2‑100 (MQ=255)
                                                                                        ggCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTCACCGCAACGCCGAACGCTATGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTg            <  2:2083915/100‑1 (MQ=255)
                                                                                        ggCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTg            <  1:213133/100‑1 (MQ=255)
                                                                                        ggCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTg            >  2:1835306/1‑100 (MQ=255)
                                                                                          cAAAAAGGCTTTTTTGCGTGGTTTAACCAGATGTTTCACCGCAACGACAAAAAAAAAAAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGa          <  2:3149059/100‑1 (MQ=255)
                                                                                             aaaGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATtg       <  1:1068500/100‑1 (MQ=255)
                                                                                              aaGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATtgt      <  2:458182/100‑1 (MQ=255)
                                                                                                 gCTTTTTTGCCAGGTTTAACCAGATGTTTAACCGCAACGACGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCGCCGTAGCCTGCGCTGGATTGTGa    <  2:3139232/99‑1 (MQ=255)
                                                                                                  cTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGTGCTGGATTGTGAtt  <  2:1235017/100‑1 (MQ=255)
                                                                                                  cTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGAtt  >  1:2547677/1‑100 (MQ=255)
                                                                                                   ttttttGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGAtt  <  1:475846/99‑1 (MQ=255)
                                                                                                   |                                                                                                  
ATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGATT  >  NC_000913/2589026‑2589223

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