Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,931,021 C→A 100% A123E (GCG→GAG)  rbsK → ribokinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,931,0210CA97.1% 739.5 / 11.3 239A123E (GCG→GAG) rbsKribokinase
Reads supporting (aligned to +/- strand):  ref base C (4/3);  new base A (143/89);  total (147/92)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 244 total aligned reads displayed.

CAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGCGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCAAAATAAGACTATCGTTGCGCTTAACCCGGCTCCGGCTC  >  CP009273/3930876‑3931164
                                                                                                                                                 |                                                                                                                                               
cagcagCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGgagc                                                                                                                                               <  1:2803181/148‑1 (MQ=255)
cagcagCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGgagc                                                                                                                                               <  1:3459576/148‑1 (MQ=255)
 agcagcTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCt                                                                                                                                              >  2:1932040/1‑148 (MQ=255)
   cagcTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCtgg                                                                                                                                            <  2:2055934/148‑1 (MQ=255)
     gcTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCtggtg                                                                                                                                          >  1:2701711/1‑148 (MQ=255)
     gcTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGGGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCtggtg                                                                                                                                          >  1:1585476/1‑148 (MQ=255)
        cGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGaa                                                                                                                                       >  2:1925765/1‑148 (MQ=255)
           cACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAgcg                                                                                                                                    <  2:3450623/148‑1 (MQ=255)
            aCTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCGGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAgcgc                                                                                                                                   >  2:958501/1‑148 (MQ=255)
            aCTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCGGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAgcgc                                                                                                                                   >  2:958492/1‑148 (MQ=255)
            aCTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAgcgc                                                                                                                                   >  2:2268093/1‑148 (MQ=255)
                aTAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACg                                                                                                                               <  1:1293448/148‑1 (MQ=255)
                aTAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACg                                                                                                                               >  2:2605110/1‑148 (MQ=255)
                  aaCATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTg                                                                                                                             >  2:3236283/1‑148 (MQ=255)
                   aCATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGa                                                                                                                            >  1:670554/1‑148 (MQ=255)
                   aCATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGa                                                                                                                            <  1:924442/148‑1 (MQ=255)
                   aCATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGa                                                                                                                            <  1:2598094/148‑1 (MQ=255)
                   aCATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGa                                                                                                                            >  1:1994653/1‑148 (MQ=255)
                    cATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAg                                                                                                                           >  1:2247564/1‑148 (MQ=255)
                    cATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAg                                                                                                                           >  2:320835/1‑148 (MQ=255)
                     aTTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGc                                                                                                                          >  1:3320605/1‑148 (MQ=255)
                      ttGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCg                                                                                                                         >  2:1099494/1‑148 (MQ=255)
                       tGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGt                                                                                                                        >  1:964558/1‑148 (MQ=255)
                        gATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTa                                                                                                                       <  1:2661576/148‑1 (MQ=255)
                         atatTACTCCGGTCAGCGTGATCAAAGGCTAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTAt                                                                                                                      >  1:938842/1‑148 (MQ=255)
                              aCTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCa                                                                                                                 <  1:4632/148‑1 (MQ=255)
                              aCTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCa                                                                                                                 <  1:467741/148‑1 (MQ=255)
                               cTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCaa                                                                                                                <  2:2212466/148‑1 (MQ=255)
                               cTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCaa                                                                                                                >  1:1138436/1‑148 (MQ=255)
                                tCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAAc                                                                                                               >  1:2544247/1‑148 (MQ=255)
                                     tCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTc                                                                                                          <  1:3034928/148‑1 (MQ=255)
                                      cAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCa                                                                                                         <  2:284532/148‑1 (MQ=255)
                                      cAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCa                                                                                                         >  1:2030217/1‑148 (MQ=255)
                                        gCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGc                                                                                                       <  2:1563614/148‑1 (MQ=255)
                                        gCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGc                                                                                                       <  2:721990/148‑1 (MQ=255)
                                         cGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCa                                                                                                      >  2:1379705/1‑148 (MQ=255)
                                          gTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCat                                                                                                     <  1:1663961/148‑1 (MQ=255)
                                            gATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTACTGGCGCAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCatta                                                                                                   <  2:443476/148‑1 (MQ=255)
                                               cAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCAttatta                                                                                                <  1:2951471/148‑1 (MQ=255)
                                                aaaGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGACTGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTaa                                                                                               <  1:1908655/148‑1 (MQ=255)
                                                    gCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCa                                                                                           >  2:2946991/1‑148 (MQ=255)
                                                    gCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCa                                                                                           >  2:881183/1‑148 (MQ=255)
                                                      gAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGc                                                                                         <  1:2138291/148‑1 (MQ=255)
                                                          cAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGa                                                                                     >  1:1025640/1‑148 (MQ=255)
                                                           aaCAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGaa                                                                                    <  1:815647/148‑1 (MQ=255)
                                                             cAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATc                                                                                  >  1:103692/1‑148 (MQ=255)
                                                               ggTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATcac                                                                                >  1:2012737/1‑148 (MQ=255)
                                                                gtgtGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATTCAGCTGGAATcacc                                                                               >  1:1867204/1‑148 (MQ=255)
                                                                gtgtGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATcacc                                                                               >  2:594335/1‑148 (MQ=255)
                                                                gtgtGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATcacc                                                                               >  1:2354972/1‑148 (MQ=255)
                                                                gtgtGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATcacc                                                                               >  1:1088602/1‑148 (MQ=255)
                                                                gtgtGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATcacc                                                                               <  1:2463090/148‑1 (MQ=255)
                                                                 tgtgGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGATGGAATcacca                                                                              <  1:1963853/148‑1 (MQ=255)
                                                                        cTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGaaa                                                                       >  1:149413/1‑148 (MQ=255)
                                                                         tGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAg                                                                      >  1:354259/1‑148 (MQ=255)
                                                                          gATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAgt                                                                     >  1:3425425/1‑148 (MQ=255)
                                                                          gATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAgt                                                                     >  1:427315/1‑148 (MQ=255)
                                                                           aTTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAgtg                                                                    >  2:3039268/1‑148 (MQ=255)
                                                                             ttttGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGtgtg                                                                  <  2:1512005/148‑1 (MQ=255)
                                                                              tttGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGa                                                                 >  1:2639795/1‑148 (MQ=255)
                                                                                tGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATg                                                               >  1:611770/1‑148 (MQ=255)
                                                                                  ttAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGc                                                             <  1:3306620/148‑1 (MQ=255)
                                                                                   tAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCa                                                            >  2:1763856/1‑148 (MQ=255)
                                                                                    aaTGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAg                                                           <  1:320835/148‑1 (MQ=255)
                                                                                     aTGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAgc                                                          >  1:1173612/1‑148 (MQ=255)
                                                                                      tGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAgcg                                                         >  2:1215616/1‑148 (MQ=255)
                                                                                      tGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAgcg                                                         >  2:3006965/1‑148 (MQ=255)
                                                                                      tGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAgcg                                                         <  2:2362188/148‑1 (MQ=255)
                                                                                        gCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAgcggc                                                       >  1:1697503/1‑148 (MQ=255)
                                                                                            aGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGaaa                                                   <  2:249752/148‑1 (MQ=255)
                                                                                            aGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGaaa                                                   <  1:3082439/148‑1 (MQ=255)
                                                                                             ggTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGaaaa                                                  <  1:420665/148‑1 (MQ=255)
                                                                                             ggTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGaaaa                                                  >  1:2352581/1‑148 (MQ=255)
                                                                                             ggTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGaaaa                                                  <  2:1648533/148‑1 (MQ=255)
                                                                                              gTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAAt                                                 >  1:243576/1‑148 (MQ=255)
                                                                                                gagaATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCg                                               >  2:3320867/1‑148 (MQ=255)
                                                                                                 agaATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGc                                              >  1:1533158/1‑148 (MQ=255)
                                                                                                   aaTGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGccc                                            <  2:3100304/148‑1 (MQ=255)
                                                                                                   aaTGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGccc                                            <  2:1164284/148‑1 (MQ=255)
                                                                                                   aaTGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGccc                                            <  2:2666838/148‑1 (MQ=255)
                                                                                                        cATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCaa                                       >  1:2225570/1‑148 (MQ=255)
                                                                                                           cGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGTATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCAAAAt                                    >  1:1059508/1‑148 (MQ=255)
                                                                                                           cGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCAAAAt                                    <  2:1194397/148‑1 (MQ=255)
                                                                                                              tATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCAAAATAAg                                 >  2:2106569/1‑148 (MQ=255)
                                                                                                               aTTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCAAAATAAGa                                >  1:2837745/1‑148 (MQ=255)
                                                                                                                 tCATGCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCAAAATAAGACt                              >  2:842154/1‑148 (MQ=255)
                                                                                                                     gCCGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCAAAATAAGACTATCg                          >  1:3356173/1‑148 (MQ=255)
                                                                                                                      ccGGCGCTAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCAAAATAAGACTATCGt                         <  2:3320604/148‑1 (MQ=255)
                                                                                                                             tAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCAAAATAAGACTATCGTTGCGCtt                  <  2:2114365/148‑1 (MQ=255)
                                                                                                                             tAATGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCAAAATAAGACTATCGTTGCGCtt                  >  2:1064910/1‑148 (MQ=255)
                                                                                                                              aaTGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCAAAATAAGACTATCGTTGCGCTTa                 <  2:3352925/148‑1 (MQ=255)
                                                                                                                              aaTGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCAAAATAAGACTATCGTTGCGCTTa                 <  2:1585476/148‑1 (MQ=255)
                                                                                                                               aTGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCAAAATAAGACTATCGTTGCGCTTaa                >  2:1839880/1‑148 (MQ=255)
                                                                                                                               aTGCTGCCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGAAGGCAGCGGCGAAAATCGCCCATCAAAATAAGACTATCGTTGCGCTTaa                >  2:811372/1‑148 (MQ=255)
                                                                                                                                 gctgcCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCAAAATAAGACTATCGTTGCGCTTAAcc              <  2:2806176/148‑1 (MQ=255)
                                                                                                                                  ctgcCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCAAAATAAGACTATCGTTGCGCTTAAccc             >  2:3563977/1‑148 (MQ=255)
                                                                                                                                  ctgcCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCAAAATAAGACTATCGTTGCGCTTAAccc             >  1:714875/1‑148 (MQ=255)
                                                                                                                                  ctgcCCTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCAAAATAAGACTATCGTTGCGCTTAAccc             >  2:2894820/1‑148 (MQ=255)
                                                                                                                                     cccTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCAAAATAAGACTATCGTTGCGCTTAACccggc          >  2:258028/1‑148 (MQ=255)
                                                                                                                                       cTTTCCCCGGAGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCAAAATAAGACTATCGTTGCGCTTAACccggctc        <  2:1728207/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                               
CAGCAGCTCGCCACTGATAACATTGATATTACTCCGGTCAGCGTGATCAAAGGCGAATCAACAGGTGTGGCGCTGATTTTTGTTAATGGCGAAGGTGAGAATGTCATCGGTATTCATGCCGGCGCTAATGCTGCCCTTTCCCCGGCGCTGGTGGAAGCGCAACGTGAGCGTATTGCCAACGCGTCAGCATTATTAATGCAGCTGGAATCACCACTCGAAAGTGTGATGGCAGCGGCGAAAATCGCCCATCAAAATAAGACTATCGTTGCGCTTAACCCGGCTCCGGCTC  >  CP009273/3930876‑3931164

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: