Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,652,812 A→T 5.6% I163F (ATC→TTC)  ynfE → putative selenate reductase, periplasmic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,652,8120AT 5.6% 904.3 / 33.4 248I163F (ATC→TTC) ynfEputative selenate reductase, periplasmic
Reads supporting (aligned to +/- strand):  ref base A (130/103);  new base T (4/10);  total (135/113)
Fisher's exact test for biased strand distribution p-value = 5.60e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.61e-01

Only 100 of 249 total aligned reads displayed.

TTACCCAATGAAACGCGTGGGCAAACGCGGCGAAGGCAAATTCGAACGGATTAGCTGGGATGAAGCCCTGGATACTATCGCCAGTAGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCATGCCCTACACCTACGGCAGTAATG  >  CP009273/1652664‑1652959
                                                                                                                                                    |                                                                                                                                                   
ttACCCAATGAAACGCGTGGGCAAACGCGGCGAAGGCAAATTCGAACGGATTAGCTGGGATGAAGCCCTGGATACTATCGCCAGTAGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGa                                                                                                                                                     <  1:2593046/149‑1 (MQ=255)
  aCCCAATGAAACGCGTGGGCAAACGCGGTGAAGGCAAATTCGAACGGATTAGCTGGGATGAAGCCCTGGATACTATCGCCAGTAGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATc                                                                                                                                                   <  2:692993/149‑1 (MQ=255)
  aCCCAATGAAACGCGTGGGCAAACGCGGCGAAGGCAAATTCGAACGGATTAGCTGGGATGAAGCCCTGGATACTATCGCCAGTAGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATc                                                                                                                                                   <  1:248850/149‑1 (MQ=255)
    ccAATGAAACGCGTGGGCAAACGCGGCGAAGGCAAATTCGAACGGATTAGCTGGGATGAAGCCCTGGATACTATCGCCAGTAGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCgg                                                                                                                                                 >  2:2875213/1‑148 (MQ=255)
     cAATGAAACGCGTGGGCAAACGCGGCGAAGGCAAATTCGAACGGATTAGCTGGGATGAAGCCCTGGATACTATCGCCAGTAGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGtt                                                                                                                                                >  1:3672634/1‑149 (MQ=255)
       aTGAAACGCGTGGGCAAACGCGGCGAAGGCAAATTCGAACGGATTAGCTGGGATGAAGCCCTGGATACTATCGCCAGTAGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTgg                                                                                                                                              <  2:185133/149‑1 (MQ=255)
       aTGAAACGCGTGGGCAAACGCGGCGAAGGCAAATTCGAACGGATTAGCTGGGATGAAGCCCTGGATACTATCGCCAGTAGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTgg                                                                                                                                              <  2:185366/149‑1 (MQ=255)
        tGAAACGCGTGGGCAAACGCGGCGAAGGCAAATTCGAACGGATTAGCTGGGATGAAGCCCTGGATACTATCGCCAGTAGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTggc                                                                                                                                             >  1:3341286/1‑149 (MQ=255)
        tGAAACGCGTGGGCAAACGCGGCGAAGGCAAATTCGAACGGATTAGCTGGGATGAAGCCCTGGATACTATCGCCAGTAGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTggc                                                                                                                                             <  2:4360515/149‑1 (MQ=255)
         gAAACGCGTGGGCAAACGCGGCGAAGGCAAATTCGAACGGATTAGCTGGGATGAAGCCCTGGATACTATCGCCAGTAGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGTTCGTTggcg                                                                                                                                            >  1:3719992/1‑149 (MQ=255)
         gAAACGCGTGGGCAAACGCGGCGAAGGCAAATTCGAACGGATTAGCTGGGATGAAGCCCTGGATACTATCGCCAGTAGCCTGAAGAAAACCGTCGAACAATATGACAATGAGGCTGTATATATTCAGTACTCTTCGGGGCTCGTTggcg                                                                                                                                            >  2:3944528/1‑149 (MQ=255)
             cgcgTGGGCAAACGCGGCGAAGGCAAATTCGAACGGATTAGCTGGGATGAAGCCCTGGATACTATCGCCAGTAGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGTTCGTTGGCGGCaa                                                                                                                                        <  2:2566245/149‑1 (MQ=255)
             cgcgTGGGCAAACGCGGCGAAGGCAAATTCGAACGGATTAGCTGGGATGAAGCCCTGGATACTATCGCCAGTAGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCaa                                                                                                                                        <  2:3547163/149‑1 (MQ=255)
              gcgTGGGCAAACGCGGCGAAGGCAAATTCGAACGGATTAGCTGGGATGAAGCCCTGGATACTATCGCCAGTAGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAat                                                                                                                                       <  1:2426477/149‑1 (MQ=255)
                 tGGGCAAACGCGGCGAAGGCAAATTCGAACGGATTAGCTGGGATGAAGCCCTGGATACTATCGCCAGTAGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATg                                                                                                                                    >  2:4157293/1‑149 (MQ=255)
                 tGGGCAAACGCGGCGAAGGCAAATTCGAACGGATTAGCTGGGATGAAGCCCTGGATACTATCGCCAGTAGCCGGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGGATAGATTCAGTACTCGTCGGGGATCGTTGGCGGCAATATg                                                                                                                                    >  1:2960062/1‑149 (MQ=255)
                  gggCAAACGCGGCGAAGGCAAATTCGAACGGATTAGCTGGGATGAAGCCCTGGATACTATCGCCAGTAGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGa                                                                                                                                   <  2:3323266/149‑1 (MQ=255)
                       aaCGCGGCGAAGGCAAATTCGAACGGATTAGCTGGGATGAAGCCCTGGATACTATCGCCAGTAGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGc                                                                                                                              >  1:3126344/1‑149 (MQ=255)
                       aaCGCGGCGAAGGCAAATTCGAACGGATTAGCTGGGATGAAGCCCTGGATACTATCGCCAGTAGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGc                                                                                                                              >  1:77525/1‑149 (MQ=255)
                        aCGCGGCGAAGGCAAATTCGAACGGATTAGCTGGGATGAAGCCCTGGATACTATCGCCAGTAGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGct                                                                                                                             <  2:1243402/149‑1 (MQ=255)
                         cgcgGCGAAGGCAAATTCGAACGGATTAGCTGGGATGAAGCCCTGGATACTATCGCCAGTAGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGctc                                                                                                                            <  1:942547/149‑1 (MQ=255)
                          gcggcgAAGGCAAATTCGAACGGATTAGCTGGGATGAAGCCCTGGATACTATCGCCAGTAGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGctct                                                                                                                           <  2:26250/149‑1 (MQ=255)
                            ggcgAAGGCAAATTCGAACGGATTAGCTGGGATGAAGCCCTGGATACTATCGCCAGTAGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTc                                                                                                                         >  1:3196228/1‑149 (MQ=255)
                                 aGGCAAATTCGAACGGATTAGCTGGGATGAAGCCCTGGATACTATCGCCAGTAGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCACTATGACCCGCTCTTCGCCAt                                                                                                                    >  1:1690611/1‑149 (MQ=255)
                                 aGGCAAATTCGAACGGATTAGCTGGGATGAAGCCCTGGATACTATCGCCAGTAGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCAt                                                                                                                    <  2:681212/149‑1 (MQ=255)
                                 aGGCAAATTCGAACGGATTAGCTGGGATGAAGCCCTGGATACTATCGCCAGTAGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCAt                                                                                                                    <  2:3021673/149‑1 (MQ=255)
                                  ggCAAATTCGAACGGATTAGCTGGGATGAAGCCCTGGATACTATCGCCAGTAGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATc                                                                                                                   >  1:1770918/1‑149 (MQ=255)
                                            aaCGGATTAGCTGGGATGAAGCCCTGGATACTATCGCCAGTAGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTcggcg                                                                                                         >  1:1214830/1‑149 (MQ=255)
                                             aCGGATTAGCTGGGATGAAGCCCTGGATACTATCGCCAGTAGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTcggcgg                                                                                                        >  1:4555801/1‑149 (MQ=255)
                                             aCGGATTAGCTGGGATGAAGCCCTGGATACTATCGCCAGTAGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTcggcgg                                                                                                        >  2:3279795/1‑149 (MQ=255)
                                              cGGATTAGCTGGGATGAAGCCCTGGATACTATCGCCAGTAGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGt                                                                                                       <  1:60829/149‑1 (MQ=255)
                                                  ttAGCTGGGATGAAGCCCTGGATACTATCGCCAGTAGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCaaa                                                                                                   <  2:987570/149‑1 (MQ=255)
                                                       tGGGATGAAGCCCTGGATACTATCGCCAGTAGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCt                                                                                              <  2:2340562/149‑1 (MQ=255)
                                                       tGGGATGAAGCCCTGGATACTATCGCCAGTAGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCt                                                                                              >  2:1186668/1‑149 (MQ=255)
                                                         ggATGAAGCCCTGGATACTCTCGCCAGTAGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCtga                                                                                            <  2:335653/149‑1 (MQ=255)
                                                          gATGAAGCCCTGGATACTATCGCCAGTAGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCtgat                                                                                           >  1:648399/1‑149 (MQ=255)
                                                          gATGAAGCCCTGGATACTATCGCCAGTAGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCtgat                                                                                           >  1:648366/1‑149 (MQ=255)
                                                             gAAGCCCTGGATACTATCGCCAGTAGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGaa                                                                                        >  1:4141527/1‑149 (MQ=255)
                                                              aaGCCCTGGATACTATCGCCAGTAGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAAc                                                                                       >  2:3216094/1‑149 (MQ=255)
                                                              aaGCCCTGGATACTATCGCCAGTAGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAAc                                                                                       <  1:1498993/149‑1 (MQ=255)
                                                               aGCCCTGGATACTATCGCCAGTAGCCTGATGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACt                                                                                      <  2:2941936/149‑1 (MQ=255)
                                                                 cccTGGATACTATCGCCAGTAGCCTGAAGAAAACCGTCGAACAATATGGCCATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGc                                                                                    >  2:839012/1‑149 (MQ=255)
                                                                  ccTGGATACTATCGCCAGTAGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCt                                                                                   <  1:4052614/149‑1 (MQ=255)
                                                                  ccTGGATACTATCGCCAGTAGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCt                                                                                   >  2:1313116/1‑149 (MQ=255)
                                                                     ggATACTATCGCCAGTAGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTAcg                                                                                >  1:2081211/1‑149 (MQ=255)
                                                                      gATACTATCGCCAGTAGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTAcgg                                                                               >  2:869861/1‑149 (MQ=255)
                                                                        tACTATCGCCAGTAGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTAcggcg                                                                             >  1:2630659/1‑149 (MQ=255)
                                                                         aCTCTCGCCAGTAGCCTGAAGAAAACCGTCGAACAATCTGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTAcggcgg                                                                            <  2:4515920/149‑1 (MQ=255)
                                                                            aTCGCCAGTAGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTc                                                                         >  1:1103129/1‑149 (MQ=255)
                                                                                ccAGTAGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTc                                                                     <  1:4630994/149‑1 (MQ=255)
                                                                                  aGTAGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCaa                                                                   >  1:1433321/1‑149 (MQ=255)
                                                                                     aGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCa                                                                <  2:3986977/149‑1 (MQ=255)
                                                                                     aGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCa                                                                <  2:2717127/149‑1 (MQ=255)
                                                                                      gCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAg                                                               <  2:2669328/149‑1 (MQ=255)
                                                                                        cTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTa                                                             >  1:2009087/1‑149 (MQ=255)
                                                                                         tGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTAt                                                            >  1:4101752/1‑149 (MQ=255)
                                                                                         tGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTAt                                                            >  2:787016/1‑149 (MQ=255)
                                                                                           aagaaAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTAAATCAACCAGTATgg                                                          >  1:4496346/1‑149 (MQ=255)
                                                                                             gaaAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCt                                                        <  1:2695003/149‑1 (MQ=255)
                                                                                                aaCCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCt                                                     >  1:3943124/1‑149 (MQ=255)
                                                                                                 aCCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTa                                                    >  2:3417531/1‑149 (MQ=255)
                                                                                                 aCCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTa                                                    <  1:1529317/149‑1 (MQ=255)
                                                                                                  ccGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTAc                                                   <  1:3359830/149‑1 (MQ=255)
                                                                                                   cGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACa                                                  <  1:2999109/149‑1 (MQ=255)
                                                                                                    gTCGAACAATATGGCCATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAg                                                 >  2:2739859/1‑149 (MQ=255)
                                                                                                      cGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCa                                               >  1:2797335/1‑149 (MQ=255)
                                                                                                            atatGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTgcgc                                         >  2:3707984/1‑149 (MQ=255)
                                                                                                            atatGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTgcgc                                         >  2:463283/1‑149 (MQ=255)
                                                                                                             tatGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCa                                        >  2:1731003/1‑149 (MQ=255)
                                                                                                                ggCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAAt                                     >  2:764014/1‑149 (MQ=255)
                                                                                                                ggCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAAt                                     <  1:1004603/149‑1 (MQ=255)
                                                                                                                ggCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAAt                                     <  2:4095156/149‑1 (MQ=255)
                                                                                                                 gCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAAtt                                    >  2:1106756/1‑149 (MQ=255)
                                                                                                                 gCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAAtt                                    <  1:1196306/149‑1 (MQ=255)
                                                                                                                 gCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAAtt                                    <  2:3128162/149‑1 (MQ=255)
                                                                                                                   aaTGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTc                                  <  2:1678340/149‑1 (MQ=255)
                                                                                                                    aTGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTcc                                 <  1:4094648/149‑1 (MQ=255)
                                                                                                                     tGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCt                                >  1:3528833/1‑149 (MQ=255)
                                                                                                                            gTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCAt                         >  1:1649387/1‑149 (MQ=255)
                                                                                                                            gTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCAt                         >  1:3625742/1‑149 (MQ=255)
                                                                                                                              atatatTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCATGc                       >  1:1632014/1‑149 (MQ=255)
                                                                                                                              atatatTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCATGc                       >  2:2671113/1‑149 (MQ=255)
                                                                                                                              atatatTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCATGc                       <  2:1565472/149‑1 (MQ=255)
                                                                                                                               tatatTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCATGcc                      <  1:106685/149‑1 (MQ=255)
                                                                                                                               tatatTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCATGcc                      >  2:433282/1‑149 (MQ=255)
                                                                                                                                atatTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCATGccc                     >  1:1166371/1‑149 (MQ=255)
                                                                                                                                 tatTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCATGCCCt                    >  1:4651137/1‑149 (MQ=255)
                                                                                                                                  atTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCATGCCCTa                   >  1:2976423/1‑149 (MQ=255)
                                                                                                                                  atTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCATGCCCTa                   <  1:4214526/149‑1 (MQ=255)
                                                                                                                                   ttCAGTACTCTTCGGGGTTCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCATGCCCTac                  <  2:93559/149‑1 (MQ=255)
                                                                                                                                   ttCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCATGCCCTac                  >  2:1361295/1‑149 (MQ=255)
                                                                                                                                    tCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCATGCCCTaca                 <  2:1463255/149‑1 (MQ=255)
                                                                                                                                      aGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCATGCCCTACAcc               >  2:295220/1‑149 (MQ=255)
                                                                                                                                         aCTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCATGCCCTACACCTAc            >  2:499183/1‑149 (MQ=255)
                                                                                                                                          ctctTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCATGCCCTACACCTACg           >  2:3758790/1‑149 (MQ=255)
                                                                                                                                           tctTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCATGCCCTACACCTACgg          >  1:2856778/1‑149 (MQ=255)
                                                                                                                                           tctTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCATGCCCTACACCTACgg          >  2:3193145/1‑149 (MQ=255)
                                                                                                                                                 gggATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCATGCCCTACACCTACGGCAGTaa    >  2:1047181/1‑149 (MQ=255)
                                                                                                                                                  ggATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCATGCCCTACACCTACGGCAGTAat   <  2:2226454/149‑1 (MQ=255)
                                                                                                                                                  ggATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCATGCCCTACACCTACGGCAGTAat   <  2:4568363/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
TTACCCAATGAAACGCGTGGGCAAACGCGGCGAAGGCAAATTCGAACGGATTAGCTGGGATGAAGCCCTGGATACTATCGCCAGTAGCCTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCATGCCCTACACCTACGGCAGTAATG  >  CP009273/1652664‑1652959

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