Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 3,873,434 T→G T201P (ACC→CCC)  glmS ← D‑fructose‑6‑phosphate amidotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6063,873,4340TG97.2% 139.4 / ‑6.7 36T201P (ACC→CCC) glmSD‑fructose‑6‑phosphate amidotransferase
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base G (18/17);  total (18/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTATCGAAGATGTTTACCGAACGGCGAGTGATTTCCGCAATATCGCCCTCTTCAAGGAAGATAAAGCGACGGGTCACCGGCAACAGCGCCAGCTGGTCAGAAGCGATAAAGTTTTCGCCCATCCCCAGGCCAATCACCAG  >  REL606/3873361‑3873500
                                                                         |                                                                  
ttATCGAAGATGTTTACCGAACGGCGAGTGATTTCCGCAATATCGCCCTCTTCAAGGAAGATAAAGCGACGGGGCa                                                                  >  1:839322/1‑76 (MQ=255)
  aTCGAAGATGTTTACCGAACGGCGAGTGATTTCCGCAATATCGCCCTCTTCAAGGAAGATAAAGCGACGGGGCAcc                                                                >  2:1095307/1‑76 (MQ=255)
  aTCGAAGATGTTTACCGAACGGCGAGTGATTTCCGCAATATCGCCCTCTTCAAGGAAGATAAAGCGACGGGGCAcc                                                                >  2:936272/1‑76 (MQ=255)
   tCGAAGATGTTTACCGAACGGCGAGTGATTTCCGCAATATCGCCCTCTTCAAGGAAGATAAAGCGACGGGGCACCg                                                               <  1:453967/76‑1 (MQ=255)
   tCGAAGATGTTTACCGAACGACNNGTGATTTCCGNNTTATCGCCCTCTTCGAGGAAGATAAAGCGACGGGGCACCg                                                               <  1:1380181/76‑1 (MQ=255)
    cGAAGATGTTTACCGAACGGCGAGTGATTTCCGCAATATCGCCCTCTTCAAGGAAGATAAAGCGACGGGGCACCgg                                                              >  1:1329668/1‑76 (MQ=255)
     gAAGATGTTTACCGAACGGCGAGTGATTTCCGCAATATCGCCCTCTTCAAGGAAGATAAAGCGACGGGGCACCGGc                                                             <  2:1014976/76‑1 (MQ=255)
     gAAGATGTTTACCGAACGGCGAGTGATTTCCGCAATATCGCCCTCTTCAAGGAAGATAAAGCGACGGGGCACCGGc                                                             <  2:1014970/76‑1 (MQ=255)
          tGTTTACCGAACGGCGAGTGATTTCCGCAATATCGCCCTCTTCAAGGAAGATAAAGCGACGGGGCACCGGCAACAg                                                        <  2:93993/76‑1 (MQ=255)
               aCCGAACGGCGAGTGCTTTCCTTAATATCGCCCTCTTCAAGGAAGATAAAGCGACGGGGCACCGGCAACAGCGCCa                                                   <  2:46138/76‑1 (MQ=255)
               aCCGAACGGCGAGTGATTTCCGCAATATCGCCCTCTTCAAGGAAGATAAAGCGACGGGGCACCGGCAACAGCGCCa                                                   <  2:653457/76‑1 (MQ=255)
                ccGAACGGCGAGTGATTTCCGCAATATCGCCCTCTTCGAGGAAGATAAAGCGACGGGGCACCGGCAACAGCGCCAg                                                  <  2:1360717/76‑1 (MQ=255)
                  gAACGGCGAGTGATTTCCGCAATATCGCCCTCTTCAAGGAAGATAAAGCGACGGGGCACCGGCAACAGCGCCAGcc                                                >  2:146903/1‑75 (MQ=255)
                   aaCGGCGAGTGATTTCCGCAATATCGCCCTCTTCAAGGAAGATAAAGCGACGGGGCACCGGCAACAGCGCCAGCTg                                               <  2:37618/76‑1 (MQ=255)
                    aCGGCGAGTGAGTTCCGCAATATCGCCCGCTTCAAGGAAGATAAAGCGTCGGGGCACCGGCAACAGCGCCAGCTgg                                              <  2:131175/76‑1 (MQ=255)
                         gAGTGATTTCCGCAATATCGCCCTCTTCAAGGAAGATAAAGCGACGGGGCACCGGCAACAGCGCCAGCTGGTCAGa                                         >  1:1115632/1‑76 (MQ=255)
                         gAGTGATTTCCGCAATATCGCCCTCTTCAAGGAAGATAAAGCGACGGGGCACCGGCAACAGCGCCAGCTGGTCAGa                                         >  1:1023666/1‑76 (MQ=255)
                          aGTGATTTCCGCAATATCGCCCTCTTCCAGGAAGATAAAGCGACGGGGCACCGGCAACAGCGCCAGCTGCTCAGaa                                        >  1:1510478/1‑76 (MQ=255)
                           gTGATTTCCGCAATATCGCCCTCTTCAAGGAAGATAAAGCGACGGGGCACCGGCAACAGCGCCAGCTGGTCAGAAg                                       <  1:892008/76‑1 (MQ=255)
                             gATTTCCGCAATATCGCCCTCTTCAAGGAAGATAAAGCGACGGGGCACCGGCAACAGCGCCAGCTGGTCAGAAGCg                                     <  2:864347/76‑1 (MQ=255)
                               tttCCGCAATATCGCCCTCTTCAAGGAAGATAAAGCGACGGGGCACCGGCAACAGCGCCAGCTGGTCAGAAGCGAt                                   >  1:1117469/1‑76 (MQ=255)
                                ttCCGCAATATCGCCCTCTTCAAGGAAGATAAAGCGACGGGGCACCGGCAACAGCGCCAGCTGGTCAGAAGCGATa                                  >  2:568170/1‑76 (MQ=255)
                                     cAATATCGCCCTCTTCAAGGAAGATAAAGCGACGGGGCACCGGCAACAGCGCCAGCTGGTCAGAAGCGATAAAGtt                             >  2:475204/1‑76 (MQ=255)
                                       atatCGCCCTCTTCAAGGAAGATAAAGCGACGGGGCACCGGCAACAGCGCCAGCTGGTCAGAAGCGATAAAGtttt                           <  2:827032/76‑1 (MQ=255)
                                           cGCCCTCTTCAAGGAAGATAAAGCGACGGGGCACCGGCAACAGCGCCAGCTGGTCAGAAGCGATAAAGTTTTCGcc                       >  2:416218/1‑76 (MQ=255)
                                            gCCCTCTTCAAGGAAGATAAAGCGACGGGGCACCGGCAACAGCGCCAGCTGGTCAGAAGCGATAAAGTTTTCGccc                      >  1:1129492/1‑76 (MQ=255)
                                             cccTCTTCAAGGAAGATAAAGCGACGGGGCACCGGCAACAGCGCCAGCTGGTCAGAAGCGATAAAGTTTTCGCCCa                     <  2:1364795/76‑1 (MQ=255)
                                              ccTCTTCAAGGAAGATAAAGCGACGGGGCACCGGCAACAGCGCCAGCTGGTCAGAAGCGATAAAGTTTTCGCCCAt                    >  1:1115385/1‑76 (MQ=255)
                                                  ttCAAGGAAGATAAAGCGACGGGGCACCGGCAACAGCGCCAGCTGGTCAGAAGCGATAAAGTTTTCGCCCATcccc                <  2:1534528/76‑1 (MQ=255)
                                                  ttCAAGGAAGATAAAGCGACGGGGCACCGGCAACAGCGCCAGCTGGTCAGAAGCGATAAAGTTTTCGCCCATcccc                >  1:650941/1‑76 (MQ=255)
                                                         aaGATAAAGCGACGGGTCACCGGCAACAGCGCCAGCTGGTCAGAAGCGATAAAGTTTTCGCCCATCCCCAGGCCaa         <  2:650684/76‑1 (MQ=255)
                                                           gATAAAGCGACGGGGCACCGGCAACAGCGCCAGCTGGTCAGAAGCGATAAAGTTTTCGCCCATCCCCAGGCCAATc       <  2:1107339/76‑1 (MQ=255)
                                                             tAAAGCGACGGGGCACCGGCAACAGCGCCAGCTGGTCAGAAGCGATAAAGTTTTCGCCCATCCCCATGCCAATCAc     >  1:54448/1‑76 (MQ=255)
                                                             tAAAGCGACGGGGCACCGGCAACAGCGCCAGCTGGTCAGAAGCGATAAAGTTTTCGCCCATCCCCAGGCCAATCAc     <  1:1219709/76‑1 (MQ=255)
                                                              aaaGCGACGGGGCACCGGCAACAGCGCCAGCTGNTCAGAAGCGATAAAGTTTTCGCNCATCCCCAAGNCAATCAcc    >  2:1276458/1‑76 (MQ=255)
                                                               aaGCGACGGGGCACCGGCAACAGCGCCAGCTGGTCAGAAGCGATAAAGTTTTCGCCCATCCCCAGGCCAATCACCa   >  2:589177/1‑76 (MQ=255)
                                                                aGCGACGGGGCACCGGCAACAGCGCCAGCTGGTCAGAAGCGATAAAGTTTTCGCCCATCCCCAGGCCAATCACCAg  >  2:385360/1‑76 (MQ=255)
                                                                         |                                                                  
TTATCGAAGATGTTTACCGAACGGCGAGTGATTTCCGCAATATCGCCCTCTTCAAGGAAGATAAAGCGACGGGTCACCGGCAACAGCGCCAGCTGGTCAGAAGCGATAAAGTTTTCGCCCATCCCCAGGCCAATCACCAG  >  REL606/3873361‑3873500

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