Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 3,873,040 C→A R332L (CGC→CTC)  glmS ← D‑fructose‑6‑phosphate amidotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6063,873,0400CA92.9% 140.1 / 3.6 42R332L (CGC→CTC) glmSD‑fructose‑6‑phosphate amidotransferase
Reads supporting (aligned to +/- strand):  ref base C (0/3);  new base A (18/21);  total (18/24)
Fisher's exact test for biased strand distribution p-value = 2.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

GGTATCCGCGGTTTCGCCAGACTGTGACAAGGTGATCATCAGGCTGTTACGACGCACGGCAGATTTGCGATAGCGGAATTCAGAGGCGATTTCGACGTCGCACGGAATACCTGCTAGCGATTCAAACCAGTAGCGGGAAACCAT  >  REL606/3872967‑3873110
                                                                         |                                                                      
ggTATCCGCGGTTTCGCCAGACTGTGACAAGGTGATCATCAGGCTGTTACGACGCACGGCAGATTTGCGATAGAgg                                                                      >  2:501030/1‑76 (MQ=255)
ggTATCCGCGGTTTCGCCAGACTGTGACAAGGTGATCATCAGGCTGTTACGACGCACGGCAGATTTGCGATAGAgg                                                                      <  2:783566/76‑1 (MQ=255)
      cgcgGTTTCGCCAGACTGTGACAAGGTGATCATCAGGCTGTTACGACGCACGGCAGATTTGCGATAGAGGAATTCa                                                                <  1:1106994/76‑1 (MQ=255)
       gcgGTTTCGCCAGACTGTGACAAGGTGATCATCAGGCTGTTACGACGCACGGCAGATTTGCGATAGAGGAATTCag                                                               <  1:583672/76‑1 (MQ=255)
            ttCGCCAGACTGTGACAAGGTGATCATCAGGCTGTTACGACGCACGGCAGATTTGCGATAGAGGAATTCAGAGGCg                                                          <  1:226841/76‑1 (MQ=255)
                      tgtgACAAGGTGATCATCAGGCTGTTACGACGCACGGCAGATTTGCGATAGAGGAATTCAGAGGCGATTTCGACGt                                                >  1:1073547/1‑76 (MQ=255)
                           cAAGGTGATTATCAGGCTGTTACGACGCACGGCAGATTTGCGATAGAGGAATTCAGAGGCGATTTCGACGTCGCAc                                           >  2:253431/1‑76 (MQ=255)
                           cAAGGTGATCATCAGGCTGTTACGACGCACGGCAGATTTGCGATAGAGGAATTCAGAGGCGATTTCGACGTCGCAc                                           <  1:652926/76‑1 (MQ=255)
                           cAAGGTGATCATCAGGCTGTTACGACGCACGGCAGATTTGCGATAGAGGAATTCAGAGGCGATTTCGACGTCGCAc                                           >  1:265182/1‑76 (MQ=255)
                                     atcaGGCTGTTACGACGCACGGCAGATTTGCGATAGAGGAATTCAGAGGCGATTTCGACGTCGCACGGAATACCTg                                 <  2:13661/76‑1 (MQ=255)
                                      tcaGGCTGTTACGACGCACGGCAGATTTGCGATAGAGGAATTCAGAGGCGATTTCGACGTCGCACGGAATACCTGc                                <  2:595559/76‑1 (MQ=255)
                                       caGGCTGTTACGACGCACGGCAGATTTGCGATAGCGGAATTCAGAGGCGATTTCGACGTCGCACGGAATACCTGCt                               <  1:240464/76‑1 (MQ=255)
                                       caGGCTGTTACGACGCACGGCAGATTTGCGATAGAGGAATTCAGAGGCGATTTCGACGTCGCACGGAATACCTGCt                               >  1:1245650/1‑76 (MQ=255)
                                         ggCTGTTACGACGCACGGCAGATTTGCGATAGAGGAATTCAGAGGCGATTTCGACGTCGCACGGAATACCTGCTAg                             <  2:472560/76‑1 (MQ=255)
                                          gCTGTTACGACGCACGGCAGATTTGCGATAGAGGAATTCAGAGGCGATTTCGACGTCGCACGGAATACCTGCTAGc                            >  1:665337/1‑76 (MQ=255)
                                            tGTTACGACGCACNGCAGATTTGCGATAGAGGAATTCAGAGGCGATTTCGACGTCGCACGGAATACCTGCTAGCGa                          >  1:853225/1‑76 (MQ=255)
                                               tACGACGCACGGCAGATTTGCGATAGAGGAATTCAGAGGCGATTTCGACGTCGCACGGAATACCTGCTAGCGATTc                       >  2:788376/1‑76 (MQ=255)
                                                acgacgCACGGCAGATTTGCGATAGAGGAATTCAGAGGCGATTTCGACGTCGCACGGAATACCTGCTAGCGATTCa                      <  1:323297/76‑1 (MQ=255)
                                                acgacgCACGGCAGATTTGCGATAGAGGAATTCAGAGGCGATTTCGACGTCGCACGGAATACCTGCTAGCGATTCa                      <  1:1149873/76‑1 (MQ=255)
                                                 cgacgCACGGCAGGTTTGCGATAGAGGAATTCAGAGGCGATTTCGACGTCGCACGCAGTACTGGCTAGCGATTCaa                     >  2:1060945/1‑76 (MQ=255)
                                                  gacgCACGGCAGATTTGCGATAGAGGAATTCAGAGGCGATTTCGACGTCGCACGGAATACCTGCTAGCGATTCaaa                    >  1:1155967/1‑76 (MQ=255)
                                                    cgCACGGCAGATTTGCGATAGAGGAATTCAGAGGCGATTTCGACGTCGCACGGAATACCTGCTAGCGATTCAAAcc                  >  2:57363/1‑76 (MQ=255)
                                                    cgCACGGCAGATTTGCGATAGAGGAATTCAGAGGCGATTTCGACGTCGCACGGAATACCTGCTAGCGATTCAAAcc                  >  1:1258503/1‑76 (MQ=255)
                                                    cgCACGGCAGATTTGCGATAGAGGAATTCAGAGGCGATTTCAACGTCGCACGGAATACCTGCTAGCGATTCAAAcc                  <  1:99790/76‑1 (MQ=255)
                                                     gCACGGCAGATTTGCGATAGAGGAATTCAGAGGCGATTTCGACGTCGCACGGAATACCTGCTAGCGATTCAAACCa                 <  2:842698/76‑1 (MQ=255)
                                                      cACGGCAGATTTGCGATAGAGGAATTCAGAGGCGATTTCGACGTCGCACGGAATACCTGCTAGCGATTCAAACCAg                <  1:629517/76‑1 (MQ=255)
                                                      cACGGCAGATTTGCGATAGAGGAATTCAGAGGCGATTTCGACGTCGCACGGAATACCTGCTAGCGATTCAAACCAg                >  2:481481/1‑76 (MQ=255)
                                                       aCGGCAGCTTTGCGATAGAGGAATTCAGAGGCGATTTCGACGTCGCACGGAATACCTGCTAGCGATTCAAACCAGt               >  1:810413/1‑76 (MQ=255)
                                                        cGGCAGATTTGCGATAGAGGAATTCAGAGGCGATTTCGACGTCGCACGGAATACCTGCTAGCGATTCAAACCAGTa              <  2:446387/76‑1 (MQ=255)
                                                         ggCAGATTTGCGATAGCGGAATTCAGAGGCGATTTCGACGTCGCACGGAATACCTGCTAGCGATTCAAACCAGTAg             <  2:364836/76‑1 (MQ=255)
                                                         ggCAGATTTGCGATAGAGGAATTCAGAGGCGATTTCGACGTCGCACGGAATACCTGCTAGCGNTTCAAACCAGTAg             <  1:857314/76‑1 (MQ=255)
                                                         ggCAGATTTGCGATAGAGGAATTCAGAGGCGATTTCGACGTCGCACGGAATACCTGCTAGCGATTCAAACCAGTAg             <  2:788882/76‑1 (MQ=255)
                                                         ggCAGATTTGCGATAGAGGAATTCAGAGGCGATTTCGACGTCGCACGGAATACCTGCTAGCGATTCAAACCAGTAg             >  1:1053908/1‑76 (MQ=255)
                                                           cAGATTTGCGATAGAGGAATTCAGAGGCGATTTCGACGTCGCACGGAATACCTGCTAGCGATTCAAACCAGTAGCg           <  2:209751/76‑1 (MQ=255)
                                                             gATTTGCGATAGAGGAATTCAGAGGCGATTTCGACGTCGCACGGAATACCTGCTAGCGATTCAAACCAGTAGCggg         <  2:865714/76‑1 (MQ=255)
                                                              aTTTGCGATAGAGGAATTCAGAGGCGATTTCGACGTCGCACGGAATACCTGCTAGCGATTCAAACCAGTAGCGGGa        >  1:737579/1‑76 (MQ=255)
                                                              aTTTGCGATAGAGGAATTCAGAGGCGATTTCGACGTCGCACGGAATACCTGCTAGCGATTCAAACCAGTAGCGGGa        >  1:565847/1‑76 (MQ=255)
                                                               tttGCGATAGAGGAATTCAGAGGCGATTTCGACGTCGCACGGAATACCTGCTAGCGATTCAAACCAGTAGCGGGan       <  1:337507/76‑2 (MQ=255)
                                                                ttGCGATAGAGGAATTCAGAGGCGATTTCGACGTCGCACGGAATACCTGCTAGCGATTCAAACCAGTAGCGGGaaa      <  1:107204/76‑1 (MQ=255)
                                                                 tGCGATAGCGGAATTCAGAGGCGATTTCGACGTCGCACGGAATACCTGCTAGCGATTCAAACCAGTAGCGGGAAAc     <  2:565025/76‑1 (MQ=255)
                                                                  gCGATAGAGGAATTCAGAGGCGATTTCGACGTCGCACGGAATACCTGCTAGCGATTCAAACCAGTAGCGGGAAAcc    >  2:145696/1‑76 (MQ=255)
                                                                    gATAGAGGAATTCAGAGGCGATTTCGACGTCGCACGGAATACCTGCTAGCGATTCAAACCAGTAGCGGGAAACCAt  <  2:1258503/76‑1 (MQ=255)
                                                                         |                                                                      
GGTATCCGCGGTTTCGCCAGACTGTGACAAGGTGATCATCAGGCTGTTACGACGCACGGCAGATTTGCGATAGCGGAATTCAGAGGCGATTTCGACGTCGCACGGAATACCTGCTAGCGATTCAAACCAGTAGCGGGAAACCAT  >  REL606/3872967‑3873110

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