Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 726,087 G→T 39.4% R350L (CGT→CTT)  rhsC → Rhs family putative polymorphic toxin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273726,0870GT39.4% 352.5 / 308.0 293R350L (CGT→CTT) rhsCRhs family putative polymorphic toxin
Reads supporting (aligned to +/- strand):  ref base G (79/98);  new base T (59/57);  total (138/155)
Fisher's exact test for biased strand distribution p-value = 3.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 293 total aligned reads displayed.

TACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCT  >  CP009273/725942‑726234
                                                                                                                                                 |                                                                                                                                                   
 aCCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTcac                                                                                                                                                 >  1:4042971/1‑149 (MQ=255)
  cccGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTcaca                                                                                                                                                <  1:4113224/149‑1 (MQ=255)
  cccGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTcaca                                                                                                                                                <  2:3096884/149‑1 (MQ=255)
   ccGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCacac                                                                                                                                               >  2:4647835/1‑149 (MQ=255)
    cGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACg                                                                                                                                              >  1:622503/1‑149 (MQ=255)
      gagaATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACggg                                                                                                                                            >  2:4026101/1‑149 (MQ=255)
        gaATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGcc                                                                                                                                          >  1:4560446/1‑149 (MQ=255)
        gaATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGcc                                                                                                                                          <  2:3155873/149‑1 (MQ=255)
        gaATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGcc                                                                                                                                          <  2:3301473/149‑1 (MQ=255)
         aaTTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCg                                                                                                                                         >  2:2479001/1‑149 (MQ=255)
                cTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGaa                                                                                                                                  >  1:1882755/1‑149 (MQ=255)
                 tGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGaaa                                                                                                                                 >  2:3055456/1‑149 (MQ=255)
                   ccgcgccgcTGGTGCGCTATGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATc                                                                                                                               <  2:1423457/149‑1 (MQ=255)
                        ccgcTGGTGCGCTATGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTa                                                                                                                          <  1:3549065/149‑1 (MQ=255)
                          gcTGGTGCGCTATGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCTGTTAcc                                                                                                                        <  1:2222907/149‑1 (MQ=255)
                            tGGTGCGCTATGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCTGTTACCGt                                                                                                                      >  2:452224/1‑149 (MQ=255)
                             ggTGCGCTATGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGtt                                                                                                                     <  1:1816043/149‑1 (MQ=255)
                             ggTGCGCTATGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGtt                                                                                                                     <  2:3470691/149‑1 (MQ=255)
                                 cgcTATGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCTGTTACCGTTACGa                                                                                                                 <  2:1545828/149‑1 (MQ=255)
                                   cTATGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACa                                                                                                               >  2:2442642/1‑149 (MQ=255)
                                   cTATGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACa                                                                                                               >  1:2106634/1‑149 (MQ=255)
                                     aTGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGc                                                                                                             <  2:2405219/149‑1 (MQ=255)
                                     aTGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGc                                                                                                             <  2:4380796/149‑1 (MQ=255)
                                      tGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCg                                                                                                            >  1:3401294/1‑149 (MQ=255)
                                      tGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCg                                                                                                            <  2:2228966/149‑1 (MQ=255)
                                       ggCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGa                                                                                                           <  1:1313406/149‑1 (MQ=255)
                                           ggACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGAcggg                                                                                                       <  1:3143457/149‑1 (MQ=255)
                                            gACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGAcgggc                                                                                                      >  2:1685751/1‑149 (MQ=255)
                                             aCGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGAcgggcg                                                                                                     <  1:3876015/149‑1 (MQ=255)
                                                 cccGCGGCGAACTTGCGGCGGTGTATGACCGAAGAAATACACAGGTGCGCAGATTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGAACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTg                                                                                                 <  1:2693619/149‑1 (MQ=17)
                                                  ccGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGa                                                                                                <  2:1124836/149‑1 (MQ=255)
                                                   cgcgGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGAc                                                                                               <  2:2593920/149‑1 (MQ=255)
                                                         gAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACa                                                                                         >  2:1477211/1‑149 (MQ=255)
                                                           aCTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGc                                                                                       >  2:636226/1‑149 (MQ=255)
                                                                 ggcggTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAAcc                                                                                 >  2:190048/1‑149 (MQ=255)
                                                                 ggcggTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAAcc                                                                                 <  2:3858201/149‑1 (MQ=255)
                                                                      tgtATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCa                                                                            <  1:280216/149‑1 (MQ=255)
                                                                         aTGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGc                                                                         <  2:342736/149‑1 (MQ=255)
                                                                         aTGACCGCAGCAAGACACAGGGGCGCAGCTTTACTTCCGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACCGCTAAACCCGGCAGGc                                                                         <  2:4270691/149‑1 (MQ=14)
                                                                          tGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCt                                                                        <  1:894119/149‑1 (MQ=255)
                                                                          tGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCt                                                                        <  2:3718713/149‑1 (MQ=255)
                                                                           gACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCtt                                                                       <  1:4388329/149‑1 (MQ=255)
                                                                            aCCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACGGAACAGCTAAACCCGGCAGGCTTa                                                                      >  1:1670874/1‑149 (MQ=255)
                                                                            aCCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTa                                                                      >  1:1945149/1‑149 (MQ=255)
                                                                             ccGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTaa                                                                     <  1:156409/149‑1 (MQ=255)
                                                                             ccGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTaa                                                                     >  2:1889027/1‑149 (MQ=255)
                                                                               gcagcaATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGt                                                                   >  2:2432742/1‑149 (MQ=255)
                                                                               gcagcaATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGt                                                                   >  1:2970543/1‑149 (MQ=255)
                                                                               gcagcaATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGt                                                                   <  1:2097786/149‑1 (MQ=255)
                                                                                 agcaATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTa                                                                 <  2:3083276/149‑1 (MQ=255)
                                                                                   caATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTaca                                                               <  2:801554/149‑1 (MQ=255)
                                                                                       acacAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTa                                                           <  2:2554544/149‑1 (MQ=38)
                                                                                          caGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCa                                                        >  2:1875499/1‑149 (MQ=37)
                                                                                             gTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTa                                                     >  2:4350303/1‑149 (MQ=37)
                                                                                               gcgcAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATg                                                   >  2:1645167/1‑149 (MQ=37)
                                                                                                   aGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGaa                                               >  1:76014/1‑149 (MQ=37)
                                                                                                   aGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGaa                                               <  2:4717564/149‑1 (MQ=37)
                                                                                                   aGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGaa                                               >  2:2642505/1‑149 (MQ=37)
                                                                                                    gCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAg                                              >  1:4209366/1‑149 (MQ=38)
                                                                                                     cTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAgg                                             >  2:3986369/1‑149 (MQ=38)
                                                                                                      tttACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGa                                            >  1:3701250/1‑149 (MQ=255)
                                                                                                      tttACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTAGCAGTATGAGAAGGa                                            >  2:2737736/1‑149 (MQ=255)
                                                                                                      tttACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCTACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGa                                            >  1:2076704/1‑149 (MQ=255)
                                                                                                      tttAATTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGTAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGa                                            <  1:3858825/149‑1 (MQ=255)
                                                                                                       ttacttacGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGAc                                           >  2:4108290/1‑149 (MQ=39)
                                                                                                       ttacttacGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGAc                                           <  1:81895/149‑1 (MQ=39)
                                                                                                       ttacttacGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGAc                                           >  2:608043/1‑149 (MQ=39)
                                                                                                         acttacGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCg                                         >  1:3970088/1‑149 (MQ=39)
                                                                                                           ttacGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGca                                       <  1:2404718/149‑1 (MQ=255)
                                                                                                             acGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGcatc                                     <  2:1206454/149‑1 (MQ=255)
                                                                                                               gatgatAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGcatcac                                   >  1:4729823/1‑149 (MQ=39)
                                                                                                               gatgatAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGcatcac                                   <  1:735033/149‑1 (MQ=39)
                                                                                                                 tgatAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGcatcacca                                 >  2:3442954/1‑149 (MQ=255)
                                                                                                                  gatAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGcatcaccat                                <  2:2214027/149‑1 (MQ=39)
                                                                                                                    tAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGcatcaccatca                              >  2:2530374/1‑149 (MQ=39)
                                                                                                                     aaaTACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGcatcaccatcac                             <  1:430396/149‑1 (MQ=39)
                                                                                                                       aTACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCg                           >  2:1572853/1‑149 (MQ=39)
                                                                                                                        tACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGa                          <  1:2980445/149‑1 (MQ=39)
                                                                                                                         aCCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGAc                         >  1:3488629/1‑149 (MQ=255)
                                                                                                                            ggggCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGGGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGc                      >  2:3077964/1‑149 (MQ=255)
                                                                                                                             gggCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGcc                     >  1:1023995/1‑149 (MQ=255)
                                                                                                                              ggCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCt                    <  1:4635755/149‑1 (MQ=39)
                                                                                                                                 cGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGaa                 <  1:1745806/149‑1 (MQ=40)
                                                                                                                                  ggATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAAc                <  2:1803441/149‑1 (MQ=255)
                                                                                                                                   gATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAAcc               >  1:3216667/1‑149 (MQ=255)
                                                                                                                                   gATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAAcc               <  1:4702548/149‑1 (MQ=255)
                                                                                                                                   gATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAAcc               >  2:127220/1‑149 (MQ=255)
                                                                                                                                     tggtggCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAAccgc             <  2:2841793/149‑1 (MQ=255)
                                                                                                                                      ggtggCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAAccgcc            >  2:1923035/1‑149 (MQ=255)
                                                                                                                                        tggCGCACCTTCACACGGGCCGACCGTAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAACCGCCGt          <  1:2003746/149‑1 (MQ=255)
                                                                                                                                        tggCGCACCGTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAACCGCCGt          <  1:2494846/149‑1 (MQ=255)
                                                                                                                                        tggCGCACCGTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAACCGCCGt          <  1:4282373/149‑1 (MQ=255)
                                                                                                                                         ggCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAACCGCCGTg         >  2:723419/1‑149 (MQ=255)
                                                                                                                                          gcgcACCGTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAACCGCCGTGa        >  1:3781427/1‑149 (MQ=255)
                                                                                                                                            tcACCGTCACACGGGCCGACGTTAATTCTTTTACGGTTACGACAGCGACGGGCGGGTTACAGAACAGCTAAACCCGGCTGGTTTAAGTTACACGTATCAGTATGAGAAGGACCGAATCACCATCACCGACAGCCTGAACCGCCGTGAAg      <  1:1220925/148‑1 (MQ=14)
                                                                                                                                            gcACCGTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAg      >  1:1415908/1‑149 (MQ=255)
                                                                                                                                             cACCGTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGt     >  2:4294825/1‑149 (MQ=255)
                                                                                                                                             cACCGTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGt     >  2:215391/1‑149 (MQ=255)
                                                                                                                                               ccGTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTcc   >  2:3170602/1‑149 (MQ=255)
                                                                                                                                                cGTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCt  <  2:1980980/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
TACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGTCCT  >  CP009273/725942‑726234

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: