Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 520,873 A→T 13.5% A719A (GCA→GCT rhsD → Rhs family putative polymorphic toxin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273520,8730AT13.5% 766.7 / 96.4 260A719A (GCA→GCTrhsDRhs family putative polymorphic toxin
Reads supporting (aligned to +/- strand):  ref base A (128/97);  new base T (22/13);  total (150/110)
Fisher's exact test for biased strand distribution p-value = 5.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.64e-01

Only 100 of 274 total aligned reads displayed.

CAGCAGGGCGGCTTTGACGGGCGGACGCAACGTTATCATTATGACCTGACCGGAAAACTCACACAGAGTGAGGATGAGGGACTTGTCATCCTCTGGTACTACGATGAATCGGACCGTATCACTCACCGCACGGTGAACGGCGAACCGGCAGAGCAGTGGCAGTATGATGGCCACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAACT  >  CP009273/520724‑521022
                                                                                                                                                     |                                                                                                                                                     
cagcagGGCGGCTTTGACGGGCGGACGCAACGTTATCATTATGACCTGACCGGAAAACTCACACAGAGTGAGGATGAGGGACTTGTCATCCTCTGGTACTACGATGAATCGGACCGTATCACTCACCGCACGGTGAACGGCGAACCGGCa                                                                                                                                                       <  2:2340550/150‑1 (MQ=255)
 agcagGGCGGCTTTGACGGGCGGACGCAACGTTATCATTATGACCTGACCGGAAAACTCACACAGAGTGAGGATGAGGGACTTGTCATCCTCTGGTACTACGATGAATCGGACCGTATCACTCACCGCACGGTGAACGGCGAACCGGCag                                                                                                                                                      >  2:209854/1‑150 (MQ=255)
  gcagGGCGGCTTTGACGGGCGGACGCAACGTTATCATTATGACCTGACCGGAAAACTCACACAGAGTGAGGATGAGGGACTTGTCATCCTCTGGTACTACGATGAATCGGACCGTATCACTCACCGCACGGTGAACGGCGAACCGGCaga                                                                                                                                                     <  1:60723/150‑1 (MQ=255)
    agGGCGGCTTTGACGGGCGGACGCAACGTTATCATTATGACCTGACCGGAAAACTCACACAGAGTGAGGATGAGGGACTTGTCATCCTCTGGTACTACGATGAATCGGACCGTATCACTCACCGCACGGTGAACGGCGAACCGGCAGAGc                                                                                                                                                   <  1:1042173/150‑1 (MQ=255)
     gggCGGCTTTGACGGGCGGACGCAACGTTATCATTATGACCTGACCGGAAAACTCACACAGAGTGAGGATGAGGGACTTGTCATCCTCTGGTACTACGATGAATCGGACCGTATCACTCACCGCACGGTGAACGGCGAACCGGCAGAGCa                                                                                                                                                  >  2:2982608/1‑150 (MQ=255)
       gcggcTTTGACGGGCGGACGCAACGTTATCATTATGACCTGACCGGAAAACTCACACAGAGTGAGGATGAGGGACTTGTCATCCTCTGGTACTACGATGAATCGGACCGTATCACTCACCGCACGGTGATCGGCGAACCGGCAGAGCAGt                                                                                                                                                <  1:2953601/150‑1 (MQ=255)
         ggcTTTGACGGGCGGACGCAACGTTATCATTATGACCTGACCGGAAAACTCACACAGAGTGAGGATGAGGGACTTGTCATCCTCTGGTACTACGATGAATCGGACCGTATCACTCACCGCACGGTGAACGGCGAACCGGCAGAGCAGTgg                                                                                                                                              >  1:2531111/1‑150 (MQ=255)
               gACGGGCGGACGCAACGTTATCATTATGACCTGACCGGAAAACTCACACAGAGTGAGGATGAGGGACTTGTCATCCTCTGGTACTACGATGAATCGGACCGTATCACTCACCGCACGGTGAACGGCGAACCGGCAGAGCAGTGGCAGTat                                                                                                                                        <  1:843773/150‑1 (MQ=255)
                 cGGGCGGACGCAACGTTATCATTATGACCTGACCGGAAAACTCACACAGAGTGAGGATGAGGGACTTGTCATCCTCTGGTACTACGATGAATCGGACCGTATCACTCACCGCACGGTGAACGGCGAACCGGCAGAGCAGTGGCAGTatga                                                                                                                                      <  1:3217159/150‑1 (MQ=255)
                 cGGGCGGACGCAACGTTATCATTATGACCTGACCGGAAAACTCACACAGAGTGAGGATGAGGGACTTGTCATCCTCTGGTACTACGATGAATCGGACCGTATCACTCACCGCACGGTGAACGGCGAACCGGCAGAGCAGTGGCAGTatga                                                                                                                                      <  1:2355158/150‑1 (MQ=255)
                         cGCAACGTTATCATTATGACCTGACCGGAAAACTCACACAGAGTGAGGATGAGGGACTTGTCATCCTCTGGTACTACGATGAATCGGACCGTATCACTCACCGCACGGTGAACGGCGAACCGGCAGAGCAGTGGCAGTATGATGGCCACg                                                                                                                              >  2:3040720/1‑150 (MQ=255)
                              cGTTATCATTATGACCTGACCGGAAAACTCACACAGAGTGAGGATGAGGGACTTGTCATCCTCTGGTACTACGATGAATCGGACCGTATCACTCACCGCACGGTGAACGGCGAACCGGCTGAGCAGTGGCAGTATGATGGCCACggctgg                                                                                                                         <  2:628695/150‑1 (MQ=255)
                                 tATCATTATGACCTGACCGGAAAACTCACACAGAGTGAGGATGAGGGACTTGTCATCCTCTGGTACTACGATGAATCGGACCGTATCACTCACCGCACGGTGAACGGCGAACCGGCAGAGCAGTGGCAGTATGATGGCCACGgctggctg                                                                                                                      >  2:1529772/1‑150 (MQ=255)
                                 tATCATTATGACCTGACCGGAAAACTCACACAGAGTGAGGATGAGGGACTTGTCATCCTCTGGTACTACGATGAATCGGACCGTATCACTCACCGCACGGTGAACGGCGAACCGGCAGAGCAGTGGCAGTATGATGGCCACGgctggctg                                                                                                                      >  2:498569/1‑150 (MQ=255)
                                   tCATTATGACCTGACCGGAAAACTCACACAGAGTGAGGATGAGGGACTTGTCATCCTCTGGTACTACGATGAATCGGACCGTATCACTCACCGCACGGTGAACGGCGAACCGGCAGAGCAGTGGCAGTATGATGGCCACGGCTGGCTgac                                                                                                                    >  2:2917147/1‑150 (MQ=255)
                                      ttATGACCTGACCGGAAAACTCACACAGAGTGAGGATGAGGGACTTGTCATCCTCTGGTACTACGATGAATCGGACCGTATCACTCACCGCACGGTGAACGGCGAACCGGCAGAGCAGTGGCAGTATGATGGCCACGGCTGGCTgacaga                                                                                                                 <  2:2352737/150‑1 (MQ=255)
                                      ttATGACCTGACCGGAAAACTCACACAGAGTGAGGATGAGGGACTTGTCATCCTCTGGTACTACGATGAATCGGACCGTATCACTCACCGCACGGTGAACGGCGAACCGGCAGAGCAGTGGCAGTATGATGGCCACGGCTGGCTgacaga                                                                                                                 <  1:100980/150‑1 (MQ=255)
                                       tATGACCTGACCGGAAAACTCACACAGAGTGAGGATGAGGGACTTGTCATCCTCTGGTACTACGATGAATCGGACCGTATCACTCACCGCACGGTGAACGGCGAACCGGCAGAGCAGTGGCAGTATGATGGCCACGGCTGGCTgacagac                                                                                                                >  1:3595320/1‑150 (MQ=255)
                                        aTGACCTGACCGGAAAACTCACACAGAGTGAGGATGAGGGACTTGTCATCCTCTGGTACTACGATGAATCGGACCGTATCACTCACCGCACGGTGAACGGCGAACCGGCAGAGCAGTGGCAGTATGATGGCCACGGCTGGCTgacagaca                                                                                                               <  1:2567163/150‑1 (MQ=255)
                                           acctgaccGGAAAACTCACACAGAGTGAGGATGAGGGACTTGTCATCCTCTGGTACTACGATGAATCGGACCGTATCACTCACCGCACGGTGAACGGCGAACCGGCAGAGCAGTGGCAGTATGATGGCCACGGCTGGCTGACAGACATCa                                                                                                            >  2:3288464/1‑150 (MQ=255)
                                            cctgaccGGAAAACTCACACAGAGTGAGGATGAGGGACTTGTCATCCTCTGGTACTACGATGAATCGGACCGTATCACTCACCGCACGGTGAACGGCGAACCGGCAGAGCAGTGGCAGTATGATGGCCACGGCTGGCTGACAGACATCAg                                                                                                           >  2:3319444/1‑150 (MQ=255)
                                                accGGAAAACTCACACAGAGTGAGGATGAGGGACTTGTCATCCTCTGGTACTACGATGAATCGGACCGTATCACTCACCGCACGGTGAACGGCGAACCGGCAGAGCAGTGGCAGTATGATGGCCACGGCTGGCTGACAGACATCAGCCAc                                                                                                       >  2:3212422/1‑150 (MQ=255)
                                                accGGAAAACTCACACAGAGTGAGGATGAGGGACTTGTCATCCTCTGGTACTACGATGAATCGGACCGTATCACTCACCGCACGGTGAACGGCGAACCGGCAGAGCAGTGGCAGTATGATGGCCACGGCTGGCTGACAGACATCAGCCAc                                                                                                       >  1:1401011/1‑150 (MQ=255)
                                                  cGGAAAACTCACACAGAGTGAGGATGAGGGACTTGTCATCCTCTGGTACTACGATGAATCGGACCGTATCACTCACCGCACGGTGAACGGCGAACCGGCAGAGCAGTGGCAGTATGATGGCCACGGCTGGCTGACAGACATCAGCCACCt                                                                                                     >  1:86693/1‑150 (MQ=255)
                                                    gAAAACTCACACAGAGTGAGGATGAGGGACTTGTCATCCTCTGGTACTACGATGAATCGGACCGTATCACTCACCGCACGGTGAACGGCGAACCGGCAGAGCAGTGGCAGTATGATGGCCACGGCTGGCTGACAGACATCAGCCACCTGa                                                                                                   <  1:102342/150‑1 (MQ=255)
                                                    gAAAACTCACACAGAGTGAGGATGAGGGACTTGTCATCCTCTGGTACTACGATGAATCGGACCGTATCACTCACCGCACGGTGAACGGCGAACCGGCAGAGCAGTGGCAGTATGATGGCCACGGCTGGCTGACAGACATCAGCCACCTGa                                                                                                   <  1:2903480/150‑1 (MQ=255)
                                                    gAAAACTCACACAGAGTGAGGATGAGGGACTTGTCATCCTCTGGTACTACGATGAATCGGACCGTATCACTCACCGCACGGTGAACGGCGAACCGGCAGAGCAGTGGCAGTATGATGGCCACGGCTGGCTGACAGACATCAGCCACCTGa                                                                                                   >  1:2790781/1‑150 (MQ=255)
                                                      aaaCTCACACAGAGTGAGGATGAGGGACTTGTCATCCTCTGGTACTACGATGAATCGGACCGTATCACTCACCGCACGGTGAACGGCGAACCGGCAGAGCAGTGGCAGTATGATGGCCACGGCTGGCTGACAGACATCAGCCACCTGAGc                                                                                                 >  2:1066606/1‑150 (MQ=255)
                                                      aaaCTCACACAGAGTGAGGATGAGGGACTTGTCATCCTCTGGTACTACGATGAATCGGACCGTATCACTCACCGCACGGTGAACGGCGAACCGGCAGAGCAGTGGCAGTATGATGGCCACGGCTGGCTGACAGACATCAGCCACCTGAGc                                                                                                 >  2:4096414/1‑150 (MQ=255)
                                                      aaaCTCACACAGAGTGAGGATGAGGGACTTGTCATCCTCTGGTACTACGATGAATCGGACCGTATCACTCACCGCACGGTGAACGGCGAACCGGCAGAGCAGTGGCAGTATGATGGCCACGGCTGGCTGACAGACATCAGCCACCTGAGc                                                                                                 >  1:1400326/1‑150 (MQ=255)
                                                        aCTCACACAGAGTGAGGATGAGGGACTTGTCATCCTCTGGTACTACGATGAATCGGACCGTATCACTCACCGCACGGTGAACGGCGAACCGGCAGAGCAGTGGCAGTATGATGGCCACGGCTGGCTGACAGACATCAGCCACCTGAGCGa                                                                                               <  1:4558155/150‑1 (MQ=255)
                                                              acaGAGTGAGGATGAGGGACTTGTCATCCTCTGGTACTACGATGAATCGGACCGTATCACTCACCGCACGGTGAACGGCGAACCGGCAGAGCAGTGGCAGTATGATGGCCACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCa                                                                                         >  1:1220580/1‑150 (MQ=255)
                                                                agagTGAGGATGAGGGACTTGTCATCCTCTGGTACTACGATGAATCGGACCGTATCACTCACCGCACGGTGAACGGCGAACCGGCAGAGCAGTGGCAGTATGATGGCCACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCAcc                                                                                       >  2:1432085/1‑150 (MQ=255)
                                                                agagTGAGGATGAGGGACTTGTCATCCTCTGGTACTACGATGAATCGGACCGTATCACTCACCGCACGGTGAACGGCGAACCGGCAGAGCAGTGGCAGGATGATGGCCACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCAcc                                                                                       >  1:805569/1‑150 (MQ=255)
                                                                 gagTGAGGATGAGGGACTTGTCATCCTCTGGTACTACGATGAATCGGACCGTATCACTCACCGCACGGTGAACGGCGAACCGGCAGAGCAGTGGCAGTATGATGGCCACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCg                                                                                      >  2:1306933/1‑150 (MQ=255)
                                                                     gAGGATGAGGGACTTGTCATCCTCTGGTACTACGATGAATCGGACCGTATCACTCACCGCACGGTGAACGGCGAACCGGCAGAGCAGTGGCAGTATGATGGCCACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGtt                                                                                  >  2:81269/1‑150 (MQ=255)
                                                                     gAGGATGAGGGACTTGTCATCCTCTGGTACTACGATGAATCGGACCGTATCACTCACCGCACGGTGAACGGCGAACCGGCAGAGCAGTGGCAGTATGATGGCCACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGtt                                                                                  <  1:1804066/150‑1 (MQ=255)
                                                                     gAGGATGAGGGACTTGTCATCCTCTGGTACTACGATGAATCGGACCGTATCACTCACCGCACGGTGAACGGCGAACCGGCAGAGCAGTGGCAGTATGATGGCCACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGtt                                                                                  <  1:1708431/150‑1 (MQ=255)
                                                                     gAGGATGAGGGACTTGTCATCCTCTGGTACTACGATGAATCGGACCGTATCACTCACCGCACGGTGAACGGCGAACCGGCAGAGCAGTGGCAGTATGATGGCCACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGtt                                                                                  <  1:1804006/150‑1 (MQ=255)
                                                                      aGGATGAGGGACTTGTCATCCTCTGGTACTACGATGAATCGGACCGTATCACTCACCGCACGGTGAACGGCGAACCGGCAGAGCAGTGGCAGTATGATGGCCACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTg                                                                                 <  1:413235/150‑1 (MQ=255)
                                                                        gATGAGGGACTTGTCATCCTCTGGTACTACGATGAATCGGACCGTATCACTCACCGCACGGTGAACGGCGAACCGGCAGAGCAGTGGCAGTATGATGGCCACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGcc                                                                               <  1:957971/150‑1 (MQ=255)
                                                                        gATGAGGGACTTGTCATCCTCTGGTACTACGATGAATCGGACCGTATCACTCACCGCACGGTGAACGGCGAACCGGCAGAGCAGTGGCAGTATGATGGCCACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGcc                                                                               >  2:4395696/1‑150 (MQ=255)
                                                                           gAGGGACTTGTCATCCTCTGGTACTACGATGAATCGGACCGTATCACTCACCGCACGGTGAACGGCGAACCGGCAGAGCAGTGGCAGTATGATGGCCACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTc                                                                            >  2:731547/1‑150 (MQ=255)
                                                                           gAGGGACTTGTCATCCTCTGGTACTACGATGAATCGGACCGTATCACTCACCGCACGGTGAACGGCGAACCGGCAGAGCAGTGGCAGTATGATGGCCACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTc                                                                            >  2:1825825/1‑150 (MQ=255)
                                                                            aGGGACTTGTCATCCTCTGGTACTACGATGAATCGGACCGTATCACTCACCGCACGGTGAACGGCGAACCGGCAGAGCAGTGGCAGTATGATGGCCACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTcc                                                                           >  1:3731733/1‑150 (MQ=255)
                                                                            aGGGACTTGTCATCCTCTGGTACTACGATGAATCGGACCGTATCACTCACCGCACGGTGAACGGCGAACCGGCAGAGCAGTGGCAGTATGATGGCCACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTcc                                                                           >  2:1640858/1‑150 (MQ=255)
                                                                              ggACTTGTCATCCTCTGGTACTACGATGAATCGGACCGTATCACTCACCGCACGGTGAACGGCGAACCGGCAGAGCAGTGGCAGTATGATGGCCACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCAc                                                                         <  1:2071986/150‑1 (MQ=255)
                                                                               gACTTGTCATCCTCTGGTACTACGATGAATCGGACCGTATCACTCACCGCACGGTGAACGGCGAACCGGCAGAGCAGTGGCAGTATGATGGCCACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACt                                                                        <  1:155952/150‑1 (MQ=255)
                                                                                aCTTGTCATCCTCTGGTACTACGATGAATCGGACCGTATCACTCACCGCACGGTGAACGGCGAACCGGCAGAGCAGTGGCAGTATGATGGCCACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTa                                                                       <  1:1274940/150‑1 (MQ=255)
                                                                                 cTTGTCATCCTCTGGTACTACGATGAATCGGACCGTATCACTCACCGCACGGTGAACGGCGAACCGGCTGAGCAGTGGCAGTATGATGGCCACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTAt                                                                      >  2:2455562/1‑150 (MQ=255)
                                                                                    gTCATCCTCTGGTACTACGATGAATCGGACCGTATCACTCACCGCACGGTGAACGGCGAACCGGCAGAGCAGTGGCAGTATGATGGCCACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGc                                                                   <  2:2613870/150‑1 (MQ=255)
                                                                                    gTCATCCTCTGGTACTACGATGAATCGGACCGTATCACTCACCGCACGGTGAACGGCGAACCGGCAGAGCAGTGGCAGTATGATGGCCACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGc                                                                   >  2:229488/1‑150 (MQ=255)
                                                                                     tCATCCTCTGGTACTACGATGAATCGGACCGTATCACTCACCGCACGGTGAACGGCGAACCGGCAGAGCAGTGGCAGTATGATGGCCACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCt                                                                  >  1:1099085/1‑150 (MQ=255)
                                                                                      cATCCTCTGGTACTACGATGAATCGGACCGTATCACTCACCGCACGGTGAACGGCGAACCGGCAGAGCAGTGGCAGTATGATGGCCACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTa                                                                 <  2:3837082/150‑1 (MQ=255)
                                                                                       aTCCTCTGGTACTACGATGAATCGGACCGTATCACTCACCGCACGGTGAACGGCGAACCGGCAGAGCAGTGGCAGTATGATGGCCACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTAt                                                                <  1:3983377/150‑1 (MQ=255)
                                                                                       aTCCTCTGGTACTACGATGAATCGGACCGTATCACTCACCGCACGGTGAACGGCGAACCGGCAGAGCAGTGGCAGTATGATGGCCACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTAt                                                                >  2:4596525/1‑150 (MQ=255)
                                                                                         ccTCTGGTACTACGATGAATCGGACCGTATCACTCACCGCACGGTGAACGGCGAACCGGCAGAGCAGTGGCAGTATGATGGCCACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGa                                                              >  1:2156331/1‑150 (MQ=255)
                                                                                               gTACTACGATGAATCGGACCGTACCACTCACCGCACGGTGAACGGCGAACCGGCAGAGCAGTGGCAGTATGATGGCCACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATaa                                                        <  2:2307818/150‑1 (MQ=255)
                                                                                                tactacGATGAATCGGACCGTATCACTCACCGCACGGTGAACGGCGAACCGGCAGAGCAGTGGCAGTATGATGGCCACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATaaa                                                       >  1:2999805/1‑150 (MQ=255)
                                                                                                 actacGATGAATCGGACCGTATCACTCACCGCACGGTGAACGGCGAACCGGCAGAGCAGTGGCAGTATGATGGCCACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAg                                                      <  1:3491514/150‑1 (MQ=255)
                                                                                                  ctacGATGAATCGGACCGTATCACTCACCGCACGGTGAACGGCGAACCGGCAGAGCAGTGGCAGTATGATGGCCACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAgg                                                     <  1:4096140/150‑1 (MQ=255)
                                                                                                      gATGAATCGGACCGTATCACTCACCGCACGGTGAACGGCGAACCGGCTGAGCAGTGGCAGTATGATGGCCACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGgccgc                                                 >  1:78040/1‑150 (MQ=255)
                                                                                                      gATGAATCGGACCGTATCACTCACCGCACGGTGAACGGCGAACCGGCAGAGCAGTGGCAGTATGATGACCACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGgccgc                                                 <  2:2681276/150‑1 (MQ=255)
                                                                                                       aTGAATCGGACCGTATCACTCACCGCACGGTGAACGGCGAACCGGCTGAGCAGTGGCAGTATGATGGCCACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGgccgcc                                                >  2:1481375/1‑150 (MQ=255)
                                                                                                       aTGAATCGGACCGTATCACTCACCGCACGGTGAACGGCGAACCGGCAGAGCAGTGGCAGTATGATGGCCACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGgccgcc                                                >  2:2733557/1‑150 (MQ=255)
                                                                                                       aTGAATCGGACCGTATCACTCACCGCACGGTGAACGGCGAACCGGCAGAGCAGTGGCAGTATGATGGCCACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGgccgcc                                                >  2:1558743/1‑150 (MQ=255)
                                                                                                         gAATCGGACCGTATCACTCACCGCACGGTGAACGGCGAACCGGCTGAGCATTGGCAGTATGATGGCCACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGTCCACAGAGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTg                                              >  2:2010148/1‑150 (MQ=37)
                                                                                                         gAATCGGACCGTATCACTCACCGCACGGTGAACGGCGAACCGGCTGAGCAGTGGCAGTATGATGGCCACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTg                                              <  1:3823396/150‑1 (MQ=255)
                                                                                                         gAATCGGACCGTATCACTCACCGCACGGTGAACGGCGAACCGGCAGAGCAGTGGCAGTATGATGGCCACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTg                                              >  2:3096470/1‑150 (MQ=255)
                                                                                                         gAATCGGAACGTATCACTCACCGCACGGTGAACGGCGAACCGGCTGAGCAGTGGCAGTAAGATGGCCACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCAATGACGATAAAGgccgccag                                              >  2:366076/1‑148 (MQ=255)
                                                                                                          aaTCGGACCGTATCACTCACCGCACGGTGAACGGCGAACCGGCAGAGCAGTGGCAGTATGATGGCCACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGa                                             <  1:3291369/150‑1 (MQ=255)
                                                                                                           aTCGGACCGTATCACTCACCGCACGGTGAACGGCGAACCGGCAGAGCAGTGGCAGTATGATGGCCACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGAc                                            >  2:4118539/1‑150 (MQ=255)
                                                                                                           aTCGGACCGTATCACTCACCGCACGGTGAACGGCGAACCGGCAGAGCAGTGGCAGTATGATGGCCACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGAc                                            >  1:1877245/1‑150 (MQ=255)
                                                                                                             cGGACCGTATCACTCACCGCACGGTGAACGGCGAACCGGCAGAGCAGTGGCAGTATGATGGCCACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCg                                          >  2:2487884/1‑150 (MQ=255)
                                                                                                              ggACCGTATCACTCACCGCACGGTGAACGGCGAACCGGCAGAGCAGTGGCAGTATGATGGCCACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCgg                                         >  1:191788/1‑150 (MQ=255)
                                                                                                                aCCGTATCACTCACCGCACGGTGAACGGCGAACCGGCAGAGCAGTGGCAGTATGATGGCCACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCg                                       <  2:704744/150‑1 (MQ=255)
                                                                                                                aCCGTATCACTCACCGCACGGTGAACGGCGAACCGGCAGAGCAGTGGCAGTATGATGGCCACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCg                                       >  1:4514774/1‑150 (MQ=255)
                                                                                                                   gTATCACTCACCGCACGGTGAACGGCGAACCGGCAGAGCAGTGGCAGTATGATGGCCACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAAt                                    <  2:1867529/150‑1 (MQ=255)
                                                                                                                     aTCACTCACCGCACGGTGAACGGCGAACCGGCTGAGCAGTGGCAGTATGATGGCCACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGc                                  >  1:1791979/1‑150 (MQ=255)
                                                                                                                     aTCACTCACCGCACGGTGAACGGCGAACCGGCAGAGCAGTGGCAGTATGATGGCCACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGc                                  >  2:3452492/1‑150 (MQ=255)
                                                                                                                     aTCACTCACAGCACGGTGAACGGCGAACCGGCAGAGCAGTGGCAGTATGATGGCCACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGc                                  >  2:3176154/1‑150 (MQ=255)
                                                                                                                         ctcacCGCACGGTGAACGGCGAACCGGCAGAGCAGTGGCAGTATGATGGCCACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGa                              >  1:3649664/1‑150 (MQ=255)
                                                                                                                          tcacCGCACGGTGAACGGCGAACCGGCAGAGCAGTGGCAGTATGATGGCCACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGAc                             <  1:489188/150‑1 (MQ=255)
                                                                                                                           cacCGCACGGTGAACGGCGAACCGGCTGAGCAGTGGCAGTATGATGGCCACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACg                            <  2:3499477/150‑1 (MQ=255)
                                                                                                                           cacCGCACGGTGAACGGCGAACCGGCTGAGCAGTGGCAGTATGATGGCCACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACg                            <  2:3684996/150‑1 (MQ=255)
                                                                                                                              cGCACGGTGAACGGCGAACCGGCAGAGCAGTGGCAGTATGATGGCCACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTg                         <  2:1128500/150‑1 (MQ=255)
                                                                                                                                  cGGTGAACGGCGAACCGGCAGAGCAGTGGCAGTATGATGGCCACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGgaga                     >  2:1168472/1‑150 (MQ=255)
                                                                                                                                   ggTGAACGGCGAACCGGCAGAGCAGTGGCAGTATGATGGCCACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGaa                    >  1:2226059/1‑150 (MQ=255)
                                                                                                                                      gAACGGCGAACCGGCAGAGCAGTGGCAGTATGATGGCCACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGTCGGTGGAGAAccc                 <  1:4204131/150‑1 (MQ=255)
                                                                                                                                        aCGGCGAACCGGCAGAGCAGTGGCAGTATGATGGCCACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCgg               <  1:1191282/150‑1 (MQ=255)
                                                                                                                                          ggCGAACCGGCTGAGCAGTGGCAGTATGATGGCCACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGgag             >  1:3393948/1‑150 (MQ=255)
                                                                                                                                           gCGAACCGGGTGAGCTGTTGCAGTTTGATGGCCTCGGCTGGCTGACAGTCATAATCCACCTGAGCGAATGCCACCGTGTTGACGTCCACTATGGATATGACGATAATGGAAGCCTGACCGGCGAATGCCTTACGGTGGAGAACCCGgaga            <  1:4010548/150‑1 (MQ=21)
                                                                                                                                            cGAACCGGCAGAGCAGTGGCAGTATGATGGCCACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGAGAc           >  2:504579/1‑150 (MQ=255)
                                                                                                                                            cGAACCGGCAGAGCAGTGGCAGTATGATGGCCACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTAGGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGAGAc           <  2:700313/150‑1 (MQ=255)
                                                                                                                                             gAACCGGCAGAGCAGTGGCAGTATGATGGCCACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACg          >  1:2228613/1‑150 (MQ=255)
                                                                                                                                                ccGGCAGAGCAGTGGCAGTATGATGGCCACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCAAGACGGTGGAGAACCCGGAGACgggg       <  1:2982608/150‑1 (MQ=255)
                                                                                                                                                 cGGCAGAGCAGTGGCAGTATGATGGCCACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGATACggggg      <  2:4248197/150‑1 (MQ=255)
                                                                                                                                                 cGGCAGAGCAGTGGCAGTATGATGGCCACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACggggg      <  2:4047369/150‑1 (MQ=255)
                                                                                                                                                  ggCAGAGCAGTGGCAGTATGATGGCCACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGa     >  1:4310921/1‑150 (MQ=255)
                                                                                                                                                     agagCAGTGGCAGTATGATGGCCACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAAct  >  2:1719276/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CAGCAGGGCGGCTTTGACGGGCGGACGCAACGTTATCATTATGACCTGACCGGAAAACTCACACAGAGTGAGGATGAGGGACTTGTCATCCTCTGGTACTACGATGAATCGGACCGTATCACTCACCGCACGGTGAACGGCGAACCGGCAGAGCAGTGGCAGTATGATGGCCACGGCTGGCTGACAGACATCAGCCACCTGAGCGAAGGCCACCGTGTTGCCGTCCACTATGGCTATGACGATAAAGGCCGCCTGACCGGCGAATGCCAGACGGTGGAGAACCCGGAGACGGGGGAACT  >  CP009273/520724‑521022

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: