Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,078,603 T→C 6.2% intergenic (‑148/+476) mdtA ← / ← PP_3586 multidrug transporter membrane fusion protein/transposase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,078,6030TC 6.2% 1101.6 / 38.2 338intergenic (‑148/+476)mdtA/PP_3586multidrug transporter membrane fusion protein/transposase
Reads supporting (aligned to +/- strand):  ref base T (171/146);  new base C (12/9);  total (183/155)
Fisher's exact test for biased strand distribution p-value = 8.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.86e-01

Only 100 of 347 total aligned reads displayed.

ATGGGTTGAACGAACCTTGTTCGGAATGACGACTTGAGGAGGATGAACGATAAGCACTGCCATCGTTTTAGCAAAGGGCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGTCTGGCCCTTTCGCGGGTGAACCCGCTTCTATCAAACAAAAGATAACCCACTGTTCTTTCTAGCGAAAAAAACAGGCGCGCGGGGT  >  NC_002947/4078454‑4078751
                                                                                                                                                     |                                                                                                                                                    
   ggTTGAACGAACCTTGTTCGGAATGACGACTTGAGGAGGATGAACGATAAGCACTGCCATCGTTTTAGCAAAGGGCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGtt                                                                                                                                                   >  1:3737338/1‑150 (MQ=255)
    gTTGAACGAACCTTGTTCGGAATGACGACTTGAGGAGGATGAACGATAAGCACTGCCATCGTTTTAGCAAAGGGCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTc                                                                                                                                                  >  1:132660/1‑150 (MQ=255)
    gTTGAACGAACCTTGTTCGGAATGACGACTTGAGGAGGATGAACGATAAGCACTGCCATCGTTTTAGCAAAGGGCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTc                                                                                                                                                  >  2:2049002/1‑150 (MQ=255)
        aacgaacCTTGTTCGGAATGACGACTTGAGGAGGATGAACGATAAGCACTGCCATCGTTTTAGCAAAGGGCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCggg                                                                                                                                              >  2:637765/1‑150 (MQ=255)
        aacgaacCTTGTTCGGAATGACGACTTGAGGAGGATGAACGATAAGCACTGCCATCGTTTTAGCAAAGGGCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCggg                                                                                                                                              >  1:2382550/1‑150 (MQ=255)
          cgaccCTTGTTCGGAATGACGACTTGAGGAGGATGAACGATAAGCACTGCCATCGTTTTAGCAAAGGGCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACACGTTCCGGGcc                                                                                                                                            <  2:397972/150‑1 (MQ=255)
           gacccTTTTTCGGAATGACGACTTGAGGAGGATGAACGATAAGCACTGCCATCGTTTTAGCAAAGGGCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCg                                                                                                                                           <  1:6187800/147‑1 (MQ=255)
               cTTGTTCGGAATGACGACTTGAGGAGGATGAACGATAAGCACTGCCATCGTTTTAGCAAAGGGCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTAc                                                                                                                                       >  1:5971948/1‑150 (MQ=255)
                ttGTTCGGAATGACGACTTGAGGAGGATGAACGATAAGCACTGCCATCGTTTTAGCAAAGGGCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACa                                                                                                                                      >  2:7206630/1‑150 (MQ=255)
                 tGTTCGTAATGACGACTTGAGGAGTATGAACGATAAGCACTGCAATCGTTTTAGCAAAGGGCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAg                                                                                                                                     <  1:3946500/150‑1 (MQ=255)
                     cGGAATGACGACTTGAGGAGGATGAACGATAAGCACTGCCATCGTTTTAGCAAAGGGCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGgtg                                                                                                                                 >  1:7144772/1‑150 (MQ=255)
                     cGGAATGACGACTTGAGGAGGATGAACGATAAGCACTGCCATCGTTTTAGCAAAGGGCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGgtg                                                                                                                                 >  1:4851225/1‑150 (MQ=255)
                        aaTGACGACTTGAGGAGGATGAACGATAAGCACTGCCATCGTTTTAGCAAAGGGCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTg                                                                                                                              >  1:3794241/1‑150 (MQ=255)
                        aaTGACGACTTGAGGAGGATGAACGATAAGCACTGCCATCGTTTTAGCAAAGGGCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTg                                                                                                                              <  2:2810626/150‑1 (MQ=255)
                        aaTGACGACTTGAGGAGGATGAACGATAAGCACTGCCATCGTTTTAGCAAAGGGCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATCTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTg                                                                                                                              >  2:890677/1‑150 (MQ=255)
                              gacTTGAGGAGGATGAACGATAAGCACTGCCATCGTTTTAGCAAAGGGCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGTAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGa                                                                                                                        <  1:5262524/150‑1 (MQ=255)
                              gacTTGAGGAGGATGAACGATAAGCACTGCCATCGTTTTAGAAAAGGGCATTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGCGGCGGATTTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGa                                                                                                                        <  1:662121/150‑1 (MQ=255)
                              gacTTGAGGAGGATGAACGATAAGCACTGCCATCGTTTTAGAAAAGGGACTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGa                                                                                                                        >  1:689965/1‑150 (MQ=255)
                               acTTGAGGAGGATGAACGATAAGCACTGCCATCGTTTTAGCAAAGGGCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAg                                                                                                                       <  1:4262522/150‑1 (MQ=255)
                                 ttGAGGAGGATGAACGATAAGCACTGCCATCGTTTTAGCAAAGGGCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCg                                                                                                                     >  1:3274373/1‑150 (MQ=255)
                                 ttGAGGAGGATGAACGATAAGCACTGCCATCGTTTTAGCAAAGGGCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCg                                                                                                                     >  2:5493424/1‑150 (MQ=255)
                                   gaggagGATGAACGATAAGCACTGCCATCGTTTTAGCAAAGGGCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGt                                                                                                                   >  2:3812794/1‑150 (MQ=255)
                                    aggaggATGAACGATAAGCACTGCCATCGTTTTAGCAAAGGGCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTg                                                                                                                  >  1:791720/1‑150 (MQ=255)
                                     ggaggaTGAACGATAAGCACTGCCATCGTTTTAGCAAAGGGCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCACGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTgg                                                                                                                 >  2:592032/1‑150 (MQ=255)
                                      gaggaTGAACGATAAGCACTGCCATCGTTTTAGCAAAGGGCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCggg                                                                                                                                              <  2:5382930/120‑1 (MQ=255)
                                       aggaTGAACGATAAGCACTGCCATCGTTTTAGCAAAGGGCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCg                                                                                                               <  1:3622200/150‑1 (MQ=255)
                                        ggaTGAACGATAAGCACTGCCATCGTTTTAGCAAAGGGCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCgg                                                                                                              <  2:2436454/150‑1 (MQ=255)
                                         gaTGAACGATAAGCACTGCCATCGTTTTAGCAAAGGGCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCggg                                                                                                             >  1:5643440/1‑150 (MQ=255)
                                           tGAACGATAAGCACTGCCATCGTTTTAGCAAAGGGCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGcc                                                                                                           >  1:340313/1‑150 (MQ=255)
                                           tGAACGATAAGCACTGCCATCGTTTTAGCAAAGGGCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGcc                                                                                                           >  1:5365337/1‑150 (MQ=255)
                                              aCGATAAGCACTGCCATCGTTTTAGCAAAGGGCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGa                                                                                                        >  2:4949399/1‑150 (MQ=255)
                                                gATAAGCACTGCCATCGTTTTAGCAACGGGCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTTAAATGAACACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAg                                                                                                      <  1:474308/150‑1 (MQ=255)
                                                gATAAGCACTGCCATCGTTTTAGCAAAGGGCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAg                                                                                                      >  2:4742192/1‑150 (MQ=255)
                                                gATAAGCACTGCCATCGTTTTAGCAAAGGGCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAg                                                                                                      <  1:3244060/150‑1 (MQ=255)
                                                  tAAGCACTGCCATCGTTTTAGCAAAGGGCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGTGGTGGCGGGCCTGAAGGt                                                                                                    >  1:703751/1‑150 (MQ=255)
                                                   aaGCACGGCCATCGTTTTAGCAACGGGCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGCGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGtt                                                                                                   <  1:242690/150‑1 (MQ=255)
                                                       aCTGCCATCGTTTTAGCAAAGGGCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGt                                                                                               >  2:7301478/1‑150 (MQ=255)
                                                           ccATCGTTTTAGCAAAGGGCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGcc                                                                                           >  1:2547338/1‑150 (MQ=255)
                                                            cATCGTTTTAGCAACGGGCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCt                                                                                          <  1:4464733/150‑1 (MQ=255)
                                                             aTCGTTTTAGCAAAGGGCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTg                                                                                         <  2:1382299/150‑1 (MQ=255)
                                                             aTCGTTTTAGCAAAGGGCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTg                                                                                         >  2:1953015/1‑150 (MQ=255)
                                                              tCGTTTTAGCAAAGGGCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTgg                                                                                        >  2:437895/1‑150 (MQ=255)
                                                                   ttAGCAAAGGGCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGTCTgg                                                                                   <  2:3243451/150‑1 (MQ=255)
                                                                    tAGCAACGGGCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGTCTGGc                                                                                  <  1:7331402/150‑1 (MQ=255)
                                                                    tAGCAAAGGGCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGTCTGGc                                                                                  >  1:1521274/1‑150 (MQ=255)
                                                                       cAACGGGCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGTCTGGCCCt                                                                               <  1:1408022/150‑1 (MQ=255)
                                                                       cAAAGGGCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGTCTGGCCCt                                                                               >  1:4392634/1‑150 (MQ=255)
                                                                       cAAAGGGCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGTCTGGCCCt                                                                               <  2:4128125/150‑1 (MQ=255)
                                                                           gggCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGTCTGGCCCTTTcg                                                                           >  2:5168418/1‑150 (MQ=255)
                                                                             gCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGTCTGGCCCTTTcgcg                                                                         <  1:5732460/150‑1 (MQ=255)
                                                                              ccTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGTCTGGCCCTTTCGCgg                                                                        <  1:365923/150‑1 (MQ=255)
                                                                                tttACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGTCTGGCCCTTTCGCGGGt                                                                      >  1:5426008/1‑150 (MQ=255)
                                                                                   aCCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGTCTGGCCCTTTCGCGGGTGaa                                                                   >  1:864413/1‑150 (MQ=255)
                                                                                     cAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGTCTGGCCCTTTCGCGGGTGAAcc                                                                 >  2:7133632/1‑150 (MQ=255)
                                                                                     cAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGTCTGGCCCTTTCGCGGGTGAAcc                                                                 >  1:2840265/1‑150 (MQ=255)
                                                                                        agGTTGACCTTTGCCTGACGTGGCCACGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGTCTGGCCCTTTCGCGGGTGAACCCGc                                                              <  2:4466972/150‑1 (MQ=255)
                                                                                         ggTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGTCTGGCCCTTTCGCGGGTGAACCCGCt                                                             <  2:2247826/150‑1 (MQ=255)
                                                                                          gTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGTCTGGCCCTTTCGCGGGTGAACCCGCtt                                                            >  1:6171154/1‑150 (MQ=255)
                                                                                          gTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGTCTGGCCCTTTCGCGGGTGAACCCGCtt                                                            >  2:5658161/1‑150 (MQ=255)
                                                                                          gTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGTCTGGCCCTTTCGCGGGTGAACCCGCtt                                                            <  2:328313/150‑1 (MQ=255)
                                                                                          gTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTTAAAGGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGTCTGGCCCTTTCGCGGGTGAACCCGCtt                                                            <  2:1433173/150‑1 (MQ=255)
                                                                                               ccTTTGCCTGACGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGTCTGGCCCTTTCGCGGGTGAACCCGCTTCTATc                                                       >  1:7406589/1‑150 (MQ=255)
                                                                                                cTTTGCCTGACGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGTCTGGCCCTTTCGCGGGTGAACCCGCTTCTATca                                                      >  1:166434/1‑150 (MQ=255)
                                                                                                cTTTGCCTGACGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGTCTGGCCCTTTCGCGGGTGAACCCGCTTCTATca                                                      >  2:5303534/1‑150 (MQ=255)
                                                                                                 tttGCCTGACGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGTCTGGCCCTTTCGCGGGTGAACCCGCTTCTATcaa                                                     <  2:6648199/150‑1 (MQ=255)
                                                                                                 tttGCCTGACGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGTCTGGCCCTTTCGCGGGTGAACCCGCTTCTATcaa                                                     >  2:3404434/1‑150 (MQ=255)
                                                                                                 tttGCCTGACGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGTCTGGCCCTTTCGCGGGTGAACCCGCTTCTATcaa                                                     >  2:7049666/1‑150 (MQ=255)
                                                                                                  ttGCCTGACGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGTCTGGCCCTTTCGCGGGTGAACCCGCTTCTATcaaa                                                    <  2:3512567/150‑1 (MQ=255)
                                                                                                  ttGCCTGACGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGTCTGGCCCTTTCGCGGGTGAACCCGCTTCTATcaaa                                                    <  2:3666615/150‑1 (MQ=255)
                                                                                                   tgtcTGACGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGTCTGGCCCTTTCGCGGGTGAACCCGCTTCTATcaaac                                                   <  2:6870526/147‑1 (MQ=255)
                                                                                                        gACGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGTCTGGCCCTTTCGCGGGTGAACCCGCTTCTATCAAACAAAAg                                              >  2:5751201/1‑150 (MQ=255)
                                                                                                           gTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGTCTGGCCCTTTCGCGGGTGAACCCGCTTCTATCAAACAAAAGATa                                           >  2:6330809/1‑150 (MQ=255)
                                                                                                            ttGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGTCTGGCCCTTTCGCGGGTGAACCCGCTTCTATCAAACAAAAGATaa                                          <  1:5175170/150‑1 (MQ=255)
                                                                                                                cAAGAGGCGTAATTTTAAATGAACACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGTCTGGCCCTTTCGCGGGTGAACCCGCTTCTATCAAACAAAAGATAACCCa                                      <  1:5064249/150‑1 (MQ=255)
                                                                                                                        ggAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGTCTgg                                                                                   <  1:1035957/97‑1 (MQ=255)
                                                                                                                        ggAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGTCTgg                                                                                   >  2:1035957/1‑97 (MQ=255)
                                                                                                                          aaTTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGTCTGGCCCTTTCGCGGGTGAACCCGCTTCTATCAAACAAAAGATAACCCACTGttctttc                            <  1:114293/150‑1 (MQ=255)
                                                                                                                           aTTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGTCTGGCCCTTTCGCGGGTGAACCCGCTTCTATCAAACAAAAGATAACCCACTGttctttct                           <  2:6225457/150‑1 (MQ=255)
                                                                                                                             tttAAATGAAAACGGCCCGGATCGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGTCTGGCCCTTTCGCGGGTGAACCCGCTTCTATCAAACAAAAGATAACCCACTGTTCTTTCTAg                         <  2:5116202/150‑1 (MQ=255)
                                                                                                                                accctgtaaCGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGTCTGGCCCTTTCGCGGGTGAACCCGCTTCTATCAAACAAAAGATAACCCACTGTTCTTTCTAGCGa                      >  2:817821/8‑150 (MQ=255)
                                                                                                                                accctgtaaCGGCCCGGAACACGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGTCTGGCCCTTTCGCGGGTGAACCCGCTTCTATCAAACAAAAGATAACCCACTGTTCTTTCTAGCGa                      >  1:2083468/8‑150 (MQ=255)
                                                                                                                                aaaTGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGCGCCTGAAGGTTCGGTTGCCTGGTCTGGCCCTTTCGCGGGTGAACCCGCTTCTATCAAACAAAAGATAACCCACTGTTCTTTCTAGCGa                      <  2:3794259/150‑1 (MQ=255)
                                                                                                                                      aaaCGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGTCTGGCCCTTTCGCGGGTGAACCCGCTTCTATCAAACAAAAGATAACCCACTGTTCTTTCTAGCGaaaaaaa                <  2:4851225/150‑1 (MQ=255)
                                                                                                                                      aaaCGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGTCTGGCCCTTTCGCGGGTGAACCCGCTTCTATCAAACAAAAGATAACCCACTGTTCTTTCTAGCGaaaaaaa                <  2:1258305/150‑1 (MQ=255)
                                                                                                                                        aCGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGTCTGGCCCTTTCGCGGGTGAACCCGCTTCTATCAAACAAAAGATAACCCACTGTTCTTTCTAGCGAAAAAAACa              <  1:5285655/150‑1 (MQ=255)
                                                                                                                                         cGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGTCTGGCCCTTTCGCGGGTGAACCCGCTTCTATCAAACAAAAGATAACCCACTGTTCTTTCTAGCGAAAAAAACAg             <  2:4291608/150‑1 (MQ=255)
                                                                                                                                          ggCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGTCTGGCCCTTTCGCGGGTGAACCCGCTTCTATCAAACAAAAGATAACCCACTGTTCTTTCTAGCGAAAAAAACAgg            <  2:1296339/150‑1 (MQ=255)
                                                                                                                                          ggCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGTCTGGCCCTTTCGCGGGTGAACCCGCTTCTATCAAACAAAAGATAACCCACTGTTCTTTCTAGCGAAAAAAACAgg            >  1:7431450/1‑150 (MQ=255)
                                                                                                                                          ggCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGTCTGGCCCTTTCGCGGGTGAACCCGCTTCTATCAAACAAAAGATAACCCACTGTTCTTTCTAGCGAAAAAAACAgg            >  2:5609439/1‑150 (MQ=255)
                                                                                                                                          ggCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGTCTGGCCCTTTCGCGGGTGAACCCGCTTCTATCAAACAAAAGATAACCCACTGTTCTTTCTAGCGAAAAAAACAgg            >  1:420613/1‑150 (MQ=255)
                                                                                                                                          ggCCCGGAACACGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGTCTGGCCCTTTCGCGGGTGAACCCGCTTCTATCAAACAAAAGATAACCCACTGTTCTTTCTAGCGAAAAAAACAg             >  1:4035018/1‑149 (MQ=255)
                                                                                                                                               ggAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGTCTGGCCCTTTCGCGGGTGAACCCGCTTCTATCAAACAAAAGATAACCCACTGTTCTTTCTAGCGAAAAAAACAGgcgcgc       >  2:1772334/1‑150 (MQ=255)
                                                                                                                                                gAACGTGTTCCGGGCCGTTACAGGGGGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGTCTGGCCCTTTCGCGGGTGAACCCGCTTCTATCAAACAAAAGATAACCCACTGTTCTTTCTAGCGAACACAACAGGcgcgcg      >  1:6509937/1‑150 (MQ=255)
                                                                                                                                                 aaCGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGTCTGGCCCTTTCGCGGGTGAACCCGCTTCTATCAAACAAAAGATAACCCACTGTTCTTTCTAGCGAAAAAAACAGGCGCGCgg     <  1:1936238/150‑1 (MQ=255)
                                                                                                                                                  aCGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGTCTGGCCCTTTCGCGGGTGAACCCGCTTCTATCAAACAAAAGATAACCCACTGTTCTTTCTAGCGAAAAAAACAGGCGCGCggg    <  2:6075366/150‑1 (MQ=255)
                                                                                                                                                  aCGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGTCTGGCCCTTTCGCGGGTGAACCCGCTTCTATCAAACAAAAGATAACCCACTGTTCTTTCTAGCGAAAAAAACAGGCGCGCggg    <  2:2864766/150‑1 (MQ=255)
                                                                                                                                                  aCGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGTCTGGCCCTTTCGCGGGTGAACCCGCTTCTATCAAACAAAAGATAACCCACTGTTCTTTCTAGCGAAAAAAACAGGCGCGCggg    >  2:1683797/1‑150 (MQ=255)
                                                                                                                                                  aCGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGTCTGGCCCTTTCGCGGGTGAACCCGCTTCTATCAAACAAAAGATAACCCACTGTTCTTTCTAGCGAAAAAAACAGGCGCGCggg    <  1:6324926/150‑1 (MQ=255)
                                                                                                                                                    gtgtTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGTCTGGCCCTTTCGCGGGTGAACCCGCTTCTATCAAACAAAAGATAACCCACTGTTCTTTCTAGCGAAAAAAACAGGCGCGCGGGGt  >  1:2296421/1‑150 (MQ=255)
                                                                                                                                                    gtgtTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGTCTGGCCCTTTCGCGGGTGAACCCGCTTCTATCAAACAAAAGATAACCCACAGTTCTTTCTAGCGAAAAAAACAGGCGCGCGGGGt  >  1:6417887/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
ATGGGTTGAACGAACCTTGTTCGGAATGACGACTTGAGGAGGATGAACGATAAGCACTGCCATCGTTTTAGCAAAGGGCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGTCTGGCCCTTTCGCGGGTGAACCCGCTTCTATCAAACAAAAGATAACCCACTGTTCTTTCTAGCGAAAAAAACAGGCGCGCGGGGT  >  NC_002947/4078454‑4078751

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