Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,782,353 C→A 15.0% 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,3530CA15.0% 635.2 / 62.4 206R218S (CGC→AGC) ppsRbifunctional regulatory protein: PEP synthase kinase and PEP synthase pyrophosphorylase
Reads supporting (aligned to +/- strand):  ref base C (76/99);  new base A (10/21);  total (86/120)
Fisher's exact test for biased strand distribution p-value = 3.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.19e-01

Only 100 of 210 total aligned reads displayed.

TCTGGCAATGCAGTTTGGTATCCGCGCGGCAAACTACCCCTTTATTGCCGACGATATGGATAATCTGGTGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGATCCCGTGGATTAACAGTACCAATTATTCGGTAGAAGAGATTGCCACCAAGATCC  >  CP009273/1782208‑1782500
                                                                                                                                                 |                                                                                                                                                   
tCTGGCAATGCAGTTTGGTATCCGCGCGGCAAACTACCCCTTTATTGCCGACGATATGGATAATCTGGTGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTcgc                                                                                                                                                   >  2:1758245/1‑148 (MQ=255)
tCTGGCAATGCAGTTTGGTATCCGCGCGGCAAACTACCCCTTTATTGCCGACGATATGGATAATCTGGTGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTcgc                                                                                                                                                   >  2:2631007/1‑148 (MQ=255)
   ggCAATGCAGTTTGGTATCCGCGCGGCAAACTACCCCTTTATTGCCGACGATATGGATAATCTGGTGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGACCGTCTGGCGGCGATTCGCGAg                                                                                                                                                >  2:1176423/1‑148 (MQ=255)
    gCAATGCAGTTTGGTATCCGCGCGGCAAACTACCCCTTTATTGCCGACGATATGGATAATCTGGTGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAgg                                                                                                                                               <  1:3545786/148‑1 (MQ=255)
     cactGCAGTTTGGTATCCGCGCGGCACACTACCCCTTTATTGCCGACGATATGGATAATCTGGTGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGa                                                                                                                                              <  2:1806922/145‑1 (MQ=255)
     cAATGCAGTTTGGTATCCGCGCGGCAAACTACCCCTTTATTGCCGACGATATGGATAATCTGGTGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGa                                                                                                                                              <  2:3462434/148‑1 (MQ=255)
     cAATGCAGTTTGGTATCCGCGCGGCAAACTACCCCTTTATTGCCGACGATATGGATAATCTGGTGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGa                                                                                                                                              <  2:1448368/148‑1 (MQ=255)
     cAATGCAGTTTGGTATCCGCGCGGCAAACTACCCCTTTATTGCCGACGATATGGATAATCTGGTGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGa                                                                                                                                              >  2:2746967/1‑148 (MQ=255)
     cAATGCAGTTTGGTATCCGCGCGGCAAACTACCCCTTTATTGCCGACGATATGGATAATCTGGTGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGa                                                                                                                                              <  2:1713740/148‑1 (MQ=255)
     cAATGCAGTTTGGTATCCGCGCGGCAAACTACCCCTTTATTGCCGACGATATGGATAATCTGGTGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTAGCGAGGa                                                                                                                                              <  1:3715376/148‑1 (MQ=255)
      aaTGCAGTTTGGTATCCGCGCGGCAAACTACCCCTTTATTGCCGACGATATGGATAATCTGGTGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGaa                                                                                                                                             <  2:3230144/148‑1 (MQ=255)
             tttGGTATCCGCGCGGCAAACTACCCCTTTATTGCCGACGATATGGATAATCTGGTGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCggg                                                                                                                                      <  1:1419512/148‑1 (MQ=255)
             tttGGTATCCGCGCGGCAAACTACCCCTTTATTGCCGACGATATGGATAATCTGGTGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCggg                                                                                                                                      >  1:2308943/1‑148 (MQ=255)
                ggTATCCGCGCGGCAAACTACCCCTTTATTGCCGACGATATGGATAATCTGGTGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGgaga                                                                                                                                   <  2:2472538/148‑1 (MQ=255)
                ggTATCCGCGCGGCAAACTACCCCTTTATTGCCGACGATATGGATAATCTGGTGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAACTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGgaga                                                                                                                                   <  1:1949127/148‑1 (MQ=255)
                  tATCCGCGCGGCAAACTACCCCTTTATTGCCGACGATATGGATAATCTGGTGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAAc                                                                                                                                 >  1:2993162/1‑148 (MQ=255)
                    tCCGCGCGGCAAACTACCCCTTTATTGCCGACGATATGGATAATCTGGTGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTAGCGAGGAACGTCGGGAGAACAg                                                                                                                               <  2:3176226/148‑1 (MQ=255)
                       gcgcgGCAAACTACCCCTTTATTGCCGACGATATGGATAATCTGGTGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAAGAACGTCGGGAGAACAGTCg                                                                                                                            <  2:333028/148‑1 (MQ=255)
                        cgcgGCAAACTACCCCTTTATTGCCGACGATATGGATAATCTGGTGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGc                                                                                                                           <  1:1184775/148‑1 (MQ=255)
                        cgcgGCAAACTACCCCTTTATTGCCGACGATATGGATAATCTGGTGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGc                                                                                                                           >  2:1247005/1‑148 (MQ=255)
                         gcgGCAAACTACCCCTTTATTGCCGACGATATGGATAATCTGGTGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCt                                                                                                                          <  1:3152998/148‑1 (MQ=255)
                         gcgGCAAACTACCCCTTTATTGCCGACGATATGGATAATCTGGTGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCt                                                                                                                          <  1:1571636/148‑1 (MQ=255)
                         gcgGCAAACTACCCCTTTATTGCCGACGATATGGATAATCTGGTGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCt                                                                                                                          <  1:3278759/148‑1 (MQ=255)
                             cAAACTACCCCTTTATTGCCGACGATATGGATAATCTGGTGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGc                                                                                                                      <  1:1638751/148‑1 (MQ=255)
                                aCTACCCCTTTATTGCCGACGATATGGATAATCTGGTGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTc                                                                                                                   >  2:3078774/1‑148 (MQ=255)
                                aCTACCCCTTTATTGCCGACGATATGGATAATCTGGTGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTc                                                                                                                   >  1:816309/1‑148 (MQ=255)
                                   aCCCCTTTATTGCCGACGATATGGATAATCTGGTGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCt                                                                                                                <  2:1071751/148‑1 (MQ=255)
                                       cTTTATTGCCGACGATATGGATAATCTGGTGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGt                                                                                                            >  2:932555/1‑148 (MQ=255)
                                        tttATTGCCGACGATATGGATAATCTGGTGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTAGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTc                                                                                                           >  2:1861027/1‑148 (MQ=255)
                                         ttATTGCCGACGATATGGATAATCTGGTGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCa                                                                                                          <  2:1291416/148‑1 (MQ=255)
                                            ttGCCGACGATATGGATAATCTGGTGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTg                                                                                                       <  1:1758245/148‑1 (MQ=255)
                                            ttGCCGACGATATGGATAATCTGGTGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTg                                                                                                       >  1:3492091/1‑148 (MQ=255)
                                              gCCGACGATATGGATAATCTGGTGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCa                                                                                                     <  2:1562208/148‑1 (MQ=255)
                                              gCCGACGATATGGATAATCTGGTGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCa                                                                                                     <  1:1476659/148‑1 (MQ=255)
                                                 gacgaTATGGATAATCTGGTGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGa                                                                                                  >  2:2241051/1‑148 (MQ=255)
                                                 gacgaTATGGATAATCTGGTGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGa                                                                                                  <  2:3535948/148‑1 (MQ=255)
                                                  acgaTATGGATAATCTGGTGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGAt                                                                                                 <  1:3287353/148‑1 (MQ=255)
                                                     atatGGATAATCTGGTGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGa                                                                                              >  1:585254/1‑148 (MQ=255)
                                                     atatGGATAATCTGGTGCGACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTAGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGa                                                                                              <  2:1208915/148‑1 (MQ=255)
                                                       atGGATAATCTGGTGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAg                                                                                            >  2:3627181/1‑148 (MQ=255)
                                                       atGGATAATCTGGTGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAg                                                                                            <  1:1891999/148‑1 (MQ=255)
                                                          gATAATCTGGTGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTcg                                                                                         >  1:2525004/1‑148 (MQ=255)
                                                           aTAATCTGGTGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTAGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTcgc                                                                                        >  2:3560799/1‑148 (MQ=255)
                                                              aTCTGGTGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGa                                                                                     >  2:650990/1‑148 (MQ=255)
                                                               tCTGGTGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGaa                                                                                    >  2:3117499/1‑148 (MQ=255)
                                                                cTGGTGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAg                                                                                   <  2:2896228/148‑1 (MQ=255)
                                                                cTGGTGCGACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAg                                                                                   <  2:2896264/148‑1 (MQ=255)
                                                                 tGGTGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGt                                                                                  <  1:1982257/148‑1 (MQ=255)
                                                                  ggTGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTg                                                                                 <  1:458905/148‑1 (MQ=255)
                                                                        aCCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGcc                                                                           <  2:121225/148‑1 (MQ=255)
                                                                          ccGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGACCGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCtt                                                                         >  1:2244353/1‑148 (MQ=255)
                                                                           cgcgTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTg                                                                        >  2:610700/1‑148 (MQ=255)
                                                                                cGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTAGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCg                                                                   >  2:1520470/1‑148 (MQ=255)
                                                                                     aacccGCGTCAGCATAAAGTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTaaaa                                                              <  1:2925745/145‑1 (MQ=255)
                                                                                     aaaCCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTaaaa                                                              <  1:1015440/148‑1 (MQ=255)
                                                                                     aaaCCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTaaaa                                                              <  1:1717444/148‑1 (MQ=255)
                                                                                      aaCCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTaaaaa                                                             <  2:3422175/148‑1 (MQ=255)
                                                                                      aaCCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTaaaaa                                                             >  1:1841805/1‑148 (MQ=255)
                                                                                      aaCCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTAGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTaaaaa                                                             <  1:666245/148‑1 (MQ=255)
                                                                                       aCCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAAt                                                            >  2:2861082/1‑148 (MQ=255)
                                                                                        ccGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATc                                                           >  2:2705516/1‑148 (MQ=255)
                                                                                         cGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGACCGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCa                                                          >  2:3808038/1‑148 (MQ=255)
                                                                                             tCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGATc                                                      >  1:2866339/1‑148 (MQ=255)
                                                                                             tCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGATc                                                      >  2:2977896/1‑148 (MQ=255)
                                                                                               aGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTAGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGATccc                                                    <  1:1158304/148‑1 (MQ=255)
                                                                                                gCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGATCCCg                                                   >  2:1960688/1‑148 (MQ=255)
                                                                                                gCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAAGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGATCCCg                                                   >  2:3549997/1‑148 (MQ=255)
                                                                                                 cATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGATACCGt                                                  >  1:9552/1‑148 (MQ=255)
                                                                                                   tAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGATCCCGTgg                                                <  2:133847/148‑1 (MQ=255)
                                                                                                   tAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGATCCCGTgg                                                >  1:1798407/1‑148 (MQ=255)
                                                                                                    acaTTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGATCCCGTGGa                                               <  2:350264/146‑1 (MQ=255)
                                                                                                    aaaTTGGTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTAGCGAGGAACGTCGGGAGTACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGACTTGTACCGTAAAAATCAGATCCCGTGGa                                               <  2:3722780/148‑1 (MQ=255)
                                                                                                     aaTTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGATCCCGTGGAt                                              <  2:4830/148‑1 (MQ=255)
                                                                                                       ttGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGATCCCGTGGATTa                                            <  2:57795/148‑1 (MQ=255)
                                                                                                        tGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGATCCCGTGGATTaa                                           <  1:1048481/148‑1 (MQ=255)
                                                                                                        tGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGATCCCGTGGATTaa                                           <  1:2788848/148‑1 (MQ=255)
                                                                                                              gCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGATCCCGTGGATTAACAGTAc                                     <  2:702876/148‑1 (MQ=255)
                                                                                                                cTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTGGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGACTTGTACCGTAAAAATCAGATCCCGTTGATTAACAGTACCa                                   <  2:936399/148‑1 (MQ=255)
                                                                                                                cTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGATCCCGTGGATTAACAGTACCa                                   <  2:584044/148‑1 (MQ=255)
                                                                                                                 tGACTATCGACCCGGAACGTCTGGCGGCGATTAGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGATCCCGTGGATTAACAGTACCaa                                  >  2:1050540/1‑148 (MQ=255)
                                                                                                                    cTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGATCCCGTGGATTAACAGTACCAatta                               <  2:2977586/148‑1 (MQ=255)
                                                                                                                     tATCGACCCGGAACGTCTGGCGGCGATTAGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGATCCCGTGGATTAACAGTACCAattat                              <  2:2760758/148‑1 (MQ=255)
                                                                                                                       tCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGATCCCGTGGATTAACAGTACCAATTATTc                            >  1:2236459/1‑148 (MQ=255)
                                                                                                                       tCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGATCCCGTGGATTAACAGTACCAATTATTc                            >  2:2367229/1‑148 (MQ=255)
                                                                                                                        cGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGATCCCGTGGATTAACAGTACCAATTATTCg                           >  1:1987325/1‑148 (MQ=255)
                                                                                                                        cGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGATCCCGTGGATTAACAGTACCAATTATTCg                           >  1:3067834/1‑148 (MQ=255)
                                                                                                                         gACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGATCCCGTGGATTAACAGTACCAATTATTCgg                          <  1:2631007/148‑1 (MQ=255)
                                                                                                                         gACCCGGAACGTCTGGCGGCGATTAGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGATCCCGTGGATTAACAGTACCAATTATTCgg                          >  2:1534365/1‑148 (MQ=255)
                                                                                                                           cccGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACTGTAAAAATCAGATCCCGTGGATTAACAGTACCAATTATTCGGTa                        >  1:438254/1‑148 (MQ=255)
                                                                                                                            ccGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGATCCCGTGGATTAACAGTACCAATTATTCGGTag                       >  2:3042911/1‑148 (MQ=255)
                                                                                                                             cGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGATCCCGTGGATTAACAGTACCAATTATTCGGTaga                      <  2:2035077/148‑1 (MQ=255)
                                                                                                                                   gTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGATCCCGTGGATTAACAGTACCAATTATTCGGTAGAAGAGAt                <  1:2461027/148‑1 (MQ=255)
                                                                                                                                    tCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGATCCCGTGGATTAACAGTACCAATTATTCGGTAGAAGAGAtt               >  1:3620655/1‑148 (MQ=255)
                                                                                                                                     cTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGATCCCGTGGATTAACTGTACCAATTATTCGGTAGAAGAGATTg              >  1:455779/1‑148 (MQ=255)
                                                                                                                                     cTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGATCCCGTGGATTAACAGTACCAATTATTCGGTAGAAGAGATTg              <  2:1541306/148‑1 (MQ=255)
                                                                                                                                      tGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGATCCCGTGGATTAACAGTACCAATTATTCGGTAGAAGAGATTGc             <  1:3331555/148‑1 (MQ=255)
                                                                                                                                           gcgATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGATCCCGTGGATTAACAGTACCAATTATTCGGTAGAAGAGATTGccacca        <  2:2525004/148‑1 (MQ=255)
                                                                                                                                           gcgATTAGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGATCCCGTGGATTAACAGTACCAATTATTCGGTAGAAGAGATTGccacca        >  2:1624713/1‑148 (MQ=255)
                                                                                                                                               ttCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTTCCGTAAAAATCAGATCCCGTGGATTAACAGTACCAATTATTCGGTAGAAGAGATTGCCACCAAgct    >  1:3193989/1‑146 (MQ=255)
                                                                                                                                                 cgcgAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGATCCCGTGGATTAACAGTACCAATTATTCGGTAGAAGAGATTGCCACCAAGATcc  <  1:2059887/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
TCTGGCAATGCAGTTTGGTATCCGCGCGGCAAACTACCCCTTTATTGCCGACGATATGGATAATCTGGTGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAGTCGCTATGCCTCGCTTCGTCAGTGCAGGATGGAAGTCGCGGAAGTGGAAGCCTTGTACCGTAAAAATCAGATCCCGTGGATTAACAGTACCAATTATTCGGTAGAAGAGATTGCCACCAAGATCC  >  CP009273/1782208‑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: