Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,589,125 T→C 17.2% 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,1250TC17.2% 203.3 / 34.7 93F511L (TTT→CTT) acrDaminoglycoside/multidrug efflux system
Reads supporting (aligned to +/- strand):  ref base T (38/38);  new base C (7/9);  total (46/47)
Fisher's exact test for biased strand distribution p-value = 7.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.02e-01

GTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGAT  >  NC_000913/2589029‑2589222
                                                                                                |                                                                                                 
gTGCTGTCAGTATTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCtttt                                                                                                >  2:1914950/1‑100 (MQ=255)
gTGCTGTCAGTAATGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCtttt                                                                                                >  2:1916790/1‑100 (MQ=255)
 tGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCttttt                                                                                               >  1:190076/1‑100 (MQ=255)
    tGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGAATAAAAAAAAAAAAGGCTTTTTTGc                                                                                            <  1:2779931/100‑1 (MQ=255)
     gTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGcc                                                                                           >  1:2595209/1‑100 (MQ=255)
     gTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGcc                                                                                           <  1:825361/100‑1 (MQ=255)
     gTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAAGCCTTTTTGcc                                                                                           <  1:106196/100‑1 (MQ=255)
     gTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAAGCCTTTTTGcc                                                                                           <  2:2877050/100‑1 (MQ=255)
     gTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAAGCCTTTTTGcc                                                                                           <  2:1591685/100‑1 (MQ=255)
       cAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGACAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTg                                                                                         >  1:1431014/1‑100 (MQ=255)
       cAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAAGCCTTTTTGccca                                                                                         <  2:1800732/100‑3 (MQ=255)
       cAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCAACATAGGAAAAAAAGCCTTTTTGccca                                                                                         <  1:55656/100‑3 (MQ=255)
        aGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTgg                                                                                        <  2:1359565/100‑1 (MQ=255)
        aGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTgg                                                                                        <  2:2162528/100‑1 (MQ=255)
          tACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCtttttttttttttt                                                                                      >  1:1665758/1‑92 (MQ=255)
          tACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGtt                                                                                      <  2:2128084/100‑1 (MQ=255)
          tACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGtt                                                                                      >  2:248345/1‑100 (MQ=255)
            cTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTa                                                                                    <  1:277098/100‑1 (MQ=255)
            cTGGTAGCGATGATCCTCACCCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAAGCCTTTTTGcccatgatg                                                                                    <  2:132083/100‑8 (MQ=38)
             tGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATAGGCAAAAAAGCCTTTTTGcccatgatgc                                                                                   <  2:231890/100‑9 (MQ=37)
              ggTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAAc                                                                                  >  1:2386442/1‑100 (MQ=255)
              ggTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTGAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAAc                                                                                  <  1:1115289/100‑1 (MQ=255)
              ggTAGCGATGATCCTCACTCCGGCTCTGTGAGCAACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGACTGGTTTAAc                                                                                  >  1:1911925/1‑100 (MQ=255)
                tAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCa                                                                                >  1:1727490/1‑100 (MQ=255)
                tAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCa                                                                                <  1:943667/100‑1 (MQ=255)
                 aGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAg                                                                               <  2:1717384/100‑1 (MQ=255)
                 aGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAg                                                                               >  2:51501/1‑100 (MQ=255)
                  gCGATAATCCTCAAACCAGCTCTGTGCGCCACACTACTAAAGCCACTAAAAAAAGGTGAGCATCATGGGCAAAAAAGCCTTTTTGcccatgatgc                                                                                   <  1:2756522/95‑9 (MQ=19)
                    gatgatCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAACAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATg                                                                            >  1:555568/1‑100 (MQ=255)
                    gatgatCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATg                                                                            <  2:1124657/100‑1 (MQ=255)
                    gatgatCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATg                                                                            <  2:540350/100‑1 (MQ=255)
                      tgatCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGtt                                                                          <  1:1201850/100‑1 (MQ=255)
                         tCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTaa                                                                       <  2:1148639/100‑1 (MQ=255)
                          ccTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAAc                                                                      >  1:339115/1‑100 (MQ=255)
                          ccTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAAc                                                                      <  1:1279342/100‑1 (MQ=255)
                             cACTCCGGCTCTGTGTGCCACACTACTTACACCACTGCAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGAGGTTTAACCGc                                                                   >  1:1916276/1‑100 (MQ=255)
                             cACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGc                                                                   >  2:373841/1‑100 (MQ=255)
                              aCTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCa                                                                  >  2:513344/1‑100 (MQ=255)
                                 ccGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACg                                                               <  2:429195/100‑1 (MQ=255)
                                    gCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTGTGCCTGGTTTAACCAGATGTTTAACCGCAACGcc                                                             <  1:2839309/99‑1 (MQ=255)
                                      tctGTGTGCGACACTACTTAAGCCACTGAAAAAAGGTGAGCGTCATGGGCAAAGAGGCGTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGaa                                                          >  1:2035667/1‑100 (MQ=255)
                                       ctGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAAc                                                         <  2:2176471/100‑1 (MQ=255)
                                       ctGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAAc                                                         <  2:2104901/100‑1 (MQ=255)
                                        tgtgtgCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACg                                                        <  1:2426196/100‑1 (MQ=255)
                                           gtgCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTa                                                     >  1:1491867/1‑100 (MQ=255)
                                            tgCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTAc                                                    <  1:1953471/100‑1 (MQ=255)
                                             gCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACg                                                   <  2:1852124/100‑1 (MQ=255)
                                              ccACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGa                                                  >  1:848915/1‑100 (MQ=255)
                                               cacaCTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGaa                                                 >  2:2315079/1‑100 (MQ=255)
                                                acacTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGaaa                                                >  2:1108129/1‑100 (MQ=255)
                                                   ctactTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAg                                             >  2:398979/1‑100 (MQ=255)
                                                    tactTAAGCCACTGAAAAAAGGTGAGCATCATGGCCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGTCGAACGCTACGAAAAAgg                                            <  1:1903620/100‑1 (MQ=255)
                                                     actTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAggg                                           >  2:770007/1‑100 (MQ=255)
                                                        tAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTg                                        >  2:2750050/1‑100 (MQ=255)
                                                        tAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTg                                        >  2:2896135/1‑100 (MQ=255)
                                                         aaGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTGTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGGAAAAGGGGTg                                        >  2:415360/1‑99 (MQ=255)
                                                          aGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGc                                      <  2:2105958/100‑1 (MQ=255)
                                                           gCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCg                                     <  2:2453989/100‑1 (MQ=255)
                                                              aCTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGaaa                                  >  1:528730/1‑100 (MQ=255)
                                                              aCTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGaaa                                  >  2:705514/1‑100 (MQ=255)
                                                                tGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAAt                                <  2:2778100/100‑1 (MQ=255)
                                                                tGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAAt                                <  2:2777903/100‑1 (MQ=255)
                                                                 gAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAAtt                               >  2:2210435/1‑100 (MQ=255)
                                                                 gAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAAtt                               >  2:569410/1‑100 (MQ=255)
                                                                  aaaaaaGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATtc                              >  1:678882/1‑100 (MQ=255)
                                                                   aaaaaGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATtct                             >  1:1883083/1‑100 (MQ=255)
                                                                   aaaaaGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATtct                             <  2:1499165/100‑1 (MQ=255)
                                                                   aaaaaGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATtct                             <  2:6472/100‑1 (MQ=255)
                                                                    aaaaGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATtctc                            >  1:1367822/1‑100 (MQ=255)
                                                                     aaaGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTcc                           >  2:1477593/1‑100 (MQ=255)
                                                                     aaaGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTcc                           <  2:2474577/100‑1 (MQ=255)
                                                                      aaGGTGAGCATCCTGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCa                          >  2:1916125/1‑100 (MQ=255)
                                                                      aaGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTcc                           >  1:1771978/1‑99 (MQ=255)
                                                                        ggTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCAcc                        >  2:2044856/1‑100 (MQ=255)
                                                                          tGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGt                      >  1:28876/1‑100 (MQ=255)
                                                                             gCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAATCGCAACGCCTAACGCTACGATACAGGGGTGGCGAAAATTCTCCACCGTAg                    >  1:683308/1‑99 (MQ=255)
                                                                                 catGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTgc               >  2:786319/1‑100 (MQ=255)
                                                                                  ccaggCAAAAAAGCCTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTgcg              >  2:2533146/4‑100 (MQ=255)
                                                                                  ccaggCAAAAAAGCCTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTgcg              >  1:2342294/4‑100 (MQ=255)
                                                                                  ccaggCAAAAAAGCCTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTgcg              >  2:1599254/4‑100 (MQ=255)
                                                                                  ccaggCAAAAAAGCCTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTgcg              >  2:238817/4‑100 (MQ=255)
                                                                                   tGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTgcg              >  2:746159/1‑99 (MQ=255)
                                                                                   tGGGCAAAAAGGCCTTTTTGGATGGTTTAACCAGATGTTTACCCGTAACGCAAAAAAAAAAAAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTgcgc             <  1:2627097/100‑1 (MQ=255)
                                                                                   caggCAAAAAAGCCTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTgcgc             >  1:2836461/3‑100 (MQ=255)
                                                                                    gggCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCt            <  2:292345/100‑1 (MQ=255)
                                                                                    gggCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCt            <  2:1194773/100‑1 (MQ=255)
                                                                                    aggCAAAAAAGCCTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCt            >  1:1480173/2‑100 (MQ=255)
                                                                                    aggCAAAAAAGCCTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCt            >  2:586342/2‑100 (MQ=255)
                                                                                     ggCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTg           <  2:293822/100‑1 (MQ=255)
                                                                                      gCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAACAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTgg          >  1:2720887/1‑100 (MQ=255)
                                                                                        aaaaaGGCTTTTTTGCCTGGTTTAACCAGATGTTTATCCGCAACGCAAAAAAAAAAAAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGAt        <  1:1641736/100‑1 (MQ=255)
                                                                                        aaaaaGGCTTTTCTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGAt        <  2:584518/100‑1 (MQ=255)
                                                                                         aaaaGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGAtt       >  2:1020471/1‑100 (MQ=255)
                                                                                          aaaGGCTTTTTTGCGTGGTTTAACCAGATGTTTAACCGCAATGCCGAACGCTACGAAAAGGGGGTGGCGAAAATTCTCCACCGTAGGCTGCGCTGGATtg      <  1:1898032/100‑1 (MQ=255)
                                                                                          aaaGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATtg      <  2:2675502/100‑1 (MQ=255)
                                                                                           aaGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATtgt     <  2:2503748/100‑1 (MQ=255)
                                                                                               cTTTATTGCCTGGTTTAACCAGATGACTAACCGTAACGCCGCACGCTACTAATAATGGCTGGCGAAAATTCTCCAGCGTATCCTGCGCTGGATTGTAAt  <  1:700604/99‑1 (MQ=255)
                                                                                                |                                                                                                 
GTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGAT  >  NC_000913/2589029‑2589222

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: