Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,589,125 T→C 12.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,1250TC12.2% 226.7 / 18.5 82F511L (TTT→CTT) acrDaminoglycoside/multidrug efflux system
Reads supporting (aligned to +/- strand):  ref base T (40/30);  new base C (5/5);  total (47/35)
Fisher's exact test for biased strand distribution p-value = 7.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

CACCGGTGCCATCTATCGCCAGTTCTCTATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGATTTATGTCCTGCTGCTTGGCGGCATGGTGTTCCTGTTCCTGCGTTTGCCGAC  >  NC_000913/2588980‑2589273
                                                                                                                                                 |                                                                                                                                                    
caccGGTGCCATCTATCGCCAGTTCTCTATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCtttt                                                                                                                                                   >  3:88812/1‑149 (MQ=255)
    ggTGCCATCTATCGCCAGTTCTCTATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCGGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGc                                                                                                                                               <  1:13392/149‑1 (MQ=255)
     gTGCCATCTATCGCCAGTTCTCTATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGcc                                                                                                                                              >  1:218676/1‑149 (MQ=255)
      tGCCATCTATCGCCAGTTCTCTATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCt                                                                                                                                             <  1:114042/149‑1 (MQ=255)
      tGCCATCTATCGCCAGTTCTCTATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAAGCCTTTTTGccc                                                                                                                                             <  2:550168/149‑2 (MQ=255)
      tGCCATCTATCGCCAGTTCTCTATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGACTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAAGCCTTTTTGccc                                                                                                                                             <  1:65467/149‑2 (MQ=255)
        ccATCTATCGCCAGTTCTCTATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTgg                                                                                                                                           >  3:184328/1‑149 (MQ=255)
        ccATCTATCGCCAGTTCTCTATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTgg                                                                                                                                           >  2:174416/1‑149 (MQ=255)
         cATCTATCGCCAGTTCTCTATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCTACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAAGCCTTTTTGccc                                                                                                                                             >  1:118798/1‑145 (MQ=255)
         cATCTATCGCCAGTTCTCTATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCTACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAAGCCTTTTTGccc                                                                                                                                             <  2:118798/146‑2 (MQ=255)
          aTCTATCGCCAGTTCTCTATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGtt                                                                                                                                         >  3:554250/1‑149 (MQ=255)
          aTCTATCGCCAGTTCTCTATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCGTTTTTGCCTGGtt                                                                                                                                         >  4:437337/1‑149 (MQ=255)
           tCTATCGCCAGTTCTCTATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGttt                                                                                                                                        <  1:170871/149‑1 (MQ=255)
              aTCGCCAGTTCTCTATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAAc                                                                                                                                     <  1:313121/149‑1 (MQ=255)
                cGCCAGTTCTCTATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAAcaa                                                                                                                                   <  4:473755/149‑3 (MQ=255)
                      ttCTCTATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGtt                                                                                                                             >  2:298021/1‑149 (MQ=255)
                         tctATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTaa                                                                                                                          >  3:59999/1‑149 (MQ=255)
                         tctATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTaa                                                                                                                          >  2:537372/1‑149 (MQ=255)
                         tctATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTaa                                                                                                                          >  3:619544/1‑149 (MQ=255)
                           tATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAAcc                                                                                                                        >  3:617154/1‑149 (MQ=255)
                                ccATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAAc                                                                                                                   <  1:174416/149‑1 (MQ=255)
                                  aTTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGt                                                                                                                                          <  2:103181/124‑1 (MQ=255)
                                  aTTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGt                                                                                                                                          >  1:103181/1‑124 (MQ=255)
                                    tgttgCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCg                                                                                                               <  4:88812/149‑1 (MQ=255)
                                         cggcgATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGc                                                                                                          >  4:566412/1‑149 (MQ=255)
                                          ggcgATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCt                                                                                                         >  2:435815/1‑149 (MQ=255)
                                          ggcgATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCt                                                                                                         <  1:121595/149‑1 (MQ=255)
                                             gATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATAGGCAAAAAAGCCTTTTTGccc                                                                                                                                             >  3:97512/1‑109 (MQ=255)
                                             gATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATAGGCAAAAAAGCCTTTTTGccc                                                                                                                                             <  4:97512/110‑2 (MQ=255)
                                                       tCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAAc                                                                                                                                     >  4:181001/1‑108 (MQ=255)
                                                       tCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAAc                                                                                                                                     <  3:181001/108‑1 (MQ=255)
                                                          gTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGc                                                                                         >  3:367368/1‑149 (MQ=255)
                                                                     gatgatCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGAt                                                                                                                                <  4:363731/99‑1 (MQ=255)
                                                                     gatgatCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGAt                                                                                                                                >  3:363731/1‑99 (MQ=255)
                                                                        gatCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAAc                                                                                                                                     <  1:207665/91‑1 (MQ=255)
                                                                        gatCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAAc                                                                                                                                     >  2:207665/1‑91 (MQ=255)
                                                                         atCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCg                                                                          >  4:587931/1‑149 (MQ=255)
                                                                          tCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGt                                                                         >  2:401732/1‑149 (MQ=255)
                                                                           ccTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCGTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTa                                                                        >  2:394802/1‑149 (MQ=255)
                                                                             tCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGc                                                                      >  3:651362/1‑149 (MQ=255)
                                                                               aCTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCt                                                                    >  3:116514/1‑149 (MQ=255)
                                                                                    ggCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCt                                                               >  2:223413/1‑149 (MQ=255)
                                                                                     gCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTg                                                              >  2:77368/1‑149 (MQ=255)
                                                                                     gCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTg                                                              >  3:358337/1‑149 (MQ=255)
                                                                                      ctctGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTgg                                                             >  2:207347/1‑149 (MQ=255)
                                                                                       tctGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGa                                                            <  4:651362/149‑1 (MQ=255)
                                                                                         tgtgtgCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGAtt                                                          <  4:358337/149‑1 (MQ=255)
                                                                                         tgtgtgCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGAtt                                                          >  4:669814/1‑149 (MQ=255)
                                                                                             tgCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGa                                                      <  4:554250/149‑1 (MQ=255)
                                                                                             tgCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGa                                                      <  1:499975/149‑1 (MQ=255)
                                                                                             tgCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGa                                                      <  2:399218/149‑1 (MQ=255)
                                                                                             tgCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGa                                                      <  1:298021/149‑1 (MQ=255)
                                                                                                cacaCTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGAttt                                                   <  4:367368/149‑1 (MQ=255)
                                                                                                cacaCTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGAttt                                                   >  1:28822/1‑149 (MQ=255)
                                                                                                     tactTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGATTTATGTc                                              <  3:566412/149‑1 (MQ=255)
                                                                                                       ctTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGATTTATGTCct                                            >  3:199330/1‑149 (MQ=255)
                                                                                                            gCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGATTTATGTCCtgctgc                                       <  4:619544/149‑1 (MQ=255)
                                                                                                              cACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGATTTATGTCCTGCTGCtt                                     >  4:588339/1‑149 (MQ=255)
                                                                                                               aCTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGATTTATGTCCTGCTGCTTg                                    >  1:238044/1‑149 (MQ=255)
                                                                                                                     aaaaGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGATTTATGTCCTGCTGCTTGGCGGCa                              <  4:577801/149‑1 (MQ=255)
                                                                                                                      aaaGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGATTTATGTCCTGCTGCTTGGCGGCAt                             <  3:108869/149‑1 (MQ=255)
                                                                                                                       aaGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGATTTATGTCCTGCTGCTTGGCGGCATg                            <  1:401732/149‑1 (MQ=255)
                                                                                                                          gTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGATTTATGTCCTGCTGCTTGGCGGCATGgtg                         >  3:599945/1‑149 (MQ=255)
                                                                                                                             aGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGATTTATGTCCTGCTGCTTGGCGGCATGGtgttc                      <  2:238044/149‑1 (MQ=255)
                                                                                                                               catcatGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGATTTATGTCCTGCTGCTTGGCGGCATGGtgttcct                    >  4:520924/1‑149 (MQ=255)
                                                                                                                                 tcatGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGCGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGATTTATGTCCTGCTGCTTGGCGGCATGGtgttcctgt                  <  4:184328/149‑1 (MQ=255)
                                                                                                                                   atGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGATTTATGTCCTGCTGCTTGGCGGCATGGtgttcctgttc                >  2:386417/1‑149 (MQ=255)
                                                                                                                                    caggCAAAAAAGCCTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGCTTGTGATTTATGTCCTGCTGCTTGGCGGCATGGtgttcctgttcc               >  4:157672/3‑149 (MQ=255)
                                                                                                                                    caggCAAAAAAGCCTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGATTTATGTCCTGCTGCTTGGCGGCATGGtgttcctgttcc               >  2:19808/3‑149 (MQ=255)
                                                                                                                                      ggCAAAAAAGCCTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGa                                                            >  1:45279/1‑102 (MQ=255)
                                                                                                                                      ggCAAAAAAGCCTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGa                                                            <  2:45279/102‑1 (MQ=255)
                                                                                                                                       gCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGATTTATGTCCTGCTGCTTGGCGGCATGGTGTTCCTGTTCCTGc            >  1:295938/1‑149 (MQ=255)
                                                                                                                                        cAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGATTTATGTCCTGCTGCTTGGCGGCATGGTGTTCCTGTTCCTGCg           >  4:551096/1‑149 (MQ=255)
                                                                                                                                         aaaaaGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGATTTATGTCCTGCTGCTTGGCGGCATg                            >  4:286997/1‑131 (MQ=255)
                                                                                                                                         aaaaaGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGATTTATGTCCTGCTGCTTGGCGGCATg                            <  3:286997/131‑1 (MQ=255)
                                                                                                                                         aaaaaGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGATTTATGTCCTGCTGCTTGGCGGCATGGTGTTCCTGTTCCTGCGt          >  2:301382/1‑149 (MQ=255)
                                                                                                                                          aaaaGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGATTTATGTCCTGCTGCTTGGCGGCATGGTGTTCCTGTTCCTGCGtt         <  4:500694/149‑1 (MQ=255)
                                                                                                                                           aaaGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGATTTATGTCCTGCTGCTTGGCGGCATGGTGTTCCTGTTCCTGCGttt        >  2:171140/1‑149 (MQ=255)
                                                                                                                                           aaaGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGATTTATGTCCTGCTGCTTGGCGGCATGGTGTTCCTGTTCCTGCGttt        >  2:430835/1‑149 (MQ=255)
                                                                                                                                             aGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGATTTATGTCCTGCTGCTTGGCGGCATGGTGTTCCTGTTCCTGCGTTTGc      <  4:617154/149‑1 (MQ=255)
                                                                                                                                             aGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGATTTATGTCCTGCTGCTTGGCGGCATGGTGTTCCTGTTCCTGCGTTTGc      <  3:551096/149‑1 (MQ=255)
                                                                                                                                               gCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGATTTATGTCCTGCTGCTTGGCGGCATGGTGTTCCTGTTCCTGCGTTTGCCg    >  4:154539/1‑149 (MQ=255)
                                                                                                                                               gCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGATTTATGTCCTGCTGCTTGGCGGCATGGTGTTCCTGTTCCTGCGTTTGCCg    >  3:26661/1‑149 (MQ=255)
                                                                                                                                                 ttttttGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGATTTATGTCCTGCTGCTTGGCGGCATGGTGTTCCTGTTCCTGCGTTTGCCGAc  <  4:59999/149‑1 (MQ=255)
                                                                                                                                                 ttttttGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGATTTATGTCCTGCTGCTTGGCGGCATGGTGTTCCTGTTCCTGCGTTTGCCGAc  >  2:465287/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
CACCGGTGCCATCTATCGCCAGTTCTCTATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGATTTATGTCCTGCTGCTTGGCGGCATGGTGTTCCTGTTCCTGCGTTTGCCGAC  >  NC_000913/2588980‑2589273

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: