Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 484,879 G→C 5.6% I247M (ATC→ATG) ‡ acrA ← multidrug efflux system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913484,8790GC 5.6% 1163.2 / 42.8 318I247M (ATC→ATG) ‡acrAmultidrug efflux system
Reads supporting (aligned to +/- strand):  ref base G (146/154);  new base C (13/5);  total (159/159)
Fisher's exact test for biased strand distribution p-value = 8.69e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 320 total aligned reads displayed.

GACGTGCGCGCACGAACATACCCGGCAGCAGAGTGTGATCCGGGTTCGGGAAGATAGCGCGTAGGGTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCAGGAAGTCGTTGCTGGACTGGGTCACATCAACGTAGATCGGATCAAGTTGCTGCACGGTTGCCAGCGCAGTCGCCTGA  >  NC_000913/484731‑485026
                                                                                                                                                    |                                                                                                                                                   
gACGTGCGCGCACGAACATACCCGGCAGCAGAGTGTGATCCGGGTTCGGGAAGATAGCGCGTAGGGTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTg                                                                                                                                                     >  3:1214974/1‑149 (MQ=255)
      cgcgcACGAACATACCCGGCAGCAGAGTGTGATCCGGGTTCGGGAAGATAGCGCGTAGGGTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGt                                                                                                                                               <  4:1760834/149‑1 (MQ=255)
       gcgcACGAACATACCCGGCAGCAGAGTGTGATCCGGGTTCGGGAAGATAGCGCGTAGGGTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTg                                                                                                                                              <  4:1307410/149‑1 (MQ=255)
             gAACATACCCGGCAGCAGAGTGTGATCCGGGTTCGGGAAGATAGCGCGTAGGGTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGt                                                                                                                                               >  2:904131/1‑142 (MQ=255)
             gAACATACCCGGCAGCAGAGTGTGATCCGGGTTCGGGAAGATAGCGCGTAGGGTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGt                                                                                                                                               >  2:404514/1‑142 (MQ=255)
             gAACATACCCGGCAGCAGAGTGTGATCCGGGTTCGGGAAGATAGCGCGTAGGGTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACtt                                                                                                                                        >  2:1490256/1‑149 (MQ=255)
              aaCATACCCGGCAGCAGAGTGTGATCCGGGTTCGGGAAGATAGCGCGTAGGGTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACttt                                                                                                                                       <  2:962719/149‑1 (MQ=255)
               aCATACCCGGCAGCAGAGTGTGATCCGGGTTCGGGAAGATAGCGCGTAGGGTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGt                                                                                                                                               >  2:229644/1‑140 (MQ=255)
                     ccGGCAGCAGAGTGTGATCCGGGTTCGGGAAGATAGCGCGTAGGGTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTg                                                                                                                                <  3:868989/149‑1 (MQ=255)
                     ccGGCAGCAGAGTGTGATCCGGGTTCGGGAAGATAGCGCGTAGGGTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTg                                                                                                                                >  3:1326013/1‑149 (MQ=255)
                      cGGCAGCAGAGTGTGATCCGGGTTCGGGAAGATAGCGCGTAGGGTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGc                                                                                                                               <  1:302915/149‑1 (MQ=255)
                      cGGCAGCAGAGTGTGATCCGGGTTCGAGAAGATAGCGCGTAGGGTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCt                                                                                                                                   <  1:1789229/145‑1 (MQ=255)
                        gcagcaGAGGGTGATCCGGGTTCGGGAAGATAGCGCGTAGGGTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACCGGTGATCAGTGACACTTTGGCTTTGCCg                                                                                                                             >  4:1629449/1‑149 (MQ=255)
                           gcagAGTGTGATCCGGGTTCGGGAAGATAGCGCGTAGGGTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTtc                                                                                                                          <  4:1114739/149‑1 (MQ=255)
                              gagTGTGATCCGGGTTCGGGAAGATAGCGCGTAGGGTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGa                                                                                                                                             <  1:1049119/127‑1 (MQ=255)
                              gagTGTGATCCGGGTTCGGGAAGATAGCGCGTAGGGTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGa                                                                                                                                             >  2:1049119/1‑127 (MQ=255)
                               agTGTGATCCGGGTTCGGGAAGATAGCGCGTAGGGTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCtt                                                                                                                      <  4:68053/149‑1 (MQ=255)
                                     aTCCGGGTTCGGGAAGATAGCGCGTAGGGTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGcc                                                                                                                              >  2:1796312/1‑135 (MQ=255)
                                     aTCCGGGTTCGGGAAGATAGCGCGTAGGGTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCa                                                                                                                >  1:614156/1‑149 (MQ=255)
                                     aTCCGGGTTCGGGAAGATAGCGCGTAGGGTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCa                                                                                                                >  3:1772712/1‑149 (MQ=255)
                                      tCCGGGTTCGGGAAGATAGCGCGTAGGGTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCa                                                                                                                                                 >  1:873820/1‑115 (MQ=255)
                                               gggAAGATAGCGCGTAGGGTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGa                                                                                                                                             <  4:1812044/110‑1 (MQ=255)
                                               gggAAGATAGCGCGTAGGGTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGa                                                                                                                                             >  3:1812044/1‑110 (MQ=255)
                                                 gAAGATAGCGCGTAGGGTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCg                                                                                                    >  1:1204597/1‑149 (MQ=255)
                                                   aGATAGCGCGTAGGGTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGcc                                                                                                  <  1:42833/149‑1 (MQ=255)
                                                              aGGGTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTtct                                                                                                                         >  4:259067/1‑115 (MQ=255)
                                                                ggTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATc                                                                                                                                                  >  4:367338/1‑88 (MQ=255)
                                                                ggTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCa                                                                                     <  1:193561/149‑1 (MQ=255)
                                                                 gTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAg                                                                                    <  4:705925/149‑1 (MQ=255)
                                                                  tGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACGCttt                                                                                                                                       >  3:1873077/1‑97 (MQ=255)
                                                                  tGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACGCttt                                                                                                                                       <  4:1873077/97‑1 (MQ=255)
                                                                    aTAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGgcg                                                                                 >  2:1597785/1‑149 (MQ=255)
                                                                    aTAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGgcg                                                                                 >  4:33507/1‑149 (MQ=255)
                                                                        aCCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCa                                                                               >  1:1348427/1‑147 (MQ=255)
                                                                        aCCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCa                                                                               <  2:1348427/147‑1 (MQ=255)
                                                                          ccAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGtt                                                                                                                   >  3:34804/1‑109 (MQ=255)
                                                                            aGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTtc                                                                                                                          <  1:399423/100‑1 (MQ=255)
                                                                              tGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCAGGAAGTCg                                                                       <  1:1136042/149‑1 (MQ=255)
                                                                               ggTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCAGGAAGTCGt                                                                      <  3:367444/149‑1 (MQ=255)
                                                                                 tCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGc                                                                                                                               <  1:1484399/90‑1 (MQ=255)
                                                                                 tCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCAGGAAGTCGTTg                                                                    >  3:215773/1‑149 (MQ=255)
                                                                                   tGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGt                                                                                                                    <  2:871203/99‑1 (MQ=255)
                                                                                   tGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGc                                                                                                   <  1:534868/116‑1 (MQ=255)
                                                                                    gATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAGCTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCAGGAAGTCGTTGCTg                                                                 <  3:1629449/149‑1 (MQ=255)
                                                                                     aTCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGcc                                                                                                                              <  3:105122/87‑1 (MQ=255)
                                                                                     aTCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCAGGAAGTCGTTGCTg                                                                 <  4:909774/148‑1 (MQ=255)
                                                                                     aTCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCAGGAAGTCGTTGCTg                                                                 >  3:1588075/1‑148 (MQ=255)
                                                                                      tCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGt                                                                                                                    <  2:1201162/96‑1 (MQ=255)
                                                                                      tCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCAGGAAGTCGTTGCTGGa                                                               >  4:1686797/1‑149 (MQ=255)
                                                                                        aaCGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGa                                                                                                                                             >  2:820022/1‑69 (MQ=255)
                                                                                        aaCGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCAGGAAGTCGTTGCTGGACt                                                             <  3:1418294/149‑1 (MQ=255)
                                                                                         aCGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTc                                                                                      >  4:470426/1‑123 (MQ=255)
                                                                                          cGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCa                                                                                     <  2:1788939/123‑1 (MQ=255)
                                                                                          cGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGAGCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCAGGAAGTCGTTGCTGGACTgg                                                           <  1:86387/149‑1 (MQ=255)
                                                                                           ggTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCAGGAAGTCGTTGCTGGACTggg                                                          >  1:1684964/1‑149 (MQ=255)
                                                                                                 gTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCAGGAAGTCGTTGCTg                                                                 <  2:111880/136‑1 (MQ=255)
                                                                                                 gTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCAGGAAGTCGTTGCTGGACTGGGTCACAt                                                    >  1:1852469/1‑149 (MQ=255)
                                                                                                  tCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCAGGAAGTCGTTGCTGGa                                                               <  1:1247051/137‑1 (MQ=255)
                                                                                                    agagAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCAGGaa                                                                           >  3:291549/1‑123 (MQ=255)
                                                                                                    agagAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCAGGAAGTCGTTg                                                                    >  2:511640/1‑130 (MQ=255)
                                                                                                    agagAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCAGGAAGTCGTTGCTGGACTGGGTCACATCaa                                                 >  1:1295280/1‑149 (MQ=255)
                                                                                                    agagAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCAGGAAGTCGTTGCTGGACTGGGTCACATCaa                                                 >  4:820816/1‑149 (MQ=255)
                                                                                                     gagaATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCAGGAAGTCGTTGCTGGACTggg                                                          >  3:689793/1‑139 (MQ=255)
                                                                                                     gagaATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCAGGAAGTCGTTGCTGGACTGGGTCACATCAAc                                                <  3:1686797/149‑1 (MQ=255)
                                                                                                     gagaATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCAGGAAGTCGTTGCTGGACTGGGTCACATCAAc                                                <  3:931942/149‑1 (MQ=255)
                                                                                                      agaATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCAGGAAGTCGTTGCTGGACTGGGTCACATCAACg                                               >  1:1774914/1‑149 (MQ=255)
                                                                                                      agaATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCAGGAAGTCGTTGCTGGACTGGGTCACATCAACg                                               >  1:1405621/1‑149 (MQ=255)
                                                                                                          ttCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCAGGAAGTCGTTGCTGGACTGGGTCACATCAACGTAGa                                           >  1:952858/1‑149 (MQ=255)
                                                                                                            ccAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCAGGAAGTCGTTGCTGGACTGGGTCACATCAACGTAGATc                                         >  2:1292745/1‑149 (MQ=255)
                                                                                                             cAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCAGGAAGTCGTTGCTGGACTGGGTCACATCAACGTAGATCg                                        >  4:1520126/1‑149 (MQ=255)
                                                                                                             cAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCAGGAAGTCGTTGCTGGACTGGGTCACATCAACGTAGATCg                                        >  2:188227/1‑149 (MQ=255)
                                                                                                             cAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCAGGAAGTCGTTGCTGGACTGGGTCACATCAACGTAGATCg                                        >  4:1327679/1‑149 (MQ=255)
                                                                                                               gCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCAGGAAGTCGTTGCTGGACTGGGTCACATCAACGTAGACCGGa                                      <  1:1597785/149‑1 (MQ=255)
                                                                                                                cGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCAGGAAGTCGTTGCTGGACTGGGTCACATCAACGTAGATCGGAt                                     <  1:1745789/149‑1 (MQ=255)
                                                                                                                  tACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCAGGAAGTCGTTGCTGGACTGGGTCACATCAACGTAGATCGGATCa                                   >  4:1422440/1‑149 (MQ=255)
                                                                                                                   aCCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCAGGAAGTCGTTGCTGGACTGGGTCACATCAACGTAGATCGGATCaa                                  <  1:1444091/149‑1 (MQ=255)
                                                                                                                   aCCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCAGGAAGTCGTTGCTGGACTGGGTCACATCAACGTAGATCGGATCaa                                  <  2:540284/149‑1 (MQ=255)
                                                                                                                   aCCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCAGGAAGTCGTTGCTGGACTGGGTCACATCAACGTAGATCGGATCaa                                  <  1:1144673/149‑1 (MQ=255)
                                                                                                                    ccGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCAGGAAGTCGTTGCTGGACTGGGTCACATCAACGTAGATCGGAt                                     >  4:1238572/1‑145 (MQ=255)
                                                                                                                      gTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCAGGAAGTCGTTGCTGGACTGGGTCACATCAACGTAGATCGGATCAAGtt                               <  3:516508/149‑1 (MQ=255)
                                                                                                                         cTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCAGGAAGTCGTTGCTGGACTGGGTCACATCAACGTAGATCGGATCAAGTtgct                            <  2:1405621/149‑1 (MQ=255)
                                                                                                                            cGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCAGGAAGTCGTTGCTGGACTGGGTcaca                                                     <  4:101364/121‑1 (MQ=255)
                                                                                                                             gggAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCAGGAAGTCGTTGCTGGACTGGGTCACATCAACGTAGATCGGATCAAGTTGCTGCAc                        >  4:418211/1‑149 (MQ=255)
                                                                                                                             gggAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCAGGAAGTCGTTGCTGGACTGGGTCACATCAACGTAGATCGAATCAAGTtgc                             <  2:1746588/144‑1 (MQ=255)
                                                                                                                                   ttAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCAGGAAg                                                                          <  2:1716057/93‑1 (MQ=255)
                                                                                                                                   ttAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCAGGAAGTCGTTGCTGGACTGGGTCACATCAACGTAGATCGGATCAAGTTGCTGCACGGTTGc                  >  3:582827/1‑149 (MQ=255)
                                                                                                                                   ttAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCAGGAAGTCGTTGCTGGACTGGGTCACATCAACGTAGATCGGATCAAGTTGCTGCACGGTTGc                  >  2:1460129/1‑149 (MQ=255)
                                                                                                                                     caaagtgTCACTGATCACCAGTGACGCTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCAGGAAGTCGTTGCTGGACTGGGTCACATCAACGTAGATCGGATCAAGTTGCTGCACGGTTGCCa                >  3:255320/7‑149 (MQ=255)
                                                                                                                                      aTGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCAGGAAGTCGTTGCTGGACTGGGTCACATCAACGTAGATCGGATCAAGTTGCTGCACGGTTGc                  >  3:1167264/1‑146 (MQ=255)
                                                                                                                                      aTGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCAGGAAGTCGTTGCTGGACTGGGTCACATCAACGTAGATCGGATCAAGTTGCTGCACGGTTGc                  <  4:1167264/146‑1 (MQ=255)
                                                                                                                                       tGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCAGGAAGTCGTTGCTGGACTGGGTCACATCAACGTAGATCGGATCAAGTTGCTGCACGGTTGCCAgc              <  4:215773/149‑1 (MQ=255)
                                                                                                                                        gCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCAGGAAGTCGTTGCTGGACTGGGTCACATCAAc                                                <  4:551296/114‑1 (MQ=255)
                                                                                                                                          tgTCACTGATCACCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCAGGAAGTCGTTGCTgg                                                                <  2:1140857/95‑1 (MQ=255)
                                                                                                                                          cGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCAGGAAGTCGTTGCTGGACTGGGTCACATCAACGTAGATCGGATCAAGTTGCTGCACGGTTGCCAGCGCa           <  3:1520126/149‑1 (MQ=255)
                                                                                                                                           gTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCAGGAAGTCGTTGCTGGACTGGGTCACATCAACGTAGATCGGATCAAGTTGCTGCACGGTTGCCAGCGCAg          >  1:782852/1‑149 (MQ=255)
                                                                                                                                           gTCACTGATCACCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCAGGAAGTCGTTGCTGGACTGGGTCACATCAACGTAGATCGGATCAAGTtg                              <  3:129378/129‑1 (MQ=255)
                                                                                                                                            tCACTGATCACCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCAGGAAGTCGTTGCTGGACTGGGTCACATCAACGTAGATCGGATCAAGTTGCTGCACGGTTGCCAGCGCAGt         >  3:785184/1‑149 (MQ=255)
                                                                                                                                            tCACTGATCACCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCAGGAAGTCGTTGCTGGACTGGGTCACATCAACGTAGATCGGATCAAGTTGCTGCACGGTTGCCAGCGCAGt         >  2:730347/1‑149 (MQ=255)
                                                                                                                                                 ggTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCAGGAAGTCGTTGCTGGACTGGGTCACATCAACGTAGATCGGATCAAGTTGCTGCACGGTTGCCAGCGCAGTCGCCt    >  3:432053/1‑149 (MQ=255)
                                                                                                                                                  gTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCAGGAAGTCGTTGCTGGACTGGGTCACATCAACGTAGATCGGATCAAGTTGCTGCACGGTTGCCAGCGCAGTCGCCTg   >  4:761652/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GACGTGCGCGCACGAACATACCCGGCAGCAGAGTGTGATCCGGGTTCGGGAAGATAGCGCGTAGGGTGATAGACCCAGTGGTCTGATCAACGGTAACGTCAGAGAATTCCAGCGTACCGTCCTGCGGGAACTTAATGCCGTCACTGGTGATCAGTGACACTTTGGCTTTGCCGTTCTCTTGTTTCAGCGTGCCATTCGCCAGTTCCTGTTTCAGGCGCAGGAAGTCGTTGCTGGACTGGGTCACATCAACGTAGATCGGATCAAGTTGCTGCACGGTTGCCAGCGCAGTCGCCTGA  >  NC_000913/484731‑485026

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