Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,782,353 C→A 5.8% R218S (CGC→AGC)  ppsR → bifunctional regulatory protein: PEP synthase kinase and PEP synthase pyrophosphorylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,782,3530CA 5.8% 774.4 / 15.9 204R218S (CGC→AGC) ppsRbifunctional regulatory protein: PEP synthase kinase and PEP synthase pyrophosphorylase
Reads supporting (aligned to +/- strand):  ref base C (96/96);  new base A (2/10);  total (98/106)
Fisher's exact test for biased strand distribution p-value = 3.48e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

Only 100 of 208 total aligned reads displayed.

CTGGCAATGCAGTTTGGTATCCGCGCGGCAAACTACCCCTTTATTGCCGACGATATGGATAATCTGGTGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGATCCCGTGGATTAACAGTACCAATTATTCGGTAGAAGAGATTGCCACCAAGATCC  >  CP009273/1782209‑1782500
                                                                                                                                                |                                                                                                                                                   
cTGGCAATGCAGTTTGGTATCCGCGCGGCAAACTACCCCTTTATTGCCGACGATATGGATAATCTGGTGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTcgcg                                                                                                                                                  >  1:3261682/1‑148 (MQ=255)
cTGGCAATGCAGTTTGGTATCCGCGCGGCAAACTACCCCTTTATTGCCGACGATATGGATAATCTGGTGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTcgcg                                                                                                                                                  >  2:2500193/1‑148 (MQ=255)
  ggCAATGCAGTTTGGTATCCGCGCGGCAAACTACCCCTTTATTGCCGACGATATGGATAATCTGGTGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAg                                                                                                                                                >  1:2255051/1‑148 (MQ=255)
        gCAGTTTGGTATCCGCGCGGCAAACTACCCCTTTATTGCCGACGATATGGATAATCTGGTGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGt                                                                                                                                          >  2:1651948/1‑148 (MQ=255)
               ggTATCCGCGCGGCAAACTACCCCTTTATTGCCGACGATATGGATAATCTGGTGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGgaga                                                                                                                                   <  2:1243553/148‑1 (MQ=255)
                   tCCGCGCGGCAAACTACCCCTTTATTGCCGACGATATGGATAATCTGGTGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAg                                                                                                                               >  1:2279986/1‑148 (MQ=255)
                       cgcgGCAAACTACCCCTTTATTGCCGACGATATGGATAATCTGGTGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGc                                                                                                                           >  1:3021522/1‑148 (MQ=255)
                           gCAAACTACCCCTTTATTGCCGACGATATGGATAATCTGGTGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATg                                                                                                                       <  1:429029/148‑1 (MQ=255)
                            cAAACTACCCCTTTATTGCCGACGATATGGATAATCTGGTGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGc                                                                                                                      >  2:423886/1‑148 (MQ=255)
                             aaaCTACCCCTTTATTGCCGACGATATGGATAATCTGGTGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGcc                                                                                                                     <  2:2228747/148‑1 (MQ=255)
                             aaaCTACCCCTTTATTGCCGACGATATGGATAATCTGGTGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTAGCGAGGAACGTCGGGAGAACAGTCGCTATGcc                                                                                                                     <  1:584440/148‑1 (MQ=255)
                              aaCTACCCCTTTATTGCCGACGATATGGATAATCTGGTGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCt                                                                                                                    >  1:1219068/1‑148 (MQ=255)
                              aaCTACCCCTTTATTGCCGACGATATGGATAATCTGGTGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCt                                                                                                                    <  1:1514455/148‑1 (MQ=255)
                              aaCTACCCCTTTATTGCCGACGATATGGATAATCTGGTGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCt                                                                                                                    >  1:3332548/1‑148 (MQ=255)
                              aaCTACCCCTTTATTGCCGACGATATGGATAATCTGGTGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCt                                                                                                                    >  1:398486/1‑148 (MQ=255)
                               aCTACCCCTTTATTGCCGACGATATGGATAATCTGGTGCTACCCGCGTCGCTAAAACCGCTTCAGCATAAATTTTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTc                                                                                                                   >  1:2144825/1‑148 (MQ=255)
                               aCTACCACTTTATTGCCGACGATATGGATAATCTGGTGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTc                                                                                                                   >  2:2698447/1‑148 (MQ=255)
                                cTACCCCTTTATTGCCGACGATATGGATAATCTGGTGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCg                                                                                                                  >  2:3131486/1‑148 (MQ=255)
                                 tACCCCTTTATTGCCGACGATATGGATAATCTGGTGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGc                                                                                                                 <  1:585339/148‑1 (MQ=255)
                                 tACCCCTTTATTGCCGACGATATGGATAATCTGGTGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTAGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGc                                                                                                                 <  1:2747532/148‑1 (MQ=255)
                                  aCCCCTTTATTGCCGACGATATGGATAATCTGGTGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCt                                                                                                                >  1:2058973/1‑148 (MQ=255)
                                  aCCCCTTTATTGCCGACGATATGGATAATCTGGTGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCt                                                                                                                <  1:2342338/148‑1 (MQ=255)
                                    cccTTTATTGCCGACGATATGGATAATCTGGTGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTc                                                                                                              >  1:115950/1‑148 (MQ=255)
                                     ccTTTATTGCCGACGATATGGATAATCTGGTGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCg                                                                                                             >  1:1921122/1‑148 (MQ=255)
                                     ccTTTATTGCCGACGATATGGATAATCTGGTGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCg                                                                                                             >  2:2939867/1‑148 (MQ=255)
                                      cTTTATTGCCGACGATATGGATAATCTGGTGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGt                                                                                                            >  2:704596/1‑148 (MQ=255)
                                        ttATTGCCGACGATATGGATAATCTGGTGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCa                                                                                                          <  1:1737870/148‑1 (MQ=255)
                                        ttATTGCCGACGATATGGATAATCTGGTGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCa                                                                                                          <  1:1636158/148‑1 (MQ=255)
                                        ttATTGCCGACGATATGGATAATCTGGTGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTAGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCa                                                                                                          <  2:2164621/148‑1 (MQ=255)
                                          aTTGCCGACGATATGGATAATCTGGTGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGt                                                                                                        >  1:935171/1‑148 (MQ=255)
                                           ttGCCGACGATATGGATAATCTGGTGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTAGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTg                                                                                                       <  1:1092391/148‑1 (MQ=255)
                                               cgacgaTATGGATAATCTGGTGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAgg                                                                                                   <  2:2623888/148‑1 (MQ=255)
                                                    atatGGATAATCTGGTGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGa                                                                                              >  1:800929/1‑148 (MQ=255)
                                                     tatGGATAATCTGGTGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGaa                                                                                             <  1:2832914/148‑1 (MQ=255)
                                                     tatGGATAATCTGGTGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCGTCGTCCGTGCAGGATGGaa                                                                                             >  1:2651952/1‑148 (MQ=255)
                                                      atGGATAATCTGGTGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGCTGGAAg                                                                                            >  2:3222721/1‑148 (MQ=255)
                                                      atGGATAATCTGGTGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAg                                                                                            <  2:2327974/148‑1 (MQ=255)
                                                       tGGATAATCTGGTGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGt                                                                                           <  1:494943/148‑1 (MQ=255)
                                                         gATAATCTGGTGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAAGTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTcg                                                                                         >  1:3307859/1‑148 (MQ=255)
                                                             aTCTGGTGCTACCCGCGTCGCTCACACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGa                                                                                     <  2:2513759/148‑1 (MQ=255)
                                                                 ggTGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTg                                                                                 <  2:3309551/148‑1 (MQ=255)
                                                                    gCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGaa                                                                              <  1:1273215/148‑1 (MQ=255)
                                                                      tACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGTAGTGGatgg                                                                            <  2:644825/148‑4 (MQ=255)
                                                                       aCCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGcc                                                                           <  2:2983326/148‑1 (MQ=255)
                                                                        cccGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTAGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCt                                                                          <  1:1621690/148‑1 (MQ=255)
                                                                         ccGCGTCGCTCAAACCGCGTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCtt                                                                         <  2:2585215/148‑1 (MQ=255)
                                                                          cgcgTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTg                                                                        <  2:1091246/148‑1 (MQ=255)
                                                                            cgTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTAGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCGTGTa                                                                      >  1:19384/1‑148 (MQ=255)
                                                                             gTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTAc                                                                     <  2:147389/148‑1 (MQ=255)
                                                                                gCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGt                                                                  >  1:2599179/1‑148 (MQ=255)
                                                                                   cAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTaaa                                                               >  1:1586616/1‑148 (MQ=255)
                                                                                    aaaCCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTaaaa                                                              <  1:2919650/148‑1 (MQ=255)
                                                                                     aaCCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTaaaaa                                                             >  1:725/1‑148 (MQ=255)
                                                                                     aaCCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTaaaaa                                                             <  2:796453/148‑1 (MQ=255)
                                                                                     aaCCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTaaaaa                                                             >  2:1967352/1‑148 (MQ=255)
                                                                                      aCCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGACGTGGAAGCCTTGTACCGTAAAAAt                                                            >  1:2988808/1‑148 (MQ=255)
                                                                                      aCCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAAt                                                            >  1:110105/1‑148 (MQ=255)
                                                                                       ccGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATc                                                           <  2:3427431/148‑1 (MQ=255)
                                                                                       ccGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATc                                                           >  1:3061359/1‑148 (MQ=255)
                                                                                           ttCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAANAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGAt                                                       >  1:1785022/1‑148 (MQ=255)
                                                                                           ttCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTAGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGAt                                                       <  1:3174148/148‑1 (MQ=255)
                                                                                            tCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGATc                                                      >  2:1691598/1‑148 (MQ=255)
                                                                                            tCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGATc                                                      <  1:2043573/148‑1 (MQ=255)
                                                                                             cAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGATcc                                                     >  2:594119/1‑148 (MQ=255)
                                                                                               gCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGATCCCg                                                   >  2:68120/1‑148 (MQ=255)
                                                                                               gCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGATCCCg                                                   >  1:1988790/1‑148 (MQ=255)
                                                                                                cATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGATCCCGt                                                  <  2:1949039/148‑1 (MQ=255)
                                                                                                cATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGATCCCGt                                                  <  2:3007644/148‑1 (MQ=255)
                                                                                                cATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGATCCCGt                                                  >  2:3020620/1‑148 (MQ=255)
                                                                                                  tAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGATCCCGTgg                                                >  2:2725321/1‑148 (MQ=255)
                                                                                                    aaTTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGATCCCGTGGAt                                              >  1:1046926/1‑148 (MQ=255)
                                                                                                        gTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGATCCCGTGGATTAAc                                          <  1:3510403/148‑1 (MQ=255)
                                                                                                         ttCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGATCCCGTGGATTAACa                                         <  1:1348024/148‑1 (MQ=255)
                                                                                                              ccTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGATCCCGTGGATTAACAGTAcc                                    >  2:2047486/1‑148 (MQ=255)
                                                                                                              ccTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGATCCCGTGGATTAACAGTAcc                                    <  2:1860166/148‑1 (MQ=255)
                                                                                                              ccTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGATCCCGTGGATTAACAGTAcc                                    >  2:275795/1‑148 (MQ=255)
                                                                                                               cTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGATCCCGTGGATTAACAGTACCa                                   <  2:2369380/148‑1 (MQ=255)
                                                                                                               cTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGATCCCGTGGATTAACAGTACCa                                   <  2:2331756/148‑1 (MQ=255)
                                                                                                               cTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGATCCCGTGGATTAACAGTACCa                                   <  2:1110610/148‑1 (MQ=255)
                                                                                                               cTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGATCCCGTGGATTAACAGTACCa                                   <  2:1033394/148‑1 (MQ=255)
                                                                                                                tGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGATCCCGTGGATTAACAGTACCaa                                  <  2:935171/148‑1 (MQ=255)
                                                                                                                tGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGATCCCGTGGATTAACAGTACCaa                                  <  2:2110214/148‑1 (MQ=255)
                                                                                                                 gACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGATCCCGTGGATTAACAGTACCAat                                 <  2:1343823/148‑1 (MQ=255)
                                                                                                                  aCTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGATCCCGTGGATTAACAGTACCAatt                                >  2:920678/1‑148 (MQ=255)
                                                                                                                    tATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGACGCCTTGTACCGTAAAAATCAGATCCCGTGGATTAACAGTACCAattat                              <  2:2823534/148‑1 (MQ=255)
                                                                                                                    tATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGATCCCGTGGATTAACAGTACCAattat                              <  2:1083426/148‑1 (MQ=255)
                                                                                                                    tATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGATCCCGTGGATTAACAGTACCAattat                              <  1:3226029/148‑1 (MQ=255)
                                                                                                                       cGACCCGGAACGTCTGGCGGCGATTAGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGATCCCGTGGATTAACAGTACCAATTATTCg                           <  2:1003605/148‑1 (MQ=255)
                                                                                                                        gACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGATCCCGTGGATTAACAGTACCAATTATTCgg                          <  1:3131486/148‑1 (MQ=255)
                                                                                                                         aCCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGATCCCGTGGATTAACAGTACCAATTATTCGGt                         <  2:401765/148‑1 (MQ=255)
                                                                                                                              gAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGATCCCGTGGATTAACAGTACCAATTATTCGGTagaag                    <  2:3240586/148‑1 (MQ=255)
                                                                                                                               aaCGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGATCCCGTGGATTAACAGTACCAATTATTCGGTagaaga                   >  1:430883/1‑148 (MQ=255)
                                                                                                                                 cGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGATCCCGTGGATTAACAGTACCAATTATTCGGTAGAAgaga                 <  1:2500194/148‑1 (MQ=255)
                                                                                                                                 cGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGATCCCGTGGATTAACAGTACCAATTATTCGGTAGAAgaga                 <  2:3168065/148‑1 (MQ=255)
                                                                                                                                 cGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGATCCCGTGGATTAACAGTACCAATTATTCGGTAGAAgaga                 >  2:863197/1‑148 (MQ=255)
                                                                                                                                       gcggcgATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGATCCCGTGGATTAACAGTACCAATTATTCGGTAGAAGAGATTGcca           <  2:443534/148‑1 (MQ=255)
                                                                                                                                        cggcgAGTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGATCCCGTGGATTAACAGTACCAATTATTCGGTAGAAGAGATTGccac          <  1:2885808/148‑1 (MQ=255)
                                                                                                                                           cgATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGATCCCGTGGATTAACAGTACCAATTATTCGGTAGAAGAGATTGCCACCaa       <  1:1002244/148‑1 (MQ=255)
                                                                                                                                           cgATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGATCCCGTGGATTAACAGTACCAATTATTCGGTAGAAGAGATTGCCACCaa       <  1:1282948/148‑1 (MQ=255)
                                                                                                                                            gATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGATCCCGTGGATTAACAGTACCAATTATTCGGTAGAAGAGATTGCCACCAAg      <  2:238137/148‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
CTGGCAATGCAGTTTGGTATCCGCGCGGCAAACTACCCCTTTATTGCCGACGATATGGATAATCTGGTGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGATCCCGTGGATTAACAGTACCAATTATTCGGTAGAAGAGATTGCCACCAAGATCC  >  CP009273/1782209‑1782500

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: