Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,589,122 G→A 15.7% 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.7% 356.9 / 50.3 147G510S (GGC→AGC) acrDaminoglycoside/multidrug efflux system
Reads supporting (aligned to +/- strand):  ref base G (60/63);  new base A (13/10);  total (73/74)
Fisher's exact test for biased strand distribution p-value = 6.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.51e-01

Only 100 of 152 total aligned reads displayed.

GCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGA  >  NC_000913/2589023‑2589221
                                                                                                   |                                                                                                   
gcgATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAg                                                                                                     <  2:2886155/100‑1 (MQ=255)
 cgATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAgg                                                                                                    <  1:508082/100‑1 (MQ=255)
   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)
          tGTCAGTACTGGTAGCGATGATCCTCACTCCGGGTCTGTGTGCCACACTACTTAAAAAAAAAAAAAAAGGTGAGCATCATGGGCAAAAAAGCCTTTTTGc                                                                                           <  2:3287674/100‑1 (MQ=255)
          tGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAGAGGCTTTTTTGc                                                                                           <  1:2506605/100‑1 (MQ=255)
          tGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGc                                                                                           >  2:2299473/1‑100 (MQ=255)
          tGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGc                                                                                           >  1:2777689/1‑100 (MQ=255)
          tGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGc                                                                                           >  2:2031173/1‑100 (MQ=255)
          tGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGc                                                                                           >  2:719685/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)
              aGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTgn                                                                                       <  1:490929/100‑2 (MQ=255)
              aGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTgg                                                                                       >  1:1311005/1‑100 (MQ=255)
              aGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAACCCTTTTTGcccat                                                                                       <  1:2566362/100‑4 (MQ=255)
               gTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGt                                                                                      >  2:321605/1‑100 (MQ=255)
                tACTGGTAGCGCTGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAAGCCTTTTTGcccatga                                                                                     <  2:2076277/100‑6 (MQ=255)
                tACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGtt                                                                                     <  1:2497533/100‑1 (MQ=255)
                 aCTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATAGGCAAAAAAGCCTTTTTGcccatgat                                                                                    <  1:101965/100‑7 (MQ=38)
                 aCTGGTAGCGATCATCCTCGGTACGGCTCTGTGTGCCACACTAACTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTCCCTGGtt                                                                                     <  2:453431/99‑1 (MQ=255)
                     gTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAAcc                                                                                <  1:2636877/100‑1 (MQ=255)
                     gTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAAcc                                                                                <  2:2805426/100‑1 (MQ=255)
                     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:1570862/1‑100 (MQ=255)
                           atgatCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGt                                                                          >  1:1570895/1‑100 (MQ=255)
                           atgatCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGt                                                                          <  1:860500/100‑1 (MQ=255)
                            tgatCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATg                                                                           <  2:1833005/98‑1 (MQ=255)
                             gatCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGttt                                                                        >  2:3428947/1‑100 (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)
                                    aCTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCa                                                                 >  2:386498/1‑100 (MQ=255)
                                          gCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCg                                                           >  1:1446382/1‑100 (MQ=255)
                                             ctGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAAc                                                        <  2:1735181/100‑1 (MQ=255)
                                              tgtgtgCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGTTTTTTTTAACAAGATGTTTAACAGCAACGCCGAACg                                                       >  2:682721/1‑100 (MQ=255)
                                              tgtgtgCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACg                                                       <  2:2283739/100‑1 (MQ=255)
                                              tgtgtgCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACg                                                       <  2:201688/100‑1 (MQ=255)
                                              tgtgtgCCACACTACTTAAGCCAAAGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTAGTTTAACAAGATGTTTAAACGCAACACCGaa                                                         >  2:3009400/1‑98 (MQ=255)
                                               gtgtgCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGc                                                      >  1:2054654/1‑100 (MQ=255)
                                                 gtgCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTa                                                    >  1:2968377/1‑100 (MQ=255)
                                                   gCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACg                                                  <  2:2866484/100‑1 (MQ=255)
                                                   gCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACg                                                  <  1:668866/100‑1 (MQ=255)
                                                   gCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACg                                                  <  1:726412/100‑1 (MQ=255)
                                                    ccACACTACTTAAGCCACTGAAAAAAGGTTAGCATCATGGTCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACg                                                  <  2:1131419/99‑1 (MQ=255)
                                                    ccACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACg                                                  >  1:3514382/1‑99 (MQ=255)
                                                     cacaCTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGaa                                                >  1:2366230/1‑100 (MQ=255)
                                                       cacTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGaaaa                                              <  2:3454755/100‑1 (MQ=255)
                                                          tacaTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAgg                                           <  2:2506206/100‑1 (MQ=255)
                                                           actTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAggg                                          >  2:3667118/1‑100 (MQ=255)
                                                              tAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTg                                       >  2:2746520/1‑100 (MQ=255)
                                                              tAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTg                                       >  1:1872437/1‑100 (MQ=255)
                                                                aGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGc                                     <  2:3785088/100‑1 (MQ=255)
                                                                aGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGc                                     <  2:3326092/100‑1 (MQ=255)
                                                                aGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGc                                     >  2:1130779/1‑100 (MQ=255)
                                                                    aCTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGaaa                                 >  2:1766514/1‑100 (MQ=255)
                                                                    aCTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGaaa                                 <  1:2377124/100‑1 (MQ=255)
                                                                     cTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTTTTTTTTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGaaaa                                >  2:175540/1‑100 (MQ=255)
                                                                     cTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGaaaa                                >  1:46286/1‑100 (MQ=255)
                                                                     cTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGaaaa                                >  1:1937004/1‑100 (MQ=255)
                                                                      ttaaaaaaGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGCTGTTTAACCTCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAAt                               <  2:2113555/98‑1 (MQ=255)
                                                                      tGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAAt                               <  1:409105/100‑1 (MQ=255)
                                                                      tGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAAt                               <  1:2181707/100‑1 (MQ=255)
                                                                       gAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAAtt                              >  1:380519/1‑100 (MQ=255)
                                                                        aaaaaaGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATtc                             >  2:2804152/1‑100 (MQ=255)
                                                                        aaaaaaGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATtc                             <  1:1631094/100‑1 (MQ=255)
                                                                        aaaaaaGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACCCTACGAAAAAGGGGTGGCGAAAATtc                             <  1:2316555/100‑1 (MQ=255)
                                                                         aaaaaGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATtct                            <  1:3388523/100‑1 (MQ=255)
                                                                           aaaGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAATGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTcc                          <  1:2475552/100‑1 (MQ=255)
                                                                              ggTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCAcc                       <  1:669939/100‑1 (MQ=255)
                                                                              ggTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCAcc                       >  1:3431622/1‑100 (MQ=255)
                                                                                 gAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTa                    >  1:455334/1‑100 (MQ=255)
                                                                                   gCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGc                  >  1:1004605/1‑100 (MQ=255)
                                                                                   gCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGc                  >  2:3797717/1‑100 (MQ=255)
                                                                                    catcatGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGcc                 <  2:1954702/100‑1 (MQ=255)
                                                                                     atcatGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCt                <  1:3189750/100‑1 (MQ=255)
                                                                                       catGGGCAAAAAGGCTTTTTTTTTTTTTTTTACCAGATGTTTAACGCCAACAGAGAACGCTACGAAGAAGGGGTGGCGAAAATTCTCCCCCGTAGCCTg               >  2:1698933/1‑99 (MQ=255)
                                                                                       accaggCAAAAAAGCCTTTTTGCCCGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTgc              >  2:3218775/5‑100 (MQ=255)
                                                                                        ccaggCAAAAAAGCCTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTgcg             >  2:426581/4‑100 (MQ=255)
                                                                                        ccaggCAAAAAAGCCTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTgcg             >  1:3124206/4‑100 (MQ=255)
                                                                                          gggCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCt           >  1:3097280/1‑100 (MQ=255)
                                                                                          aggCAAAAAAGCCTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCt           >  2:1400301/2‑100 (MQ=255)
                                                                                          aggCAAAAAAGCCTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCt           >  1:3282444/2‑100 (MQ=255)
                                                                                          aggCAAAAAAGCCTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCt           >  2:2778475/2‑100 (MQ=255)
                                                                                           ggCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTg          >  2:1835306/1‑100 (MQ=255)
                                                                                           ggCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTg          <  1:213133/100‑1 (MQ=255)
                                                                                                aaaGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATtg     <  1:1068500/100‑1 (MQ=255)
                                                                                                 aaGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATtgt    <  2:458182/100‑1 (MQ=255)
                                                                                                 aaGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATtgt    <  1:3649467/100‑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 < 32 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: