Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,967,504 A→C 10.2% K179N (AAA→AAC rffH → glucose‑1‑phosphate thymidylyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,967,5040AC10.2% 475.6 / 50.8 187K179N (AAA→AACrffHglucose‑1‑phosphate thymidylyltransferase
Reads supporting (aligned to +/- strand):  ref base A (94/74);  new base C (8/11);  total (102/85)
Fisher's exact test for biased strand distribution p-value = 3.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.42e-01

Only 100 of 192 total aligned reads displayed.

GACGGTTTTTGGCTATCAGGTGATGGACCCGGAACGCTTTGGCGTGGTGGAGTTTGACGACAATTTCCGCGCTATCTCGCTGGAAGAAAAGCCAAAACAGCCGAAGTCAAACTGGGCGGTGACCGGGCTTTATTTCTACGACAGTAAAGTCGTGGAGTACGCAAAGCAGGTGAAGCCGTCGGAGCGTGGTGAACTGGAGATTACCTCCATCAACCAGATGTACCTCGAGGCGGGCAACCTGACCGTTGAACTGCTCGGGCGCGGATTTGCCTGGCTGGATACCGGCACTCACGAC  >  CP009273/3967357‑3967651
                                                                                                                                                   |                                                                                                                                                   
 aCGGTTTTTGGCTATCAGGTGATGGACCCGGAACGCTTTGGCGTGGTGGAGTTTGACGACAATTTCCGCGCTATCTCGCTGGAAGAAAAGCCAAAACAGCCGAAGTCAAACTGGGCTGTGACCGGGCTTTATTTCTACGACAGTAAAg                                                                                                                                                    <  1:3538210/148‑1 (MQ=255)
 aCGGTTTTTGGCTATCAGGTGATGGACCCGGAACGCTTTGGCGTGGTGGAGTTTGACGACAATTTCCGCGCTATCTCGCTGGAAGAAAAGCCAAAACAGCCGAAGTCAAACTGGGCTGTGACCGGGCTTTATTTCTACGACAGTAAAg                                                                                                                                                    <  1:1776222/148‑1 (MQ=255)
  cGGTTTTTGGCTATCAGGTGATGGACCCGGAACGCTTTGGCGTGGTGGAGTTTGACGACAATTTCCGCGCTATCTCGCTGGAAGAAAAGCCAAAACAGCCGAAGTCAAACTGGGCTGTGACCGGGCTTTATTTCTACGACAGTAAAGt                                                                                                                                                   <  1:2052845/148‑1 (MQ=255)
  cGGTTTTTGGCTATCAGGTGATGGACCCGGAACGCTTTGGCGTGGTGGAGTTTGACGACAATTTCCGCGCTATCTCGCTGGAAGAAAAGCCAAAACAGCCGAAGTCAAACTGGGCTGTGACCGGGCTTTATTTCTACGACAGTAAAGt                                                                                                                                                   <  2:3506566/148‑1 (MQ=255)
  cGGTTTTTGGCTATCAGGTGATGGACCCGGAACGCTTTGGCGTGGTGGAGTTTGACGACAATTTCCGCGCTATCTCGCTGGAAGAAAAGCCAAAACAGCCGAAGTCAAACTGGGCTGTGACCGGGCTTTATTTCTACGACAGTAAAGt                                                                                                                                                   <  2:3506593/148‑1 (MQ=255)
    gTTTTTGGCTATCAGGTGATGGACCCGGAACGCTTTGGCGTGGTGGAGTTTGACGACAATTTCCGCGCTATCTCGCTGGAAGAAAAGCCAAAACAGCCGAAGTCAAACTGGGCTGTGACCGGGCTTTATTTCTACGACAGTAAAGTCg                                                                                                                                                 <  1:1684837/148‑1 (MQ=255)
               tCAGGTGATGGACCCGGAACGCTTTGGCGTGGTGGAGTTTGACGACAATTTCCGCGCTATCTCGCTGGAAGAAAAGCCAAAACAGCCGAAGTCAAACTGGGCTGTGACCGGGCTTTATTTCTACGACAGTAAAGTCGTGGAGTACGCa                                                                                                                                      >  1:2039255/1‑148 (MQ=255)
                  ggTGATGGACCCGGAACGCTTTGGCGTGGTGGCGTTTGACGACAATTTCCGCGCTATCTCGCTGGAAGAAAAGCCAAAACAGCCGAAGTCAAACTGGGCTGTGACCGGGCTTTATTTCTACGACAGTAAAGTCGTGGAGTACGCAAAg                                                                                                                                   <  2:1592645/148‑1 (MQ=255)
                      aTGGACCCGGAACGCTTGGGCGTGGTGGAGTTTGACGACAATTTCCGCGCTATCTCGCTGGAAGAAAAGCCAAAACAGCCGAAGTCAAACTGGGCTGTGACCGGGCTTTATTTCTACGACAGTAAAGTCGTGGAGTACGCAAAGCAgg                                                                                                                               <  2:1547057/148‑1 (MQ=255)
                          aCCCGGAACGCTTTGGCGTGGTGGAGTTTGACGACAATTTCCGCGCTATCTCGCTGGAAGAAAAGCCAAAACAGCCGAAGTCAAACTGGGCTGTGACCGGGCTTTATTTCTACGACAGTAAAGTCGTGGAGTACGCAAAGCAGGTGaa                                                                                                                           >  1:2120527/1‑148 (MQ=255)
                           cccGGAACGCTTTGGCGTGGTGGAGTTTGACGACAATTTCCGCGCTATCTCGCTGGAAGAAAAGCCAAAACAGCCGAAGTCAAACTGGGCTGTGACCGGGCTTTATTTCTACGACAGTAAAGTCGTGGAGTACGCAAAGCAGGTGAAg                                                                                                                          <  2:2966201/148‑1 (MQ=255)
                               gAACGCTTTGGCGTGGTGGAGTTTGACGACAATTTCCGCGCTATCTCGCTGGAAGAAAAGCCAAAACAGCCGAAGTCAAACTGGGCTGTGACCGGGCTTTATTTCTACGACAGTAACGTCGTGGAGTACGCAAAGCAGGTGAAGCCGt                                                                                                                      <  2:1914152/148‑1 (MQ=255)
                               gAACGCTTTGGCGTGGTGGAGTTTGACGACAATTTCCGCGCTATCTCGCTGGAAGAAAAGCCAAAACAGCCGAAGTCAAACTGGGCTGTGACCGGGCTTTATTTCTACGACAGTAAAGTCGTGGAGTACGCAAAGCAGGTGAAGCCGt                                                                                                                      <  1:3224605/148‑1 (MQ=255)
                                aaCGCTTTGGCGTGGTGGAGTTTGACGACAATTTCCGCGCTATCTCGCTGGAAGAAAAGCCAAAACAGCCGAAGTCAAACTGGGCTGTGACCGGGCTTTATTTCTACGACAGTAAAGTCGTGGAGTACGCAAAGCAGGTGAAGCCGTc                                                                                                                     >  1:2346104/1‑148 (MQ=255)
                                  cGCTTTGGCGTGGTGGAGTTTGACGACAATTTCCGCGCTATCTCGCTGGAAGAAAAGCCAAAACAGCCGAAGTCAAACTGGGCTGTGACCGGGCTTTATTTCTACGACAGTAAAGTCGTGGAGTACGCAAAGCAGGTGAAGCCGTCgg                                                                                                                   <  2:2976436/148‑1 (MQ=255)
                                    cTTTGGCGTGGTGGAGTTTGACGACAATTTCCGCGCTATCTCGCTGGAAGAAAAGCCAAAACAGCCGAAGTCAAACTGGGCTGTGACCGGGCTTTATTTCTACGACAGTAACGTCGTGGAGTACGCAAAGCAGGTGAAGCCGTCGGAg                                                                                                                 >  1:1396128/1‑148 (MQ=255)
                                    cTTTGGCGTGGTGGAGTTTGACGACAATTTCCGCGCTATCTCGCTGGAAGAAAAGCCAAAACAGCCGAAGTCAAACTGGGCTGTGACCGGGCTTTATTTCTACGACAGTAAAGTCGTGGAGTACGCAAAGCAGGTGAAGCCGTCGGAg                                                                                                                 >  2:1929522/1‑148 (MQ=255)
                                    cTTTGGCGTGGTGGAGTTTGACGACAATTTCCGCGCTATCTCGCTGGAAGAAAAGCCAAAACAGCCGAAGTCAAACTGGGCTGTGACCGGGCTTTATTTCTACGACAGTAAAGTCGTGGAGTACGCAAAGCAGGTGAAGCCGTCGGAg                                                                                                                 <  2:497164/148‑1 (MQ=255)
                                    cTTTGGCGTGGTGGAGTTTGACGACAATTTCCGCGCTATCTCGCTGGAAGAAAAGCCAAAACAGCCGAAGTCAAACTGGGCTGTGACCGGGCTTTATTTCTACGACAGTAAAGTCGTGGAGTACGCAAAGCAGGTGAAGCCGTCGGAg                                                                                                                 >  2:256891/1‑148 (MQ=255)
                                     tttGGCGTGGTGGAGTTTGACGACAATTTCCGCGCTATCTCGCTGGAAGAAAAGCCAAAACAGCCGAAGTCAAACTGGGCTGTGACCGGGCTTTATTTCTACGACAGTAAAGTCGTGGAGTACGCAAAGCAGGTGAAGCCGTCGGAGc                                                                                                                >  1:2205054/1‑148 (MQ=255)
                                              gtggAGTTTGACGACAATTTCCGCGCTATCTCGCTGGAAGAAAAGCCAAAACAGCCGAAGTCAAACTGGGCTGTGACCGGGCTTTATTTCTACGACAGTAAAGTCGTGGAGTACGCAAAGCAGGTGAAGCCGTCGGAGCGTGGTGAAc                                                                                                       <  2:310724/148‑1 (MQ=255)
                                              gtggAGTTTGACGACAATTTCCGCGCTATCTCGCTGGAAGAAAAGCCAAAACAGCCGAAGTCAAACTGGGCTGTGACCGGGCTTTATTTCTACGACAGTAAAGTCGTGGAGTACGCAAAGCAGGTGAAGCCGTCGGAGCGTGGTGAAc                                                                                                       <  2:3182189/148‑1 (MQ=255)
                                               tggAGTTTGACGACAATTTCCGCGCTATCTCGCTGGAAGAAAAGCCAAAACAGCCGAAGTCAAACTGGGCTGTGACCGGGCTTTATTTCTACGACAGTAAAGTCGTGGAGTACGCAAAGCAGGTGAAGCCGTCGGAGCGTGGTGAACt                                                                                                      >  2:2790117/1‑148 (MQ=255)
                                               tggAGTTTGACGACAATTTCCGCGCTATCTCGCTGGAAGAAAAGCCAAAACAGCCGAAGTCAAACTGGGCTGTGACCGGGCTTTATTTCTACGACAGTAAAGTCGTGGAGTACGCAAAGCAGGTGAAGCCGTCGGAGCGTGGTGAACt                                                                                                      >  2:3130726/1‑148 (MQ=255)
                                               tggAGTTTGACGACAATTTCCGCGCTATCTCGCTGGAAGAAAAGCCAAAACAGCCGAAGTCAAACTGGGCGGTGACCGGGCTTTATTTCTACGACAGTAAAGTCGTGGAGTACGCAAAGCAGGTGAAGCCGTCGGAGCGTGGTGAACt                                                                                                      <  2:761709/148‑1 (MQ=255)
                                                ggAGTTTGACGACAATTTCCGCGCTATCTCGCTGGAAGAAAAGCCAAAACAGCCGAAGTCAAACTGGGCTGTGACCGGGCTTTATTTCTACGACAGTAAAGTCGTGGAGTACGCAAAGCAGGTGAAGCCGTCGGAGCGTGGTGAACTg                                                                                                     >  1:3384272/1‑148 (MQ=255)
                                                 gAGTTTGACGACAATTTCCGCGCTATCTCGCTGGAAGAAAAGCCAAAACAGCCGAAGTCAAACTGGGCTGTGACCGGGCTTTATTTCTACGACAGTAAAGTCGTGGAGTACGCAAAGCAGGTGAAGCCGTCGGAGCGTGGTGAACTgg                                                                                                    <  2:2296161/148‑1 (MQ=255)
                                                 gAGTTTGACGACAATTTCCGCGCTATCTCGCTGGAAGAAAAGCCAAAACAGCCGAAGTCAAACTGGGCTGTGACCGGGCTTTATTTCTACGACAGTAAAGTCGTGGAGTACGCAAAGCAGGTGAAGCCGTCGGAGCGTGGTGAACTgg                                                                                                    >  2:3530832/1‑148 (MQ=255)
                                                   gTTTGACGACAATTTCCGCGCTATCTCGCTGGAAGAAAAGCCAAAACAGCCGAAGTCAAACTGGGCTGTGACCGGGCTTTATTTCTACGACAGTAAAGTCGTGGAGTACGCAAAGCAGGTGAAGCCGTCGGAGCGTGGTGAACTGgag                                                                                                  <  1:187624/148‑1 (MQ=255)
                                                    tttGACGACAATTTCCGCGCTATCTCGCTGGAAGAAAAGCCAAAACAGCCGAAGTCAAACTGGGCTGTGACCGGGCTTTATTTCTACGACAGTAACGTCGTGGAGTACGCAAAGCAGGTGAAGCCGTCGGAGCGTGGTGAACTGgaga                                                                                                 >  1:1044182/1‑148 (MQ=255)
                                                    tttGACGACAATTTCCGCGCTATCTCGCTGGAAGAAAAGCCAAAACAGCCGAAGTCAAACTGGGCTGTGACCGGGCTTTATTTCTACGACAGTAACGTCGTGGAGTACGCAAAGCAGGTGAAGCCGTCGGAGCGTGGTGAACTGgaga                                                                                                 >  1:131082/1‑148 (MQ=255)
                                                    tttGACGACAATTTCCGCGCTATCTCGCTGGAAGAAAAGCCAAAACAGCCGAAGTCAAACTGGGCTGTGACCGGGCTTTATTTCTACGACAGTAACGTCGTGGAGTACGCAAAGCAGGTGAAGCCGTCGGAGCGTGGTGAACTGgaga                                                                                                 >  1:1929478/1‑148 (MQ=255)
                                                    tttGACGACAATTTCCGCGCTATCTCGCTGGAAGAAAAGCCAAAACAGCCGAAGTCAAACTGGGCTGTGACCGGGCTTTATTTCTACGACAGTAAAGTCGTGGAGTACGCAAAGCAGGTGAAGCCGTCGGAGCGTGGTGAACTGgaga                                                                                                 >  1:3424005/1‑148 (MQ=255)
                                                     ttGACGACAATTTCCGCGCTATCTCGCTGGAAGAAAAGCCAAAACAGCCGAAGTCAAACTGGGCTGTGACCGGGCTTTATTTCTACGACAGTAAAGTCGTGGAGTACGCAAAGCAGGTGAAGCCGTCGGAGCGTGGTGAACTGGAGAt                                                                                                >  2:764434/1‑148 (MQ=255)
                                                     ttGACGACAATTTCCGCGCTATCTCGCTGGAAGAAAAGCCAAAACAGCCGAAGTCAAACTGGGCTGTGACCGGGCTTTATTTCTACGACAGTAAAGTCGTGGAGTACGCAAAGCAGGTGAAGCCGTCGGAGCGTGGTGAACTGGAGAt                                                                                                >  2:2551590/1‑148 (MQ=255)
                                                       gacgacAATTTCCGCGCTATCTCGCTGGAAGAAAAGCCAAAACAGCCGAAGTCAAACTGGGCTGTGACCGGGCTTTATTTCTACGACAGTAAAGTCGTGGAGTACGCAAAGCAGGTGAAGCCGTCGGAGCGTGGTGAACTGGAGATTa                                                                                              <  1:207497/148‑1 (MQ=255)
                                                           acAATTTCCGCGCTATCTCGCTGGAAGAAAAGCCAAAACAGCCGAAGTCAAACTGGGCTGTGACCGGGCTTTATTTCTACGACAGTAAAGTCGTGGAGTACGCAAAGCAGGTGAAGCCGTCGGAGCGTGGTGAACTGGAGATTACCTc                                                                                          <  2:1227374/148‑1 (MQ=255)
                                                            cAATTTCCGCGCTATCTCGCTGGAAGAAAAGCCAAAACAGCCGAAGTCAAACTGGGCTGTGACCGGGCTTTATTTCTACGACAGTAAAGTCGTGGAGTACGCAAAGCAGGTGAAGCCGTCGGAGCGTGGTGAACTGGAGATTACCTcc                                                                                         >  1:1943749/1‑148 (MQ=255)
                                                              aTTTCCGCGCTATCTCGCTGGAAGAAAAGCCAAAACAGCCGAAGTCAAACTGGGCTGTGACCGGGCTTTATTTCTACGACAGTAAAGTCGTGGAGTACGCAAAGCAGGTGAAGCCGTCGGAGCGTGGTGAACTGGAGATTACCTCCAt                                                                                       >  2:2147400/1‑148 (MQ=255)
                                                                ttCCGCGCTATCTCGCTGGAAGAAAAGCCAAAACAGCCGAAGTCAAACTGGGCTGTGACCGGGCTTTATTTCTACGACAGTAACGTCGTGGAGTACGCAAAGCAGGTGAAGCCGTCGGAGCGTGGTGAACTGGAGATTACCTCCATCa                                                                                     >  2:3323538/1‑148 (MQ=255)
                                                                   cgcgCTATCTCGCTGGAAGAAAAGCCAAAACAGCCGAAGTCAAACTGGGCTGTGACCGGGCTTTATTTCTACGACAGTAAAGTCGTGGAGTACGCAAAGCAGGTGAAGCCGTCGGAGCGTGGTGAACTGGAGATTACCTCCATCAAcc                                                                                  <  1:1594066/148‑1 (MQ=255)
                                                                     cgcTATCTCGCTGGAAGAAAAGCCAAAACAGCCGAAGTCAAACTGGGCTGTGACCGGGCTTTATTTCTACGACAGTAAAGTCGTGGAGTACGCAAAGCAGGTGAAGCCGTCGGAGCGTGGTGAACTGGAGATTACCTCCATCAACCAg                                                                                >  2:1247144/1‑148 (MQ=255)
                                                                     cgcTATCTCGCTGGAAGAAAAGCCAAAACAGCCGAAGTCAAACTGGGCGGTGACCGGGCTTTATTTCTACGACAGTAAAGTCGTGGAGTACGCAAAGCAGGTGAAGCCGTCGGAGCGTGGTGAACTGGAGATTACCTCCATCAACCAg                                                                                >  1:3436640/1‑148 (MQ=255)
                                                                      gcTATCTCGCTGGAAGAAAAGCCAAAACAGCCGAAGTCAAACTGGGCTGTGACCGGGCTTTATTTCTACGACAGTAACGTCGTGGAGTACGCAAAGCAGGTGAAGCCGTCGGAGCGTGGTGAACTGGAGATTACCTCCATCAACCAGa                                                                               <  1:3491666/148‑1 (MQ=255)
                                                                      gcTATCTCGCTGGAAGAAAAGCCAAAACAGCCGAAGTCAAACTGGGCTGTGACCGGGCTTTATTTCTACGACAGTAAAGTCGTGGAGTACGCAAAGCAGGTGAAGCCGTCGGAGCGTGGTGAACTGGAGATTACCTCCATCAACCAGa                                                                               <  1:1860209/148‑1 (MQ=255)
                                                                       cTATCTCGCTGGAAGAAAAGCCAAAACAGCCGAAGTCAAACTGGGCTGTGACCGGGCTTTATTTCTACGACAGTAAAGTCGTGGAGTACGCAAAGCAGGTGAAGCCGTCGGAGCGTGGTGAACTGGAGATTACCTCCATCAACCAGAt                                                                              >  2:269759/1‑148 (MQ=255)
                                                                        tATCTCGCTGGAAGAAAAGCCAAAACAGCCGAAGTCAAACTGGGCTGTGACCGGGCTTTATTTCTACGACAGTAAAGTCGTGGAGTACGCAAAGCAGGTGAAGCCGTCGGAGCGTGGTGAACTGGAGATTACCTCCATCAACCAGATg                                                                             >  2:290968/1‑148 (MQ=255)
                                                                        tATCTCGCTGGAAGAAAAGCCAAAACAGCCGAAGTCAAACTGGGCGGTGACCGGGCTTTATTTCTACGACAGTAAAGTCGTGGAGTACGCAAAGCAGGTGAAGCCGTCGGAGCGTGGTGAACTGGAGATTACCTCCATCAACCAGATg                                                                             <  2:2078099/148‑1 (MQ=255)
                                                                          tctcGCTGGAAGAAAAGCCAAAACAGCCGAAGTCAAACTGGGCTGTGACCGGGCTTTATTTCTACGACAGTAAAGTCGTGGAGTACGCAAAGCAGGTGAAGCCGTCGGAGCGTGGTGAACTGGAGATTACCTCCATCAACCAGATGTa                                                                           <  1:261666/148‑1 (MQ=255)
                                                                          tctcGCTGGAAGAAAAGCCAAAACAGCCGAAGTCAAACTGGGCTGTGACCGGGCTTTATTTCTACGACAGTAAAGTCGTGGAGTACGCAAAGCAGGTGAAGCCGTCGGAGCGTGGTGAACTGGAGATTACCTCCATCAACCAGATGTa                                                                           <  1:577334/148‑1 (MQ=255)
                                                                            tcGCTGGAAGAAAAGCCAAAACAGCCGAAGTCAAACTGGGCTGTGACCGGGCTTTATTTCTACGACAGTAAAGTCGTGGAGTACGCAAAGCAGGTGAAGCCGTCGGAGCGTGGTGAACTGGAGATTACCTCCATCAACCAGATGTAcc                                                                         >  1:3494220/1‑148 (MQ=255)
                                                                            tcGCTGGAAGAAAAGCCAAAACAGCCGAAGTCAAACTGGGCTGTGACCGGGCTTTATTTCTACGACAGTAAAGTCGTGGAGTACGCAAAGCAGGTGAAGCCGTCGGAGCGTGGTGAACTGGAGATTACCTCCATCAACCAGATGTAcc                                                                         <  2:3242262/148‑1 (MQ=255)
                                                                             cGCTGGAAGAAAAGCCAAAACAGCCGAAGTCAAACTGGGCTGTGACCGGGCTTTATTTCTACGACAGTAAAGTCGTGGAGTACGCAAAGCAGGTGAAGCCGTCGGAGCGTGGTGAACTGGAGATTACCTCCATCAACCAGATGTACCt                                                                        >  1:678084/1‑148 (MQ=255)
                                                                               cTGGAAGAAAAGCCAAAACAGCCGAAGTCAAACTGGGCTGTGACCGGGCTTTATTTCTACGACAGTAAAGTCGTGGAGTACGCAAAGCAGGTGAAGCCGTCGGAGCGTGGTGAACTGGAGATTACCTCCATCAACCAGATGTACCTCg                                                                      <  2:2827626/148‑1 (MQ=255)
                                                                                 ggAAGAAAAGCCAAAACAGCCGACGTCAAACTGGGCTGTGACCGGGCTTTATTTCTACGACAGTAAAGTCGTGGAGTACGCAAAGCAGGTGAAGCCGTCGGAGCGTGGTGAACTGGAGATTACCTCCATCAACCAGATGTACCTCGAg                                                                    <  2:89843/148‑1 (MQ=255)
                                                                                     gaaAAGCCAAAACAGCCGAAGTCAAACTGGGCTGTGACCGGGCTTTATTTCTACGACAGTAAAGTCGTGGAGTACGCAAAGCAGGTGAAGCCGTCGGAGCGTGGTGAACTGGAGATTACCTCCATCAACCAGATGTACCTCGAGGCgg                                                                >  2:3074898/1‑148 (MQ=255)
                                                                                      aaaaGCCAAAACAGCCGAAGTCAAACTGGGCTGTGACCGGGCTTTATTTCTACGACAGTAAAGTCGTGGAGTACGCAAAGCAGGTGAAGCCGTCGGAGCGTGGTGAACTGGAGATTACCTCCATCAACCAGATGTACCTCGAGGCggg                                                               >  1:1743980/1‑148 (MQ=255)
                                                                                      aaaaGCCAAAACAGCCGAAGTCAAACTGGGCTGTGACCGGGCTTTATTTCTACGACAGTAAAGTCGTGGAGTACGCAAAGCAGGTGAAGCCGTCGGAGCGTGGTGAACTGGAGATTACCTCCATCAACCAGATGTACCTCGAGGCggg                                                               <  1:998224/148‑1 (MQ=255)
                                                                                        aaGCCAAAACAGCCGAAGTCAAACTGGGCTGTGACCGGGCTTTATTTCTACGACAGTAAAGTCGTGGAGTACGCAAAGCAGGTGAAGCCGTCGGAGCGTGGTGAACTGGAGATTACCTCCATCAACCAGATGTACCTCGAGGCGGGCa                                                             >  2:3554093/1‑148 (MQ=255)
                                                                                        aaGCCAAAACAGCCGAAGTCAAACTGGGCTGTGACCGGGCTTTATTTCTACGACAGTAAAGTCGTGGAGTACGCAAAGCAGGTGAAGCCGTCGGAGCGTGGTGAACTGGAGATTACCTCCATCAACCAGATGTACCTCGAGGCGGGCa                                                             <  1:2685111/148‑1 (MQ=255)
                                                                                         aGCCAAAACAGCCGAAGTCAAACTGGGCTGTGACCGGGCTTTATTTCTACGACAGTAAAGTCGTGGAGTACGCAAAGCAGGTGAAGCCGTCGGAGCGTGGTGAACTGGAGATTACCTCCATCAACCAGATGTACCTCGAGGCGGGCaa                                                            <  1:876552/148‑1 (MQ=255)
                                                                                            cAAAACAGCCGAAGTCAAACTGGGCTGTGACCGGGCTTTATTTCTACGACAGTAAAGTCGTGGAGTACGCAAAGCAGGTGAAGCCGTCGGAGCGTGGTGAACTGGAGATTACCTCCATCaa                                                                                    >  2:163587/1‑121 (MQ=255)
                                                                                               aaCAGCCGAAGTCAAACTGGGCTGTGACCGGGCTTTATTTCTACGACAGTAAAGTCGTGGAGTACGCAAAGCAGGTGAAGCCGTCGGAGCGTGGTGAACTGGAGATTACCTCCATCAACCAGATGTACCTCGAGGCGGGCAACCTGAc                                                      >  2:301083/1‑148 (MQ=255)
                                                                                                   gCCGAAGTCAAACTGGGCTGTGACCGGGCTTTATTTCTACGACAGTAAAGTCGTGGAGTACGCAAAGCAGGTGAAGCCGTCGGAGCGTGGTGAACTGGAGATTACCTCCATCAACCAGATGTACCTCGAGGCGGGCAACCTGACCGtt                                                  <  1:1172884/148‑1 (MQ=255)
                                                                                                      gAAGTCAAACTGGGCTGTGACCGGGCTTTATTTCTACGACAGTAAAGTCGTGGAGTACGCAAAGCAGGTGAAGCCGTCGGAGCGTGGTGAACTGGAGATTACCTCCATCAACCAGATGTACCTCGAGGCGGGCAACCTGACCGTTGaa                                               <  1:744922/148‑1 (MQ=255)
                                                                                                       aaGTCAAACTGGGCTGTGACCGGGCTTTATTTCTACGACAGTAAAGTCGTGGAGTACGCAAAGCAGGTGAAGCCGTCGGAGCGTGGTGAACTGGAGATTACCTCCATCAACCAGATGTACCTCGAGGCGGGCAACCTGACCGTTGAAc                                              <  1:1029222/148‑1 (MQ=255)
                                                                                                         gTCAAACTGGGCTGTGACCGGGCTTTATTTCTACGACAGTAAAGTCGTGGAGTACGCAAAGCAGGTGAAGCCGTCGGAGCGTGGTGAACTGGAGATTACCTCCATCAACCAGATGTACCTCGAGGCGGGCAACCTGACCGTTGAACTg                                            <  2:2284339/148‑1 (MQ=255)
                                                                                                         gTCAAACTGGGCGGTGACCGGGCTTTATTTCTACGACAGTAAAGTCGTGGAGTACGCAAAGCAGGTGAAGCCGTCGGAGCGTGGTGAACTGGAGATTACCTCCATCAACCAGATGTACCTCGAGGCGGGCAACCTGACCGTTGAACTg                                            >  2:579774/1‑148 (MQ=255)
                                                                                                          tCAAACTGGGCTGTGACCGGGCTTTATTTCTACGACAGTAAAGTCGTGGAGTACGCAAAGCAGGTGAAGCCGTCGGAGCGTGGTGAACTGGAGATTACCTCCATCAACCAGATGTACCTCGAGGCGGGCAACCTGACCGTTGAACTGc                                           >  2:537300/1‑148 (MQ=255)
                                                                                                               cTGGGCTGTGACCGGGCTTTATTTCTACGACAGTAAAGTCGTGGAGTACGCAAAGCAGGTGAAGCCGTCGGAGCGTGGTGAACTGGAGATTACCTCCATCAACCAGATGTACCTCGAGGCGGACAACCTGACCGTTGAACTGCTCggg                                      >  1:3069072/1‑148 (MQ=255)
                                                                                                                 gggCTGTGACCGGGCTTTATTTCTACGACAGTAAAGTCGTGGAGTACGCAAAGCAGGTGAAGCCGTCGGAGCGTGGTGAACTGGAGATTACCTCCATCAACCAGATGTACCTCGAGGCGGGCAACCTGACCGTTGAACTGCTCGGgcg                                    >  2:731199/1‑148 (MQ=255)
                                                                                                                  ggCTGTGACCGGGCTTTATTTCTACGACAGTAAAGTCGTGGAGTACGCAAAGCAGGTGAAGCCGTCGGAGCGTGGTGAACTGGAGATTACCTCCATCAACCAGATGTACCTCGAGGCGGGCAACCTGACCGTTGAACTGCTCGGgcgc                                   >  2:2356811/1‑148 (MQ=255)
                                                                                                                    ctgTGACCGGGCTTTATTTCTACGACAGTAAAGTCGTGGAGTACGCAAAGCAGGTGAAGCCGTCGGAGCGTGGTGAACTGGAGATTACCTCCATCAACCAGATGTACCTCGAGGCGGGCAACCTGACCGTTGAACTGCTCGGGCGCgg                                 <  2:2347148/146‑1 (MQ=255)
                                                                                                                      gTGACCGGGCTTTATTTCTACGACAGTAAAGTCGTGGAGTACGCAAAGCAGGTGAAGCCGTCGGAGCGTGGTGAACTGGAGATTACCTCCATCAACCAGATGTACCTCGAGGCGGGCAACCTGACCGTTGAACTGCTCGGGCGCGGAt                               >  1:1293355/1‑148 (MQ=255)
                                                                                                                       tGACCGGGCTTTATTTCTACGACAGTAACGTCGTGGAGTACGCAAAGCAGGTGAAGCCGTCGGAGCGTGGTGAACTGGAGATTACCTCCATCAACCAGATGTACCTCGAGGCGGGCAACCTGACCGTTGAACTGCTCGGGCGCGGAtt                              <  1:3453305/148‑1 (MQ=255)
                                                                                                                            gggCTTTATTTCTACGACAGTAAAGTCGTGGAGTACGCAAAGCAGGTGAAGCCGTCGGAGCGTGGTGAACTGGAGATTACCTCCATCAACCAGATGTACCTCGAGGCGGGCAACCTGACCGTTGAACTGCTCGGGCGCGGATTTGCct                         <  2:1732445/148‑1 (MQ=255)
                                                                                                                             ggCTTTATTTCTACGACAGTAAAGTCGTGGAGTACGCAAAGCAGGTGAGGCCGTCGGAGCGTGGTGAACTGGAGATTACCTCCATCAACCAGATGTACCTCGAGGCGGGCAACCTGACCGTTGAACTGCTCGGGCGCGGATTTGCctg                        <  1:2361686/148‑1 (MQ=255)
                                                                                                                             ggCTTTATTTCTACGACAGTAAAGTCGTGGAGTACGCAAAGCAGGTGAAGCCGTCGGAGCGTGGTGAACTGGAGATTACCTCCATCAACCAGATGTACCTCGAGGCGGGCAACCTGACCGTTGAACTGCTCGGGCGCGGATTTGCctg                        >  1:1525738/1‑148 (MQ=255)
                                                                                                                              gCTTTATTTCTACGACAGTAAAGTCGTGGAGTACGCAAAGCAGGTGAAGCCGTCGGAGCGTGGTGAACTGGAGATTACCTCCATCAACCAGATGTACCTCGAGGCGGGCAACCTGACCGTTGAACTGCTCGGGCGCGGATTTGCctgg                       >  1:2339790/1‑148 (MQ=255)
                                                                                                                              gCTTTATTTCTACGACAGTAAAGTCGTGGAGTACGCAAAGCAGGTGAAGCCGTCGGAGCGTGGTGAACTGGAGATTACCTCCATCAACCAGATGTACCTCGAGGCGGGCAACCTGACCGTTGAACTGCTCGGGCGCGGATTTGCctgg                       <  2:320701/148‑1 (MQ=255)
                                                                                                                                tttATTTCTACGACAGTAAAGTCGTGGAGTACGCAAAGCAGGTGAAGCCGTCGGAGCGTGGTGAACTGGAGATTACCTCCATCAACCAGATGTACCTCGAGGCGGGCAACCTGACCGTTGAACTGCTCGGGCGCGGATTTGCctggct                     >  2:1104119/1‑148 (MQ=255)
                                                                                                                                  tATTTCTACGACAGTAAAGTCGTGGAGTACGCAAAGCAGGTGAAGCCGTCGGAGCGTGGTGAACTGGAGAGTACCTCCATCAACCAGATGTACCTCGAGGCGGGCAACCTGACCGTTGAACTGCTCGGGCGCGGATTTGCctggctgg                   <  2:653199/148‑1 (MQ=255)
                                                                                                                                   aTTTCTACGACAGTAAAGTCGTGGAGTACGCAAAGCAGGTGAAGCCGTCGGAGCGTGGTGAACTGGAGATTACCTCCATCAACCAGATGTACCTCGAGGCGGGCAACCTGACCGTTGAACTGCTCGGGCGCGGATTTGCCTGGCTGGa                  >  1:2791012/1‑148 (MQ=255)
                                                                                                                                     ttCTACGACAGTAAAGTCGTGGAGTACGCAAAGCAGGTGAAGCCGTCGGAGCGTGGTGAACTGGAGATTACCTCCATCAACCAGATGTACCTCGAGGCGGGCAACCTGACCGTTGAACTGCTCGGGCGCGGATTTGCCTGGCTGGATa                >  2:800799/1‑148 (MQ=255)
                                                                                                                                     ttCTACGACAGTAAAGTCGTGGAGTACGCAAAGCAGGTGAAGCCGTCGGAGCGTGGTGAACTGGAGATTACCTCCATCAACCAGATGTACCTCGAGGCGGGCAACCTGACCGTTGAACTGCTCGGGCGCGGATTTGCCTGGCTGGATa                >  1:237037/1‑148 (MQ=255)
                                                                                                                                         aCGACAGTAAAGTCGTGGAGTACGCAAAGCAGGTGAAGCCGTCGGAGCGTGGTGAACTGGAGATTACCTCCATCAACCAGATGTACCTCGAGGCGGGCAACCTGACCGTTGAACTGCTCGGGCGCGGATTTGCCTGGCTGGATACCgg            >  1:231884/1‑148 (MQ=255)
                                                                                                                                         aCGACAGTAAAGTCGTGGAGTACGCAAAGCAGGTGAAGCCGTCGGAGCGTGGTGAACTGGAGATTACCTCCATCAACCAGATGTACCTCGAGGCGGGCAACCTGACCGTTGAACTGCTCGGGCGCGGATTTGCCTGGCTGGATACCgg            >  1:1451129/1‑148 (MQ=255)
                                                                                                                                           gACAGTAACGTCGTGGAGTACGCAAAGCAGGTGAAGCCGTCGGAGCGTGGTGAACTGGAGATTACCTCCATCAACCAGATGTACCTCGAGGCGGGCAACCTGACCGTTGAACTGCTCGGGCGCGGATTTGCCTGGCTGGATACCGGCa          <  2:1396128/148‑1 (MQ=255)
                                                                                                                                           gACAGTAACGTCGTGGAGTACGCAAAGCAGGTGAAGCCGTCGGAGCGTGGTGAACTGGAGATTACCTCCATCAACCAGATGTACCTCGAGGCGGGCAACCTGACCGTTGAACTGCTCGGGCGCGGATTTGCCTGGCTGGATACCGGCa          <  1:3323538/148‑1 (MQ=255)
                                                                                                                                           gACAGTAAAGTCGTGGAGTACGCAAAGCAGGTGAAGCCGTCGGAGCGTGGTGAACTGGAGATTACCTCCATCAACCAGATGTACCTCGAGGCGGGCAACCTGACCGTTGAACTGCTCGGGCGCGGATTTGCCTGGCTGGATACCGGCa          <  2:2960212/148‑1 (MQ=255)
                                                                                                                                            aCAGTAAAGTCGTGGAGTACGCAAAGCAGGTGAAGCCGTCGGAGCGTGGTGAACTGGAGATTACCTCCATCAACCAGATGTACCTCGAGGCGGGCAACCTGACCGTTGAACTGCTCGGGCGCGGATTTGCCTGGCTGGATACCGGCAc         >  1:2255721/1‑148 (MQ=255)
                                                                                                                                             cAGTAAAGTCGTGGAGTACGCAAAGCAGGTGAAGCCGTCGGAGCGTGGTGAACTGGAGATTACCTCCATCAACCAGATGTACCTCGAGGCGGGCAACCTGACCGTTGAACTGCTCGGGCGCGGATTTGCCTGGCTGGATACCGGCACt        >  1:2191902/1‑148 (MQ=255)
                                                                                                                                              aGTAACGTCGTGGAGTACGCAAAGCAGGTGAAGCCGTCGGAGCGTGGTGAACTGGAGATTACCTCCATCAACCAGATGTACCTCGAGGCGGGCAACCTGACCGTTGAACTGCTCGGGCGCGGATTTGCCTGGCTGGATACCGGCACTc       <  2:2720559/148‑1 (MQ=255)
                                                                                                                                              aGTAAAGTCGTGGAGTACGCAAAGCAGGTGAAGCCGTCGGAGCGTGGTGAACTGGAGATTACCTCCATCAACCAGATGTACCTCGAGGCGGGCAACCTGACCGTTGAACTGCTCGGGCGCGGATTTGCCTGGCTGGATACCGGCACTc       <  1:2023578/148‑1 (MQ=255)
                                                                                                                                               gTAAAGTCGTGGAGTACGCAAAGCAGGTGAAGCCGTCGGAGCGTGGTGAACTGGAGATTACCTCCATCAACCAGATGTACCTCGAGGCGGGCAACCTGACCGTTGAACTGCTCGGGCGCGGATTTGCCTGGCTGGATACCGGCACTCa      <  2:3388849/148‑1 (MQ=255)
                                                                                                                                               gTAAAGTCGTGGAGTACGCAAAGCAGGTGAAGCCGTCGGAGCGTGGTGAACTGGAGATTACCTCCATCAACCAGATGTACCTCGAGGCGGGCAACCTGACCGTTGAACTGCTCGGGCGCGGATTTGCCTGGCTGGATACCGGCACTCa      >  1:247771/1‑148 (MQ=255)
                                                                                                                                                tAAAGTCGTGGAGTACGCAAAGCAGGTGAAGCCGTCGGAGCGTGGTGAACTGGAGATTACCTCCATCAACCAGATGTACCTCGAGGCGGGCAACCTGACCGTTGAACTGCTCGGGCGCGGATTTGCCTGGCTGGATACCGGCACTCAc     >  2:1169512/1‑148 (MQ=255)
                                                                                                                                                 aaaGTCGTGGAGTACGCAAAGCAGGTGAAGCCGTCGGAGCGTGGTGAACTGGAGATTACCTCCATCAACCAGATGTACCTCGAGGCGGGCAACCTGACCGTTGAACTGCTCGGGCGCGGATTTGCCTGGCTGGATACCGGCACTCACg    >  2:1140576/1‑148 (MQ=255)
                                                                                                                                                  aaGTCGTGGAGTACGCAAAGCAGGTGAAGCCGTCGGAGCGTGGTGAACTGGAGATTACCTCCATCAACCAGATGTACCTCGAGGCGGGCAACCTGACCGTTGAACTGCTCGGGCGCGGATTTGCCTGGCTGGATACCGGCACTCACGa   >  2:1040110/1‑148 (MQ=255)
                                                                                                                                                   aGTCGTGGAGTACGCAAAGCAGGTGAAGCCGTCGGAGCGTGGTGAACTGGAGATTACCTCCATCAACCAGATGTACCTCGAGGCGGGCAACCTGACCGTTGAACTGCTCGGGCGCGGATTTGCCTGGCTGGATACCGGCACTCACGAc  >  1:2569835/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GACGGTTTTTGGCTATCAGGTGATGGACCCGGAACGCTTTGGCGTGGTGGAGTTTGACGACAATTTCCGCGCTATCTCGCTGGAAGAAAAGCCAAAACAGCCGAAGTCAAACTGGGCGGTGACCGGGCTTTATTTCTACGACAGTAAAGTCGTGGAGTACGCAAAGCAGGTGAAGCCGTCGGAGCGTGGTGAACTGGAGATTACCTCCATCAACCAGATGTACCTCGAGGCGGGCAACCTGACCGTTGAACTGCTCGGGCGCGGATTTGCCTGGCTGGATACCGGCACTCACGAC  >  CP009273/3967357‑3967651

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: