Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,116,216 G→T 5.8% R317S (CGT→AGT)  glcB ← malate synthase G

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,116,2160GT 5.8% 1321.4 / 39.4 346R317S (CGT→AGT) glcBmalate synthase G
Reads supporting (aligned to +/- strand):  ref base G (147/178);  new base T (15/6);  total (162/184)
Fisher's exact test for biased strand distribution p-value = 2.38e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 353 total aligned reads displayed.

CATGACGCCATCAAGAATGCCTTCCGGGATTTCATTGCCTTCGCTGTCCCAAATCACAGGAATGGTCATCAAATGACCCACGTTGCGGATAAACAGCAGCGAGCGTCCGTGCAGAGAAATTTCAGAGCCATCGGCGGCGGTGTAATGACGATCGTCATTCAGTTTACGCACGATTTGCCGACCGTTTTTCTCCATTTTCTCTTGCAGAGTCCCCTGCATCAGGCCCAGCAGGTTGCGGTACAGCAGGATTTTATCTTCCGCATCAACCGCCGCGACCGAATCTTCGCAGTCGAGAATGG  >  CP009273/3116067‑3116365
                                                                                                                                                     |                                                                                                                                                     
cATGACGCCATCAAGAATGCCTTCCGGGATTTCATTGCCTTCGCTGTCCCAAATCACAGGAATGGTCATCAAATGACCCACGTTGCGGATAAACAGCAGCGAGCGTCCGTGCAGAGAAATTTCAGAGCCATCGGCGGCGGTGTAATGACg                                                                                                                                                       <  2:3509875/150‑1 (MQ=255)
  tGACGCCATCAAGAATGCCTTCCGGGATTTCATTGCCTTCGCTGTCCCAAATCACAGGAATGGTCATCAAATGACCCACGTTGCGGATAAACAGCAGCGAGCGTCCGTGCAGAGAAATTTCAGAGCCATCGGCGGCGGTGTAATGACGAt                                                                                                                                                     <  2:79428/150‑1 (MQ=255)
   gACGCCATCAAGAATGCCTTCCGGGATTTCATTGCCTTCGCTGTCCCAAATCACAGGAATGGTCATCAAATGACCCACGTTGCGGATAAACAGCAGCGAGCGTCCGTGCAGAGAAATTTCAGAGCCATCGGCGGCGGTGTAATGACGATc                                                                                                                                                    >  1:3877227/1‑150 (MQ=255)
    aCGCCATCAAGAATGCCTTCCGGGATTTCATTGCCTTCGCTGTCCCAAATCACAGGAATGGTCATCAAATGACCCACGTTGCGGATAAACAGCAGCGAGCGTCCGTGCAGAGAAATTTCAGAGCCATCGGCGGCGGTGTAATGACGATCg                                                                                                                                                   <  1:1904909/150‑1 (MQ=255)
    aCGCCATCAAGAATGCCTTCCGGGATTTCATTGCCTTCGCTGTCCCAAATCACAGGAATGGTCATCAAATGACCCACGTTGCGGATAAACAGCAGCGAGCGTCCGTGCAGAGAAATTTCAGAGCCATCGGCGGCGGTGTAATGACGAGCg                                                                                                                                                   >  2:4328318/1‑150 (MQ=255)
      gCCATCAAGAATGCCTTCCGGGATTTCATTGCCTTCGCTGTCCCAAATCACAGGAATGGTCATCAAATGACCCACGTTGCGGATAAACAGCAGCGAGCGTCCGTGCAGAGAAATTTCAGAGCCATCGGCGGCGGTGTAATGACGATCGTc                                                                                                                                                 <  2:4616428/150‑1 (MQ=255)
       ccATCAAGAATGCCTTCCGGGATTTCATTGCCTTCGCTGTCCCAAATCACAGGAATGGTCATCAAATGACCCACGTTGCGGATAAACAGCAGCGAGCGTCCGTGCAGAGAAATTTCAGAGCCATCGGCGGCGGTGTAATGACGATCGTCa                                                                                                                                                >  2:1487696/1‑150 (MQ=255)
        cATCAAGAATGCCTTCCGGGATTTCATTGCCTTCGCTGTCCCAAATCACAGGAATGGTCATCAAATGACCCACGTTGCGGATAAACAGCAGCGAGTGTCCGTGCAGAGAAATTTCAGAGCCATCGGCGGCGGTGTAATGACGATCGTCAt                                                                                                                                               >  2:2818617/1‑150 (MQ=255)
         aTCAAGAATGCCTTCCGGGATTTCATTGCCTTCGCTGTCCCAAATCACAGGAATGGTCATCAAATGACCCACGTTGCGGATAAACAGCAGCGAGCGTCCGTGCAGAGAAATTTCAGAGCCATCGGCGGCGGTGTAATGACGATCGTCAtt                                                                                                                                              >  1:1811632/1‑150 (MQ=255)
         aTCAAGAATGCCTTCCGGGATTTCATTGCCTTCGCTGTCCCAAATCACAGGAATGGTCATCAAATGACCCACGTTGCGGATAAACAGCAGCGAGCGTCCGTGCAGAGAAATTTCAGAGCCATCGGCGGCGGTGTAATGACGATCGTCAtt                                                                                                                                              >  1:5652426/1‑150 (MQ=255)
         aTCAAGAATGCCTTCCGGGATTTCATTGCCTTCGCTGTCCCAAATCACAGGAATGGTCATCAAATGACCCACGTTGCGGATAAACAGCAGCGAGCGTCCGTGCAGAGAAATTTCAGAGCCATCGGCGGCGGTGTAATGACGATCGTCAtt                                                                                                                                              <  1:101446/150‑1 (MQ=255)
         aTCAAGAATGCCTTCCGGGATTTCATTGCCTTCGCTGTCCCAAATCACAGGAATGGTCATCAAATGACCCACGTTGCGGATAAACAGCAGCGAGCGTCCGTGCAGAGAAATTTCAGAGCCATCGGCGGCGGTGTAATGACGATCGTCAtt                                                                                                                                              >  1:2850314/1‑150 (MQ=255)
            aaGAATGCCTTCCGGGATTTCATTGCCTTCGCTGTCCCAAATCACAGGAATGGTCATCAAATGACCCACGTTGCGGATAAACAGCAGCGAGCGTCCGTGCAGAGAAATTTCAGAGCCATCGGCGGCGGTGTAATGACGATCGTCATTCAg                                                                                                                                           <  1:3982371/150‑1 (MQ=255)
               aaTGCCTTCCGGGATTTCATTGCCTTCGCTGTCCCAAATCACAGGAATGGTCATCAAATGACCCACGTTGCGGATAAACAGCAGCGAGCGTCCGTGCAGAGAAATTTCAGAGCCATCGGCGGCGGTGTAATGACGATCGTCATTCAGttt                                                                                                                                        <  2:406122/150‑1 (MQ=255)
                 tGCCTTCCGGGATTTCATTGCCTTCGCTGTCCCAAATCACAGGAATGGTCATCAAATGACCCACGTTGCGGATAAACAGCAGCGAGCGTCCGTGCAGAGAAATTTCAGAGCCATCGGCGGCGGTGTAATGACGATCGTCATTCAGTTTAc                                                                                                                                      <  1:1368690/150‑1 (MQ=255)
                 tGCCTTCCGGGATTTCATTGCCTTCGCTGTCCCAAATCACAGGAATGGTCATCAAATGACCCACGTTGCGGATAAACAGCAGCGAGCGTCCGTGCAGAGAAATTTCAGAGCCATCGGCGGCGGTGTAATGACGATCGTCATTCAGTTTAc                                                                                                                                      <  1:1238847/150‑1 (MQ=255)
                     ttCCGGGATTTCATTGCCTTCGCTGTCCCAAATCACAGGAATGGTCATCAAATGACCCACGTTGCGGATAAACAGCAGCGAGCGTCCGTGCAGAGAAATTTCAGAGCCATCGGCGGCGGTGTAATGACGATCGTCATTCAGTTTACGCAc                                                                                                                                  <  2:5661586/150‑1 (MQ=255)
                     ttCCGGGATTTCATTGCCTTCGCTGTCCCAAATCACAGGAATGGTCATCAAATGACCCACGTTGCGGATAAACAGCAGCGAGCGTCCGTGCAGAGAAATTTCAGAGCCATCGGCGGCGGTGTAATGACGATCGTCATTCAGTTTACGCAc                                                                                                                                  >  1:5070715/1‑150 (MQ=255)
                     ttCCGGGATTTCATTGCCTTCGCTGTCCCAAATCACAGGAATGGTCATCAAATGACCCACGTTGCGGATAAACAGCAGCGAGCGTCCGTGCAGAGAAATTTCAGAGCCATCGGCGGCGGTGTAATGACGATCGTCATTCAGTTTACGCAc                                                                                                                                  <  1:3191179/150‑1 (MQ=255)
                        cGGGATTTCATTGCCTTCGCTGTCCCAAATCACAGGAATGGTCATCAAATGACCCACGTTGCGGATAAACAGCAGCGAGCGTCCGTGCAGAGAAATTTCAGAGCCATCGGCGGCGGTGTAATGACGATCGTCATTCAGTTTACGCACGAt                                                                                                                               <  1:1153400/150‑1 (MQ=255)
                        cGGGATTTCATTGCCTTCGCTGTCCCAAATCACAGGAATGGTCATCAAATGACCCACGTTGCGGATAAACAGCAGCGAGCGTCCGTGCAGAGAAATTTCAGAGCCATCGGCGGCGGTGTAATGACGATCGTCATTCAGTTTACGCACGAt                                                                                                                               <  1:4026310/150‑1 (MQ=255)
                        cGGGATTTCATTGCCTTCGCTGTCCCAAATCACAGGAATGGTCATCAAATGACCCACGTTGCGGATAAACAGCAGCGAGCGTCCGTGCAGAGAAATTTCAGAGCCATCGGCGGCGGTGTAATGACGATCGTCATTCAGTTTACGCACGAt                                                                                                                               <  1:1531740/150‑1 (MQ=255)
                               tCATTGCCTTCGCTGTCCCAAATCACAGGAATGGTCATCAAATGACCCACGTTGCGGATAAACAGCAGCGAGCGTCCGTGCAGAGAAATTTCAGAGCCATCGGCGGCGGTGTAATGACGATCGTCATTCAGTTTACGCACGATTTGCCGa                                                                                                                        <  1:4553453/150‑1 (MQ=255)
                                 aTTGCCTTCGCTGTCCCAAATCACAGGAATGGTCATCAAATGACCCATGTTGCGGATAAACAGCAGCGAGCGTCCGTGCAGAGAAATTTCAGAGCCATCGGCGGCGGTGTAATGACGATCGTCATTCAGTTTACGCACGATTTGCCGAcc                                                                                                                      <  2:4253319/150‑1 (MQ=255)
                                  ttGCCTTCGCTGTCCCAAATCACAGGAATGGTCATCAAATGACCCACGTTGCGGATAAACAGCAGCGAGCGTCCGTGCAGAGAAATTTCAGAGCCATCGGCGGCGGTGTAATGACGATCGTCATTCAGTTTACGCACGATTTGCCGACCg                                                                                                                     >  1:5233294/1‑150 (MQ=255)
                                   tGCCTTCGCTGTCCCAAATCACAGGAATGGTCATCAAATGACCCACGTTGCGGATAAACAGCAGCGAGCGTCCGTGCAGAGAAATTTCAGAGCCATCGGCGGCGGTGTAATGACGATCGTCATTCAGTTTACGCACGATTTGCCGACCGt                                                                                                                    <  2:2147159/150‑1 (MQ=255)
                                   tGCCTTCGCTGTCCCAAATCACAGGAATGGTCATCAAATGACCCACGTTGCGGATAAACAGCAGCGAGCGTCCGTGCAGAGAAATTTCAGAGCCATCGGCGGCGGTGTAATGACGATCGTCATTCAGTTTACGCACGATTTGCCGACCGt                                                                                                                    >  1:2295421/1‑150 (MQ=255)
                                    gCCTTCGCTGTCCCAAATCACAGGAATGGTCATCAAATGACCCACGTTGCGGATAAACAGCAGCGAGCGTCCGTGCAGAGAAATTTCAGAGCCATCGGCGGCGGTGTAATGACGATCGTCATTCAGTTTACGCACGATTTGCCGACCGtt                                                                                                                   >  1:2392334/1‑150 (MQ=255)
                                     ccTTCGCTGTCCCAAATCACAGGAATGGTCATCAAATGACCCACGTTGCGGATAAACAGCAGCGAGCGTCCGTGCAGAGAAATTTCAGAGCCATCGGCGGCGGTGTAATGACGATCGTCATTCAGTTTACGCACGATTTGCCGACCGttt                                                                                                                  <  2:1160250/150‑1 (MQ=255)
                                      cTTCGCTGTCCCAAATCACAGGAATGGTCATCAAATGACCCACGTTGCGGATAAACAGCAGCGAGCGTCCGTGCAGAGAAATTTCAGAGCCATCGGCGGCGGTGTAATGACGATCGTCATTCAGTTTACGCACGATTTGCCGACCGtttt                                                                                                                 <  1:3678597/150‑1 (MQ=255)
                                        tctcTGTCCCAAAACAAAGGAAAGGTCAACAAATGACCCACGTTTCGGATAAACAGCAGCGAGCGTCCGTGAAGAGAAATTTCAGAGCCATAGGCGGCGTAGTAATGACGATCGAAAATCAGTTTACTCACGATATGCCGACCGTTTTtc                                                                                                               >  2:1904343/4‑150 (MQ=255)
                                        tCGCTGTCCCAAATCACAGGAATGGTCATCAAATGACCCACGTTGCGGATAAACAGCAGCGAGCGTCCGTGCAGAGAAATTTCAGAGCCATCGGCGGCGGTGTAATGACGATCGTCATTCAGTTTACGCACGATTTGCCGACCGTTTTtc                                                                                                               >  2:5573179/1‑150 (MQ=255)
                                         cGCTGTCCCAAATCACAGGAATGGTCATCAAATGACCCACGTTGCGGATAAACAGCAGCGAGCGTCCGTGCAGAGAAATTTCAGAGCCATCGGCGGCGGTGTAATGACGATCGTCATTCAGTTTACGCACGATTTGCCGACCGTTTTtct                                                                                                              <  2:1679150/150‑1 (MQ=255)
                                             gTCCCAAATCACAGGAATGGTCATCAAATGACCCACGTTGCGGATAAACAGCAGCGAGCGTCCGTGCAGAGAAATTTCAGAGCCATCGGCGGCGGTGTAATGACGATCGTCATTCAGTTTACGCACGATTTGCCGACCGTTTTTCTCCAt                                                                                                          >  1:4407403/1‑150 (MQ=255)
                                              tCCCAAATCACAGGAATGGTCATCAAATGACCCACGTTGCGGATAAACAGCAGCGAGCGTCCGTGCAGAGAAATTTCAGAGCCATCGGCGGCGGTGTAATGACGATCGTCATTCAGTTTACGCACGATTTGCCGACCGTTTTTCTCCAtt                                                                                                         <  1:3309578/150‑1 (MQ=255)
                                               cccAAATCACAGGAATGGTCATCAATTGACCCACGTTGCGGATAAACAGCAGCGAGCGTCCGTGCAGAGAAATTTCAGAGCCATCGGCGGCGGTGTAATGTCGATCGTCATTCAGTTTACGCACGATTTGCAGACCGTTTTTCTCCAttt                                                                                                        <  2:3308874/150‑1 (MQ=255)
                                                ccAAATCACAGGAATGGTCATCAAATGACCCACGTTGCGGATAAACAGCAGCGAGCGTCCGTGCAGAGAAATTTCAGAGCCATCGGCGGCGGTGTAATGACGATCGTCATTCAGTTTACGCACGATTTGCCGACCGTTTTTCTCCAtttt                                                                                                       <  1:5423267/150‑1 (MQ=255)
                                                ccAAATCACAGGAATGGTCATCAAATGACCCACGTTGCGGATAAACAGCAGCGAGCGTCCGTGCAGAGAAATTTCAGAGCCATCGGCGGCGGTGTAATGACGATCGTCATTCAGTTTACGCACGATTTGCCGACCGTTTTTCTCCAtttt                                                                                                       >  1:5128836/1‑150 (MQ=255)
                                                 cAAATCACAGGAATGGTCATCAAATGACCCACGTTGCGGATAAACAGCAGCGAGCGTCCGTGCAGAGAAATTTCAGAGCCATCGGCGGCGGTGTAATGACGATCGTCATTCAGTTTACGCACGATTTGCCGACCGTTTTTCTCCATTTtc                                                                                                      <  1:2966618/150‑1 (MQ=255)
                                                 cAAATCACAGGAATGGTCATCAAATGACCCACGTTGCGGATAAACAGCAGCGAGCGTCCGTGCAGAGAAATTTCAGAGCCATCGGCGGCGGTGTAATGACGATCGTCATTCAGTTTACGCACGATTTGCCGACCGTTTTTCTCCAATTtc                                                                                                      >  2:4525324/1‑150 (MQ=255)
                                                  aaaTCACAGGAATGGTCATCAAATGACCCACGTTGCGGATAAACAGCAGCGAGCGTCCGTGCAGAGAAATTTCAGAGCCATCGGCGGCGGTGTAATGACGATCGTCATTCAGTTTACGCACGATTTGCCGACCGTTTTTCTCCATTTtct                                                                                                     >  2:1275850/1‑150 (MQ=255)
                                                    aTCACAGGAATGGTCATCAAATGACCCACGTTGCGGATAAACAGCAGCGAGCGTCCGTGCAGAGAAATTTCAGAGCCATCGGCGGCGGTGTAATGACGATCGTCATTCAGTTTACGCACGATTTGCCGACCGTTTTTCTCCATTTTctct                                                                                                   <  1:5205962/150‑1 (MQ=255)
                                                     tCACAGGAATGGTCATCAAATGACCCACGTTGCGGATAAACAGCAGCGAGCGTCCGTGCAGAGAAATTTCAGAGCCATCGGCGGCGGTGTAATGACGATCGTCATTCAGTTTACGCACGATTTGCCGACCGTTTTTCTCCATTTTCTCtt                                                                                                  <  2:4031693/150‑1 (MQ=255)
                                                     tCACAGGAATGGTCATCAAATGACCCACGTTGCGGATAAACAGCAGCGAGCGTCCGTGCAGAGAAATTTCAGAGCCATCGGCGGCGGTGTAATGACGATCGTCATTCAGTTTACGCACGATTTGCCGACCGTTTTTCTCCATTTTCTCtt                                                                                                  <  1:2245298/150‑1 (MQ=255)
                                                        caGGAATGGTCATCAAATGACCCACGTTGCGGATAAACAGCAGCGAGCGTCCGTGCAGAGAAATTTCAGAGCCATCGGCGGCGGTGTAATGACGATCGTCATTCAGTTTACGCACGATTTGCCGACCGTTTTTCTCCATTTTCTCTTGCa                                                                                               <  2:734554/150‑1 (MQ=255)
                                                         aGGAATGGTCATCAAATGACCCACGTTGCGGATAAACAGCAGCGAGCGTCCGTGCAGAGAAATTTCAGAGCCATCGGCGGCGGTGTAATGACGATCGTCATTCAGTTTACGCACGATTTGCCGACCGTTTTTCTCCATTTTCTCTTGCag                                                                                              >  1:3717662/1‑150 (MQ=255)
                                                          ggAATGGTCATCAAATGACCCACGTTGCGGATAAACAGCAGCGAGCGTCCGTGCAGAGAAATTTCAGAGCCATCGGCGGCGGTGTAATGACGATCGTCATTCAGTTTACGCACGATTTGCCGACCGTTTTTCTCCATTTTCTCTTGCaga                                                                                             <  1:1252295/150‑1 (MQ=255)
                                                             aTGGTCATCAAATGACCCACGTTGCGGATAAACAGCAGCGAGCGTCCGTGCAGAGAAATTTCAGAGCCATCGGCGGCGGTGTAATGACTATCGTCATTCAGTTTACGCACGATTTGCCGACCGTTTTTCTCCATTTTCTCTTGCAGAGTc                                                                                          <  1:1721016/150‑1 (MQ=255)
                                                                 tcatcaAATGACCCACGTTGCGGATAAACAGCAGCGAGCGTCCGTGCAGAGAAATTTCAGAGCCATCGGCGGCGGTGTAATGACTATCGTCATTCAGTTTACGCACGATTTGCCGACCGTTTTTCTCCATTTTCTCTTGCAGAGTCCCCt                                                                                      >  2:2892887/1‑150 (MQ=255)
                                                                 tcatcaAATGACCCACGTTGCGGATAAACAGCAGCGAGCGTCCGTGCAGAGAAATTTCAGAGCCATCGGCGGCGGTGTAATGACTATCGTCATTCAGTTTACGCACGATTTGCCGACCGTTTTTCTCCATTTTCTCTTGCAGAGTCCCCt                                                                                      >  2:1387242/1‑150 (MQ=255)
                                                                  catcaAATGACCCACGTTGCGGATAAACAGCAGCGAGCGTCCGTGCAGAGAAATTTCAGAGCCATCGGCGGCGGTGTAATGACGATCGTCATTCAGTTTACGCACGATTTGCCGACCGTTTTTCTCCATTTTCTCTTGCAGAGTCCCCTg                                                                                     >  2:1540423/1‑150 (MQ=255)
                                                                   atcaAATGACCCACGTTGCGGATAAACAGCAGCGAGCGTCCGTGCAGAGAAATTTCAGAGCCATCGGCGGCGGTGTAATGACGATCGTCATTCAGTTTACGCACGATTTGCCGACCGTTTTTCTCCATTTTCTCTTGCAGAGTCCCCTGc                                                                                    <  2:3066964/150‑1 (MQ=255)
                                                                     caAATGACCCACGTTGCGGATAAACAGCAGCGAGCGTCCGTGCAGAGAAATTTCAGAGCCATCGGCGGCGGTGTAATGACGATCGTCATTCAGTTTACGCACGATTTGCCGACCGTTTTTCTCCATTTTCTCTTGCAGAGTCCCCTGCAt                                                                                  >  2:2338299/1‑150 (MQ=255)
                                                                     caAATGACCCACGTTGCGGATAAACAGCAGCGAGCGTCCGTGCAGAGAAATTTCAGAGCCATCGGCGGCGGTGTAATGACGATCGTCATTCAGTTTACGCACGATTTGCCGACCGTTTTTCTCCATTTTCTCTTGCAGAGTCCCCTGCAt                                                                                  <  1:736453/150‑1 (MQ=255)
                                                                      aaaTGACCCACGTTGCGGATAAACAGCAGCGAGCGTCCGTGCAGAGAAATTTCAGAGCCATCGGCGGCGGTGTAATGACGATCGTCATTCAGTTTACGCACGATTTGCCGACCGTTTTTCTCCATATTCTCTTGAAGAGTCCCCTGCATc                                                                                 >  1:4954351/1‑150 (MQ=255)
                                                                        aTGACCCACGTTGCGGATAAACAGCAGCGAGTGTCCGTGCAGAGAAATTTCAGAGCCATCGGCGGCGGTGTAATGACGATCGTCATTCAGTTTACGCACGATTTGCCGACCGTTTTTCTCCATTTTCTCTTGCAGAGTCCCCTGCATCAg                                                                               >  1:5451969/1‑150 (MQ=255)
                                                                         tGACCCACGTTGCGGATAAACAGCAGCGAGCGTCCGTGCAGAGAAATTTCAGAGCCATCGGCGGCGGTGTAATGACGATCGTCATTCAGTTTACGCACGATTTGCCGACCGTTTTTCTCCATTTTCTCTTGCAGAGTCCCCTGCATCAgg                                                                              >  1:2728529/1‑150 (MQ=255)
                                                                             ccACGTTGCGGATAAACAGCAGCGAGCGTCCGTGCAGAGAAATTTCAGAGCCATCGGCGGCGGTGTAATGACTATCGTCATTCAGTTTACGCACGATTTGCCGACCGTTTTTCTCCATTTTCTCTTGCAGAGTCCCCTGCATCAGGCCca                                                                          >  2:874663/1‑150 (MQ=255)
                                                                              cACGTTGCGGATAAACAGCAGCGAGCGTCCGTGCAGAGAAATTTCAGAGCCATCGGCGGCGGTGTAATGACGATCGTCATTCAGTTTACGCACGATTTGCCGACCGTTTTTCTCCATTTTCTCTTGCAGAGTCCCCTGCATCAGGCCcag                                                                         <  1:949034/150‑1 (MQ=255)
                                                                              cACGTTGCGGATAAACAGCAGCGAGCGTCCGTGCAGAGAAATTTCAGAGCCATCGGCGGCGGTGTAATGACGATCGTCATTCAGTTTACGCACGATTTGCCGACCGTTTTTCTCCATTTTCTCTTGCAGAGTCCCCTGCATCAGGCCcag                                                                         <  2:3877227/150‑1 (MQ=255)
                                                                                  ttGCGGATAAACAGCAGCGAGCGTCCGTGCAGAGAAATTTCAGAGCCATCGGCGGCGGTGTAATGACGATCGTCATTCAGTTTACGCACGATTTGCCGACCGTTTTTCTCCATTTTCTCTTGCAGAGTCCCCTGCATCAGGCCCAGCAgg                                                                     <  1:2553873/150‑1 (MQ=255)
                                                                                   tGCGGATAAACAGCAGCGAGCGTCCGTGCAGAGAAATTTCAGAGCCATCGGCGGCGGTGTAATGACGATCGTCATTCAGTTTACGCACGATTTGCCGACCGTTTTTCTCCATTTTCTCTTGCAGAGTCCCCTGCATCAGGCCCAGCAGGt                                                                    <  1:2014172/150‑1 (MQ=255)
                                                                                    gCGGATAAACAGCAGCGAGCGTCCGTGCAGAGAAATTTCAGAGCCATCGGCGGCGGTGTAATGACGATCGTCATTCAGTTTACGCACGATTTGCCGACCGTTTTTCTCCATTTTATCTTGCAGAGTCCCCTGCATCAGGCCCAGCAGGtt                                                                   <  2:2219938/150‑1 (MQ=255)
                                                                                      ggATAAACAGCAGCGAGCGTCCGTGCAGAGAAATTTCAGAGCCATCGGCGGCGGTGTAATGACTATCGTCATTCAGTTTACGCACGATTTGCCGACCGTTTTTCTCCATTTTCTCTTGCAGAGTCCCCTGCATCAGGCCCAGCAGGTTGc                                                                 >  2:3572633/1‑150 (MQ=255)
                                                                                         tAAACAGCAGCGAGCGTCCGTGCAGAGAAATTTCAGAGCCATCGGCGGCGGTGTAATGACGATCGTCATTCAGTTTACGCACGATTTGCCGACCGTTTTTCTCCATTTTCTCTTGCAGAGTCCCCTGCATCAGGCCCAGCAGGTTGCGGt                                                              <  2:3675610/150‑1 (MQ=255)
                                                                                            aCAGCAGCGAGCGTCCGTGCAGAGAAATTTCAGAGCCATCGGCGGCGGTGTAATGACTATCGTCATTCAGTTTACGCACGATTTGCCGACCGTTTTTCTCCATTTTCTCTTGCAGAGTCCCCTGCATCAGGCCCAGCAGGTTGCGGTAca                                                           <  1:3124359/150‑1 (MQ=255)
                                                                                              agcagcGAGCGTCCGTGCAGAGAAATTTCAGAGCCATCGGCGGCGGTGTAATGACGATCGTCATTCAGTTTACGCACGATTTGCCGACCGTTTTTCTCCATTTTCTCTTGCAGAGTCCCCTGCATCAGGCCCAGCAGGTTGCGGTAcagc                                                         >  1:5351945/1‑150 (MQ=255)
                                                                                                cagcGAGCGTCCGTGCAGAGAAATTTCAGAGCCATCGGCGGCGGTGTAATGACGATCGTCATTCAGTTTACGCACGATTTGCCGACCGTTTTTCTCCATTTTCTCTTGCAGAGTCCCCTGCATCAGGCCCAGCAGGTTGCGGTAcagcag                                                       <  2:4657992/150‑1 (MQ=255)
                                                                                                 agcgagcgTCCGTGCAGAGAAATTTCAGAGCCATCGGCGGCGGTGTAATGACGATCGTCATTCAGTTTACGCACGATTTGCCGACCGTTTTTCTCCATTTTCTCTTGCAGAGTCCCCTGCATCAGGCCCAGCAGGTTGCGGTACAGCAgg                                                      <  1:4273464/150‑1 (MQ=255)
                                                                                                    gagcgTCCGTGCAGAGAAATTTCAGAGCCATCGGCGGCGGTGTAATGACGATCGTCATTCAGTTTACGCACGATTTGCCGACCGTTTTTCTCCATTTTCTCTTGCAGAGTCCCCTGCATCAGGCCCAGCAGGTTGCGGTACAGCAGGAtt                                                   <  1:2038738/150‑1 (MQ=255)
                                                                                                    gagcgTCCGTGCAGAGAAATTTCAGAGCCATCGGCGGCGGTGTAATGACGATCGTCATTCAGTTTACGCACGATTTGCCGACCGTTTTTCTCCATTTTCTCTTGCAGAGTCCCCTGCATCAGGCCCAGCAGGTTGCGGTACAGCAGGAtt                                                   <  1:4867333/150‑1 (MQ=255)
                                                                                                    gagcgTCCGTGCAGAGAAATTTCAGAGCCATCGGCGGCGGTGTAATGACGATCGTCATTCAGTTTACGCACGATTTGCCGACCGTTTTTCTCCATTTTCTCTTGCAGAGTCCCCTGCATCAGGCCCAGCAGGTTGCGGTACAGCAGGAtt                                                   >  1:3393854/1‑150 (MQ=255)
                                                                                                                 gagaAATTTCAGAGCCATCGGCGGCGGTGTAATGACGATCGTCATTCAGTTTACGCACGATTTGCCGACCGTTTTTCTCCATTTTCTCTTGCAGAGTCCCCTGCATCAGGCCCAGCAGGTTGCGGTACAGCAGGATTTTATCTTCCGCAt                                      <  2:1990806/150‑1 (MQ=255)
                                                                                                                     aaTTTCAGAGCCATCGGCGGCGGTGTAATGACGATCGTCATTCAGTTTACGCACGATTTGCCGACCGTTTTTCTCCATTTTCTCTTGCAGAGTCCCCTGCATCAGGCCCAGCAGGTTGCGGTACAGCAGGATTTTATCTTCCGCATCAAc                                  >  2:5332641/1‑150 (MQ=255)
                                                                                                                         tCAGAGCCATCGGCGGCGGTGTAATGACGATCGTCATTCAGTTTACGCACGATTTGCCGACCGTTTTTCTCCATTTTCTCTTGCAGAGTCCCCTGCATCAGGCCCAGCAGGTTGCGGTACAGCAGGATTTTATCTTCCGCATCAAccgcc                              <  2:1736220/150‑1 (MQ=255)
                                                                                                                          cAGAGCCATCGGCGGCGGTGTAATGACGATCGTCATTCAGTTTACGCACGATTTGCCGACCGTTTTTCTCCATTTTCTCTTGCAGAGTCCCCTGCATCAGGCCCAGCAGGTTGCGGTACAGCAGGATTTTATCTTCCGCATCAAccgccg                             >  1:3434609/1‑150 (MQ=255)
                                                                                                                           agagCCATCGGCGGCGGTGTAATGACGATCGTCATTCAGTTTACGCACGATTTGCCGACCGTTTTTCTCCATTTTCTCTTGCAGAGTCCCCTGCATCAGGCCCAGCAGGTTGCGGTACAGCAGGATTTTATCTTCCGCATCAACcgccgc                            >  2:5303905/1‑150 (MQ=255)
                                                                                                                            gagCCATCGGCGGCGGTGTAATGACGATCGTCATTCAGTTTACGCACGATTTGCCGACCGTTTTTCTCCATTTTCTCTTGCAGAGTCCCCTGCATCAGGCCCAGCAGGTTGCGGTACAGCAGGATTTTATCTTCCGCATCAACCGCcgcg                           >  2:1320583/1‑150 (MQ=255)
                                                                                                                               ccATCGGCGGCGGTGTAATGACGATCGTCATTCAGTTTACGCACGATTTGCCGACCGTTTTTCTCCATTTTCTCTTGCAGAGTCCCCTGCATCAGGCCCAGCAGGTTGCGGTACAGCAGGATTTTATCTTCCGCATCAACCGCCGCGAcc                        >  2:5043613/1‑150 (MQ=255)
                                                                                                                                 aTCGGCGGCGGTGTAATGACGATCGTCATTCAGTTTACGCACGATTTGCCGACCGTTTTTCTCCATTTTCTCTTGCAGAGTCCCCTGCATCAGGCCCAGCAGGTTGCGGTACAGCAGGATTTTATCTTCCGCATCAACCGCCGCGACCGa                      >  1:1843374/1‑150 (MQ=255)
                                                                                                                                  tCGGCGGCGGTGTAATGACGATCGTCATTCAGTTTACGCACGATTTGCCGACCGTTTTTCTCCATTTTCTCTTGCAGAGTCCCCTGCATCAGGCCCAGCAGGTTGCGGTACAGCAGGATTTTATCTTCCGCATCAACCGCCGCGACCGaa                     <  2:4789564/150‑1 (MQ=255)
                                                                                                                                    ggcggcggTGTAATGACGATCGTCATTCAGTTTACGCACGATTTGCCGACCGTTTTTCTCCATTTTCTCTTGCAGAGTCCCCTGCATCAGGCCCAGCAGGTTGCGGTACAGCAGGATTTTATCTTCCGCATCAACCGCCGCGACCGAATc                   >  2:3136955/1‑150 (MQ=255)
                                                                                                                                     gcggcggTGTAATGACGATCGTCATTCAGTTTACGCACGATTTGCCGACCGTTTTTCTCCATTTTCTCTTGCAGAGTCCCCTGCATCAGGCCCAGCAGGTTGCGGTACAGCAGGATTTTATCTTCCGCATCAACCGCCGCGACCGAATCt                  <  2:1071478/150‑1 (MQ=255)
                                                                                                                                      cggcggTGTAATGACGATCGTCATTCAGTTTACGCACGATTTGCCGACCGTTTTTCTCCATTTTCTCTTGCAGAGTCCCCTGCATCAGGCCCAGCAGGTTGCGGTACAGCAGGATTTTATCTTCCGCATCAACCGCCGCGACCGAATCtt                 <  1:3606517/150‑1 (MQ=255)
                                                                                                                                      cggcggTGTAATGACGATCGTCATTCAGTTTACGCACGATTTGCCGACCGTTTTTCTCCATTTTCTCTTGCAGAGTCCCCTGCATCAGGCCCAGCAGGTTGCGGTACAGCAGGATTTTATCTTCCGCATCAACCGCCGCGACCGAATCtt                 >  1:5032055/1‑150 (MQ=255)
                                                                                                                                      cggcggTGTAATGACGATCGTCATTCAGTTTACGCACGATTTGCCGACCGTTTTTCTCCATTTTCTCTTGCAGAGTCCCCTGCATCAGGCCCAGCAGGTTGCGGTACAGCAGGATTTTATCTTCCGCATCAACCGCCGCGACCGAATCtt                 <  2:2392334/150‑1 (MQ=255)
                                                                                                                                       ggcggTGTAATGACGATCGTCATTCAGTTTACGCACGATTTGCCGACCGTTTTTCTCCATTTTCTCTTGCAGAGTCCCCTGCATCAGGCCCAGCAGGTTGCGGTACAGCAGGATTTTATCTTCCGCATCAACCGCCGCGACCGAATCTTc                <  1:2162958/150‑1 (MQ=255)
                                                                                                                                       ggcggTGTAATGACGATCGTCATTCAGTTTACGCACGATTTGCCGACCGTTTTTCTCCATTTTCTCTTGCAGAGTCCCCTGCATCAGGCCCAGCAGGTTGCGGTACAGCAGGATTTTATCTTCCGCATCAACCGCCGCGACCGAATCTTc                >  2:5191361/1‑150 (MQ=255)
                                                                                                                                        gcggTGTAATGACGATCGTCATTCAGTTTACGCACGATTTGCCGACCGTTTTTCTCCATTTTCTCTTGCAGAGTCCCCTGCATCAGGCCCAGCAGGTTGCGGTACAGCAGGATTTTATCTTCCGCATCAACCGCCGCGACCGAATCTTCg               <  1:4614735/150‑1 (MQ=255)
                                                                                                                                        gcggTGTAATGACGATCGTCATTCAGTTTACGCACGATTTGCCGACCGTTTTTCTCCATTTTCTCTTGCAGAGTCCCCTGCATCAGGCCCAGCAGGTTGCGGTACAGCAGGATTTTATCTTCCGCATCAACCGCCGCGACCGAATCTTCg               <  2:1834726/150‑1 (MQ=255)
                                                                                                                                          ggTGTAATGACGATCGTCATTCAGTTTACGCACGATTTGCCGACCGTTTTTCTCCATTTTCTCTTGCAGAGTCCCCTGCATCAGGCCCAGCAGGTTGCGGTACAGCAGGATTTTATCTTCCGCATCAACCGCCGCGACCGAATCTTCGCa             >  1:4801373/1‑150 (MQ=255)
                                                                                                                                          ggTGTAATGACGATCGTCATTCAGTTTACGCACGATTTGCCGACCGTTTTTCTCCATTTTCTCTTGCAGAGTCCCCTGCATCAGGCCCAGCAGGTTGCGGTACAGCAGGATTTTATCTTCCGCATCAACCGCCGCGACCGAATCTTCGCa             >  2:1916444/1‑150 (MQ=255)
                                                                                                                                          ggTGTAATGACGATCGTCATTCAGTTTACGCACGATTTGCCGACCGTTTTTCTCCATTTTCTCTTGCAGAGTCCCCTGCATCAGGCCCAGCAGGTTGCGGTACAGCAGGATTTTATCTTCCGCATCAACCGCCGCGACCGAATCTTCGCa             <  2:2643513/150‑1 (MQ=255)
                                                                                                                                          ggTGTAATGACGATCGTCATTCAGTTTACGCACGATTTGCCGACCGTTTTTCTCCATTTTCTCTTGCAGAGTCCCCTGCATCAGGCCCAGCAGGTTGCGGTACAGCAGGATTTTATCTTCCGCATCAACCGCCGCGACCGAATCTTCGCa             <  2:2558366/150‑1 (MQ=255)
                                                                                                                                           gtgtAATGACGATCGTCATTCAGTTTACGCACGATTTGCCGTCCGTTTTTCTCCATTTTCTCTTGCAGAGTCCCCTGCATCAGGCCCAGCAGGTTGCGGTACAGCAGGATTTTATCTTCCGCATCAACCGCCGCGACCGAATCTTCGCAg            <  2:3360135/150‑1 (MQ=255)
                                                                                                                                             gtAATGACGATCGTCATTCAGTTTACGCACGATTTGCCGACCGTTTTTCTCCATTTTCTCTTGCAGAGTCCCCTGCATCAGGCCCAGCAGGTTGCGGTACAGCAGGATTTTATCTTCCGCATCAACCGCCGCGACCGAATCTTCGCAGTc          <  1:2817761/150‑1 (MQ=255)
                                                                                                                                               aaTGACGATCGTCATTCAGTTTACGCACGATTTGCCGACCGTTTTTCTCCATTTTCTCTTGCAGAGTCCCCTGCATCAGGCCCAGCAGGTTGCGGTACAGCAGGATTTTATCTTCCGCATCAACCGCCGCGACCGAATCTTCGCAGTCga        >  1:805710/1‑150 (MQ=255)
                                                                                                                                                aTGACGATCGTCATTCAGTTTACGCACGATTTGCCGACCGTTTTTCTCCATTTTCTCTTGCAGAGTCCCCTGCATCAGGCCCAGCAGGTTGCGGTACAGCAGGATTTTATCTTCCGCATCAACCGCCGCGACCGAATCTTCGCAGTCgag       <  1:3365610/150‑1 (MQ=255)
                                                                                                                                                 tGACGATCGTCATTCAGTTTACGCACGATTTGCCGACCGTTTTTCTCCATTTTCTCTTGCAGAGTCCCCTGCATCAGGCCCAGCAGGTTGCGGTACAGCAGGATTTTATCTTCCGCATCAACCGCCGCGACCGAATCTTCGCAGTCgaga      <  1:4132408/150‑1 (MQ=255)
                                                                                                                                                    cGATCGTCATTCAGTTTACGCACGATTTGCCGACCGTTTTTCTCCATTTTCTCTTGCAGAGTCCCCTGCATCAGGCCCAGCAGGTTGCGGTACAGCAGGATTTTATCTTCCGCATCAACCGCCGCGACCGAATCTTCGCAGTCGAGAATg   >  2:4180987/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CATGACGCCATCAAGAATGCCTTCCGGGATTTCATTGCCTTCGCTGTCCCAAATCACAGGAATGGTCATCAAATGACCCACGTTGCGGATAAACAGCAGCGAGCGTCCGTGCAGAGAAATTTCAGAGCCATCGGCGGCGGTGTAATGACGATCGTCATTCAGTTTACGCACGATTTGCCGACCGTTTTTCTCCATTTTCTCTTGCAGAGTCCCCTGCATCAGGCCCAGCAGGTTGCGGTACAGCAGGATTTTATCTTCCGCATCAACCGCCGCGACCGAATCTTCGCAGTCGAGAATGG  >  CP009273/3116067‑3116365

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: