Predicted mutation
evidence seq id position mutation freq annotation gene description
RA BS168_reference_genome 2,358,944 A→G 7.4% R127R (CGT→CGC)  mgsA ← methylglyoxal synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*BS168_reference_genome2,358,9440AG 7.4% 160.5 / 4.3 54R127R (CGT→CGC) mgsAmethylglyoxal synthase
Reads supporting (aligned to +/- strand):  ref base A (24/26);  new base G (2/2);  total (26/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.99e-01

TCGCATATCATTACTTTTTTTTCCTGTTTGACAAACTTCGCTATTGTGCCGCCCATTCCGATCTCGACATCATCACTGTGGGCGCCAAAAGCAAGAACGTCAGCATTATACATTCGGCTCTTCTCCCCGAAGAAGGTCACGGAATTCGAAAACACCTTCATCAAGTGTGCGCACAAGAATTTCCGCAGTACCCATATTTGTGGCGAGCGGAATGGAATACACATCACATAAACGGATTAATGCCGAGACATCCGGTTCATGCGGCTGCGCGGTC  >  BS168_reference_genome/2358806‑2359079
                                                                                                                                          |                                                                                                                                       
tCGCATATCATTACTTTTTTTTCCTGTTTGACAAACTTCGCTATTGTGCCGCCCATTCCGATCTCGACATCATCACTGTGGGCGCCAAAAGCAAGAACGTCAGCATTATACATTCGGCTCTTCTCCCCGAAGAAGGTCa                                                                                                                                         <  3:215853/139‑1 (MQ=255)
tCGCATATCATTACTTTTTTTTCCTGTTTGACAAACTTCGCTATTGTGCCGCCCATTCCGATCTCGACATCATCACTGTGGGCGCCAAAAGCAAGAACGTCAGCATTATACATTCGGCTCTTCTCCCCGAAGAAGGTCa                                                                                                                                         <  2:220698/139‑1 (MQ=255)
  gCATATCATTACTTTTTTTTCCTGTTTGACAAACTTCGCTATTGTGCCGCCCATTCCGATCTCGACATCATCACTGTGGGCGCCAAAAGCAAGAACGTCAGCATTATACATTCGGCTCTTCTCCCCGAAGAAGGTCACg                                                                                                                                       <  4:424045/139‑1 (MQ=255)
   cATATCATTACTTTTTTTTCCTGTTTGACAAACTTCGCTATTGTGCCGCCCCTTCCGATCTCGACATCATCACTGGGGGCGCCAAAAGCAAGAACGTCAGCATTATACATTCGGCTCTTCTCCCCGAAGAAGGTCACgg                                                                                                                                      <  6:170114/139‑1 (MQ=255)
             cTTTTTTTTCCTGTTTGACAAACTTCGCTATTGTGCCGCCCATTCCGATCTCGACATCATCACTGTGGGCGCCAAAAGCAAGAACGTCAGCATTATACATTCGGCTCTTCTCCCCGAAGAAGGTCACGGAATTCGaaaa                                                                                                                            >  2:122882/1‑139 (MQ=255)
                ttttttCCTGTTTGACAAACTTCGCTATTGTGCCGCCCATTCCGATCTCGACATCATCACTGTGGGCGCCAAAAGCAAGAACGTCAGCATTATACATTCGGCTCTTCTCCCCGAAGAAGGTCACGGAATTCGAAAaca                                                                                                                          <  3:309092/138‑1 (MQ=255)
                  ttttCCTGTTTGACAAACTTCGCTATTGTGCCGCCCATTCCGATCTCGACATCATCACTGTGGGCGCCAAAAGCAAGAACGTCAGCATTATACATTCGGCTCTTCTCCCCGAAGAAGGTCACGGAATTCGAAAACAcc                                                                                                                        <  6:384065/138‑1 (MQ=255)
                  ttttCCTGTTTGACAAACTTCGCTATTGTGCCGCCCATTCCGATCTCGACATCATCACTGTGGGCGCCAAAAGCAAGAACGTCAGCATTATACATTCGGCTCTTCTCCCCGAAGAAGGTCACGGAATTCGAAAACAcc                                                                                                                        <  1:122882/138‑1 (MQ=255)
                  ttttCCTGTTTGACAAACTTCGCTATTGTGCCGCCCATTCCGATCTCGACATCATCACTGTGGGCGCCAAAAGCAAGAACGTCAGCATTATACATTCGGCTCTTCTCCCCGAAGAAGGTCACGGAATTCGAAAACAcc                                                                                                                        <  7:202400/138‑1 (MQ=255)
                  ttttCCTGTTTGACAAACTTCGCTATTGTGCCGCCCATTCCGATCTCGACATCATCACTGTGGGCGCCAAAAGCAAGAACGTCAGCATTATACATTCGGCTCTTCTCCCCGAAGAAGGTCACGGAATTCGAAAACAcc                                                                                                                        <  8:36716/138‑1 (MQ=255)
                  ttttCCTGTTTGACAAACTTCGCTATTGTGCCGCCCATTCCGATCTCGACATCATCACTGTGGGCGCCAAAAGCAAGAACGTCAGCATTATACATTCGGCTCTTCTCCCCGAAGAAGGTCACGGAATTCGAAAACAcc                                                                                                                        <  7:241030/138‑1 (MQ=255)
                            tGACAAACTTCGCTATTGTGCCGCCCATTCCGATCTCGACATCATCACTGTGGGCGCCAAAAGCAAGAACGTCAGCATTATACATTCGGCTCTTCTCCCCGAAGAAGGTCACGGAATTCGAAAACACCTTCATCAAgtg                                                                                                             >  1:346370/1‑139 (MQ=255)
                                      cgcgATTGTGCCGCCCATTCCGATCTCGACATCATCACTATGGGCGCCAAAAGCAAGAACGTCAGTATTATACATTCGGCTCTTCTCCCCGAAGAAGGTCGCGGAATTCGAAAACACCTTCATCAAGTGTGCGgaca                                                                                                     <  7:183120/137‑1 (MQ=255)
                                              tgCCGCCCATTCCGATCTCGACATCATCACTGTGGGCGCCAAAAGCAAGAACGTCAGCATTATACATTCGGCTCTTCTCCCCGAAGAAGGTCACGGAATTCGAAAACACCTTCATCAAGTGTGCGCACAAGAATTTCCg                                                                                           <  4:191046/139‑1 (MQ=255)
                                              tgCCGCCCATTCCGATCTCGACATCATCACTGTGGGCGCCAAAAGCAAGAACGTCAGCATTATACATTCGGCTCTTCTCCCCGAAGAAGGTCACGGAATTCGAAAACACCTTCATCAAGTGTGCGCACAAGAATTTCCg                                                                                           <  2:346370/139‑1 (MQ=255)
                                                ccgccCATTCCTATCTCGACATCATCACTGTGGGCGACAAAAGCAAGAACGTCAGCATTATACATTCGGCTCTTCTCCCCGAAGAAGGGCACTTAATTCTAAAACACCTTAATCAAGTGTGAGCACAATAATTTCCGCa                                                                                         <  5:188058/139‑1 (MQ=255)
                                                ccgccCATTCCGATCTCGACATCATCACTGTGGGCGCCAAAAGCAAGAACGTCAGCATTATACATTCGGCTCTTCTCCCCGAAGAAGGTCACGGAATTCGAAAACACCTTCATCAAGTGTGCGCACAAGAATTTCCGCa                                                                                         <  2:247030/139‑1 (MQ=255)
                                                ccgccCATTCCGATCTCGACATCATCAATGTGGGCGACAAAATCAAGAACGTAATCATTATACATTCGGCTATGCTCCCCGAAGAATGTCACGGAATTAGAAAACACCTTACTCAAGTGTGAGCCCAAGAATTTcctcc                                                                                         >  6:188058/1‑136 (MQ=255)
                                                   cccATTCCGATCTCGACATCATCACTGTGGGCGCCAAAAGCAAGAACGTCAGCATTATACATTCGGCTCTTCTCCCCGAAGAAGGTCACGGAATTCGAAAaca                                                                                                                          >  5:224182/1‑103 (MQ=255)
                                                   cccATTCCGATCTCGACATCATCACTGTGGGCGCCAAAAGCAAGAACGTCAGCATTATACATTCGGCTCTTCTCCCCGAAGAAGGTCACGGAATTCGAAAaca                                                                                                                          <  6:224182/103‑1 (MQ=255)
                                                                cGACATCATCACTATGGGCGCCAAAAGCAAGAACGTCAGTATTATACATTCGGCTCTTCTCCCCGAAGAAGGTCGCGGAATTCGAAAACACCTTCATCAAGTGTGCGGACAAGAATCTCCGCAGTACCCATATTtgtag                                                                         >  3:366494/1‑137 (MQ=255)
                                                                      catcaCTGTGGGCGCCAAAAGCAAGAACGTCAGCATTATACATTCGGCTCTTCTCCCCGAAGAAGGTCACGGAATTCGAAAACACCTTCATCAAGTGTGCGCACAAGAATTTCCGCa                                                                                         <  2:861/117‑1 (MQ=255)
                                                                      catcaCTGTGGGCGCCAAAAGCAAGAACGTCAGCATTATACATTCGGCTCTTCTCCCCGAAGAAGGTCACGGAATTCGAAAACACCTTCATCAAGTGTGCGCACAAGAATTTCCGCa                                                                                         >  1:861/1‑117 (MQ=255)
                                                                      catcaCTGTGGGCGCCAAAAGCAAGAACGTCAGCATTATACATTCGGCTCTTCTCCCCGAAGAAGGTCACGGAATTCGAAAACACCTTCATCAAGTGTGCGCACAAGAATTTCCGCAGTACCCATATTTGTGGCGAGCg                                                                   >  7:192415/1‑139 (MQ=255)
                                                                       atcaCTGTGGGCGCCAAAAGCAAGAACGTCAGCATTATACATTCGGCTCTTCTCCCCGAAGAAGGTCACGGAATTCGAAAACACCTTCATCAAGTGTGCGCACAAGAATTTCCGCAGTACCCATATTTGTGGCGAGCgg                                                                  <  3:183801/139‑1 (MQ=255)
                                                                       atcaCTATGGGCGCCAAAAGCAAGAACGTCAGTATTATACATTCGGCTCTTCTCCCCGAAGAAGGTCGCGGAATTCGAAAACACCTTCATCAAGTGTGCGGACAAGAATCTCCGCAGTACCCATATTTGTagcgagcg                                                                   >  1:262381/1‑138 (MQ=255)
                                                                        tcaCTGTGGGCGCCAAAAGCAAGAACGTCAGCATTATACATTCGGCTCTTCTCCCCGAAGAAGGTCACGGAATTCGAAAACACCTTCATCAAGTGTGCGCACAAGAATTTCCGCAGTACCCATATTTGTGGCGAGCgga                                                                 <  7:205643/139‑1 (MQ=255)
                                                                         caCTGTGGGCGCCAAAAGCAAGAACGTCAGCATTATACATTCGGCTCTTCTCCCCGAAGAAGGTCACGGAATTCGAAAACACCTTCATCAAGTGTGCGCACAAGAATTTCCGCAGTACCCATATTTGTGGCGAGCggaa                                                                >  1:127960/1‑139 (MQ=255)
                                                                              tgGGCGCCAAAAGCAAGAACGTCAGCATTATACATTCGGCTCTTCTCCCCGAAGAAGGTCACGGAATTCGAAAACACCTTCATCAAGTGTGCGCACAAGAATTTCCGCAGTACCCATATTTGTGGCGAGCggaatggaa                                                           >  6:424544/1‑139 (MQ=255)
                                                                              tgGGCGCCAAAAGCAAGAACGTCAGCATTATACATTCGGCTCTTCTCCCCGAAGAAGGTCACGGAATTCGAAAACACCTTCATCAAGTGTGCGCACAAGAATTTCCGCAGTACCCATATTTGTGGCGAGCg                                                                   <  6:129853/131‑1 (MQ=255)
                                                                              tgGGCGCCAAAAGCAAGAACGTCAGCATTATACATTCGGCTCTTCTCCCCGAAGAAGGTCACGGAATTCGAAAACACCTTCATCAAGTGTGCGCACAAGAATTTCCGCAGTACCCATATTTGTGGCGAGCg                                                                   >  5:129853/1‑131 (MQ=255)
                                                                                  cgcCAAAAGCAAGAACGTCAGTATTATACATTCGGCTCTTCTCCCCGAAGAAGGTCGCGGAATTCGAAAACACCTTCATCAAGTGTGCGGACAAGAATCTCCGCAGTACCCATATTTGTAGCGAGCAGAATGGAATaca                                                       <  4:366494/139‑1 (MQ=255)
                                                                                  cgcCAAAAGCAAGAACGTCAGCATTATACATTCGGCTCTTCTCCCCGAAGAAGGTCACGGAATTCGAAAACACCTTCATCAAGTGTGCGCACAAGAATTTCCGCAGTACCCATATTTGTGGCGAGCGGAATGGAATaca                                                       <  1:51056/139‑1 (MQ=255)
                                                                                             aGAACGTCAGCATTATACATTCGGCTCTTCTCCCCGAAGAAGGTCACGGAATTCGAAAACACCTTCATCAAGTGTGCGCACAAGAATTTCCGCAGTACCCATATTTGTGGCGAGCGGAATGGAATACACATCACATaaa                                            <  3:210037/139‑1 (MQ=255)
                                                                                                    cAGCATTATACATTCGGCTCTTCTCCCCGAAGAAGGTCACGGAATTCGAAAACACCTTCATCAAGTGTGCGCACAAGAATTTCCGCAGTAcc                                                                                    >  7:206598/1‑92 (MQ=255)
                                                                                                    cAGCATTATACATTCGGCTCTTCTCCCCGAAGAAGGTCACGGAATTCGAAAACACCTTCATCAAGTGTGCGCACAAGAATTTCCGCAGTAcc                                                                                    <  8:206598/92‑1 (MQ=255)
                                                                                                    cAGCATTATACATTCGGCTCTTCTCCCCGAAGAAGGTCACGGAATTCGAAAACACCTTCATCAAGTGTGCGCACAAGAATTTCCGCAGTACCCATATTTGTGGCGAGCGGAATGGAATACACATCACATaaa                                            >  2:91116/1‑132 (MQ=255)
                                                                                                    cAGCATTATACATTCGGCTCTTCTCCCCGAAGAAGGTCACGGAATTCGAAAACACCTTCATCAAGTGTGCGCACAAGAATTTCCGCAGTACCCATATTTGTGGCGAGCGGAATGGAATACACATCACATaaa                                            <  1:91116/132‑1 (MQ=255)
                                                                                                            taCATTCGGCTCTTCTCCCCGAAGAAGGTCACGGAATTCGAAAACACCTTCATCAAGTGTGCGCACAAGAATTTCCGCAGTACCCATATTTGTGGCGAGCggaa                                                                <  3:387335/104‑1 (MQ=255)
                                                                                                            taCATTCGGCTCTTCTCCCCGAAGAAGGTCACGGAATTCGAAAACACCTTCATCAAGTGTGCGCACAAGAATTTCCGCAGTACCCATATTTGTGGCGAGCggaa                                                                >  4:387335/1‑104 (MQ=255)
                                                                                                            taCATTCGGCTCTTCTCCCCGAAGAAGGTCACGGAATTCGAAAACACCTTCATCAAGTGTGCGCACAAGAATTTCCGCAGTACCCATATTTGTGGCGAGCGGAATGGAATACACATCACATAAACGGATTAATGcc                                >  5:136025/1‑136 (MQ=255)
                                                                                                            taCATTCGGCTCTTCTCCCCGAAGAAGGTCACGGAATTCGAAAACACCTTCATCAAGTGTGCGCACAAGAATTTCCGCAGTACCCATATTTGTGGCGAGCGGAATGGAATACACATCACATAAACGGATTAATGcc                                >  6:392216/1‑136 (MQ=255)
                                                                                                               aTTCGGCTCTTCTCCCCGAAGAAGGTCACGGAATTCGAAAACACCTTCATCAAGTGTGCGCACAAGAATTTCCGCAGTACCCATATTTGTGGCGAGCGGAATGGAATACACATCACATAAACGGATTAATGCCGAGACa                          >  5:199240/1‑139 (MQ=255)
                                                                                                                 tCGGCTCTTCTCCCCGAAGAAGGTCACGGAATTCGAAAACACCTTCATCAAGTGTGCGCACAAGAATTTCCGCAGTACCCATATTTGTGGCGAGCggaa                                                                >  1:227416/1‑99 (MQ=255)
                                                                                                                 tCGGCTCTTCTCCCCGAAGAAGGTCACGGAATTCGAAAACACCTTCATCAAGTGTGCGCACAAGAATTTCCGCAGTACCCATATTTGTGGCGAGCggaa                                                                <  2:227416/99‑1 (MQ=255)
                                                                                                                 tCGGCTCTTCTCCCCGAAGAAGGTCACGGAATTCGAAAACACCTTCATCAAGTGTGCGCACAAGAATTTCCGCAGTACCCATATTTGTGGCGAGCGGAATGGAATACACATCACATAAACGGATTAATGCCGAGACATc                        >  8:246852/1‑139 (MQ=255)
                                                                                                                 tCGGCTCTTCTCCCCGAAGAAGGTCACGGAATTCGAAAACACCTTCATCAAGTGTGCGCACAAGAATTTCCGCAGTACCCATATTTGTGGCGAGCGGAATGGAATACACATCACATAAACGGATTAATGCCGAGACATc                        >  3:384594/1‑139 (MQ=255)
                                                                                                                 tCGGCTCTTCTCCCCGAAGAAGGTCACGGAATTCGAAAACACCTTCATCAAGTGTGCGCACAAGAATTTCCGCAGTACCCATATTTGTGGCGAGCGGAATGGAATACACATCACATAAACGGATTAATGCCGAGACATc                        >  4:81931/1‑139 (MQ=255)
                                                                                                                 tCGGCTCTTCTCACCGAAGAAGGTCACGGAATTCGAAAACACCTTCATCAAGTGTGCTCACAAGAATTTCCGCATTACCCATATTTGTTGCGAGCGGAATGGATTACACATCACATAAACTGATTAATGCCTAGACagc                        >  5:394246/1‑137 (MQ=255)
                                                                                                                   ggCTCTTCTCCCCGAAGAAGGTCACGGAATTCGAAAACACCTTCATCAAGTGTGCGCACAAGAATTTCCGCAGTACCCATATTTGTGGCGAGCGGAATGGAATACACATCACATAAACGGATTAATGCCGAGACATCCg                      <  6:199240/139‑1 (MQ=255)
                                                                                                                    gCTCTTCTCCCCGAAGAAGGTCACGGAATTCGAAAACACCTTCATCAAGTGTGCGCACAAGAATTTCCGCAGTACCCATATTTGTGGCGAGCGGAATGGAATACACATCACATAAACGGATTAATGCCGAGACATCCgg                     >  8:324862/1‑139 (MQ=255)
                                                                                                                                  agaagGTCACGGAATTCGAAAACACCTTCATCAAGTGTGCGCACAAGAATTTCCGCAGTACCCATATTTGTGGCGAGCGGAATGGAATACACATCACATAAACGGATTAATGCCGAGACATCCGGTTCATGCGGCTgcg       <  5:430077/139‑1 (MQ=255)
                                                                                                                                       gTCACGGAATTCGAAAACACCTTCATCAAGTGTGCGCACAAGAATTTCCGCAGTACCCATATTTGTGGCGAGCGGAATGGAATACACATCACATAACCGGATTAATGCCGAGACATCCGGTGCCTGCGGCTGCGCGGTc  >  2:303401/1‑139 (MQ=255)
                                                                                                                                       gTCACGGAATTCGAAAACACCTTCATCAAGTGTGCGCACAAGAATTTCCGCAGTACCCATATTTGTGGCGAGCGGAATGGAATACACATCACATAAACGGATTAATGCCGAGACATCCGGTTCATGCGGCTGCGCgg    <  2:407054/137‑1 (MQ=255)
                                                                                                                                       gTCACGGAATTCGAAAACACCTTCATCAAGTGTGCGCACAAGAATTTCCGCAGTACCCATATTTGTGGCGAGCGGAATGGAATACACATCACATAAACGGATTAATGCCGAGACATCCGGTTCATGCGGCTGCGCgg    >  1:407054/1‑137 (MQ=255)
                                                                                                                                       gTCACGGAATTCGAAAACACCTTCATCAAGTGTGCGCACAAGAATTTCCGCAGTACCCATATTTGTGGCGAGCGGAATGGAATACACATCACATAAACGGATTAATGCCGAGACATCCGGTTCATGCGGCTGCGCGGTc  >  1:180469/1‑139 (MQ=255)
                                                                                                                                          |                                                                                                                                       
TCGCATATCATTACTTTTTTTTCCTGTTTGACAAACTTCGCTATTGTGCCGCCCATTCCGATCTCGACATCATCACTGTGGGCGCCAAAAGCAAGAACGTCAGCATTATACATTCGGCTCTTCTCCCCGAAGAAGGTCACGGAATTCGAAAACACCTTCATCAAGTGTGCGCACAAGAATTTCCGCAGTACCCATATTTGTGGCGAGCGGAATGGAATACACATCACATAAACGGATTAATGCCGAGACATCCGGTTCATGCGGCTGCGCGGTC  >  BS168_reference_genome/2358806‑2359079

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