Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 726,087 G→T 93.7% 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,0870GT93.7% 491.2 / 23.1 178R350L (CGT→CTT) rhsCRhs family putative polymorphic toxin
Reads supporting (aligned to +/- strand):  ref base G (9/2);  new base T (87/79);  total (97/81)
Fisher's exact test for biased strand distribution p-value = 6.76e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01

Only 100 of 178 total aligned reads displayed.

TACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGT  >  CP009273/725942‑726231
                                                                                                                                                 |                                                                                                                                                
 aCCCGGAGAATTTACCTGCCGCTCCGCTGGTGCGCTATGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTca                                                                                                                                               <  2:2630747/148‑1 (MQ=255)
 aCCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTca                                                                                                                                               >  1:436799/1‑148 (MQ=255)
    cGGAGAATTTCCCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCacac                                                                                                                                            <  2:3311508/148‑1 (MQ=255)
    cGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCacac                                                                                                                                            >  1:2550628/1‑148 (MQ=255)
    cGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCacac                                                                                                                                            >  2:1809440/1‑148 (MQ=255)
    cGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCacac                                                                                                                                            >  2:1809462/1‑148 (MQ=255)
      gagaATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACgg                                                                                                                                          >  1:1237176/1‑148 (MQ=255)
      gagaATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACgg                                                                                                                                          >  2:1176674/1‑148 (MQ=255)
      gagaATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACgg                                                                                                                                          >  1:3006483/1‑148 (MQ=255)
       agaATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACggg                                                                                                                                         >  2:752451/1‑148 (MQ=255)
       agaATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACggg                                                                                                                                         >  2:2534881/1‑148 (MQ=255)
       agaATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACggg                                                                                                                                         >  2:2001071/1‑148 (MQ=255)
        gaATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGc                                                                                                                                        <  2:2478568/148‑1 (MQ=255)
         aaTTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGcc                                                                                                                                       >  1:806895/1‑148 (MQ=255)
                   ccgcgccgcTGGTGCGCTATGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAAt                                                                                                                             <  2:3135939/148‑1 (MQ=255)
                       gccgcTGGTGCGCTATGTCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACCCGGGCCGACCGGAAATCTGt                                                                                                                         <  2:615046/148‑1 (MQ=255)
                        ccgcTGGTGCGCTATGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGtt                                                                                                                        <  1:2087743/148‑1 (MQ=255)
                              gTGCGCTATGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGtt                                                                                                                  <  2:2855421/148‑1 (MQ=255)
                                              cGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTCCGACAGCGAcgggcg                                                                                                  >  2:731171/1‑148 (MQ=255)
                                              cGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGAcgggcg                                                                                                  >  2:2897024/1‑148 (MQ=255)
                                                ccccGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGAcgggcggg                                                                                                <  1:2792670/148‑1 (MQ=255)
                                                  ccGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTg                                                                                              >  1:1282871/1‑148 (MQ=255)
                                                        cgAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGaa                                                                                        >  1:1044807/1‑148 (MQ=255)
                                                        cgAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGaa                                                                                        <  1:1769919/148‑1 (MQ=255)
                                                          aaCTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTCCCGTTACGACAGCGACGGGCGGGTGACAGAACa                                                                                      >  2:515009/1‑148 (MQ=255)
                                                           aCTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAg                                                                                     <  1:2787135/148‑1 (MQ=255)
                                                               gcggcgGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTaa                                                                                 >  1:2233007/1‑148 (MQ=255)
                                                               gcggcgGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTaa                                                                                 >  1:2232981/1‑148 (MQ=255)
                                                                  gcggTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAAcc                                                                              <  2:2940726/148‑1 (MQ=255)
                                                                       gtATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCa                                                                         <  1:876136/148‑1 (MQ=255)
                                                                       gtATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCa                                                                         <  1:875883/148‑1 (MQ=255)
                                                                       gtATGACCGCAGCAATACACAGGTGCGCAGCTTTAATTACGTTGATAAATACCGGGGCCGGTTGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCa                                                                         <  1:1627091/148‑1 (MQ=255)
                                                                        tATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAg                                                                        >  2:932171/1‑148 (MQ=255)
                                                                        tATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAg                                                                        >  1:2279565/1‑148 (MQ=255)
                                                                        tATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAg                                                                        >  1:3175670/1‑148 (MQ=255)
                                                                         aTGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAgg                                                                       <  2:1160696/148‑1 (MQ=255)
                                                                          tGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGc                                                                      >  2:1080000/1‑148 (MQ=255)
                                                                            aCCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATGTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCtt                                                                    >  1:1819955/1‑148 (MQ=255)
                                                                            aCCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACAGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCtt                                                                    >  1:372994/1‑148 (MQ=255)
                                                                             ccGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTa                                                                   >  1:1184518/1‑148 (MQ=255)
                                                                               gcagcaATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAg                                                                 <  1:2534881/148‑1 (MQ=255)
                                                                               gcagcaATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAg                                                                 <  1:752451/148‑1 (MQ=255)
                                                                               gcagcaATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAg                                                                 <  1:50740/148‑1 (MQ=255)
                                                                               gcagcaATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAg                                                                 <  1:1147986/148‑1 (MQ=255)
                                                                                  gcaATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTa                                                              >  2:855292/1‑148 (MQ=255)
                                                                                  gcaATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTa                                                              >  2:2930995/1‑148 (MQ=255)
                                                                                   caATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTac                                                             <  1:664545/148‑1 (MQ=255)
                                                                                   caATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTac                                                             >  1:3204192/1‑148 (MQ=255)
                                                                                    aaTACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTaca                                                            <  2:3200805/148‑1 (MQ=255)
                                                                                       acacAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGt                                                         <  1:3330761/148‑1 (MQ=38)
                                                                                        cacaGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTa                                                        <  2:442302/148‑1 (MQ=38)
                                                                                        cacaGGTGCGCAGCTTTACTTACGATGATAAAGACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTa                                                        <  2:768846/148‑1 (MQ=255)
                                                                                            ggTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAg                                                    >  1:14382/1‑148 (MQ=37)
                                                                                            ggTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAg                                                    >  2:294387/1‑148 (MQ=37)
                                                                                             gTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGt                                                   <  1:1809440/148‑1 (MQ=37)
                                                                                             gTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGAATCAGt                                                   >  2:1735994/1‑148 (MQ=255)
                                                                                                cgcAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACAGGGGCCGACCGGAATTATGTTACCGTTACGACATCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATg                                                <  2:3266997/148‑1 (MQ=11)
                                                                                                cgcAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATg                                                >  1:2068614/1‑148 (MQ=37)
                                                                                                cgcAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATg                                                <  1:3405671/148‑1 (MQ=37)
                                                                                                cgcAGCTTTACTTACGATGAGAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATg                                                <  2:1038006/148‑1 (MQ=255)
                                                                                                  cAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATgag                                              >  2:582590/1‑148 (MQ=37)
                                                                                                    gCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGaa                                            <  1:3060271/148‑1 (MQ=37)
                                                                                                      tttACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAgg                                          <  1:1430250/148‑1 (MQ=38)
                                                                                                      tttACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAgg                                          <  1:2925272/148‑1 (MQ=38)
                                                                                                      tttACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAgg                                          <  1:678523/148‑1 (MQ=38)
                                                                                                       ttacttacGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGa                                         <  2:1529876/148‑1 (MQ=255)
                                                                                                       ttacttacGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGa                                         >  2:345878/1‑148 (MQ=255)
                                                                                                       ttacttacGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGa                                         >  2:2958774/1‑148 (MQ=39)
                                                                                                       ttacttacGATGATAAATAACGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTCAGTTACACGGATCAGTATGAGAAGGa                                         >  2:2371276/1‑148 (MQ=255)
                                                                                                        tacttacGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGAc                                        >  2:11637/1‑148 (MQ=255)
                                                                                                           ttacGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGc                                     >  2:234600/1‑148 (MQ=255)
                                                                                                            tacGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGca                                    >  2:1662043/1‑148 (MQ=255)
                                                                                                             acGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGcat                                   >  1:373631/1‑148 (MQ=255)
                                                                                                               gatgatAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACTACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGcatca                                 >  1:1160883/1‑148 (MQ=255)
                                                                                                                 tgatAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGcatcacc                               >  2:2105672/1‑148 (MQ=255)
                                                                                                                 tgatAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGcatcacc                               >  2:561331/1‑148 (MQ=255)
                                                                                                                  gatAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGcatcacca                              <  2:643857/148‑1 (MQ=255)
                                                                                                                   atAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGcatcaccat                             >  1:1021542/1‑148 (MQ=255)
                                                                                                                   atAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGcatcaccat                             >  1:680442/1‑148 (MQ=255)
                                                                                                                     aaaTACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGcatcaccatca                           >  1:3189147/1‑148 (MQ=255)
                                                                                                                     aaaTACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGcatcaccatca                           >  1:3237602/1‑148 (MQ=255)
                                                                                                                      aaTACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGcatcaccatcac                          >  2:1905581/1‑148 (MQ=39)
                                                                                                                         aCCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGa                       >  2:3520471/1‑148 (MQ=255)
                                                                                                                         aCCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGa                       <  2:3036919/148‑1 (MQ=255)
                                                                                                                          ccGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGAc                      <  1:2600701/148‑1 (MQ=255)
                                                                                                                             gggCCGGATGGTGGCGCACCCTCACACGGGCCGACCGGAAACCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGc                   >  2:1184393/1‑148 (MQ=255)
                                                                                                                              ggCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGcc                  <  2:82374/148‑1 (MQ=255)
                                                                                                                              ggCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGcc                  >  1:588127/1‑148 (MQ=39)
                                                                                                                              ggCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGcc                  >  1:3482703/1‑148 (MQ=39)
                                                                                                                                   gATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAAc             >  2:2533598/1‑148 (MQ=255)
                                                                                                                                   gATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAAc             >  2:216155/1‑148 (MQ=255)
                                                                                                                                     tggtggCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAAccg           >  2:1285721/1‑148 (MQ=255)
                                                                                                                                      ggtggCGCACCGTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAAccgc          <  1:920783/148‑1 (MQ=255)
                                                                                                                                             cACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAg   >  2:2553353/1‑148 (MQ=255)
                                                                                                                                             cACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAg   >  2:2795180/1‑148 (MQ=255)
                                                                                                                                             cACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAg   <  2:373631/148‑1 (MQ=255)
                                                                                                                                             cACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAg   <  2:373630/148‑1 (MQ=255)
                                                                                                                                             cACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAg   >  2:120977/1‑148 (MQ=255)
                                                                                                                                              aCCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGt  <  1:294288/148‑1 (MQ=255)
                                                                                                                                              aCCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGt  <  1:450649/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
TACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAGT  >  CP009273/725942‑726231

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: