Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,340,038 C→A 5.2% I34I (ATC→ATA)  ptsN → sugar‑specific enzyme IIA component of PTS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,340,0380CA 5.2% 713.6 / 9.0 190I34I (ATC→ATA) ptsNsugar‑specific enzyme IIA component of PTS
Reads supporting (aligned to +/- strand):  ref base C (95/82);  new base A (8/3);  total (104/86)
Fisher's exact test for biased strand distribution p-value = 3.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 202 total aligned reads displayed.

CCGAGCCTGTTCCACTGTTTGAGTGGGCAGGTTCTTAGGTGAAATTATGACAAATAATGATACAACTCTACAGCTTAGCAGTGTTCTTAACAGGGAATGTACGCGAAGCCGCGTCCACTGTCAGAGCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTGACGCGGGAAAAAATGGGCAGTACCGGTATCGGCAATGGTATTGCCATTCCGCATGGCAAACTGGAAGAAGATACTCTGCGCGC  >  CP009273/3339891‑3340184
                                                                                                                                                   |                                                                                                                                                  
   aGCCTGTTCCACTGTTTGAGTGGGCAGGTTCTTAGGTGAAATTATGACAAATAATGATACAACTCTACAGCTTAGCAGTGTTCTTAACAGGGAATGTACGCGAAGCCGCGTCCACTGTCAGAGCAAAAAACGCGCCCTGGAAatcatc                                                                                                                                                 <  1:138806/148‑1 (MQ=255)
      cTGTTCCACTTTTTGAGTGGGCAGGTTCTTAGGTGAAATTATGACAAATAATTATACAACTCTACAGCTTAGCAGTGTTCTTAACAGGGTATGTACTCGAAGCCGCGTCCACTGTCAGAGCAAAAAACGCGCCCTGGAAATCATCAGc                                                                                                                                              <  1:3114608/148‑1 (MQ=255)
       tGTTCCACTGTTTGAGTGGGCAGGTTCTTAGGTGAAATTATGACAAATAATGATACAACTCTACAGCTTAGCAGTGTTCTTAACAGGGAATGTACGCGAAGCCGCGTCCACTGTCAGAGCAAAAAACGCGCCCTGGAAATCATCAGCg                                                                                                                                             <  2:1801240/148‑1 (MQ=255)
         ttCCACTGTTTGAGTGGGCAGGTTCTTAGGTGAAATTATGACAAATAATGATACAACTCTACAGCTTAGCAGTGTTCTTAACAGGGAATGTACGCGAAGCCGCGTCCACTGTCAGAGCAAAAAACGCGCCCTGGAAATCATCAGCGAg                                                                                                                                           <  2:502682/148‑1 (MQ=255)
         ttCCACTGTTTGAGTGGGCAGGTTCTTAGGTGAAATTATGACAAATAATGATACAACTCTACAGCTTAGCAGTGTTCTTAACAGGGAATGTACGCGAAGCCGCGTCCACTGTCAGAGCAAAAAACGCGCCCTGGAAATCATCAGCGAg                                                                                                                                           >  2:1719878/1‑148 (MQ=255)
           ccACTGTTTGAGTGGGCAGGTTCTTAGGTGAAATTATGACAAATAATGATACAACTCTACAGCTTAGCAGTGTTCTTAACAGGGAATGTACGCGAAGCCGCGTCCACTGTCAGAGCAAAAAACGCGCCCTGGAAATCATCAGCGAGCt                                                                                                                                         >  2:692377/1‑148 (MQ=255)
            cACTGTTTGAGTGGGCAGGTTCTTAGGTGAAATTATGACAAATAATGATACAACTCTACAGCTTAGCAGTGTTCTTAACAGGGAATGTACGCGAAGCCGCGTCCACTGTCAGAGCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTg                                                                                                                                        >  2:2737035/1‑148 (MQ=255)
            cACTGTTTGAGTGGGCAGGTTCTTAGGTGAAATTATGACAAATAATGATACAACTCTACAGCTTAGCAGTGTTCTTAACAGGGAATGTACGCGAAGCCGCGTCCACTGTCAGAGCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTg                                                                                                                                        >  2:588813/1‑148 (MQ=255)
             aCTGTTTGAGTGGGCAGGTTCTTAGGTGAAATTATGACAAATAATGATACAACTCTACAGCTTAGCAGTGTTCTTAACAGGGAATGTACGCGAAGCCGCGTCCACTGTCAGAGCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTgg                                                                                                                                       <  1:1157739/148‑1 (MQ=255)
             aCTGTTTGAGTGGGCAGGTTCTTAGGTGAAATTATGACAAATAATGATACAACTCTACAGCTTAGCAGTGTTCTTAACAGGGAATGTACGCGAAGCCGCGTCCACTGTCAGAGCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTgg                                                                                                                                       >  2:3474646/1‑148 (MQ=255)
             aCTGTTTGAGTGGGCAGGTTCTTAGGTGAAATTATGACAAATAATGATACAACTCTACAGCTTAGCAGTGTTCTTAACAGGGAATGTACGCGAAGCCGCGTCCACTGTCAGAGCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTgg                                                                                                                                       >  2:873780/1‑148 (MQ=255)
             aCTGTTTGAGTGGGCAGGTTCTTAGGTGAAATTATGACAAATAATGATACAACTCTACAGCTTAGCAGTGTTCTTAACAGGGAATGTACGCGAAGCCGCGTCCACTGTCAGAGCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTgg                                                                                                                                       >  1:3374861/1‑148 (MQ=255)
               tGTTTGAGTGGGCAGGTTCTTAGGTGAAATTATGACAAATAATGATACAACTCTACAGCTTAGCAGTGTTCTTAACAGGGAATGTACGCGAAGCCGCGTCCACTGTCAGAGCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTggcg                                                                                                                                     >  2:1153029/1‑148 (MQ=255)
                       tGGGCAGGTTCTTAGGTGAAATTATGACAAATAATGATACAACTCTACAGCTTAGCAGTGTTCTTAACAGGGAATGTACGCGAAGCCGCGTCCACTGTCAGAGCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAaaca                                                                                                                             <  1:1840166/148‑1 (MQ=255)
                          gCAGGTTCTTAGGTGAAATTATGACAAATAATGATACAACTCTACAGCTTAGCAGTGTTCTTAACAGGGAATGTACGCGAAGCCGCGTCCACTGTCAGAGCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACt                                                                                                                          >  1:858782/1‑148 (MQ=255)
                          gCAGGTTCTTAGGTGAAATTATGACAAATAATGATACAACTCTACAGCTTAGCAGTGTTCTTAACAGGGAATGTACGCGAAGCCGCGTCCACTGTCAGAGCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACt                                                                                                                          >  1:858631/1‑148 (MQ=255)
                             ggTTCTTAGGTGAAATTATGACAAATAATGATACAACTCTACAGCTTAGCAGTGTTCTTAACAGGGAATGTACGCGAAGCCGCGTCCACTGTCAGAGCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAg                                                                                                                       >  1:1405535/1‑148 (MQ=255)
                              gTTCTTAGGTGAAATTATGACAAATAATGATACAACTCTACAGCTTAGCAGTGTTCTTAACAGGGAATGTACGCGAAGCCGCGTCCACTGTCAGAGCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGc                                                                                                                      >  2:1836441/1‑148 (MQ=255)
                               ttCTTAGGTGAAATTATGACAAATAATGATACAACTCTACAGCTTAGCAGTGTTCTTAACAGGGAATGTACGCGAAGCCGCGTCCACTGTCAGAGCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGcc                                                                                                                     <  1:1536104/148‑1 (MQ=255)
                               ttCTTAGGTGAAATTATGACAAATAATGATACAACTCTACAGCTTAGCAGTGTTCTTAACAGGGAATGTACGCGAAGCCGCGTCCACTGTCAGAGCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGcc                                                                                                                     <  2:2335078/148‑1 (MQ=255)
                               ttCTTAGGTGAAATTATGACAAATAATGATACAACTCTACAGCTTAGCAGTGTTCTTAACAGGGAATGTACGCGAAGCCGCGTCCACTGTCAGAGCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGcc                                                                                                                     <  2:377162/148‑1 (MQ=255)
                               ttCTTAGGTGAAATTATGACAAATAATGATACAACTCTACAGCTTAGCAGTGTTCTTAACAGGGAATGTACGCGAAGCCGCGTCCACTGTCAGAGCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGcc                                                                                                                     <  2:48880/148‑1 (MQ=255)
                               ttCTTAGGTGAAATTATGACAAATAATGATACAACTCTACAGCTTAGCAGTGTTCTTAACAGGGAATGTACGCGAAGCCGCGTCCACTGTCAGAGCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGcc                                                                                                                     >  2:2879675/1‑148 (MQ=255)
                                tCTTAGGTGAAATTATGACAAATAATGATACAACTCTACAGCTTAGCAGTGTTCTTAACAGGGAATGTACGCGGAGCCGCGTCCACTGTCAGAGCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCt                                                                                                                    >  1:2642419/1‑148 (MQ=255)
                                     ggTGAAATTATGACAAATAATGATACAACTCTACAGCTTAGCAGTGTTCTTAACAGGGAATGTACGCGAAGCCGCGTCCACTGTCAGAGCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTgccgc                                                                                                               >  2:3151875/1‑148 (MQ=255)
                                      gTGAAATTATGACAAATAATGATACAACTCTACAGCTTAGCAGTGTTCTTAACAGGGAATGTACGCGAAGCCGCGTCCACTGTCAGAGCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTgccgcc                                                                                                              >  2:2623369/1‑148 (MQ=255)
                                       tGAAATTATGACAAATAATGATACAACTCTACAGCTTAGCAGTGTTCTTAACAGGGAATGTACGCGAAGCCGCGTCCACTGTCAGAGCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCt                                                                                                             <  1:2485106/148‑1 (MQ=255)
                                         aaaTTATGACAAATAATGATACAACTCTACAGCTTAGCAGTGTTCTTAACAGGGAATGTACGCGAAGCCGCGTCCACTGTCAGAGCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCa                                                                                                           <  2:1389461/148‑1 (MQ=255)
                                         aaaTTATGACAAATAATGATACAACTCTACAGCTTAGCAGTGTTCTTAACAGGGAATGTACGCGAAGCCGCGTCCACTGTCAGAGCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCa                                                                                                           <  2:2478671/148‑1 (MQ=255)
                                          aaTTATGACAAATAATGATACAACTCTACAGCTTAGCAGTGTTCTTAACAGGGAATGTACGCGAAGCCGCGTCCACTGTCAGAGCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAg                                                                                                          <  2:3250589/148‑1 (MQ=255)
                                          aaTTATGACAAATAATGATACAACTCTACAGCTTAGCAGTGTTCTTAACAGGGAATGTACGCGAAGCCGCGTCCACTGTCAGAGCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAg                                                                                                          >  1:3205675/1‑148 (MQ=255)
                                           aTTATGACAAATAATGATACAACTCTACAGCTTAGCAGTGTTCTTAACAGGGAATGTACGCGAAGCCGCGTCCACTGTCAGAGCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAgg                                                                                                         >  2:1788243/1‑148 (MQ=255)
                                            ttATGACAAATAATGATACAACTCTACAGCTTAGCAGTGTTCTTAACAGGGAATGTACGCGAAGCCGCGTCCACTGTCAGAGCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAggt                                                                                                        <  1:293650/148‑1 (MQ=255)
                                            ttATGACAAATAATGATACAACTCTACAGCTTAGCAGTGTTCTTAACAGGGAATGTACGCGAAGCCGCGTCCACTGTCAGAGCAAAAAACGCGCCCTGGAAATAATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAggt                                                                                                        <  2:2645915/148‑1 (MQ=255)
                                             tATGACAAATAATGATACAACTCTACAGCTTAGCAGTGTTCTTAACAGGGAATGTACGCGAAGCCGCGTCCACTGTCAGAGCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCTggtg                                                                                                       >  2:3314367/1‑148 (MQ=255)
                                             tATGACAAATAATGATACAACTCTACAGCTTAGCAGTGTTCTTAACAGGGAATGTACGCGAAGCCGCGTCCACTGTCAGAGCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAggtg                                                                                                       >  2:3329857/1‑148 (MQ=255)
                                              aTGACAAATAATGATACAACTCTACAGCTTAGCAGTGTTCTTAACAGGGAATGTACGCGAAGCCGCGTCCACTGTCAGAGCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAggtgg                                                                                                      <  1:2362702/148‑1 (MQ=255)
                                                    aataatGATACAACTCTACAGCTTAGCAGTGTTCTTAACAGGGAATGTACGCGAAGCCGCGTCCACTGTCAGAGCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTg                                                                                                >  2:892986/1‑148 (MQ=255)
                                                      taatGATACAACTCTACAGCTTAGCAGTGTTCTTAACAGGGAATGTACGCGAAGCCGCGTCCACTGTCAGAGCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGaa                                                                                              >  1:1637290/1‑148 (MQ=255)
                                                       aatGATACAACTCTACAGCTTAGCAGTGTTCTTAACAGGGAATGTACGCGAAGCCGCGTCCACTGTCAGAGCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAg                                                                                             <  1:3096267/148‑1 (MQ=255)
                                                       aatGATACAACTCTACAGCTTAGCAGTGTTCTTAACAGGGAATGTACGCGAAGCCGCGTCCACTGTCAGAGCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAg                                                                                             >  2:1553250/1‑148 (MQ=255)
                                                       aatGATACAACTCTACAGCTTAGCAGTGTTCTTAACAGGGAATGTACGCGAAGCCGCGTCCACTGTCAGAGCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAg                                                                                             >  1:2574726/1‑148 (MQ=255)
                                                        atGATACAACTCTACAGCTTAGCAGTGTTCTTAACAGGGAATGTACGCGAAGCCGCGTCCACTGTCAGAGCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGc                                                                                            <  2:2790164/148‑1 (MQ=255)
                                                           aTACAACTCTACAGCTTAGCAGTGTTCTTAACAGGGAATGTACGCGAAGCCGCGTCCACTGTCAGAGCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTAt                                                                                         >  1:431608/1‑148 (MQ=255)
                                                           aTACAACTCTACAGCTTAGCAGTGTTCTTAACAGGGAATGTACGCGAAGCCGCGTCCACTGTCAGAGCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTAt                                                                                         >  2:51028/1‑148 (MQ=255)
                                                             aCAACTCTACAGCTTAGCAGTGTTCTTAACAGGGAATGTACGCGAAGCCGCGTCCACTGTCAGAGCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATcc                                                                                       >  1:2603124/1‑148 (MQ=255)
                                                             aCAACTCTACAGCTTAGCAGTGTTCTTAACAGGGAATGTACGCGAAGCCGCGTCCACTGTCAGAGCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATcc                                                                                       >  1:2300410/1‑148 (MQ=255)
                                                             aCAACTCTACAGCTTAGCAGTGTTCTTAACAGGGAATGTACGCGAAGCCGCGTCCACTGTCAGAGCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATcc                                                                                       >  1:73257/1‑148 (MQ=255)
                                                             aCAACTCTACAGCTTAGCAGTGTTCTTAACAGGGAATGTACGCGAAGCCGCGACCACTGTCAGAGCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATcc                                                                                       >  1:155720/1‑148 (MQ=255)
                                                               aaCTCTACAGCTTAGCAGTGTTCTTGACAGGGAATGTACGCGAAGCCGCGTCCACTGTCAGAGCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTg                                                                                     <  2:1910100/148‑1 (MQ=255)
                                                               aaCTCTACAGCTTAGCAGTGTTCTTAACAGGGAATGTACGCGAAGCCGCGTCCACTGTCAGAGCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTg                                                                                     >  1:895039/1‑148 (MQ=255)
                                                               aaCTCTACAGCTTAGCAGTGTTCTTAACAGGGAATGTACGCGAAGCCGCGTCCACTGTCAGAGCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTg                                                                                     <  2:1289563/148‑1 (MQ=255)
                                                                 ctctACAGCTTAGCAGTGTTCTTAACAGGGAATGTACGCGAAGCCGCGTCCACTGTCAGAGCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTGAc                                                                                   >  1:650900/1‑148 (MQ=255)
                                                                  tctACAGCTTAGCAGTGTTCTTAACAGGGAATGTACGCGAAGCCGCGTCCACTGTCAGAGCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTGAcg                                                                                  >  2:2089055/1‑148 (MQ=255)
                                                                     aCAGCTTAGCAGTGTTCTTAACAGGGAATGTACGCGAAGCCGCGTCCACTGTCAGAGCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTGACGCgg                                                                               <  1:2306844/148‑1 (MQ=255)
                                                                      cAGCTTAGCAGTGTTCTTAACAGGGAATGTACGCGAAGCCGCGTCCACTGTCAGAGCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTGACGCggg                                                                              <  2:3245528/148‑1 (MQ=255)
                                                                         cTTAGCAGTGTTCTTAACAGGGAATGTACGCGAAGCCGCGTCCACTGTCAGAGCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTGACGCGGGaaa                                                                           >  1:2208130/1‑148 (MQ=255)
                                                                         cTTAGCAGTGTTCTTAACAGGGAATGTACGCGAAGCCGCGTCCACTGTCAGAGCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTGACGCGGGaaa                                                                           >  2:2177163/1‑148 (MQ=255)
                                                                          ttAGCAGTGTTCTTAACAGGGAATGTACGCGAAGCCGCGTCCACTGTCAGAGCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTGACGCGGGaaaa                                                                          <  2:1492132/148‑1 (MQ=255)
                                                                            aGCAGTGTTCTTAACAGGGAATGTACGCGAAGCCGCGTCCACTGTCAGAGCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTGACGCGGGaaaaaa                                                                        <  2:23701/148‑1 (MQ=255)
                                                                               aGTGTTCTTAACAGGGAATGTACGCGAAGCCGCGTCCACTGTCAGAGCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTGACGCGGGAAAAAATgg                                                                     <  1:2684362/148‑1 (MQ=255)
                                                                                gtgtTCTTAACAGGGAATGTACGCGAAGCCGCGTCCACTGTCAGAGCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTGACGCGGGAAAAAATggg                                                                    >  1:2016678/1‑148 (MQ=255)
                                                                                 tgtTCTTAACAGGGAATGTACGCGAAGCCGCGTCCACTGTCAGAGCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTGACGCGGGAAAAAATGGGc                                                                   <  1:715874/148‑1 (MQ=255)
                                                                                   ttCTTAACAGGGAATGTACGCGAAGCCGCGTCCACTATCAGAGCAAAAAACGCGCACTGGAAAAAATCAGCGAGCAGGCGGCCAAACAACTAAGCCAGCCGCCACAGGAAATTTATGAAGCAATCCTGACTACGGAAAAAAAGCGCAg                                                                 >  2:514886/1‑148 (MQ=255)
                                                                                    tCTTAACAGGGAATGTACGCGAAGCCGCGTCCACTGTCAGAGCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTGACGCGGGAAAAAATGGGCAGt                                                                <  1:1719283/148‑1 (MQ=255)
                                                                                     cTTAACTGGGAATGTACGCGAAGCCGAGTCCACTGTCAGAGCAAAAACCGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGTTTTTTGAATCTTTCCTGACGCGGGAAAAAATGGGCAGTa                                                               <  1:974189/148‑1 (MQ=255)
                                                                                     cTTAACAGGGAATGTACGCGAAGCCGCGTCCACTGTCAGAGCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTGACGCGGGAAAAAATGGGCAGTa                                                               <  2:921117/148‑1 (MQ=255)
                                                                                      ttAACAGGGAATGTACGCGAAGCCGCGTCCACTGTCAGAGCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTGACGCGGGAAAAAATGGGCAGTAc                                                              >  1:2140515/1‑148 (MQ=255)
                                                                                             ggAATGTACGCGAAGCCGCGTCCACTGTCAGAGCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTGACGCGGGAAAAAATGGGCAGTACCGGTATc                                                       >  1:3160964/1‑148 (MQ=255)
                                                                                             ggAATGTACGCGAAGCCGCGTCCACTGTCAGAGCAAAAAACGCGCCCTGGAAATAATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTGACGCGGGAAAAAATGGGCAGTACCGGTATc                                                       >  1:2186163/1‑148 (MQ=255)
                                                                                              gAATGTACGCGAAGCCGCGTCCACTGTCAGAGCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTGACGCGGGAAAAAATGGGCAGTACCGGTATCg                                                      >  1:562993/1‑148 (MQ=255)
                                                                                                 tGTACGCGAAGCCGCGTCCACTGTCAGAGCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTGACGCGGGAAAAAATGGGCAGTACCGGTATCGGCa                                                   >  1:3361342/1‑148 (MQ=255)
                                                                                                      gcgAAGCCGCGTCCACTGTCAGAGCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTGACGCGGGAAAAAATGGGCAGTACCGGTATCGGCAATGGt                                              <  2:3555398/148‑1 (MQ=255)
                                                                                                       cgAAGCCGCGTCCACTGTCAGAGCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTGACGCGGGAAAAAATGGGCAGTACCGGTATCGGCAATGGTa                                             <  1:3474646/148‑1 (MQ=255)
                                                                                                        gAAGCCGCGTCCACTGTCAGAGCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTGACGCGGGAAAAAATGGGCAGTACCGGTATCGGCAATGGTAt                                            <  1:2705779/148‑1 (MQ=255)
                                                                                                              gcgTCCACTGTCAGAGCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTGACGCGGGAAAAAATGGGCAGTACCGGTATCGGCAATGGTATTGCCAt                                      >  2:2096128/1‑148 (MQ=255)
                                                                                                              gcgTCCACTGTCAGAGCAAAAAACGCGCCCTGGAAATAATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTGACGCGGGAAAAAATGGGCAGTACCGGTATCGGCAATGGTATTGCCAt                                      >  1:626604/1‑148 (MQ=255)
                                                                                                              gcgTCCACTGTCAGAGCAAAAAACGCGCCCTGGAAATAATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTGACGCGGGAAAAAATGGGCAGTACCGGTATCGGCAATGGTATTGCCAt                                      >  1:502736/1‑148 (MQ=255)
                                                                                                                     cTGTCAGAGCAAAAAACGCGCCCTGGAAATAATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTGACGCGGGAAAAAATGGGCAGTACCGGTATCGGCAATGGTATTGCCATTCCGCAt                               <  1:1011049/148‑1 (MQ=255)
                                                                                                                          agagCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTGACGCGGGAAAAAATGGGCAGTACCGGTATCGGCAATGGTATTGCCATTCCGCATGGCaa                          <  1:2737035/148‑1 (MQ=255)
                                                                                                                            agCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTGACGCGGGAAAAAATGGGCAGTACCGGTATCGGCAATGGTATTGCCATTCCGCATGGCAAAc                        >  2:2284299/1‑148 (MQ=255)
                                                                                                                            agCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTGACGCGGGAAAAAATGGGCAGTACCGGTATCGGCAATGGTATTGCCATTCCGCATGGCAAAc                        >  2:241074/1‑148 (MQ=255)
                                                                                                                                aaaaaCGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTGACGCGGGAAAAAATGGGCAGTACCGGTATCGGCAATGGTATTGCCATTCCGCATGGCAAACTGga                    >  1:1984737/1‑148 (MQ=255)
                                                                                                                                 aaaaCGCGCCCTGGAAATAATCAGCGAGCTGGAGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTGACGCGGGAAAAAATGGGCAGAACCGGTATCGTCAATGGTATTGCAATTCCGCATGGCAAACTGgaa                   >  2:1982646/1‑148 (MQ=255)
                                                                                                                                  aaaCGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTGACGCGGGAAAAAATGGGCAGTACCGGTATCGGCAATGGTATTGCCATTCCGCATGGCAAACTGgaag                  <  2:207969/148‑1 (MQ=255)
                                                                                                                                   aaCGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTGACGCGGGAAAAAATGGGCAGTACCGGTATCGGCAATGGTATTGCCATTCCGCATGGCAAACTGgaaga                 <  1:3062788/148‑1 (MQ=255)
                                                                                                                                     cgcgCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTGACGCGGGAAAAAATGGGCAGTACCGGTATCGGCAATGGTATTGCCATTCCGCATGGCAAACTGGaagaag               <  2:3118344/148‑1 (MQ=255)
                                                                                                                                       cgcCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTGACGCGGGAAAAAATGGGCAGTACCGGTATCGGCAATGGTATTGCCATTCCGCATGGCAAACTGGAAGAAGAt             <  2:1033671/148‑1 (MQ=255)
                                                                                                                                         cccTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTGACGCGGGAAAAAATGGGCAGTACCGGTATCGGCAATGGTATTGCCATTCCGCATGGCAAACTGGAAGAAGATAc           >  2:1671456/1‑148 (MQ=255)
                                                                                                                                         cccTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTGACGCGGGAAAAAATGGGCAGTACCGGTATCGGCAATGGTATTGCCATTCCGCATGGCAAACTGGAAGAAGATAc           >  2:1253077/1‑148 (MQ=255)
                                                                                                                                         cccTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTGACGCGGGAAAAAATGGGCAGTACCGGTATCGGCAATGGTATTGCCATTCCGCATGGCAAACTGGAAGAAGATAc           <  1:1911468/148‑1 (MQ=255)
                                                                                                                                         cccTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTGACGCGGGAAAAAATGGGCAGTACCGGTATCGGCAATGGTATTGCCATTCCGCATGGCAAACTGGAAGAAGATAc           <  1:1344351/148‑1 (MQ=255)
                                                                                                                                          ccTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTGACGCGGGAAAAAATGGGCAGTACCGGTATCGGCAATGGTATTGCCATTCCGCATGGCAAACTGGAAGAAGATAct          <  1:3202789/148‑1 (MQ=255)
                                                                                                                                             ggAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTGACGCGGGAAAAAATGGGCAGTACCGGTATCGGCAATGGTATTGCCATTCCGCATGGCAAACTGGAAGAAGATACTCTg       >  1:1796192/1‑148 (MQ=255)
                                                                                                                                             ggAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTGACGCGGGAAAAAATGGGCAGTACCGGTATCGGCAATGGTATTGCCATTCCGCATGGCAAACTGGAAGAAGATACTCTg       >  1:1790195/1‑148 (MQ=255)
                                                                                                                                               aaaTCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTGACGCGGGAAAAAATGGGCAGTACCGGTATCGGCAATGGTATTGCCATTCCGCATGGCAAACTGGAAGAAGATACTCTgcg     >  1:3312896/1‑148 (MQ=255)
                                                                                                                                                aaTCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTGACGCGGGAAAAAATGGGCAGTACCGGTATCGGCAATGGTATTGCCATTCCGCATGGCAAACTGGAAGAAGATACTCTgcgc    >  2:2005680/1‑148 (MQ=255)
                                                                                                                                                 atcatcAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTGACGCGGGAAAAAATGGGCAGTACCGGTATCGGCAATGGTATTGCCATTCCGCATGGCAAACTGGAAGAAGATACTCTgcgcg   >  2:368722/1‑148 (MQ=255)
                                                                                                                                                  tcatcaGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTGACGCGGGAAAAAATGGGCAGTACCGGTATCGGCAATGGTATTGCCATTCCGCATGGCAAACTGGAAGAAGATACTCTgcgcgc  <  2:3205686/148‑1 (MQ=255)
                                                                                                                                                  tcatcaGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTGACGCGGGAAAAAATGGGCAGTACCGGTATCGGCAATGGTATTGCCATTCCGCATGGCAAACTGGAAGAAGATACTCTgcgcgc  <  2:3205671/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CCGAGCCTGTTCCACTGTTTGAGTGGGCAGGTTCTTAGGTGAAATTATGACAAATAATGATACAACTCTACAGCTTAGCAGTGTTCTTAACAGGGAATGTACGCGAAGCCGCGTCCACTGTCAGAGCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTGACGCGGGAAAAAATGGGCAGTACCGGTATCGGCAATGGTATTGCCATTCCGCATGGCAAACTGGAAGAAGATACTCTGCGCGC  >  CP009273/3339891‑3340184

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: ‑