Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 709,827 G→T 12.1% D272Y (GAT→TAT)  pgm → phosphoglucomutase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273709,8270GT12.1% 1014.1 / 72.7 286D272Y (GAT→TAT) pgmphosphoglucomutase
Reads supporting (aligned to +/- strand):  ref base G (122/129);  new base T (20/15);  total (142/144)
Fisher's exact test for biased strand distribution p-value = 3.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 287 total aligned reads displayed.

GCGATTCAGAAAGCGGGCCTGACGCTGGGCGTTGATCCGCTGGGCGGTTCCGGTATCGAATACTGGAAGCGTATTGGCGAGTATTACAACCTCAACCTGACTATCGTTAACGATCAGGTCGATCAAACCTTCCGCTTTATGCACCTTGATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAGTGTGCGATGGCGGGCCTGCTGGCACTGCGTGATAAGTTCGATCTGGCGTTTGCTAACGACCCGGATTATGACCGTCACGGTATCGTCACTCCGGCAGGTTTGATGAATCCG  >  CP009273/709680‑709976
                                                                                                                                                   |                                                                                                                                                     
  gATTCAGAAAGCGGGCCTGACGCTGGGCGTTGATCCGCTGGGCGGTTCCGGTATCGAATACTGGAAGCGTATTGGCGAGTATTACAACCTCAACCTGACTATCGTTAACGATCAGGTCGATCAAACCTTCCGCTTTATGCACCTTGATaa                                                                                                                                                   >  2:960258/1‑150 (MQ=255)
    ttCAGAAAGCGGGCCTGACGCTGGGCGTTGATCCGCTGGGCGGTTCCGGTATCGAATACTGGAAGCGTATTGGCGAGTATTACAACCTCAACCTGACTATCGTTAACGATCAGGTCGATCAAACCTTCCGCTTTATGCACCTTGATAAAg                                                                                                                                                 <  1:1075179/150‑1 (MQ=255)
      cAGAAAGCGGGCCTGACGCTGGGCGTTGATCCGCTGGGCGGTTCCGGTATCGAATACTGGAAGCGTATTGGCGAGTATTACAACCTCAACCTGACTATCGTTAACGATCAGGTCGATCAAACCTTCCGCTTTATGCACCTTGATAAAGAc                                                                                                                                               >  2:3015939/1‑150 (MQ=255)
      cAGAAAGCGGGCCTGACGCTGGGCGTTGATCCGCTGGGCGGTTCCGGTATCGAATACTGGAAGCGTATTGGCGAGTATTACAACCTCAACCTGACTATCGTTAACGATCAGGTCGATCAAACCTTCCGCTTTATGCACCTTGATAAAGAc                                                                                                                                               <  2:3648683/150‑1 (MQ=255)
       aGAAAGCGGGCCTGACGCTGGGCGTTGATCCGCTGGGCGGTTCCGGTATCGAATACTGGAAGCGTATTGGCGAGTATTACAACCTCAACCTGACTATCGTTAACGATCAGGTCGATCAAACCTTCCGCTTTATGCACCTTGATAAAGACg                                                                                                                                              <  2:2503422/150‑1 (MQ=255)
       aGAAAGCGGGCCTGACGCTGGGCGTTGATCCGCTGGGCGGTTCCGGTATCGAATACTGGAAGCGTATTGGCGAGTATTACAACCTCAACCTGACTATCGTTAACGATCAGGTCGATCAAACCTTCCGCTTTATGCACCTTGATAAAGACg                                                                                                                                              <  2:2797270/150‑1 (MQ=255)
       aGAAAGCGGGCCTGACGCTGGGCGTTGATCCGCTGGGCGGTTCCGGTATCGAATACTGGAAGCGTATTGGCGAGTATTACAACCTCAACCTGACTATCGTTAACGATCAGGTCGATCAAACCTTCCGCTTTATGCACCTTGATAAAGACg                                                                                                                                              >  1:25823/1‑150 (MQ=255)
         aaaGCGGGCCTGACGCTGGGCGTTGATCCGCTGGGCGGTTCCGGTATCGAATACTGGAAGCGTATTGGCGAGTATTACAACCTCAACCTGACTATCGTTAACGATCAGGTCGATCAAACCTTCCGCTTTATGCACCTTGATAAAGACGgc                                                                                                                                            <  2:284639/150‑1 (MQ=255)
           aGCGGGCCTGACGCTGGGCGTTGATCCGCTGGGCGGTTCCGGTATCGAATACTGGAAGCGTATTGGCGAGTATTACAACCTCAACCTGACTATCGTTAACGATCAGGTCGATCAAACCTTCCGCTTTATGCACCTTGATAAAGACGgcgc                                                                                                                                          >  1:3047768/1‑150 (MQ=255)
            gCGGGCCTGACGCTGGGCGTTGATCCGCTGGGCGGTTCCGGTATCGAATACTGGAAGCGTATTGGCGAGTATTACAACCTCAACCTGACTATCGTTAACGATCAGGTCGATCAAACCTTCCGCTTTATGCACCTTTATAAAGACGGcgcg                                                                                                                                         <  2:2623429/150‑1 (MQ=255)
            gCGGGCCTGACGCTGGGCGTTGATCCGCTGGGCGGTTCCGGTATCGAATACTGGAAGCGTATTGGCGAGTATTACAACCTCAACCTGACTATCGTTAACGATCAGGTCGATCAAACCTTCCGCTTTATGCACCTTGATAAAGACGGcgcg                                                                                                                                         <  1:2447485/150‑1 (MQ=255)
             cGGGCCTGACGCTGGGCGTTGATCCGCTGGGCGGTTCCGGTATCGAATACTGGAAGCGTATTGGCGAGTATTACAACCTCAACCTGACTATCGTTAACGATCAGGTCGATCAAACCTTCCGCTTTATGCACCTTGATAAAGACGGCGCGa                                                                                                                                        >  1:9259/1‑150 (MQ=255)
              gggCCTGACGCTGGGCGTTGATCCGCTGGGCGGTTCCGGTATCGAATACTGGAAGCGTATTGGCGAGTATTACAACCTCAACCTGACTATCGTTAACGATCAGGTCGATCAAACCTTCCGCTTTATGCACCTTGATAAAGACGGCGCGAt                                                                                                                                       >  2:2197064/1‑150 (MQ=255)
               ggCCTGACGCTGGGCGTTGATCCGCTGGGCGGTTCCGGTATCGAATACTGGAAGCGTATTGGCGAGTATTACAACCTCAACCTGACTATCGTTAACGATCAGGTCGATCAAACCTTCCGCTTTATGCACCTTTATAAAGACGGCGCGATc                                                                                                                                      <  2:816781/150‑1 (MQ=255)
                  cTGACGCTGGGCGTTGATCCGCTGGGCGGTTCCGGTATCGAATACTGGAAGCGTATTGGCGAGTATTACAACCTCAACCTGACTATCGTTAACGATCAGGTCGATCAAACCTTCCGCTTTATGCACCTTGATAAAGACGGCGCGATCCGt                                                                                                                                   <  1:667208/150‑1 (MQ=255)
                   tGACGCTGGGCGTTGATCCGCTGGGCGGTTCCGGTATCGAATACTGGAAGCGTATTGGCGAGTATTACAACCTCAACCTGACTATCGTTAACGATCAGGTCGATCAAACCTTCCGCTTTATGCACCTTGATAAAGACGGCGCGATCCGTa                                                                                                                                  >  2:2559926/1‑150 (MQ=255)
                        cTGGGCGTTGATCCGCTGGGCGGTTCCGGTATCGAATACTGGAAGCGTATTGGCGAGTATTACAACCTCAACCTGACTATCGTTAACGATCAGGTCGATCAAACCTTCCGCTTTATGCACCTTGATAAAGACGGCGCGATCCGTATGGAc                                                                                                                             >  1:3943549/1‑150 (MQ=255)
                         tGGGCGTTGATCCGCTGGGCGGTTCCGGTATCGAATACTGGAAGCGTATTGGCGAGTATTACAACCTCAACCTGACTATCGTTAACGATCAGGTCGATCAAACCTTCCGCTTTATGCACCTTGATAAAGACGGCGCGATCCGTATGGACt                                                                                                                            >  2:1373315/1‑150 (MQ=255)
                         tGGGCGTTGATCCGCTGGGCGGTTCCGGTATCGAATACTGGAAGCGTATTGGCGAGTATTACAACCTCAACCTGACTATCGTTAACGATCAGGTCGATCAAACCTTCCGCTTTATGCACCTTGATAAAGACGGCGCGATCCGTATGGACt                                                                                                                            >  2:4192092/1‑150 (MQ=255)
                                 gATCCGCTGGGCGGTTCCGGTATCGAATACTGGAAGCGTATTGGCGAGTATTACAACCTCAACCTGACTATCGTTAACGATCAGGTCGATCAAACCTTCCGCTTTATGCACCTTGATAAAGACGGCGCGATCCGTATGGACTGCtcctcc                                                                                                                    <  2:1812624/150‑1 (MQ=255)
                                 gATCCGCTGGGCGGTTCCGGTATCGAATACTGGAAGCGTATTGGCGAGTATTACAACCTCAACCTGACTATCGTTAACGATCAGGTCGATCAAACCTTCCGCTTTATGCACCTTGATAAAGACGGCGCGATCCGTATGGACTGCtcctcc                                                                                                                    >  2:34952/1‑150 (MQ=255)
                                 gATCCGCTGGGCGGTTCCGGTATCGAATACTGGAAGCGTATTGGCGAGTATTACAACCTCAACCTGACTATCGTTAACGATCAGGTCGATCAAACCTTCCGCTTTATGCACCTTGATAAAGACGGCGCGATCCGTATGGACTGCtcctcc                                                                                                                    >  1:210759/1‑150 (MQ=255)
                                  aTCCGCTGGGCGGTTCCGGTATCGAATACTGGAAGCGTATTGGCGAGTATTACAACCTCAACCTGACTATCGTTAACGATCAGGTCGATCAAACCTTCCGCTTTATGCACCTTGATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCg                                                                                                                   <  2:4339453/150‑1 (MQ=255)
                                   tCCGCTGGGCGGTTCCGGTATCGAATACTGGAAGCGTATTGGCGAGTATTACAACCTCAACCTGACTATCGTTAACGATCAGGTCGATCAAACCTTCCGCTTTATGCACCTTGATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGa                                                                                                                  >  2:89057/1‑150 (MQ=255)
                                   tCCGCTGGGCGGTTCCGGTATCGAATACTGGAAGCGTATTGGCGAGTATTACAACCTCAACCTGACTATCGTTAACGATCAGGTCGATCAAACCTTCCGCTTTATGCAACTTGATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGa                                                                                                                  >  2:2402570/1‑150 (MQ=255)
                                          ggCGGTTCCGGTATCGAATACTGGAAGCGTATTGGCGAGTATTACAACCTCAACCTGACTATCGTTAACGATCAGGTCGATCAAACCTTCCGCTTTATGCACCTTGATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAGTGTGCg                                                                                                           <  2:185284/150‑1 (MQ=255)
                                                  cGGTATCGAATACTGGAAGCGTATTGGCGAGTATTACAACCTCAACCTGACTATCGTTAACGATCAGGTCGATCAAACCTTCCGCTTTATGCACCTTGATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAGTGTGCGATGGCggg                                                                                                   >  1:4463197/1‑150 (MQ=255)
                                                   ggTATCGAATACTGGAAGCGTATTGGCGAGTATTACAACCTCAACCTGACTATCGTTAACGATCAGGTCGATCAAACCTTCCGCTTTATGCACCTTGATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAGTGTGCGATGGCGGGc                                                                                                  >  1:3468924/1‑150 (MQ=255)
                                                     tATCGAATACTGGAAGCGTATTGGCGAGTATTACAACCTCAACCTGACTATCGTTAACGATCAGGTCGATCAAACCTTCCGCTTTATGCACCTTGATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAGTGTGCGATGGCGGGCct                                                                                                <  2:3053068/150‑1 (MQ=255)
                                                      aTCGAATACTGGAAGCGTATTGGCGAGTATTACAACCTCAACCTGACTATCGTTAACGATCAGGTCGATCAAACCTTCCGCTTTATGCACCTTGATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAGTGTGCGATGGCGGGCctg                                                                                               <  2:1612800/150‑1 (MQ=255)
                                                        cGAATTCTGGAAGCGTATTGGCGAGTATTACAACCTCAACCTGACTTTCGTTAACGATCAGGTCGATCAAACCTTCCGCTTTATGCACCTTGATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAGTGTGCGATGGCGGGCctgct                                                                                             <  1:1981659/150‑1 (MQ=255)
                                                        cGAATACTGGAAGCGTATTGGCGAGTATTACAACCTCAACCTGACTATCGTTAACGATCAGGTCGATCAAACCTTCCGCTTTATGCACCTTTATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAGTGTGCGATGGCGGGCctgct                                                                                             <  1:3608390/150‑1 (MQ=255)
                                                        cGAATACTGGAAGCGTATTGGCGAGTATTACAACCTCAACCTGACTATCGTTAACGATCAGGTCGATCAAACCTTCCGCTTTATGCACCTTGATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAGTGTGCGATGGCGGGCctgct                                                                                             <  1:2973096/150‑1 (MQ=255)
                                                         gAATACTGGAAGCGTATTGGCGAGTATTACAACCTCAACCCGACTATCGTTAACGATCAGGTCGATCAAACCTTCCGCTTTATGCACCTTGATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAGTGTGCGATGGCGGGCctgctg                                                                                            >  2:1214384/1‑150 (MQ=255)
                                                              cTGGAAGCGTATTGGCGAGTATTACAACCTCAACCTGACTATCGTTAACGATCAGGTCGATCAAACCTTCCGCTTTATGCACCTTGATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAGTGTGCGATGGCGGGCCTGCTGGCACt                                                                                       <  2:2528303/150‑1 (MQ=255)
                                                               tGGAAGCGTATTGGCGAGTATTACAACCTCAACCTGACTATCGTTAACGATCAGGTCGATCAAACCTTCCGCTTTATGCACCTTGATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAGTGTGCGATGGCGGGCCTGCTGGCACTg                                                                                      <  2:682705/150‑1 (MQ=255)
                                                               tGGAAGCGTATTGGCGAGTATTACAACCTCAACCTGACTATCGTTAACGATCAGGTCGATCAAACCGTCCGCTTTATGCACCTTTATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAGTGTGCGATGGCGGGCCTGCTGGCACTg                                                                                      >  2:1999038/1‑150 (MQ=255)
                                                                ggAAGCGTATTGGCGAGTATTACAACCTCAACCTGACTATCGTTAACGATCAGGTCGATCAAACCTTCCGCTTTATGCACCTTGATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAGTGTGCGATGGCGGGCCTGCTGGCACTGc                                                                                     <  1:3880502/150‑1 (MQ=255)
                                                                ggAAGCGTATTGGCGAGTATTACAACCTCAACCTGACTATCGTTAACGATCAGGTCGATCAAACCTTCCGCTTTATGCACCTTGATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAGTGTGCGATGGCGGGCCTGCTGGCACTGc                                                                                     <  1:4439893/150‑1 (MQ=255)
                                                                 gAAGCGTATTGGCGAGTATTACAACCTCAACCTGACTATCGTTAACGATCAGGTCGATCAAACCTTCCGCTTTATGCACCTTGATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAGTGTGCGATGGCGGGCCTGCTGGCACTGCg                                                                                    <  1:3291796/150‑1 (MQ=255)
                                                                 gAAGCGTATTGGCGAGTATTACAACCTCAACCTGACTATCGTTAACGATCAGGTCGATCAAACCTTCCGCTTTATGCACCTTGATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAGTGTGCGATGGCGGGCCTGCTGGCACTGCg                                                                                    >  2:1389040/1‑150 (MQ=255)
                                                                  aaGCGTATTGGCGAGTATTACAACCTCAACCTGACTATCGTTAACGATCAGGTCGATCAAACCTTCCGCTTTATGCACCTTGATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAGTGTGCGATGGCGGGCCTGCTGGCACTGCGt                                                                                   <  1:825405/150‑1 (MQ=255)
                                                                    gcatATTGGCGAGTATTACAACCTCAACCTGACTATCGTTAACGATCAGGTCGATCAAACCTTCCGCTTTATGCACCTTGATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAGTGTGCGATGGCGGGCCTGCTGGCACTGCGTGa                                                                                 <  1:2518215/147‑1 (MQ=255)
                                                                      gTATTGGCGAGTATTACAACCTCAACCTGACTATCGTTAACGATCAGGTCGATCAAACCTTCCGCTTTATGCACCTTGATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAGTGTGCGATGGCGGGCCTGCTGGCACTGCGTGATa                                                                               >  2:3156539/1‑150 (MQ=255)
                                                                           ggCGAGTATTACAACCTCAACCTGACTATCGTTAACGATCAGGTCGATCAAACCTTCCGCTTTATGCACCTTGATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAGTGTGCGATGGCGGGCCTGCTGGCACTGCGTGATAAGTTc                                                                          <  1:1830271/150‑1 (MQ=255)
                                                                             cGAGTATTACAACCTCAACCTGACTATCGTTAACGATCAGGTCGATCAAACCTTCCGCTTTATGCACCTTGATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAGTGTGCGATGGCGGGCCTGCTGGCACTGCGTGATAAGTTCGa                                                                        <  2:3559255/150‑1 (MQ=255)
                                                                              gAGTATTACAACCTCAACCTGACTATCGTTAACGATCAGGTCGATCAAACCTTCCGCTTTATGCACCTTGATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAGTGTGCGATGGCGGGCCTGCTGGCACTGCGTGATAAGTTCGAt                                                                       >  1:3327869/1‑150 (MQ=255)
                                                                              gAGTATTACAACCTCAACCTGACTATCGTTAACGATCAGGTCGATCAAACCTTCCGCTTTATGCACCTTGATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAGTGTGCGATGGCGGGCCTGCTGGCACTGCGTGATAAGTTCGAt                                                                       >  2:1552807/1‑150 (MQ=255)
                                                                                 tATTACAACCTCAACCTGACTATCGTTAACGATCAGGTCGATCAAACCGTCCGCTTTATGCACCTTGATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAGTGTGCGATGGCGGGCCTGCTGGCACTGCGTGATAAGTTCGATCTg                                                                    <  1:2429750/150‑1 (MQ=255)
                                                                                  aTTACAACCTCAACCTGACTATCGTTAACGATCAGGTCGATCAAACCTTCCGCTTTATGCACCTTTATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAGTGTGCGATGGCGGGCCTGCTGGCACTGCGTGATAAGTTCGATCTgg                                                                   >  1:3367798/1‑150 (MQ=255)
                                                                                    tACAACCTCAACCTGACTATCGTTAACGATCAGGTCGATCAAACCTTCCGCTTTATGCACCTTGATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAGTGTGCGATGGCGGGCCTGCTGGCACTGCGTGATAAGTTCGATCTGGCg                                                                 <  2:2510560/150‑1 (MQ=255)
                                                                                     aCAACCTCAACCTGACTATCGTTAACGATCAGGTCGATCAAACCTTCCGCTTTATGCACCTTGATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAGTGTGCGATGGCGGGCCTGCTGGCACTGCGTGATAAGTTCGATCTGGCGt                                                                >  1:2910212/1‑150 (MQ=255)
                                                                                      caacctcaacctGACTATCGTTAACGATCAGGTCGATCAAACCTTCCGCTTTATGCACCTTGATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAGTGTGCGATGGCGGGCCTGCTGGCACTGCGTGATAAGTTCGATCTGGCGtt                                                               <  2:4508865/150‑1 (MQ=255)
                                                                                       aacctcaacctGACTATCGTTAACGATCAGGTCGATCAAACCTTCCGCTTTATGCACCTTGATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAGTGTGCGATGGCGGGCCTGCTGGCACTGCGTGATAAGTTCGATCTGGCGttt                                                              >  2:1900492/1‑150 (MQ=255)
                                                                                        acctcaacctGACTATCGTTAACGATCAGGTCGATCAAACCTTCCGCTTTATGCACCTTGATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAGTGTGCGATGGCGGGCCTGCTGGCACTGCGTGATAAGTTCGATCTGGCGTTTg                                                             >  2:589477/1‑150 (MQ=255)
                                                                                            caacctGACTATCGTTAACGATCAGGTCGATCAAACCTTCCGCTTTATGCACCTTGATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAGTGTGCGATGGCGGGCCTGCTGGCACTGCGTGATAAGTTCGATCTGGCGTTTGCTaa                                                         <  1:4343582/150‑1 (MQ=255)
                                                                                             aacctGACTATCGTTAACGATCAGGTCGATCAAACCTTCCGCTTTATGCACCTTTATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAGTGTGCGAGGGCGGGCCTGCTGGCACTGCGTGATAAGTTCGATCTGGCGTTTGCTAAc                                                        >  1:3437875/1‑150 (MQ=255)
                                                                                              acctGACTATCGTTAACGATCAGGTCGATCAAACCTTCCGCTTTATGCACCTTGATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAGTGTGCGATGGCGGGCCTGCTGGCACTGCGTGATAAGTTCGATCTGGCGTTTGCTAACg                                                       >  1:2945511/1‑150 (MQ=255)
                                                                                                ctGACTATCGTTAACGATCAGGTCGATCAAACCTTCCGCTTTATGCACCTTGATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAGTGTGCGATGGCGGGCCTGCTGGCACTGCGTGATAAGTTCGATCTGGCGTTTGCTAACGAc                                                     <  2:3120035/150‑1 (MQ=255)
                                                                                                ctGACTATCGTTAACGATCAGGTCGATCAAACCTTCCGCTTTATGCACCTTGATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAGTGTGCGATGGCGGGCCTGCTGGCACTGCGTGATAAGTTCGATCTGGCGTTTGCTAACGAc                                                     >  1:3102564/1‑150 (MQ=255)
                                                                                                ctGACTATCGTTAACGATCAGGTCGATCAAACCTTCCGCTTTATGCACCTTGATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAGTGTGCGATGGCGGGCCTGCTGGCACTGCGTGATAAGTTCGATCTGGCGTTTGCTAACGAc                                                     >  1:1484090/1‑150 (MQ=255)
                                                                                                 tGACTATCGTTAACGATCAGGTCGATCAAACCTTCCGCTTTATGCACCTTGATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAGTGTGCGATGGCGGGCCTGCTGGCACTGCGTGATAAGTTCGATCTGGCGTTTGCTAACGAcc                                                    >  1:1365848/1‑150 (MQ=255)
                                                                                                  gACTATCGTTAACGATCAGGTCGATCAAACCTTCCGCTTTATGCACCTTGATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAGTGTGCGATGGCGGGCCTGCTGGCACTGCGTGATAAGTTCGATCTGGCGTTTGCTAACGAccc                                                   >  1:2542862/1‑150 (MQ=255)
                                                                                                  gACTATCGTTAACGATCAGGTCGATCAAACCTTCCGCTTTATGCACCTTGATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAGTGTGCGATGGCGGGCCTGCTGGCACTGCGTGATAAGTTCGATCTGGCGTTTGCTAACGAccc                                                   <  2:2008986/150‑1 (MQ=255)
                                                                                                  gACTATCGTTAACGATCAGGTCGATCAAACCTTCCGCTTTATGCACCTTGATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAGTGTGCGATGGCGGGCCTGCTGGCACTGCGTGATAAGTTCGATCTGGCGTTTGCTAACGAccc                                                   <  2:2430510/150‑1 (MQ=255)
                                                                                                    cTATCGTTAACGATCAGGTCGATCAAACCTTCCGCTTTATGCACCTTGATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAGTGTGCGATGGCGGGCCTGCTGGCACTGCGTGATAAGTTCGATCTGGCGTTTGCTAACGACCCgg                                                 >  2:515977/1‑150 (MQ=255)
                                                                                                      atagTTAACTATCAGGTCGATCAACCCTTCCTCTTTATGCACCTTGATAAAGACGGCGCGATCCGTATGGACTGCTCCTCAGAGTGTGCGATGGCGGGCCTGCTGGCACTGCGTGATAAGTTCGATCTGGCGTTTGCTAACGACCCGGAt                                               <  2:1058947/147‑1 (MQ=255)
                                                                                                      aTCGTTAACGATCAGGTCGATCAAACCTTCCGCTTTATGCACCTTGATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAGTGTGCGATGGCGGGCCTGCTGGCACTGCGTGATAAGTTCGATCTGGCGTTTGCTAACGACCCGGAt                                               >  2:1189468/1‑150 (MQ=255)
                                                                                                        cGTTAACGATCAGGTCGATCAAACCTTCCGCTTTATGCACCTTGATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAGTGTGCGATGGCGGGCCTGCTGGCACTGCGTGATAAGTTCGATCTGGCGTTTGCTAACGACCCGGATTa                                             >  2:68740/1‑150 (MQ=255)
                                                                                                        cGTTAACGATCAGGTCGATCAAACCTTCCGCTTTATGCACCTTGATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAGTGTGCGATGGCGGGCCTGCTGGCACTGCGTGATAAGTTCGATCTGGCGTTTGCTAACGACCCGGATTa                                             >  1:2285102/1‑150 (MQ=255)
                                                                                                          ttAACGATCAGGTCGATCAAACCTTCCGCTTTATGCACCTTGATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAGTGTGCGATGGCGGGCCTGCTGGCACTGCGTGATAAGTTCGATCTGGCGTTTGCTAACGACCCGGATTATg                                           <  1:1142312/150‑1 (MQ=255)
                                                                                                           tAACGATCAGGTCGATCAAACCTTCCGCTTTATGCACCTTGATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAGTGTGCGATGGCGGGCCTGCTGGCACTGCGTGATAAGTTCGATCTGGCGTTTGCTAACGACCCGGATTATGa                                          <  1:3506010/150‑1 (MQ=255)
                                                                                                            aaCGATCAGGTCGATCAAACCTTCCGCTTTATGCACCTTGATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAGTGTGCGATGGCGGGCCTGCTGGCACTGCGTGATAAGTTCGATCTGGCGTTTGCTAACGACCCGGATTATGAc                                         >  1:516881/1‑150 (MQ=255)
                                                                                                             aCGATCAGGTCGATCAAACCTTCCGCTTTATGCACCTTGATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAGTGTGCGATGGCGGGCCTGCTGGCACTGCGTGATAAGTTCGATCTGGCGTTTGCTAACGACCCGGATTATGAcc                                        <  2:3436394/150‑1 (MQ=255)
                                                                                                              cGATCAGGTCGATCAAACCTTCCGCTTTATGCACCTTTATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAGTGTGCGATGGCGGGCCTGCTGGCACTGCGTGATAAGTTCGATCTGGCGTTTGCTAACGACCCGGATTATGACCg                                       >  2:1263089/1‑150 (MQ=255)
                                                                                                                aTCAGGTCGATCAAACCTTCCGCTTTATGCACCTTGATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAGTGTGCGATGGCGGGCCTGCTGGCACTGCGTGATAAGTTCGATCTGGCGTTTGCTAACGACCCGGATTATGACCGTc                                     >  1:2604971/1‑150 (MQ=255)
                                                                                                                 tCAGGTCGATCAAACCTTCCGCTTTATGCACCTTGATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAGTGTGCGATGGCGGGCCTGCTGGCACTGCGTGATAAGTTCGATCTGGCGTTTGCTAACGACCCGGATTATGACCGTCa                                    <  2:2680785/150‑1 (MQ=255)
                                                                                                                     gTCGATCAAACCTTCCGCTTTATGCACCTTGATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAGTGTGCGATGGCGGGCCTGCTGGCACTGCGTGATAAGTTCGATCTGGCGTTTGCTAACGACCCGGATTATGACCGTCACGGt                                <  1:984181/150‑1 (MQ=255)
                                                                                                                       cGATCAAACCTTCCGCTTTTTGCACCTTGATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAGTGTGCGATGGCGGGCCTGCTGGCACTGCGTGATAAGTTCGATCTGGCGTTTGCTAACGACCCGGATTATGACCGTCACGGTAt                              <  1:4345163/150‑1 (MQ=255)
                                                                                                                       cGATCAAACCTTCCGCTTTATGCACCTTTATAAAGATGGCGCGATCCGTATGGACTGCTCCTCCGAGTGTGCGATGGCGGGCCTGCTGGCACTGCGTGATAAGTTCGATCTGGCGTTTGCTAACGACCCGGATTATGACCGTCACGGTAt                              <  1:1873778/150‑1 (MQ=255)
                                                                                                                       cGATCAAACCTTCCGCTTTATGCACCTTGATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAGTGTGCGATGGCGGGCCTGCTGGCACTGCGTGATAAGTTCGATCTGGCGTTTGCTAACGACCCGGATTATGACCGTCACGGTAt                              <  2:1508828/150‑1 (MQ=255)
                                                                                                                       cGATCAAACCTTCCGCTTTATGCACCTTGATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAGTGTGCGATGGCGGGCCTGCTGGCACTGCGTGATAAGTTCGATCTGGCGTTTGCTAACGACCCGGATTATGACCGTCACGGTAt                              <  1:4344833/150‑1 (MQ=255)
                                                                                                                        gATCAAACCTTCCGCTTTATGCACCTTTATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAGTGTGCGATGGCGGGCCTGCTGGCACTGCGTGATAAGTTCGATCTGGCGTTTGCTAACGACCCGGATTATGACCGTCACGGTATc                             >  1:1306062/1‑150 (MQ=255)
                                                                                                                        gATCAAACCTTCCGCTTTATGCACCTTGATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAGTGTGCGATGGCGGGCCTGCTGGCACTGCGTGATAAGTTCGATCTGGCGTTTGCTAACGACCCGGATTATGACCGTCACGGTATc                             >  1:3036920/1‑150 (MQ=255)
                                                                                                                           cAAACCTTCCGCTTTATGCACCTTGATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAGTGTGCGATGGCGGGCCTGCTGGCACTGCGTGATAAGTTCGATCTGGCGTTTGCTAACGACCCGGATTATGACCGTCACGGTATCGTc                          >  1:4248564/1‑150 (MQ=255)
                                                                                                                            aaaCCTTCCGCTTTATGCACCTTGATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAGTGTGCGATGGCGGGCCTGCTGGCACTGCGTGATAAGTTCGATCTGGCGTTTGCTAACGACCCGGATTATGACCGTCACGGTATCGTCa                         >  2:1014237/1‑150 (MQ=255)
                                                                                                                            aaaCCTTCCGCTTTATGCACCTTGATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAGTGTGCGATGGCGGGCCTGCTGGCACTGCGTGATAAGTTCGATCTGGCGTTTGCTAACGACCCGGATTATGACCGTCACGGTATCGTCa                         >  2:1800757/1‑150 (MQ=255)
                                                                                                                             aaCCTTCCGCTTTATGCACCTTGATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAGTGTGCGATGGCGGGCCTGCTGGCACTGCGTGATAAGTTCGATCTGGCGTTTGCTAACGACCCGGATTATGACCGTCACGGTATCGTCAc                        <  2:3764244/150‑1 (MQ=255)
                                                                                                                             aaCCTTCCGCTTTATGCACCTTGATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAGTGTGCGATGGCGGGCCTGCTGGCACTGCGTGATAAGTTCGATCTGGCGTTTGCTAACGACCCGGATTATGACCGTCACGGTATCGTCAc                        >  2:3993094/1‑150 (MQ=255)
                                                                                                                             aaCCTTCCGCTTTATGCACCTTGATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAGTGTGCGATGGCGGGCCTGCTGGCACTGCGTGATAAGTTCGATCTGGCGTTTGCTAACGACCCGGATTATGACCGTCACGGTATCGTCAc                        >  1:1398291/1‑150 (MQ=255)
                                                                                                                                cTTCCGCTTTATGCACCTTTATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAGTGTGCGATGGCGGGCCTGCTGGCACTGCGTGATAAGTTCGATCTGGCGTTTGCTAACGACCCGGATTATGACCGTCACGGTATCGTCACTcc                     >  1:276831/1‑150 (MQ=255)
                                                                                                                                     gCTTTATGCACCTTGATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAGTGTGCGATGGCGGGCCTGCTGGCACTGCGTGATAAGTTCGATCTGGCGTTTGCTAACGACCCGGATTATGACCGTCACGGTATCGTCACTCCGGCAg                <  1:739009/150‑1 (MQ=255)
                                                                                                                                      cTTTATGCACCTTTATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAGTGTGCGATGGCGGGCCTGCTGGCACTGCGTGATAAGTTCGATCTGGCGTTTGCTAACGACCCGGATTATGACCGTCACGGTATCGTCACTCCGGCAgg               >  2:1144062/1‑150 (MQ=255)
                                                                                                                                          aTGCACCTTGATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAGTGTGCGATGGCGGGCCTGCTGGCACTGCGTGATAAGTTCGATCTGGCGTTTGCTAACGACCCGGATTATGACCGTCACGGTATCGTCACTCCGGCAGGTTtg           >  2:3205841/1‑150 (MQ=255)
                                                                                                                                             cACCTTGATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAGTGTGCGATGGCGGGCCTGCTGGCACTGCGTGATAAGTTCGATCTGGCGTTTGCTAACGACCCGGATTATGACCGTCACGGTATCGTCACTCCGGCAGGTTtgatg        <  1:461590/150‑1 (MQ=255)
                                                                                                                                             cACCTTGATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAGTGTGCGATGGCGGGCCTGCTGGCACTGCGTGATAAGTTCGATCTGGCGTTTGCTAACGACCCGGATTATGACCGTCACGGTATCGTCACTCCGGCAGGTTtgatg        >  1:1761246/1‑150 (MQ=255)
                                                                                                                                                cTTGATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAGTGTGCGATGGCGGGCCTGCTGGCACTGCGTGATAAGTTCGATCTGGCGTTTGCTAACGACCCGGATTATGACCGTCACGGTATCGTCACTCCGGCAGGTTTGATGAAt     <  1:2197056/150‑1 (MQ=255)
                                                                                                                                                 ttGATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAGTGTGCGATGGCGGGCCTGCTGGCACTGCGTGATAAGTTCGATCTGGCGTTTGCTAACGACCCGGATTATGACCGTCACGGTATCGTCACTCCGGCAGGTTTGATGAATc    >  2:4159490/1‑150 (MQ=255)
                                                                                                                                                 ttGATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAGTGTGCGATGGCGGGCCTGCTGGCACTGCGTGATAAGTTCGATCTGGCGTTTGCTAACGACCCGGATTATGACCGTCACGGTATCGTCACTCCGGCAGGTTTGATGAATc    <  1:4559730/150‑1 (MQ=255)
                                                                                                                                                   gATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAGTGTGCGATGGCGGGCCTGCTGGCACTGCGTGATAAGTTCGATCTGGCGTTTGCTAACGACCCGGATTATGACCGTCACGGTATCGTCACTCCGGCAGGTTTGATGAATCCg  >  1:4065017/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GCGATTCAGAAAGCGGGCCTGACGCTGGGCGTTGATCCGCTGGGCGGTTCCGGTATCGAATACTGGAAGCGTATTGGCGAGTATTACAACCTCAACCTGACTATCGTTAACGATCAGGTCGATCAAACCTTCCGCTTTATGCACCTTGATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAGTGTGCGATGGCGGGCCTGCTGGCACTGCGTGATAAGTTCGATCTGGCGTTTGCTAACGACCCGGATTATGACCGTCACGGTATCGTCACTCCGGCAGGTTTGATGAATCCG  >  CP009273/709680‑709976

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: