Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 3,259,207 G→A A33V (GCG→GTG)  glmM ← phosphoglucosamine mutase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6063,259,2070GA96.3% 198.0 / ‑2.8 54A33V (GCG→GTG) glmMphosphoglucosamine mutase
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base A (27/25);  total (28/26)
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.94e-01

ATATAGCCAGAAATACGCGTGTCTTTACCAATAATGATCTTACGGGAGCCGTGGCGCGCCAGCACTTTACCCGCGGCCCAACCCAGCTTAAGCACAAAATCAGGTGTGATCGGCGCATCCCCTACACGACCACGAATCCCATC  >  REL606/3259135‑3259277
                                                                        |                                                                      
atatAGCCAGAAATACGCGTGTCTTTACCAATAATGATCTTACGGGAGCCGTGGCGCGCCAGCACTTTACCcacg                                                                      >  1:392352/1‑75 (MQ=255)
atatAGCCAGAAATACGCGTGTCTTTACCAATAATGATCTTACGGGAGCCGTGGCGCGCCAGCACTTTACCcacg                                                                      >  2:628039/1‑75 (MQ=255)
 tataGCCAGAAATACGCGTGTCTTTACCAATAATGATCTTACGGGAGCCGTGGCGCGCCAGCACTTTACCCACgg                                                                     <  2:1128952/75‑1 (MQ=255)
 tataGCCAGAAATACGCGTGTCTTTACCAATAATGATCTTACGGGAGCCGTGGCGCGCCAGCACTTTACCCACgg                                                                     >  2:1171209/1‑75 (MQ=255)
 tataGCCAGAAATACGCGTGTCTTTACCAATAATGATCTTACGGGAGCCGTGGCGCGCCAGCACTTTACCCACgg                                                                     >  2:1462269/1‑75 (MQ=255)
 tataGCCAGAAATACGCGTGTCTTTACCAATAATGATCTTACGGGAGCCGTGGCGCGCCAGCACTTTACCCACgg                                                                     <  2:820490/75‑1 (MQ=255)
 tataGCCAGAAATACGCGTGTCTTTACCAATAATGATCTTACGGGAGCCGTGGCGCGCCAGCACTTTACCCACgg                                                                     >  1:742042/1‑75 (MQ=255)
    aGCCAGAAATACGCGTGTCTTTACCAATAATGATCTTACGGGAGCCGTGGCGCGCCAGCACTTTACCCACGGccc                                                                  >  1:1159877/1‑75 (MQ=255)
    aGCCAGAAATACGCGTGTCTTTACCAATAATGATCTTACGGGAGCCGTGGCGCGCCAGCACTTTACCCACGGccc                                                                  >  1:439850/1‑75 (MQ=255)
     gCCAGAAATACGCGTGTCTTTACCAATAATGATCTTACGGGAGCCGTGGCGCGCCAGCACTTTACCCGCGGCCCa                                                                 >  1:1357850/1‑75 (MQ=255)
      ccAGAAATCCGCGTGTCTTTCCAAATAATGTTCTTCCGGGAGCCGTGGCGCGCCAGCACTTTACCCACGGCCCaa                                                                <  1:893568/75‑1 (MQ=255)
        aGAAATACGCGTGTCTTTACCAATAATGATCTTACGGGAGCCGTGGCGCGCCAGCACTTTACCCACGGCCCAAcc                                                              >  2:14982/1‑75 (MQ=255)
        aGAAATACGCGTGTCTTTACCAATAATGATCTTACGGGAGCCGTGGCGCGCCAGCACTTTACCCACGGCCCAAcc                                                              >  1:1242394/1‑75 (MQ=255)
        aGAAATACGCGTGTCTTTACCAATAATGATCTTACGGGAGCCGTGGCGCGCCAGCACTTTACCCACGGCCCAAcc                                                              >  1:1220719/1‑75 (MQ=255)
        aGAAATACGCGTGTCTTTACCAATAATGATCTTACGGGAGCCGTGGCGCGCCAGCACTTTACCCACGGCCCAAcc                                                              >  1:1164048/1‑75 (MQ=255)
            aTACGCGTGTCTTTACCAATAATGATCTTACGGGAGCCGTGGCGCGCCAGCACTTTACCCACGGCCCAACCCAGc                                                          >  1:332161/1‑75 (MQ=255)
             tACGCGTGTCTTTACCAATAATGATCTTACGGGAGCCGTGGCGCGCCAGCACTTTACCCACGGCCCAACCCAGCt                                                         >  2:338682/1‑75 (MQ=255)
              aCGCGTGTCTTTACCAATAATGATCTTACGGGAGCCGTGGCGCGCCAGCACTTTACCCACGGCCCAACCCCGCtt                                                        >  2:1155131/1‑75 (MQ=255)
              aCGCGTGTCTTTACCAATAATGATCTTACGGGAGCCGTGGCGCGCCAGCACTTTACCCACGGCCCAACCCACCtt                                                        >  2:1507915/1‑75 (MQ=255)
                gcgTGTCTTTACCAATAATGATCTTACGGGAGCCGTGGCGCGCCAGCACTTTACCCACGGCCCAACCCAGCTTaa                                                      <  1:881521/75‑1 (MQ=255)
                 cgTGTCTTTACCAATAATGATCTTACGGGAGCCGTGGCGCGCCAGCACTTTACCCACGGCCCAACCCAGCTTAAg                                                     <  2:462579/75‑1 (MQ=255)
                  gtgtCTTTACCAATAATGATCTTACGGGAGCCGTGGCGCGCCAGCACTTTACCCACGGCCCAACCCAGCTTAAGc                                                    <  2:852087/75‑1 (MQ=255)
                   tgtCTTTACCAATAATGATCTTACGGGAGCCGTGGCGCGCCAGCACTTTACCCACGGCCCAACCCAGCTTAAGca                                                   >  1:587720/1‑75 (MQ=255)
                     tCTTTACCAATAATGATCTTACGGGAGCCGTGGCGCGCCAGCACTTTACCCACGGCCCAACCCAGCTTAAGcaca                                                 >  1:1131089/1‑75 (MQ=255)
                     tCTTTACCAATAATGATCTTACGGGAGCCGTGGCGCGCCAGCACTTTACCCACGGCCCAACCCAGCTTAAGcaca                                                 >  2:526773/1‑75 (MQ=255)
                      cTTTACCAATAATGATCTTACGGGAGCCGTGGCGCGCCAGCACTTTACCCACGGCCCAACCCAGCTTAAGCACaa                                                <  1:1287611/75‑1 (MQ=255)
                        ttACCAATAATGATCTTACGGGAGCCGTGGCGCGCCAGCACTTTACCCACGGCCCAACCCAGCTTAAGCACaaaa                                              >  1:1345345/1‑75 (MQ=255)
                            cAATAATGATCTTACGGGAGCCGTGGCGCGCCAGCACTTTACCCACGGCCCAACCCAGCTTAAGCACAAAATCAg                                          <  2:1360641/75‑1 (MQ=255)
                             aataatGATCTTACGGGAGCCGTGGCGCGCCAGCACTTTACCCACGGCCCAACCCAGCTTAAGCACAAAATCAgg                                         <  1:1462270/75‑1 (MQ=255)
                             aataatGATCTTACGGGAGCCGTGGCGCGCCAGCACTTTACCCACGGCCCAACCCAGCTTAAGCACAAAATCAgg                                         >  1:1582133/1‑75 (MQ=255)
                              ataatGATCTTACGGGAGCCGTGGCGCGCCAGCACTTTAGCCACGGCCCAACCCAGCtt                                                        >  2:304768/1‑59 (MQ=255)
                              ataatGATCTTACGGGAGCCGTGGCGCGCCAGCACTTTACCCACGGCCCAACCCAGCtt                                                        <  1:304768/59‑1 (MQ=255)
                                aatGATCTTACGGGAGCCGTGGCGCGCCAGCACTTTACCCGCGGCCCAACCCAGCTTAAGCACAAAATCAGgtgt                                      <  1:645924/75‑1 (MQ=255)
                                aatGATCTTACGGGAGCCGTGGCGCGCCAGCACTTTACCCACGGCCCAACCCAGCTTAAGCACAAAATCAGgtgt                                      <  1:1333844/75‑1 (MQ=255)
                                aatGATCTTACGGGAGCCGTGGCGCGCCAGCACTTTACCCACGGCCCAACCCAGCTTAAGCACAAAATCAGgtgt                                      >  2:1573082/1‑75 (MQ=255)
                                 atGATCTTACGGGAGCCGTGGCGCGCCAGCACTTTACCCACGGCCCAACCCAGCTTAAGCACAAAATCAGGtgtg                                     <  1:515217/75‑1 (MQ=255)
                                  tGATCTTACGGGAGCCGTGGCGCGCCAGCACTTTACCCACGGCCCAACCCAGCTTAAGCACAAAATCAGGTGTGa                                    <  1:338682/75‑1 (MQ=255)
                                   gATCTTACGGGAGCCGTGGCGCGCCAGCACTTTACCCACGGCCCAACCCAGCTTAAGCACAAAATCAGGTGTGAt                                   >  1:908327/1‑75 (MQ=255)
                                     tCTTACGGGAGCCGTGGCGCGCCAGCACTTTACCCACGGTCCAACCCAGCTTAAGCACAAAATCAGGTGTGATCg                                 <  2:1323282/75‑1 (MQ=255)
                                      cTTACGGGAGCCGTGGCGCGCCAGCACTTTACCCACGGCCCAACCCAGCTTAAGCACAAAATCAGGTGTGATCgg                                >  1:1299228/1‑75 (MQ=255)
                                        tACGGGAGCCGTGGCGCGCCAGCACTTTACCCACGGCCCAACCCAGCTTAAGCACAAAATCAGGTGTGATCGgcg                              >  1:1474165/1‑75 (MQ=255)
                                           gggAGCCGTGGCGCGCCAGCACTTTACCCACGGCCCAACCCAGCTTAAGCACAAAATCAGGTGTGATCGGCGCAt                           >  1:487222/1‑75 (MQ=255)
                                            ggAGCCGTGGCGCGCCAGCACTTTACCCACGGCCCAACCCAGCTTAAGCACAAAATCAGGTGTGATCGGCGCATc                          >  2:524476/1‑75 (MQ=255)
                                             gAGCCGTGGCGCGCCAGCACTTTACCCACGGCCCAACCCAGCTTAAGCACAAAATCAGGTGTGATCGGCGCATcc                         <  2:1035980/75‑1 (MQ=255)
                                               gCCGTGGCGCGCCAGCACTTTACCCACGGCCCAACCCAGCTTAAGCACAAAATCAGGTGTGATCGGCGCATcccc                       <  1:1398066/75‑1 (MQ=255)
                                               gCCGTGGCGCGCCAGCACTTTACCCACGGCCCAACCCAGCTTAAGCACAAAATCAGGTGTGATCGGCGCATcccc                       >  2:1241956/1‑75 (MQ=255)
                                                  gTGGCGCGCCAGCACTTTACCCACGGCCCAACCCAGCTTAAGCACAAAATCAGGTGTGATCGGCGCATCCCCTac                    <  2:112418/75‑1 (MQ=255)
                                                  gTGGCGCGCCAGCACTTTACCCACGGCCCAACCCAGCTTAAGCACAAAATCAGGTGTGATCGGCGCATCCCCTac                    >  1:316547/1‑75 (MQ=255)
                                                    ggCGCGCCAGCACTTTACCCACGGCCCAACCCAGCTTAAGCACAAAATCAGGTGTGATCGGCGCATCCCCTacac                  <  1:1495635/75‑1 (MQ=255)
                                                      cgcgcCAGCACTTTACCCACGGCCCAACCCAGCTTAAGCACAAAATCAGGTGTGATCGGCGCATCCCCTACACGa                <  1:188162/75‑1 (MQ=255)
                                                          ccAGCACTTTACCCACGGCCCAACCCAGCTTAAGCACAAAATCAGGTGTGATCGGCGCATCCCCTACACGACCAc            <  2:1535295/75‑1 (MQ=255)
                                                          ccAGCACTTTACCCACGGCCCAACCCAGCTTAAGCACAAAATCAGGTGTGATCGGCGCATCCCCTACACGACCAc            <  2:1158987/75‑1 (MQ=255)
                                                               aCTTTACCCACGGCCCAACCCAGCTTAAGCACAAAATCAGGTGTGATCGGCGCATCCCCTACACGACCACGAATc       <  1:1608809/75‑1 (MQ=255)
                                                                cTTTACCCACGGCCCAACCCAGCTTAAGCACAAAATCAGGTGTGATCGGCGCCTCCCCTACACGACCACGAATcc      <  2:1315879/75‑1 (MQ=255)
                                                                cTTTACCCACGGCCCAACCCAGCTTAAGCACAAAATCAGGTGTGATCGGCGCATCCCCTACACGACCACGAATcc      <  1:100556/75‑1 (MQ=255)
                                                                  ttACCCACGGCCCAACCCAGCTTAAGCACAAAATCAGGTGTGATCGGCGCATCCCCTACACGACCACGAATCCCa    <  2:12435/75‑1 (MQ=255)
                                                                   tACCCACGGCCCAACCCAGCTTAAGCACAAAATCAGGTGTGATCGGCGCATCCCCTACACGACCACGAATCCCAt   <  1:429282/75‑1 (MQ=255)
                                                                   tACCCACGGCCCAACCCAGCTTAAGCACAAAATCAGGTGTGATCGGCGCATCCCCTACACGACCACGAATCCCAt   <  1:1405231/75‑1 (MQ=255)
                                                                   tACCCACGGCCCAACCCAGCTTAAGCACAAAATCAGGTGTGATCGGCGCATCCCCTACACGACCACGAATCCCAt   >  1:167262/1‑75 (MQ=255)
                                                                    cccccCGGCCCAACCCAGCTTAAGCACAAAATCAGGTGTGATCGGCGCATCCCCTACACGACCACGAATCCCATc  <  1:698254/74‑1 (MQ=255)
                                                                        |                                                                      
ATATAGCCAGAAATACGCGTGTCTTTACCAATAATGATCTTACGGGAGCCGTGGCGCGCCAGCACTTTACCCGCGGCCCAACCCAGCTTAAGCACAAAATCAGGTGTGATCGGCGCATCCCCTACACGACCACGAATCCCATC  >  REL606/3259135‑3259277

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: