Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 843,375 A→G 8.9% A80A (GCA→GCG galU → glucose‑1‑phosphate uridylyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE843,3750AG 8.9% 223.5 / 11.8 79A80A (GCA→GCGgalUglucose‑1‑phosphate uridylyltransferase
Reads supporting (aligned to +/- strand):  ref base A (19/53);  new base G (4/3);  total (23/56)
Fisher's exact test for biased strand distribution p-value = 1.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.40e-01

GTGCTGGTTACACACTCATCTAAAAACTCTATTGAAAACCACTTTGATACCAGTTTTGAACTGGAAGCAATGCTGGAAAAACGTGTAAAACGTCAACTGCTTGATGAAGTGCAGTCTATTTGTCCACCGCA  >  minE/843307‑843437
                                                                    |                                                              
gtgCTGGTTACACACTCATCTAAAAACTCTATTGAAAACCACTTTGATACCAGTTTTGAACTGGAAGCAAt                                                              >  1:431206/1‑71 (MQ=255)
        tACACACTCATCTAAAAACTCTATTGAAAACCACTTTGATACCAGTTTTGAACTGGAAGCAATGCTGGaa                                                       >  1:3451898/1‑70 (MQ=255)
         acacacTCATCTAAAAACTCTATTGAAAACCACTTTGATACCAGTTTTGAACTGGAAGCAATGCTGGaaa                                                      <  1:2354200/70‑1 (MQ=255)
         acacacTCATCTAAAAACTCTATTGAAAACCACTTTGATACCAGTTTTGAACTGGAAGCAATGCTGGaaa                                                      <  1:3392178/70‑1 (MQ=255)
         acacacTCATCTAAAAACTCTATTGAAAACCACTTTGATACCAGTTTTGAACTGGAACCAATGCTGGaaa                                                      <  1:3730025/70‑1 (MQ=255)
          cacacTCATCTAAAAACTCTATTGAAAACCACTTTGATACCAGTTTTGAACTGGAAGCAATGCTGGaaaa                                                     <  1:1546445/70‑1 (MQ=255)
          cacacTCATCTAAAAACTCTATTGAAAACCACTTTGATACCAGTTTTGAACTGGAAGCAATGCTGGaaaa                                                     <  1:2621307/70‑1 (MQ=255)
               tCATCTAAAAACTCTATTGAAAACCACTTTGATACCAGTTTTGAACTGGAAGCAATGCTGGAAAAACgtg                                                >  1:70131/1‑70 (MQ=255)
                cATCTAAAAACTCTATTGAAAACCACTTTGATACCAGTTTTGAACTGGAAGCAATGCTGGAAAAACgtgt                                               >  1:1898037/1‑70 (MQ=255)
                cATCTAAAAACTCTATTGAAAACCACTTTGATACCAGTTTTGAACTGGAAGCAATGCTGGAAAAACgtgt                                               >  1:1584228/1‑70 (MQ=255)
                cATCTAAAAACTCTATTGAAAACCACTTTGATACCAGTTTTGAACTGGAAGCAATGCTGGAAAAACgtgt                                               >  1:3584288/1‑70 (MQ=255)
                cATCTAAAAACTCTATTGAAAACCACTTTGATACCAGTTTTGAACTGGAAGCAATGCTGGAAAAACgtgt                                               >  1:297797/1‑70 (MQ=255)
                 aTCTAAAAACTCTATTGAAAACCACTTTGATACCAGTTTTGAACTGGAAGCAATGCTGGAAAAACGTGTaa                                             <  1:3664131/71‑1 (MQ=255)
                 aTCTAAAAACTCTATTGAAAACCACTTTGATACCAGTTTTGAACTGGAAGCAATGCTGGAAAAACGTGTaa                                             <  1:83389/71‑1 (MQ=255)
                  tCTAAAAACTCTATTGAAAACCACTTTGATACCAGTTTTGAACTGGAAGCAATGCTTGaaa                                                      <  1:2252033/61‑1 (MQ=255)
                  tCTAAAAACTCTATTGAAAACCACTTTGATACCAGTTTTGAACTGGAAGCAATGCTGGaaa                                                      <  1:247140/61‑1 (MQ=255)
                  tCTAAAAACTCTATTGAAAACCACTTTGATACCAGTTTTGAACTGGAAGCAATGCTGGaaa                                                      <  1:3099747/61‑1 (MQ=255)
                  tCTAAAAACTCTATTGAAAACCACTTTGATACCAGTTTTGAACTGGAAGCAATGCTGGaaa                                                      <  1:2349497/61‑1 (MQ=255)
                   cTAAAAACTCTATTGAAAACCACTTTGATACCAGTTTTGAACTGGAAGCAATGCTGGaaa                                                      <  1:720358/60‑1 (MQ=255)
                    tAAAAACTCTATTGAAAACCACTTTGATACCAGTTTTGAACTGGAAGCAATGCTGGAAAAACGTGTAAAAc                                          >  1:3996312/1‑71 (MQ=255)
                    tAAAAACTCTATTGAAAACCACTTTGATACCAGTTTTGAACTGGAAGCAATGCTGGAAAAACGTGTAAAAc                                          >  1:465419/1‑71 (MQ=255)
                    tAAAAACTCTATTGAAAACCACTTTGATACCAGTTTTGAACTGGAAGCAATGCTGGAAAAACGTGTAAAAc                                          >  1:3665563/1‑71 (MQ=255)
                      aaaaCTCTATTGAAAACCACTTTGATACCAGTTTTGAACTGGAAGCAATGCTGGAAAAACGTGTAAAACGt                                        >  1:1538186/1‑71 (MQ=255)
                      aaaaCTCTATTGAAAACCACTTTGATACCAGTTTTGAACTGGAAGCAATGCTGGAAAAACGTGTAAAACGt                                        >  1:2393483/1‑71 (MQ=255)
                         aCTCTATTGAAAACCACTTTGATACCAGTTTTGAACTGGAAGCAATGCTGGAAAAAc                                                   >  1:303252/1‑57 (MQ=255)
                           tctATTGAAAACCACTTTGATACCAGTTTTGAACTGGAAGCaa                                                               >  1:3173670/1‑43 (MQ=255)
                            ctATTGAAAACCACTTTGATACCAGTTTTGAACTGGAAGCGATGCTGGAAAAACGTGTAAAACGTCAACTg                                  >  1:4107494/1‑71 (MQ=255)
                            ctATTGAAAACCACTTTGATACCAGTTTTGAACTGGAAGCGATGCTGGAAAAACGTGTAAAACGTCAACTg                                  >  1:1106985/1‑71 (MQ=255)
                            ctATTGAAAACCACTTTGATACCAGTTTTGAACTGGAAGCGATGCTGGAAAAACGTGTAAAACGTCAACTg                                  >  1:3581307/1‑71 (MQ=255)
                            ctATTGAAAACCACTTTGATACCAGTTTTGAACTGGAAGCGATGCTGGAAAAACGTGTAAAACGTCAACTg                                  >  1:3450354/1‑71 (MQ=255)
                             tATTGAAAACCACTTTGATACCAGTTTTGAACTGGAAGCAATGCTGGaaa                                                      <  1:2420326/50‑1 (MQ=255)
                             tATTGAAAACCACTTTGATACCAGTTTTGAACTGGAAGCAATGCTGGaaa                                                      <  1:1326835/50‑1 (MQ=255)
                                tGAAAACCACTTTGATACCAGTTTTGAACTGGAAGCAATGCTGGAAAAAGGTGTAAAACGTCAACTGCtt                               >  1:3281755/1‑70 (MQ=255)
                                tGAAAACCACTTTGATACCAGTTTTGAACTGGAAGCAATGCTGGAAAAACGTGTAAAACGTCAACTGCtt                               >  1:4106329/1‑70 (MQ=255)
                                 gAAAACCACTTTGATACCAGTTTTGAACTGGAAGAAATGCTGGAAAAACGTGTAAAACGTCAACTGCTtga                             >  1:3422388/1‑71 (MQ=255)
                                   aaaCCACTTTGATACCAGTTTTGAACTGGAAGCGATGCTGGAAAAACGTGTAAAACGTCAACTGCTtgatg                           <  1:621677/71‑1 (MQ=255)
                                   aaaCCACTTTGATACCAGTTTTGAACTGGAAGCAATGCTGGAACAACGTGTAAAACGTCAACTGCTtgatg                           <  1:1101615/71‑1 (MQ=255)
                                   aaaCCACTTTGATACCAGTTTTGAACTGGAAGCAATGCTGGAAAAACGTGTAAAACGTCAACTGCTtgatg                           <  1:2781185/71‑1 (MQ=255)
                                   aaaCCACTTTGATACCAGTTTTGAACTGGAAGCAATGCTGGAAAAACGTGTAAAACGTCAACTGCTtgatg                           <  1:4100367/71‑1 (MQ=255)
                                   aaaCCACTTTGATACCAGTTTTGAACTGGAAGCAATGCTGGAAAAACGTGTAAAACGTCAACTGCTtgatg                           <  1:2403270/71‑1 (MQ=255)
                                   aaaCCACTTTGATACCAGTTTTGAACTGGAAGCAATGCTGGAAAAACGTGTAAAACGTCAACTGCTtgatg                           <  1:1589393/71‑1 (MQ=255)
                                   aaaCCACTTTGATACCAGTTTTGAACTGGAAGCAATGCTGGAAAAACGTGTAAAACGTCAACTGCTtgatg                           <  1:1106036/71‑1 (MQ=255)
                                    aaCCACTTTGATACCAGTTTTGAACTGGAAGCAATGCTGGAAAAACGTGTAAAACGTCAACTGCTtgatga                          <  1:971783/71‑1 (MQ=255)
                                    aaCCACTTTGATACCAGTTTTGAACTGGAAGCAATGCTGGAAAAACGTGTAAAACGTCAACTGCTtgatga                          <  1:1364330/71‑1 (MQ=255)
                                     aCCACTTTGATACCAGTTTTGAACTGGAAGCAATGCTGGAAAAACGTGTAAAACGTCAACTGCTtgatga                          <  1:843633/70‑1 (MQ=255)
                                        aCTTTGATACCAGTTTTGAACTGGAAGCAATGCTGGAAAAACGTGTAAAACGTCAACTGCTTGATGAAGTg                      <  1:731133/71‑1 (MQ=255)
                                        aCTTTGATACCAGTTTTGAACTGGAAGCAATGCTGGAAAAACGTGTAAAACGTCAACTGCTTGATGAAGTg                      <  1:1401556/71‑1 (MQ=255)
                                         cTTTGATACCAGTTTTGAACTGGAAGCAATGCTGGAAAAACGTGTAAAACGTCAACTGCTtgatga                          <  1:3513569/66‑1 (MQ=255)
                                         cTTTGATACCAGTTTTGAACTGGAAGCAATGCTGGAAAAACGTGTAAAACGTCAACTGCTTGATGAAGTg                      <  1:2967145/70‑1 (MQ=255)
                                               tACCAGTTTTGAACTGGAAGCGATGCTGGAAAAACGTGTAAATCGTCAACTGCTTGATGAAGTGCAGTCTa               <  1:3836433/71‑1 (MQ=255)
                                               tACCAGTTTTGAACTGGAAGCGATGCTGGAAAAACGTGTAAAACGTCAACTGCTTGATGAAGTGCAGTCTa               <  1:342370/71‑1 (MQ=255)
                                                 ccAGTTTTGAACTGGAAGCAATGCTGGAAAAACGTGTAAAACGTCAACTGCTTGATGAAGTGCAGTCTAt              <  1:3017526/70‑1 (MQ=255)
                                                 ccAGTTTTGAACTGGAAGCAATGCTGGAAAAACGTGTAAAACGTCAACTGCTTGATGAAGTGCAGTCTAt              <  1:791292/70‑1 (MQ=255)
                                                 ccAGTTTTGAACTGGAAGCAATGCTGGAAAAACGTGTAAAACGTCAACTGCTTGATGAAGTGCAGTCTAt              <  1:2653594/70‑1 (MQ=255)
                                                 ccAGTTTTGAACTGGAAGCAATGCTGGAAAAACGTGTAAAACGTCAACTGCTTGATGAAGTGCAGTCTAt              <  1:2840771/70‑1 (MQ=255)
                                                   aGTTTTGAACTGGAAGCAATGCTGGAAAAACGTGTAAAACGTCAACTGCTTGATGAAGTGCAGTCTATTTg           <  1:1447278/71‑1 (MQ=255)
                                                   aGTTTTGAACTGGAAGCAATGCTGGAAAAACGTGTAAAACGTCAACTGCTTGATGAAGTGCAGTCTATTTg           <  1:2686272/71‑1 (MQ=255)
                                                   aGTTTTGAACTGGAAGCAATGCTGGAAAAACGTGTAAAACGTCAACTGCTTGATGAAGTGCAGTCTATTTg           <  1:3314617/71‑1 (MQ=255)
                                                   aGTTTTGAACTGGAAGCAATGCTGGAAAAACGTGTAAAACGTCAACTGCTTGATGAAGTGCAGTCTATTTg           <  1:2991011/71‑1 (MQ=255)
                                                     ttttGAACTGGAAGCAATGCTGGAAAAACGTGTAAAACGTCAACTGCTTGATGAAGTGCAGTCTATTTGTc         <  1:2686510/71‑1 (MQ=255)
                                                       ttGAACTGGAAGCAATGCTGGAAAAACGTGTAAAACGTCAACTGCTTGATGAAGTGCAGTCTATTTGTCCa       >  1:2851339/1‑71 (MQ=255)
                                                       ttGAACTGGAAGCAATGCTGGAAAAACGTGTAAAACGTCAACTGCTTGATGAAGTGCAGTCTATTTGTCCa       >  1:2564792/1‑71 (MQ=255)
                                                       ttGAACTGGAAGCAATGCTGGAAAAACGTGTAAAACGTCAACTGCTTGATGAAGTGCAGTCTATTTGTCCa       >  1:336650/1‑71 (MQ=255)
                                                           aCTGGAAGCAATGCTGGAAAAACGTGTAAAACGTCAACTGCTTGATGAAGTGCAGTCTATTTGTCCACCGc   <  1:3167805/71‑1 (MQ=255)
                                                           aCTGGAAGCAATGCTGGAAAAACGTGTAAAACGTCAACTGCTTGATGAAGTGCAGTCTATTTGTCCACCGc   <  1:2148758/71‑1 (MQ=255)
                                                            cTGGAAGCAATGCTGGAAAAACGTGTAAAACGTCAACTGCTtgatg                           <  1:1242984/46‑1 (MQ=255)
                                                            cTGGAAGCAATGCTGGAAAAACGTGTAAAACGTCAACTGCTtgatg                           <  1:3150843/46‑1 (MQ=255)
                                                            cTGGAAGCAATGCTGGAAAAACGTGTAAAACGTCAACTGCTtgatg                           <  1:1189503/46‑1 (MQ=255)
                                                            cTGGAAGCAATGCTGGAAAAACGTGTAAAACGTCAACTGCTtgatg                           <  1:159121/46‑1 (MQ=255)
                                                            cTGGAAGCAATGCTGGAAAAACGTGTAAAACGTCAACTGCTtgatg                           <  1:4113189/46‑1 (MQ=255)
                                                            cTGGAAGCAATGCTGGAAAAACGTGTAAAACGTCAACTGCTTGATGAAGTGCAGTCTATTTGTCCACCg    <  1:2920010/69‑1 (MQ=255)
                                                            cTGGAAGCAATGCTGGAAAAACGTGTAAAACGTCAACTGCTTGATGAAGTGCAGTCTATTTGTCCACCGCa  <  1:1868228/71‑1 (MQ=255)
                                                            cTGGAAGCAATGCTGGAAAAACGTGTAAAACGTCAACTGCTTGATGAAGTGCAGTCTATTTGTCCACCGCa  <  1:1521556/71‑1 (MQ=255)
                                                             tGGAAGCAATGCTGGAAAAACGTGTAAAACGTCAACTGCTTGATGAAGTGCAGTCTATTTGTCCACCGCa  <  1:3431219/70‑1 (MQ=255)
                                                                  gCAATGCTGGAAAAACGTGTAAAACGTCAACTGCTtgatg                           <  1:4084853/40‑1 (MQ=255)
                                                                  gCAATGCTGGAAAAACGTGTAAAACGTCAACTGCTtgatg                           <  1:753995/40‑1 (MQ=255)
                                                                  gCAATGCTGGAAAAACGTGTAAAACGTCAACTGCTtgatg                           <  1:3759364/40‑1 (MQ=255)
                                                                   cAATGCTGGAAAAACGTGTAAAACGTCAACTGCtt                               <  1:1932205/35‑1 (MQ=255)
                                                                   cAATGCTGGAAAAACGTGTAAAACGTCAACTGCtt                               <  1:3815667/35‑1 (MQ=255)
                                                                   cAATGCTGGAAAAACGTGTAAAACGTCAACTGCtt                               <  1:2378435/35‑1 (MQ=255)
                                                                    aaTGCTGGAAAAACGTGTAAAACGTCAACTGCTTGATGAAg                        >  1:1871109/1‑41 (MQ=255)
                                                                    |                                                              
GTGCTGGTTACACACTCATCTAAAAACTCTATTGAAAACCACTTTGATACCAGTTTTGAACTGGAAGCAATGCTGGAAAAACGTGTAAAACGTCAACTGCTTGATGAAGTGCAGTCTATTTGTCCACCGCA  >  minE/843307‑843437

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: