Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,561,258 C→A 35.4% D46Y (GAT→TAT)  glgA ← glycogen synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,561,2580CA35.4% 457.2 / 269.2 279D46Y (GAT→TAT) glgAglycogen synthase
Reads supporting (aligned to +/- strand):  ref base C (78/102);  new base A (55/44);  total (133/146)
Fisher's exact test for biased strand distribution p-value = 6.03e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 283 total aligned reads displayed.

ATACGGGCTTCCCGGACGATCATAGAGATGCGGCGCGTCAATCAGGTAAATGCCAACCCCGTTGTAATGACCGAACAACAGCGTGATATGTCCGGCGAAGGTATCACGACGGGATACTACCTGCGCATCGGTCACGCCACGGCGAATATCGGGAAATGCAGGCAACAGTACGCGAGCGTCAACGCCGTCTGCGATTTGTGCTGCGGGTAATGCCCCAATAACATCAGCCAGACCGCCGGTTTTAAGCAGCGGGAACATCTCTGAACATACATGTAAAACCTGCATTATCGCTCCTGTTT  >  CP009273/3561109‑3561407
                                                                                                                                                     |                                                                                                                                                     
  aCGGGCTTCCCGGACGATCATAGAGATGCGGCGCGTCAATCAGGTAAATGCCAACCCCGTTGTAATGACCGAACAACAGCGTGATATGTCCGGCGAAGGTATCACGACGGGATACTACCTGCGCATCGGTCACGCCACGGCGAATATAgg                                                                                                                                                     <  1:2783357/150‑1 (MQ=255)
    gggCTTCCCGGACGATCATAGAGATGCGGCGCGTCAATCAGGTAAATGCCAACCCCGTTGTAATGACCGAACAACAGCGTGATATGTCCGGCGAAGGTATCACGACGGGATACTACCTGCGCATCGGTCACGCCACGGCGAATATCGGGa                                                                                                                                                   >  2:2953529/1‑150 (MQ=255)
         tCCCGGACGATCATAGAGATGCGGCGCGTCAATCAGGTAAATGCCAACCCCGTTGTAATGACCGAACAACAGCGTGATATGTCCGGCGAAGGTATCACGACGGGATACTACCTGCGCATCGGTCACGCCACGGCGAATATAGGGAAATGc                                                                                                                                              <  2:3015233/150‑1 (MQ=255)
            cGGACGATCATAGAGATGCGGCGCGTCAATCAGGTAAATGCCAACCCCGTTGTAATGACCGAACAACAGCGTGATATGTCCGGCGAAGGTATCACGACGGGATACTACCTGCGCATCGGTCACGCCACGGCGAATATCGGGAAATGCAgg                                                                                                                                           >  2:1927702/1‑150 (MQ=255)
            cGGACGATCATAGAGATGCGGCGAGTCAATCAGGTATATGCCAACCCCGTTGTAATGACCGAACAACAGCGTGATATGTCCGGCGAAGGTATCACGACGGGATACTACCTGCGCATCGGTCACGCCACGGCGAATATAGGGAAATGCAgg                                                                                                                                           <  1:3944608/150‑1 (MQ=255)
              gACGATCATAGAGATGCGGCGCGTCAATCAGGTAAATGCCAACCCCGTTGTAATGACCGAACAACAGCGTGATATGTCCGGCGAAGGTATCACGACGGGATACTACCTGCGCATCGGTCACGCCACGGCGAATATCGGGAAATGCAGGCa                                                                                                                                         <  1:28204/150‑1 (MQ=255)
                 gATCATAGAGATGCGGCGCGTCAATCAGGTAAATGCCAACCCCGTTGTAATGACCGAACAACAGCGTGATATGTCCGGCGAAGGTATCACGACGGGATACTACCTGCGCATCGGTCACGCCACGGCGAATATCGGGAAATGCAGGCAACa                                                                                                                                      <  2:2863461/150‑1 (MQ=255)
                   tCATAGAGATGCGGCGCGTCAATCAGGTAAATGCCAACCCCGTTGTAATGACCGAACAACAGCGTGATATGTCCGGCGAAGGTATCACGACGGGATACTACCTGCGCATCGGTCACGCCACGGCGAATATCGGGAAATGCAGGCAACAGt                                                                                                                                    >  1:1139192/1‑150 (MQ=255)
                     aTAGAGATGCGGCGCGTCAATCAGGTAAATGCCAACCCCGTTGTAATGACCGAACAACAGCGTGATATGTCCGGCGAAGGTATCACGACGGGATACTACCTGCGCATCGGTCACGCCACGGCGAATATCGGGAAATGCAGGCAACAGTAc                                                                                                                                  >  2:83023/1‑150 (MQ=255)
                     aTAGAGATGCGGCGCGTCAATCAGGTAAATGCCAACCCCGTTGTAATGACCGAACAACAGCGTGATATGTCCGGCGAAGGTATCACGACGGGATACTACCTGCGCATCGGTCACGCCACGGCGAATATCGGGAAATGCAGGCAACAGTAc                                                                                                                                  >  2:3716966/1‑150 (MQ=255)
                     aTAGAGATGCGGCGCGTCAATCAGGTAAATGCCAACCCCGTTGTAATGACCGAACAACAGCGTGATATGTCCGGCGAAGGTATCACGACGGGATACTACCTGCGCATCGGTCACGCCACGGCGAATATAGGGAAATGCAGGCAACAGTAc                                                                                                                                  >  1:2151160/1‑150 (MQ=255)
                      tAGAGATGCGGCGCGTCAATCAGGTAAATGCCAACCCCGTTGTAATGACCGAACAACAGCGTGATATGTCCGGCGAAGGTATCACGACGGGATACTACCTGCGCATCGGTCACGCCACGGCGAATATCGGGAAATGCAGGCAACAGTAcg                                                                                                                                 <  1:4144719/150‑1 (MQ=255)
                      tAGAGATGCGGCGCGTCAATCAGGTAAATGCCAACCCCGTTGTAATGACCGAACAACAGCGTGATATGTCCGGCGAAGGTATCACGACGGGATACTACCTGCGCATCGGTCACGCCACGGCGAATATAGGGAAATGCAGGCAACAGTAcg                                                                                                                                 <  2:4605929/150‑1 (MQ=255)
                       agagATGCGGCGCGTCAATCAGGTAAATGCCAACCCCGTTGTAATGACCGAACAACAGCGTGATATGTCCGGCGAAGGTATCACGACGGGATACTACCTGCGCATCGGTCACGCCACGGCGAATATAGGGAAATGCAGGCAACAGTAcgc                                                                                                                                >  1:4623116/1‑150 (MQ=255)
                        gagaTGCGGCGCGTCAATCAGGTAAATGCCAACCCCGTTGTAATGACCGAACAACAGCGTGATATGTCCGGCGAAGGTATCACGACGGGATACTACCTGCGCATCGGTCACGCCACGGCGAATATCGGGAAATGCAGGCAACAGTAcgcg                                                                                                                               >  2:1177374/1‑150 (MQ=255)
                         agaTGCGGCGCGTCAATCAGGTAAATGCCAACCCCGTTGTAATGACCGAACAACAGCGTGATATGTCCGGCGAAGGTATCACGACGGGATACTACCTGCGCATCGGTCACGCCACGGCGAATATCGGGAAATGCAGGCAACAGTACGCGa                                                                                                                              >  1:1470641/1‑150 (MQ=255)
                          gaTGCGGCGCGTCAATCAGGTAAATGCCAACCCCGTTGTAATGACCGAACAACAGCGTGATATGTCCGGCGAAGGTATCACGACGGGATACTACCTGCGCATCGGTCACGCCACGGCGAATATCGGGAAATGCAGGCAACAGTACGCGAg                                                                                                                             <  1:3397214/150‑1 (MQ=255)
                          gaTGCGGCGCGTCAATCAGGTAAATGCCAACCCCGTTGTAATGACCGAACAACAGCGTGATATGTCCGGCGAAGGTATCACGACGGGATACTACCTGCGCATCGGTCACGCCACGGCGAATATCGGGAAATGCAGGCAACAGTACGCGAg                                                                                                                             <  1:1964695/150‑1 (MQ=255)
                              cggcgCGTCAATCAGGTAAATGCCAACCCCGTTGTAATGACCGAACAACAGCGTGATATGTCCGGCGAAGGTATCACGACGGGATACTACCTGCGCATCGGTCACGCCACGGCGAATATCGGGAAATGCAGGCAACAGTACGCGAGCGTc                                                                                                                         <  1:3119614/150‑1 (MQ=255)
                               ggcgCGTCAATCAGGTAAATGCCAACCCCGTTGTAATGACCGAACAACAGCGTGATATGTCCGGCGAAGGTATCACGACGGGATACTACCTGCGCATCGGTCACGCCACGGCGAATATCGGGAAATGCAGGCAACAGTACGCGAGCGTCa                                                                                                                        >  2:1374658/1‑150 (MQ=255)
                                 cgcgTCAATCAGGTAAATGCCAACCCCGTTGTAATGACCGAACAACAGCGTGATATGTCCGGCGAAGGTATCACGACGGGATACTACCTGCGCATCGGTCACGCCACGGCGAATATCGGGAAATGCAGGCAACAGTACGCGAGCGTCAAc                                                                                                                      >  1:1552633/1‑150 (MQ=255)
                                 cgcgTCAATCAGGTAAATGCCAACCCCGTTGTAATGACCGAACAACAGCGTGATATGTCCGGCGAAGGTATCACGACGGGATACTACCTGCGCATCGGTCACGCCACGGCGAATATAGGGAAATGTAGGCAACAGTACGCGAGCGTCAAc                                                                                                                      >  1:930402/1‑150 (MQ=255)
                                     tCAATCAGGTAAATGCCAACCCCGTTGTAATGACCGAACAACAGCGTGATATGTCCGGCGAAGGTATCACGACGGGATACTACCTGCGCATCGGTCACGCCACGGCGAATATCGGGAAATGCAGGCAACAGTACGCGAGCGTCAACGCCg                                                                                                                  <  1:401287/150‑1 (MQ=255)
                                     tCAATCAGGTAAATGCCAACCCCGTTGTAATGACCGAACAACAGCGTGATATGTCCGGCGAAGGTATCACGACGGGATACTACCTGCGCATCGGTCACGCCACGGCGAATATCGGGAAATGCAGGCAACAGTACGCGAGCGTCAACGCCg                                                                                                                  >  1:2360845/1‑150 (MQ=255)
                                     tCAATCAGGTAAATGCCAACCCCGTTGTAATGACCGAACAACAGCGTGATATGTCCGGCGAAGGTATCACGACGGGATACTACCTGCGCATCGGTCACGCCACGGCGAATATAGGGAAATGCAGGCAACAGTACGCGAGCGTCAACGCCg                                                                                                                  >  1:205993/1‑150 (MQ=255)
                                       aaTCAGGTAAATGCCAACCCCGTTGTAATGACCGAACAACAGCGTGATATGTCCGGCGAAGGTATCACGACGGGATACTACCTGCGCATCGGTCACGCCACGGCGAATATAGGGAAATGCAGGCAACAGTACGCGAGCGTCAACGCCGTc                                                                                                                >  1:953126/1‑150 (MQ=255)
                                       aaTCAGGTAAATGCCAACCCCGTTGTAATGACCGAACAACAGCGTGATATGTCCGGCGAAGGTATCACGACGGGATACTACCTGCGCATCGGTCACGCCACGGCGAATATAGGGAAATGCAGGCAACAGTACGCGAGCGTCAACGCCGTc                                                                                                                >  2:384880/1‑150 (MQ=255)
                                        aTCAGGTAAATGCCAACCCCGTTGTAATGACCGAACAACAGCGTGATATGTCCGGCGAAGGTATCACGACGGGATACTACCTGCGCATCGGTCACGCCACGGCGAATATCGGGAAAAGCAGGCAACAGTACGCGAGCGTCAACGCCGTCt                                                                                                               >  1:1202592/1‑150 (MQ=255)
                                         tCAGGTAAATGCCAACCCCGTTGTAATGACCGAACAACAGCGTGATATGTCCGGCGAAGGTATCACGACGGGATACTACCTGCGCATCGGTCACGCCACGGCGAATATCGGGAAATGCAGGCAACAGTACGCGAGCGTCAACGCCGTCTg                                                                                                              <  2:2269668/150‑1 (MQ=255)
                                          cAGGTAAATGCCAACCCCGTTGTAATGACCGAACAACAGCGTGATATGTCCGGCGAAGGTATCACGACGGGATACTACCTGCGCATCGGTCACGCCACGGCGAATATAGGGAAATGCAGGCAACAGTACGCGAGCGTCAACGCCGTCTGc                                                                                                             <  1:3666310/150‑1 (MQ=255)
                                           aGGTAAATGCCAACCCCGTTGTAATGACCGAACAACAGCGTGATATGTCCGGCGAAGGTATCACGACGGGATACTACCTTCGCATCGGTCACGCCACGGCGAATATAGGGAAATGCAGGCAACAGTACGCGAGCGTCAACGCCGTCTGCg                                                                                                            >  1:3730198/1‑150 (MQ=255)
                                           aGGTAAATGCCAACCCCGTTGTAATGACCGAACAACAGCGTGATATGTCCGGCGAAGGTATCACGACGGGATACTACCTGCGCATCGGTCACGCCACGGCGAATATCGGGAAATGCAGGCAACAGTACGCGAGCGTCAACGCCGTCTGCg                                                                                                            <  2:125515/150‑1 (MQ=255)
                                            ggTAAATGCCAACCCCGTTGTAATGACCGAACAACAGCGTGATATGTCCGGCGAAGGTATCACGACGGGATACTACCTGCGCATCGGTCACGCCACGGCGAATATAGGGAATGCAGGCAACAGTACGCGAGCGTCAACGCCGTCTGCGAt                                                                                                          >  1:1557155/1‑150 (MQ=255)
                                             gTAAATGCCAACCCCGTTGTAATGACCGAACAACAGCGTGATATGTCCGGCGAAGGTATCACGACGGGATACTACCTGCGCATCGGTCACGCCACGGCGAATATCGGGAAATGCAGGCAACAGTACGCGAGCGTCAACGCCGTCTGCGAt                                                                                                          <  1:3139028/150‑1 (MQ=255)
                                             gTAAATGCCAACCCCGTTGTAATGACCGAACAACAGCGTGATATGTCCGGCGAAGGTATCACGACGGGATACTACCTGCGCATCGGTCACGCCACGGCGAATATCGGGAAATGCAGGCAACAGTACGCGAGCGTCAACGCCGTCTGCGAt                                                                                                          <  2:3142327/150‑1 (MQ=255)
                                               aaaTGCCAACCCCGTTGTAATGACCGAACAACAGCGTGATATGTCCGGCGAAGGTATCACGACGGGATACTACCTGCGCATCGGTCACGCCACGGCGAATATCGGGAAATGCAGGCAACAGTACGCGAGCGTCAACGCCGTCTGCGAttt                                                                                                        <  2:2412777/150‑1 (MQ=255)
                                                aaTGCCAACCCCGTTGTAATGACCGAACAACAGCGTGATATGTCCGGCGAAGGTATCACGACGGGATACTACCTGCGCATCGGTCACGCCACGGCGAATATAGGGAAATGCAGGCAACAGTACGCGAGCGTCAACGCCGTCTGCGATTtg                                                                                                       >  1:3281286/1‑150 (MQ=255)
                                                 aTGCCAACCCCGTTGTAATGACCGAACAACAGCGTGATATGTCCGGCGAAGGTATCACGACGGGATACTACCTGCGCATCGGTCACGCCACGGCGAATATAGGGAAATGCAGGCAACAGTACGCGAGCGTCAACGCCGTCTGCGATTtgt                                                                                                      >  1:3681176/1‑150 (MQ=255)
                                                  tGCCAACCCCGTTGTAATGACCGAACAACAGCGTGATATGTCCGGCGAAGGTATCACGACGGGATACTACCTGCGCATCGGTCACGCCACGGCGAATATCGGGAAATGCAGGCAACAGTACGCGAGCGTCAACGCCGTCTGCGATTtgtg                                                                                                     <  2:266643/150‑1 (MQ=255)
                                                    ccAACCCCGTTGTAATGACCGAACAACAGCGTGATATGTCCGGCGAAGGTATCACGACGGGATACTACCTGCGCATCGGTCACGCCACGGCGAATATAGGGAAATGCAGGCAACAGTACGCGAGCGTCAACGCCGTCTGCGATTTGtgct                                                                                                   <  1:3056455/150‑1 (MQ=255)
                                                     cAACCCCGTTGTAATGACCGAACAACAGCGTGATATGTCCGGCGAAGGTATCACGACGGGATACTACCTGCGCATCGGTCACGCCACGGCGAATATCGGGAAATGCAGGCAACAGTACGCGAGCGTCAACGCCGTCTGCGATTTGtgctg                                                                                                  >  2:1294623/1‑150 (MQ=255)
                                                     cAACCCCGTTGTAATGACCGAACAACAGCGTGATATGTCCGGCGAAGGTATCACGACGGGATACTACCTGCGCATCGGTCACGCCACGGCGAATATCGGGAAATGCAGGCAACAGTACGCGAGCGTCAACGCCGTCTGCGATTTGtgctg                                                                                                  <  2:3853199/150‑1 (MQ=255)
                                                      aaCCCCGTTGTAATGACCGAACAACAGCGTGATATGTCCGGCGAAGGTATCACGACGGGATACTACCTGCGCATCGGTCACGCCACGGCGAATATAGGGAAATGCAGGCAACAGTACGCGAGCGTCAACGCCGTCTGCGATTTGtgctgc                                                                                                 >  2:4059669/1‑150 (MQ=255)
                                                         cccGTTGTAATGACCGAACAACAGCGTGATATGTCCGGCGAAGGTATCACGACGGGATACTACCTGCGCATCGGTCACGCCACGGCGAATATAGGGAAATGCAGGCAACAGTACGCGAGCGTCAACGCCGTCTGCGATTTGTGCTGCggg                                                                                              >  2:357317/1‑150 (MQ=255)
                                                          ccGTTGTAATGACCGAACAACAGCGTGATATGTCCGGCGAAGGTATCACGACGGGATACTACCTGCGCATCGGTCACGCCACGGCGAATATCGGGAAATGCAGGCAACAGTACGCGAGCGTCAACGCCGTCTGCGATTTGTGCTGCGGGt                                                                                             <  1:4110315/150‑1 (MQ=255)
                                                            gTTGTAATGACCGAACAACAGCGTGATATGTCCGGCGAAGGTATCACGACGGGATACTACCTGCGCATCGGTCACGCCACGGCGAATATCGGGAAATGCAGGCAACAGTACGCGAGCGTCAACGCCGTCTGCGATTTGTGCTGCGGGTaa                                                                                           <  1:3654743/150‑1 (MQ=255)
                                                             ttGTAATGACCGAACAACAGCGTGATATGTCCGGCGAAGGTATCACGACGGGATACTACCTGCGCATCGGTCACGCCACGGCGAATATAGGGAAATGCAGGCAACAGTACGCGAGCGTCAACGCCGTCTGCGATTTGTGCTGCGGGTAAt                                                                                          <  1:2254449/150‑1 (MQ=255)
                                                              tGTAATGACCGAACAACAGCGTGATATGTCCGGCGAAGGTATCACGACGGGATACTACCTGCGCATCGGTCACGCCACGGCGAATATAGGGAAATGCAGGCAACAGTACGCGAGCGTCAACGCCGTCTGCGATTTGTGCTGCGGGTAATg                                                                                         >  2:3087700/1‑150 (MQ=255)
                                                                   tGACCGAACAACAGCGTGATATGTCCGGCGAAGGTATCACGACGGGATACTACCTGCGCATCGGTCACGCCACGGCGAATATCGGGAAATGCAGGCAACAGTACGCGAGCGTCAACGCCGTCTGCGATTTGTGCTGCGGGTAATGCCCCa                                                                                    <  1:986189/150‑1 (MQ=255)
                                                                   tGACCGAACAACAGCGTGATATGTCCGGCGAAGGTATCACGACGGGATACTACCTGCGCATCGGTCACGCCACGGCGAATATAGGGAAATGCAGGCAACAGTACGCGAGCGTCAACGCCGTCTGCGATTTGTGCTGCGGGTAATGCCCCa                                                                                    >  2:4438596/1‑150 (MQ=255)
                                                                      ccGAACAACAGCGTGATATGTCCGGCGAAGGTATCACGACGGGATACTACCTGCGCATCGGTCACGCCACGGCGAATATCGGGAAATGCAGGCAACAGTACGCGAGCGTCAACGCCGTCTGCGATTTGTGCTGCGGGTAATGCCCCAATa                                                                                 <  1:106490/150‑1 (MQ=255)
                                                                       cGAACAACAGCGTGATATGTCCGGCGAAGGTATCACGACGGGATACTACCTGCGCATCGGTCACGCCACGGCGTATATCGGGAAATGCAGGCAACTGTACGCGAGCGTCAACGCCGTCTGCGATTTGTGCTGCGGGTAATGCCCCAATaa                                                                                <  1:650865/150‑1 (MQ=255)
                                                                        gAACAACAGCGTGATATGTCCGGCGAAGGTATCACGACGGGATACTACCTGCGCATCGGTCACGCCACGGCGAATATCGGGAAATGCAGGCAACAGTACGCGAGCGTCAACGCCGTCTGCGATTTGTGCTGCGGGTAATGCCCCAATAAc                                                                               >  2:3352842/1‑150 (MQ=255)
                                                                        gAACAACAGCGTGATATGTCCGGCGAAGGTATCACGACGGGATACTACCTGCGCATCGGTCACGCCACGGCGAATATCGGGAAATGCAGGCAACAGTACGCGAGCGTCAACGCCGTCTGCGATTTGTGCTGCGGGTAATGCCCCAATAAc                                                                               >  2:2094248/1‑150 (MQ=255)
                                                                          acaacaGCGTGATATGTCCGGCGAAGGTATCACGACGGGATACTACCTGCGCATCGGTCACGCCACGGCGAATATAGGGAAATGCAGGCAACAGTACGCGAGCGTCAACGCCGTCTGCGATTTGTGCTGCGGGTAATGCCCCAATAACAt                                                                             >  1:3896125/1‑150 (MQ=255)
                                                                            aacaGCGTGATATGTCCGGCGAAGGTATCACGACGGGATACTACCTGCGCATCGGTCACGCCACGGCGAATATCGGGAAATGCAGGCAACAGTACGCGAGCGTCAACGCCGTCTGCGATTTGTGCTGCGGGTAATGCCCCAATAACATCa                                                                           <  1:1691/150‑1 (MQ=255)
                                                                             acaGCGTGATATGTCCGGCGAAGGTATCACGACGGGATACTACCTGCGCATCGGTCACGCCACGGCGAATATCGGGAAATGCAGGCAACAGTACGCGAGCGTCAACGCCGTCTGCGATTTGTGCTGCGGGTAATGCCCCAATAACATCAg                                                                          >  1:2173527/1‑150 (MQ=255)
                                                                              caGCGTGATATGTCCGGCGAAGGTATCACGACGGGATACTACCTGCGCATCGGTCACGCCACGGCGAATATAGGGAAATGCAGGCAACAGTACGCGAGCGTCAACGCCGTCTGCGATTTGTGCTGCGGGTAATGCCCCAATAACATCAGc                                                                         >  1:4265528/1‑150 (MQ=255)
                                                                               aGCGTGATATGTCCGGCGAAGGTATCACGACGGGATACTACCTGCGCATCGGTCACGCCACGGCGAATATAGGGAAATGCAGGCAACAGTACGCGAGCGTCAACGCCGTCTGCGATTTGTGCTGCGGGTAATGCCCCAATAACATCAGcc                                                                        <  2:1743329/150‑1 (MQ=255)
                                                                                     atatGTCCGGCGAAGGTATCACGACGGGATACTACCTGCGCATCGGTCACGCCACGGCGAATATCGGGAAATGCAGGCAACAGTACGCGAGCGTCAACGCCGTCTGCGATTTGTGCTGCGGGTAATGCCCCAATAACATCAGCCAGAccg                                                                  <  1:465763/150‑1 (MQ=255)
                                                                                       atGTCCGGCGAAGGTATCACGACGGGATACTACCTGCGCATCGGTCACGCCACGGCGAATATCGGGAAATGCAGGCAACAGTACGCGAGCGTCAACGCCGTCTGCGATTTGTGCTGCGGGTAATGCCCCAATAACATCAGCCAGAccgcc                                                                >  1:2674890/1‑150 (MQ=255)
                                                                                        tGTCCGGCGAAGGTATCACGACGGGATACTACCTGCGCATCGGTCACGCCACGGCGAATATCGGGAAATGCAGGCAACAGTACGCGAGCGTCAACGCCGTCTGCGATTTGTGCTGCGGGTAATGCCCCAATAACATCAGCCAGAccgccg                                                               <  2:4290797/150‑1 (MQ=255)
                                                                                        tGTCCGGCGAAGGTATCACGACGGGATACTACCTGCGCATCGGTCACGCCACGGCGAATATAGGGAAATGCAGGCAACAGTACGCGAGCGTCAACGCCGTCTGCGATTTGTGCTGCGGGTAATGCCCCAATAACATCAGCCAGAccgccg                                                               >  2:3968928/1‑150 (MQ=255)
                                                                                             ggCGAAGGTATCACGACGGGATACTACCTGCGCATCGGTCACGCCACGGCGAATATCGGGAAATGCAGGCAACAGTACGCGAGCGTCAACGCCGTCTGCGATTTGTGCTGCGGGTAATGCCCCAATAACATCAGCCAGACCGCCGGtttt                                                          <  2:184193/150‑1 (MQ=255)
                                                                                             ggCGAAGGTATCACGACGGGATACTACCTGCGCATCGGTCACGCCACGGCGAATATCGGGAAATGCAGGCAACAGTACGCGAGCGTCAACGCCGTCTGCGATTTGTGCTGCGGGTAATGCCCCAATAACATCAGCCAGACCGCCGGtttt                                                          >  2:2510696/1‑150 (MQ=255)
                                                                                             ggCGAAGGTATCACGACGGGATACTACCTGCGCATCGGTCACGCCACGGCGAATATCGGGAAATGCAGGCAACAGTACGCGAGCGTCAACGCCGTCTGCGATTTGTGCTGCGGGTAATGCCCCAATAACATCAGCCAGACCGCCGGtttt                                                          <  1:1801098/150‑1 (MQ=255)
                                                                                             ggCGAAGGTATCACGACGGGATACTACCTGCGCATCGGTCACGCCACGGCGAATATAGGGAAATGCAGGCAACAGTACGCGAGCGTCAACGCCGTCTGAGATTTGTGCTGCGGGTAATGCCCCAATAACATCAGACATACCGCCGGtttt                                                          >  2:1654343/1‑150 (MQ=255)
                                                                                               cGAAGGTATCACGACGGGATACTACCTGCGCATCGGTCACGCCACGGCGAATATCGGGAAATGCAGGCAACAGTACGCGAGCGTCAACGCCGTCTGCGATTTGTGCTGCGGGTAATGCCCCAATAACATCAGCCAGACCGCCGGTTTTaa                                                        <  1:638061/150‑1 (MQ=255)
                                                                                               cGAAGGTATCACGACGGGATACTACCTGCGCATCGGTCACGCCACGGCGAATATCGGGAAATGCAGGCAACAGTACGCGAGCGTCAACGCCGTCTGCGATTTGTGCTGCGGGTAATGCCCCAATAACATCAGCCAGACCGCCGGTTTTaa                                                        >  2:4163113/1‑150 (MQ=255)
                                                                                                gAAGGTATCACGACGGGATACTACCTGCGCATCGGTCACGCCACGGCGAATATCGGGAAATGCAGGCAACAGTACGCGAGCGTCAACGCCGTCTGCGATTTGTGCTGCGGGTAATGCCCCAATAACATCAGCCAGACCGCCGGTTTTAag                                                       <  1:1260861/150‑1 (MQ=255)
                                                                                                  aGGTATCACGACGGGATACTACCTGCGCATCGGTCACGCCACGGCGAATATCGGGAAATGCAGGCAACAGTACGCGAGCGTCAACGCCGTCTGCGATTTGTGCTGCGGGTACTGCCCCAATAACATCAGCCAGACCGCCGGTTTTAagca                                                     >  2:1992055/1‑150 (MQ=255)
                                                                                                     tATCACGACGGGATACTACCTGCGCATCGGTCACGCCACGGCGAATATCGGGAAATGCAGGCAACAGTACGCGAGCGTCAACGCCGTCTGCGATTTGTGCTGCGGGTAATGCCCCAATAACATCAGCCAGACCGCCGGTTTTAAGCAGCg                                                  >  1:2456325/1‑150 (MQ=255)
                                                                                                      aTCACGACGGGATACTACCTGCGCATCGGTCACGCCACGGCGAATATCGGGAAATGCAGGCAACAGTACGCGAGCGTCAACGCCGTCTGCGATTTGTGCTGCGGGTAATGCCCCAATAACATCAGCCAGACCGCCGGTTTTAAGCAGCgg                                                 >  1:4410098/1‑150 (MQ=255)
                                                                                                      aTCACGACGGGATACTACCTGCGCATCGGTCACGCCACGGCGAATATCGGGAAATGCAGGCAACAGTACGCGAGCGTCAACGCCGTCTGCGATTTGTGCTGCGGGTAATGCCCCAATAACATCAGCCAGACCGCCGGTTTTAAGCAGCgg                                                 >  1:587025/1‑150 (MQ=255)
                                                                                                        cACGACGGGATACTACCTGCGCATCGGTCACGCCACGGCGAATATCGGGAAATGCAGGCAACAGTACGCGAGCGTCAACGCCGTCTGCGATTTGTGCTGCGGGTAATGCCCCAATAACATCAGCCAGACCGCCGGTTTTAAGCAGCGGGa                                               <  1:2450709/150‑1 (MQ=255)
                                                                                                          cgacgGGATACTACCTGCGCATCGGTCACGCCACGGCGAATATCGGGAAATGCAGGCAACAGTACGCGAGCGTCAACGCCGTCTGCGATTTGTGCTGCGGGTAATGCCCCAATAACATCAGCCAGACCGCCGGTTTTAAGCAGCGGGAAc                                             <  1:520101/150‑1 (MQ=255)
                                                                                                          cgacgGGATACTACCTGCGCATCGGTCACGCCACGGCGAATATAGGGAAATGCAGGCAACAGTACGCGAGCGTCAACGCCGTCTGCGATTTGTGCTGCGGGTAATGCCCCAATAACATCAGCCAGACCGCCGGTTTTAAGCAGCGGGAAc                                             <  1:2605053/150‑1 (MQ=255)
                                                                                                          cgacgGGATACTACCTGCGCATCGGTCACGCCACGGCGAATATAGGGAAATGCAGGCAACAGTACGCGAGCGTCAACGCCGTCTGCGATTTGTGCTGCGGGTAATGCCCCAATAACATCAGCCAGACCGCCGGTTTTAAGCAGCGGGAAc                                             <  1:2164379/150‑1 (MQ=255)
                                                                                                              gggATACTACCTGCGCATCGGTCACGCCACGGCGAATATAGGGAAATGCAGGCAACAGTACGCGAGCGTCAACGCCGTCTGCGATTTGTGCTGCGGGTAATGCCCCAATAACATCAGCCAGACCGCCGGTTTTAAGCAGCGGGAACAtct                                         >  2:4572466/1‑150 (MQ=255)
                                                                                                                gATACTACCTGCGCATCGGTCACGCCACGGCGAATATAGGGAAATGCAGGCAACAGTACGCGAGCGTCAACGCCGTCTGCGATTTGTGCTGCGGGTAATGCCCCAATAACATCAGCCAGACCGCCGGTTTTAAGCAGCGGGAACATctct                                       <  1:726446/150‑1 (MQ=255)
                                                                                                                 aTACTACCTGCGCATCGGTCACGCCACGGCGAATATCGGGAAATGCAGGCAACAGTACGCGAGCGTCAACGCCGTCTGCGATTTGTGCTGCGGGTAATGCCCCAATAACATCAGCCAGACCGCCGGTTTTAAGCAGCGGGAACATCTCTg                                      <  2:350715/150‑1 (MQ=255)
                                                                                                                   actacCTGCGCATCGGTCACGCCACGGCGAATATAGGGAAATGCAGGCAACAGTACGCGAGCGTCAACGCCGTCTGCGATTTGTGCTGCGGGTAATGCCCCAATAACATCAGCCAGACCGCCGGTTTTAAGCAGCGGGAACATCTCTGaa                                    >  2:3433417/1‑150 (MQ=255)
                                                                                                                     tacCTGCGCATCGGTCACGCCACGGCGAATATCGGGAAATGCAGGCAACAGTACGCGAGCGTCAACGCCGTCTGCGATTTGTGCTGCGGGTAATGCCCCAATAACATCAGCCAGACCGCCGGTTTTAAGCAGCGGGAACATCTCTGAaca                                  <  2:1132348/150‑1 (MQ=255)
                                                                                                                     tacCTGCGCATCGGTCACGCCACGGCGAATATAGGGAAATGCAGGCAACAGTACGCGAGCGTCAACGCCGTCTGCGATTTGTGCTGCGGGTAATGCCCCAATAACATCAGCCAGACCGCCGGTTTTAAGCAGCGGGAACATCTCTGAaca                                  >  1:1599503/1‑150 (MQ=255)
                                                                                                                      acCTGCGCATCGGTCACGCCACGGCGAATATCGGGAAATGCAGGCAACAGTACGCGAGCGTCAACGCCGTCTGCGATTTGTGCTGCGGGTAATGCCCCAATAACATCAGCCAGACCGCCGGTTTTAAGCAGCGGGAACATCTCTGAacat                                 <  2:4559545/150‑1 (MQ=255)
                                                                                                                       ccTGCGCATCGGTCACGCCACGGCGAATATCGGGAAATGCAGGCAACAGTACGCGAGCGTCAACGCCGTCTGCGATTTGTGCTGCGGGTAATGCCCCAATAACATCAGCCAGACCGCCGGTTTTAAGCAGCGGGAACATCTCTGAacata                                >  2:3827458/1‑150 (MQ=255)
                                                                                                                       ccTGCGCATCGGTCACGCCACGGCGAATATCGGGAAATGCAGGCAACAGTACGCGAGCGTCAACGCCGTCTGCGATTTGTGCTGCGGGTAATGCCCCAATAACATCAGCCAGACCGCCGGTTTTAAGCAGCGGGAACATCTCTGAacata                                >  2:47101/1‑150 (MQ=255)
                                                                                                                         tGCGCATCGGTCACGCCACGGCGAATATAGGGAAATGCAGGCAACAGTACGCGAGCGTCAACGCCGTCTGCGATTTGTGCTGCGGGTAATGCCCCAATAACATCAGCCAGACCGCCGGTTTTAAGCAGCGGGAACATCTCTGAacataca                              >  2:4188320/1‑150 (MQ=255)
                                                                                                                           cgcATCGGTCACGCCACGGCGAATATCGGGAAATGCAGGCAACAGTACGCGAGCGTCAACGCCGTCTGCGATTTGTGCTGCGGGTAATGCCCCAATAACATCAGCCAGACCGCCGGTTTTAAGCAGCGGGAACATCTCTGAACATACATg                            <  2:2072069/150‑1 (MQ=255)
                                                                                                                                cGGTCACGCCACGGCGAATATCGGGAAATGCAGGCAACAGTACGCGAGCGTCAACGCCGTCTGCGATTTGTGCTGCGGGTAATGCCCCAATAACATCAGCCAGACCGCCGGTTTTAAGCAGCGGGAACATCTCTGAACATACATGTaaaa                       >  2:64068/1‑150 (MQ=255)
                                                                                                                                  gTCACGCCACGGCGAATATCGGGAAATGCAGGCAACAGTACGCGAGCGTCAACGCCGTCTGCGATTTGTGCTGCGGGTAATGCCCCAATAACATCAGCCAGACCGCCGGTTTTAAGCAGCGGGAACATCTCTGAACATACATGTAAAAcc                     <  1:83023/150‑1 (MQ=255)
                                                                                                                                       gCCACGGCGAATATAGGGAAATGCAGGCAACAGTACGCGAGCGTCAACGCCGTCTGCGATTTGTGCTGCGGGTAATGCCCCAATAACATCAGCCAGACCGCCGGTTTTAAGCAGCGGGAACATCTCTGAACATACCTGTAAAACCTGCAt                <  2:4498189/150‑1 (MQ=255)
                                                                                                                                       gCCACGGCGAATATAGGGAAATGCAGGCAACAGTACGCGAGCGTCAACGCCGTCTGCGATTTGTGCTGCGGGTAATGCCCCAATAACATCAGCCAGACCGCCGGTTTTAAGCAGCGGGAACATCTCTGAACATACATGTAAAACCTGCAt                >  1:4462840/1‑150 (MQ=255)
                                                                                                                                         cACGGCGAATATAGGGAAATGCAGGCAACAGTACGCGAGCGTCAACGCCGTCTGCGATTTGTGCTGCGGGTAATGCCCCAATAACATCAGCCAGACCGCCGGTTTTAAGCAGCGGGAACATCTCTGAACATACATGTAAAACCTGCATTa              >  2:2136517/1‑150 (MQ=255)
                                                                                                                                              cGAATATAGGGAAATGCAGGCAACAGTACGCGAGCGTCAACGCCGTCTGCGATTTGTGCTGCGGGTAATGCCCCAATAACATCAGCCAGACCGCCGGTTTTAAGCAGCGGGAACATCTCTGAACATACATGTAAAACCTGCATTATCGCt         <  2:110325/150‑1 (MQ=255)
                                                                                                                                                 atatAGGGAAATGCAGGCAACAGTACGCGAGCGTCAACGCCGTCTGCGATTTGTGCTGCGGGTAATGCCCCAATAACATCAGCCAGACCGCCGGTTTTAAGCAGCGGGAACATCTCTGAACATACATGTAAAACCTGCATTATCGCTCCt      >  1:4059920/1‑150 (MQ=255)
                                                                                                                                                   atCGGGAAATGCAGGCAACAGTTCGCGAGCGTCAACGCCGTCTGCGATTTGTGCTGCGGGTAATGCCCCTATTACATCAGGCAGACCGCCGGTTTTCAGCAGCGGGAACATCTCTGTTCATACATGTAAAACCTGCTTTAGCGCTGCTGt    <  2:1605033/150‑1 (MQ=255)
                                                                                                                                                    tCGGGAAATGCAGGCAACAGTACGCGAGCGTCAACGCCGTCTGCGATTTGTGCTGCGGGTAATGCCCCAATAACATCAGCCAGACCGCCGGTTTTAAGCAGCGGGAACATCTCTGAACATACATGTAAAACCTGCATTATCGCTCCTGtt   <  2:3648922/150‑1 (MQ=255)
                                                                                                                                                    tCGGGAAATGCAGGCAACAGTACGCGAGCGTCAACGCCGTCTGCGATTTGTGCTGCGGGTAATGCCCCAATAACATCAGCCAGACCGCCGGTTTTAAGCAGCGGGAACATCTCTGAACATACATGTAAAACCTGCATTATCGCTCCTGtt   <  2:1552633/150‑1 (MQ=255)
                                                                                                                                                    tCGGGAAATGCAGGCAACAGTACGCGAGCGTCAACGCCGTCTGCGATTTGTGCTGCGGGTAATGCCCCAAGAACATCAGCCAGACCGCCGGTTTTAAGCAGCGGGAACATCTCTGAACATACATGTAAAACCTGCATTATCGCTCCTGtt   <  2:1552452/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ATACGGGCTTCCCGGACGATCATAGAGATGCGGCGCGTCAATCAGGTAAATGCCAACCCCGTTGTAATGACCGAACAACAGCGTGATATGTCCGGCGAAGGTATCACGACGGGATACTACCTGCGCATCGGTCACGCCACGGCGAATATCGGGAAATGCAGGCAACAGTACGCGAGCGTCAACGCCGTCTGCGATTTGTGCTGCGGGTAATGCCCCAATAACATCAGCCAGACCGCCGGTTTTAAGCAGCGGGAACATCTCTGAACATACATGTAAAACCTGCATTATCGCTCCTGTTT  >  CP009273/3561109‑3561407

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: