Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 159,666 G→T 100% G359C (GGT→TGT)  hrpB → putative ATP‑dependent helicase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273159,6660GT97.0% 564.4 / 6.7 193G359C (GGT→TGT) hrpBputative ATP‑dependent helicase
Reads supporting (aligned to +/- strand):  ref base G (4/3);  new base T (97/89);  total (101/92)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 193 total aligned reads displayed.

AACGCGTTAGCCAGGCATCCATGACGCAGCGTGCCGGGCGCGCCGGGCGTCTGGAGCCGGGTATCAGCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCGGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGCCGCGAAACGTCTGTTACAAATGCTGGGGGCACTGGAGGGTGAACGGCTTAGTGCGCA  >  CP009273/159520‑159811
                                                                                                                                                  |                                                                                                                                                 
aaCGCGTTAGCCAGGCATCCATGACGCAGCGTGCCGGGCGCGCCGGGCGTCTGGAGCCGGGTATCAGCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGt                                                                                                                                                 <  2:3100712/149‑1 (MQ=255)
 aCGCGTTAGCCAGGCATCCATGACGCAGCGTGCCGGGCGCGCCGGGCGTCTGGAGCCGGGTATCAGCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGtt                                                                                                                                                <  1:724467/149‑1 (MQ=255)
        aGCCAGGCATCCATGACGCAGCGTGCCGGGCGCGCCGGGCGTCTGGAGCCGGGTATCAGCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGGTGAt                                                                                                                                         >  1:2215344/1‑149 (MQ=255)
        aGCCAGGCATCCATGACGCAGCGTGCCGGGCGCGCCGGGCGTCTGGAGCCGGGTATCAGCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGAt                                                                                                                                         <  2:1991684/149‑1 (MQ=255)
         gCCAGGCATCCATGACGCAGCGTGCCGGGCGCGCCGGGCGTCTGGAGCCGGGTATCAGCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATg                                                                                                                                        >  1:383345/1‑149 (MQ=255)
           cAGGCATCCATGACGCAGCGTGCCGGGCGCGCCGGGCGTCTGGAGCCGGGTATCAGCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGa                                                                                                                                      >  1:3551947/1‑149 (MQ=255)
           cAGGCATCCATGACGCAGCGTGCCGGGCGCGCCGGGCGTCTGGAGCCGGGTATCAGCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGa                                                                                                                                      >  2:1186230/1‑149 (MQ=255)
             ggCATCCATGACGCAGCGTGCCGGGCGCGCCGGGCGTCTGGAGCCGGGTATCAGCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAAt                                                                                                                                    <  2:1370524/149‑1 (MQ=255)
              gCATCCATGACGCAGCGTGCCGGGCGCGCCGGGCGTCTGGAGCCGGGTATCAGCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAAtt                                                                                                                                   >  2:3478267/1‑149 (MQ=255)
                aTCCATGACGCAGCGTGCCGGGCGCGCCGGGCGTCTGGAGCCGGGTATCAGCCTGCATTTAATCGCCACAGAACAAGCAGCTCGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCGGTGTGCGGATGGTATTag                                                                                                                                 >  2:2119114/1‑148 (MQ=255)
                 tCCATGACGCAGCGTGCCGGGCGCGCCGGGCGTCTGGAGCCGGGTATCAGCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACt                                                                                                                                >  2:1912645/1‑149 (MQ=255)
                 tCCATGACGCAGCGTGCCGGGCGCGCCGGGCGTCTGGAGCCGGGTATCAGCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCGGTTTGCTGATGGAATTACt                                                                                                                                >  1:868046/1‑149 (MQ=255)
                    aTGACGCAGCGTGCCGGGCGCGCCGGGCGTCTGGAGCCGGGTATCAGCCTGCATTTAATCGCCCAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCa                                                                                                                             >  2:3541726/1‑149 (MQ=255)
                     tGACGCAGCGTGCCGGGCGCGCCGGGCGTCTGGAGCCGGGTATCAGCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCaa                                                                                                                            >  2:3248390/1‑149 (MQ=255)
                      gACGCAGCGTGCCGGGCGCGCCGGGCGTCTGGAGCCGGGTATCAGCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAAt                                                                                                                           >  2:2315046/1‑149 (MQ=255)
                       aCGCAGCGTGCCGGGCGCGCCGGGCGTCTGGAGCCGGGTATCAGCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATg                                                                                                                          <  2:4085142/149‑1 (MQ=255)
                       aCGCAGCGTGCCGGGCGCGCCGGGCGTCTGGAGCCGGGTATCAGCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATg                                                                                                                          <  1:3455382/149‑1 (MQ=255)
                        cGCAGCGTGCCGGGCGCGCCGGGCGTCTGGAGCCGGGTATCAGCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATgg                                                                                                                         <  2:2261970/149‑1 (MQ=255)
                         gCAGCGTGCCGGGCGCGCCGGGCGTCTGGAGCCGGGTATCAGCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAAgggg                                                                                                                        >  1:2680047/1‑149 (MQ=255)
                          cAGCGTGCCGGGCGCGCCGGGCGTCTGGAGCCGGGTATCAGCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATgggg                                                                                                                       <  1:1203685/149‑1 (MQ=255)
                           aGCGTGCCGGGCGCGCCGGGCGTCTGGAGCCGGGTATCAGCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGa                                                                                                                      >  1:752754/1‑149 (MQ=255)
                            gCGTGCCGGGCGCGCCGGGCGTCTGGAGCCGGGTATCAGCCTGCATTTAATCGCCAAAGACCAAGCAGAAAGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGATGATGGAATTACTGCAATGGGGAt                                                                                                                     <  1:1705539/149‑1 (MQ=255)
                            gCGTGCCGGGCGCGCCGGGCGTCTGGAGCCGGGTATCAGCCTGCATTTAAGCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCGTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGAt                                                                                                                     <  1:4077455/149‑1 (MQ=255)
                             cGTGCCGGGCGCGCCGGGCGTCTGGAGCCGGGTATCAGCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATg                                                                                                                    <  1:2551369/149‑1 (MQ=255)
                              gTGCCGGGCGCGCCGGGCGTCTGGAGCCGGGTATCAGCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGc                                                                                                                   <  2:2750070/149‑1 (MQ=255)
                                   gggCGCGCCGGGCGTCTGGAGCCGGGTATCAGCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGa                                                                                                              <  1:3970968/149‑1 (MQ=255)
                                   gggCGCGCCGGGCGTCTGGAGCCGGGTATCAGCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGa                                                                                                              >  1:2979981/1‑149 (MQ=255)
                                    ggCGCGCCGGGCGTCTGGAGCCGGGTATCAGCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGCTGCAGCGAt                                                                                                             >  1:2441059/1‑149 (MQ=255)
                                    ggCGCGCCGGGCGTCTGGAGCCGGGTATCAGCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCCCCGCGGCGCAAAGTGAACCGGATATATTACAAAGCGATCATTCATGTTTGCTGATTGAATTACTGCAATGGGGATGCAGCGaa                                                                                                             >  1:369809/1‑148 (MQ=255)
                                      cgcgcCGGGCGTCTGGAGCCGGGTATCAGCCTGCATTTCATCGCCAACGAACACGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATcc                                                                                                           <  1:2071836/149‑1 (MQ=255)
                                         gcCGGGCGTCTGGAGCCGGGTATCAGCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGgc                                                                                                        >  2:112891/1‑149 (MQ=255)
                                            gggCGTCTGGAGCCGGGTATCAGCCGGCATTTAATCGCCAAAGACCAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCa                                                                                                     <  2:193547/149‑1 (MQ=255)
                                             ggCGTCTGGAGCCGGGTATCAGCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCCATGGGGATGCAGCGATCCGGCGCAg                                                                                                    >  1:2715783/1‑149 (MQ=255)
                                             ggCGTCTGGAGCCGGGTATCAGCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAg                                                                                                    >  2:2409029/1‑149 (MQ=255)
                                              gCGTCTGGAGCCGGGTATCAGCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACAGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGa                                                                                                   <  2:1626178/149‑1 (MQ=255)
                                                gTCTGGAGCCGGGTATCAGCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATg                                                                                                 <  2:46674/149‑1 (MQ=255)
                                                 tCTGGAGCCGGGTATCAGCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGa                                                                                                >  1:2199966/1‑149 (MQ=255)
                                                 tCTGGAGCCGGGTATCAGCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGa                                                                                                <  2:212690/149‑1 (MQ=255)
                                                  cTGGAGCCGGGTATCAGCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAg                                                                                               <  2:839730/149‑1 (MQ=255)
                                                   tGGAGCCGGGTATCAGCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAgc                                                                                              >  2:3223246/1‑149 (MQ=255)
                                                           ggTATCAGCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGa                                                                                      >  1:2974644/1‑149 (MQ=255)
                                                            gTATCAGCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGAt                                                                                     >  1:3530475/1‑149 (MQ=255)
                                                               tCAGCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTAAAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCACTGGGGCTGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCaa                                                                                  >  1:1925601/1‑149 (MQ=255)
                                                               tCAGCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTAAAAAGCGATCTTTCCTGTGTGCTGATGGAATTAATGCAATGGGGATGCAGCGATCAGGCGCAGATGAGCTGGCTGGATCaa                                                                                  >  1:3961343/1‑149 (MQ=255)
                                                                 aGCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAAcc                                                                                >  2:3260457/1‑149 (MQ=255)
                                                                  gCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTTCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCg                                                                               >  1:1264480/1‑149 (MQ=255)
                                                                  gCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCg                                                                               >  1:888478/1‑149 (MQ=255)
                                                                      gCATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCag                                                                           <  1:3902161/149‑1 (MQ=255)
                                                                         tttAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCagtag                                                                        >  2:666240/1‑149 (MQ=255)
                                                                            aaTCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGa                                                                     >  2:3361771/1‑149 (MQ=255)
                                                                              tCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAAt                                                                   <  2:337643/149‑1 (MQ=255)
                                                                               cGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATc                                                                  <  1:699058/149‑1 (MQ=255)
                                                                                      gAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGc                                                           <  1:694144/149‑1 (MQ=255)
                                                                                       aaCAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGcc                                                          >  2:2202682/1‑149 (MQ=255)
                                                                                        aCAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGCcg                                                         >  2:529391/1‑149 (MQ=255)
                                                                                        aCAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGCcg                                                         >  2:1344984/1‑149 (MQ=255)
                                                                                         cAAGCAGAACGCTCTGCCGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCGGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTAATGGCcgc                                                        <  1:2683951/149‑1 (MQ=255)
                                                                                              aGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGTAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGCCGCGAAAc                                                   <  1:2766100/149‑1 (MQ=255)
                                                                                               gAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGCCGCGAAACg                                                  >  1:596926/1‑149 (MQ=255)
                                                                                                aaCGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGACGCGAAACGt                                                 >  2:2616953/1‑149 (MQ=255)
                                                                                                 aCGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGCCGCGAAACGTc                                                <  1:2605128/149‑1 (MQ=255)
                                                                                                 aCGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGCCGCGAAACGTc                                                <  1:1513568/149‑1 (MQ=255)
                                                                                                 aCGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGCCGCGAAACGTc                                                <  2:4022016/149‑1 (MQ=255)
                                                                                                 aCGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGCCGCGAAACGTc                                                >  2:3054222/1‑149 (MQ=255)
                                                                                                   gcgcCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGACTCTACTGGCCGCGAAACGTCTg                                              >  2:2785126/1‑149 (MQ=255)
                                                                                                          ggcgCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGCCGCGAAACGTCTGTTACaaa                                       <  1:4053266/149‑1 (MQ=255)
                                                                                                          ggcgCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGCCGCGAAACGTCTGTTACaaa                                       <  2:1333717/149‑1 (MQ=255)
                                                                                                           gcgcAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGCCGCGAAACGTCTGTTACAAAt                                      <  2:216142/149‑1 (MQ=255)
                                                                                                                aaGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGCCGCGAAACGTCTGTTACAAATGCTgg                                 >  2:1438164/1‑149 (MQ=255)
                                                                                                                aaGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGCCGCGAAACGTCTGTTACAAATGCTgg                                 >  2:3797021/1‑149 (MQ=255)
                                                                                                                 aGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGCCGCGAAACGTCTGTTACAAATGCTggg                                >  2:570872/1‑149 (MQ=255)
                                                                                                                 aGTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGCCGCGAAACGTCTGTTACAAATGCTggg                                <  1:3654345/149‑1 (MQ=255)
                                                                                                                  gTGAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGCCGCGAAACGTCTGTTACAAATGCTgggg                               >  1:3611569/1‑149 (MQ=255)
                                                                                                                    gAACCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGCCGCGAAACGTCTGTTACAAATGCTGGGGGc                             <  1:1354995/149‑1 (MQ=255)
                                                                                                                     aaCCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGCCGCGAAACGTCTGTTACAAATGCTGGGGGCa                            >  1:1252654/1‑149 (MQ=255)
                                                                                                                     aaCCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGCCGCGAAACGTCTGTTACAAATGCTGGGGGCa                            >  2:3002249/1‑149 (MQ=255)
                                                                                                                      aCCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGCCGCGAAACGTCTGTTACAAATGCTGGGGGCAc                           <  2:122439/149‑1 (MQ=255)
                                                                                                                      aCCGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGCCGCGAAACGTCTGTTACAAATGCTGGGGGCAc                           <  1:2409029/149‑1 (MQ=255)
                                                                                                                        cGGAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGCCGCGAAACGTCTGTTACAAATGCTGGGGGCACTg                         >  1:3423384/1‑149 (MQ=255)
                                                                                                                         ggAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGCCGCGAAACGTCTGTTACAAATGCTGGGGGCACTgg                        >  1:4039929/1‑149 (MQ=255)
                                                                                                                         ggAGATCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGCCGCGAAACGTCTGTTACAAATGCTGGGGGCACTgg                        >  1:3121148/1‑149 (MQ=255)
                                                                                                                            gaTCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGCCGCGAAACGTCTGTTACAAATGCTGGGGGCACTGGAgg                     >  1:503958/1‑149 (MQ=255)
                                                                                                                            gaTCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGCCGCGAAACGTCTGTTACAAATGCTGGGGGCACTGGAgg                     <  2:585084/149‑1 (MQ=255)
                                                                                                                             aTCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGCCGCGAAACGTCTGTTACAAATGCTGGGGGCACTGGAggg                    >  1:2024617/1‑149 (MQ=255)
                                                                                                                              tCTTACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGCCGCGAAACGTCTGTTACAAATGCTGGGGGCACTGGAGGGt                   >  1:1227761/1‑149 (MQ=255)
                                                                                                                                ttACAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGCCGCGAAACGTCTGTTACAAATGCTGGGGGCACTGGAGGGTGa                 <  1:3248390/149‑1 (MQ=255)
                                                                                                                                 tacaccgCGATCTGTCCTGTTTGCTGAGGGAAGTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGCCGCGAACCGTCTGTTACAAATGCTGGGGGCACTGGAGGGGGaa                <  2:2990453/143‑1 (MQ=255)
                                                                                                                                  aCAAAGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGCCGCGAAACGTCTGTTACAAATGCTGGGGGCACTGGAGGGTGAAc               >  1:1483613/1‑149 (MQ=255)
                                                                                                                                    aaaGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGCCGCGAAACGTCTGTTACAAATGCTGGGGGCACTGGAGGGTGAACgg             <  1:2651261/149‑1 (MQ=255)
                                                                                                                                      aGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGCCGCGAAACGTCTGTTACAAATGCTGGGGGCACTGGAGGGTGAACGGCt           <  2:1004132/149‑1 (MQ=255)
                                                                                                                                      aGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGCCGCGAAACGTCTGTTACAAATGCTGGGGGCACTGGAGGGTGAACGGCt           <  1:1804147/149‑1 (MQ=255)
                                                                                                                                      aGCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGCCGAGAAACGTCTGTTACAAATGCTGTGGGCACTGGAGGGTGAACGGCt           <  1:1180348/149‑1 (MQ=255)
                                                                                                                                       gCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGCCGCGAAACGTCTGTTACAAATGCTGGGGGCACTGGAGGGTGAACGGCtt          >  2:3909405/1‑149 (MQ=255)
                                                                                                                                       gCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGCCGCGAAACGTCTGTTACAAATGCTGGGGGCACTGGAGGGTGAACGGCtt          <  1:185798/149‑1 (MQ=255)
                                                                                                                                       gCGATCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGCCGCGAAACGTCTGTTACAAATGCTGGGGGCACTGGAGGGTGAACGGCtt          >  1:1886038/1‑149 (MQ=255)
                                                                                                                                          aTCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGCCGCGAAACGTCTGTTACAAATGCTGGGGGCACTGGAGGGTGAACGGCTTAGt       <  1:3109500/149‑1 (MQ=255)
                                                                                                                                           tCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGCCGCGAAACGTCTGTTACAAATGCTGGGGGCACTGGAGGGTGAACGGCTTAGTg      >  2:1694281/1‑149 (MQ=255)
                                                                                                                                           tCTTTCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGCCGCGAAACGTCTGTTACAAATGCTGGGGGCACTGGAGGGTGAACGGCTTAGTg      >  2:443617/1‑149 (MQ=255)
                                                                                                                                             tttCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGCCGCGAAACGTCTGTTACAAATGCTGGGGGCACTGGAGGGTGAACGGCTTAGTgcg    <  2:1886577/149‑1 (MQ=255)
                                                                                                                                               tCCTGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGCCGCGAAACGTCTGTTACAAATGCTGGGGGCACTGGAGGGTGAACGGCTTAGTGCGCa  >  2:1841665/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
AACGCGTTAGCCAGGCATCCATGACGCAGCGTGCCGGGCGCGCCGGGCGTCTGGAGCCGGGTATCAGCCTGCATTTAATCGCCAAAGAACAAGCAGAACGCGCCGCGGCGCAAAGTGAACCGGAGATCTTACAAAGCGATCTTTCCGGTTTGCTGATGGAATTACTGCAATGGGGATGCAGCGATCCGGCGCAGATGAGCTGGCTGGATCAACCGCCAGTAGTGAATCTACTGGCCGCGAAACGTCTGTTACAAATGCTGGGGGCACTGGAGGGTGAACGGCTTAGTGCGCA  >  CP009273/159520‑159811

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