Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 2,872,245 G→C 13.8% A166P (GCC→CCC)  ulaF → L‑ribulose 5‑phosphate 4‑epimerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported2,872,2450GC13.8% 157.6 / 16.3 65A166P (GCC→CCC) ulaFL‑ribulose 5‑phosphate 4‑epimerase
Reads supporting (aligned to +/- strand):  ref base G (31/25);  new base C (6/3);  total (37/28)
Fisher's exact test for biased strand distribution p-value = 7.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.48e-01

TGATTATCGAAACGCTGGGCAACGCCGAGCCGCTGCATACGCCGGGAATTGTGGTGTATCAGCACGGGCCGTTCGCCTGGGGGAAAGATGCTCACGATGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGAAAA  >  Exported/2872178‑2872315
                                                                   |                                                                      
tGATTATCGAAACGCTGGGCAACGCCGAGCCGCTGCATACGCCGGGAATTGTGGTGTATCAGCACGGGCCg                                                                     <  1:4367174/71‑1 (MQ=255)
  aTTATCGAAACGCTGGGCAACGCCGAGCCGCTGCATACGCCGGGAATTGTGGTGTATCAGCACAGGCCGt                                                                    <  1:2309922/70‑1 (MQ=255)
   ttATCGAAACGCTGGGCAACGCCGAGCCGCTGCATACGCCGGGAATTGTGGTGTATCAGCACGGGCCGTTc                                                                  <  1:4236091/71‑1 (MQ=255)
    tATCGAAACGCTGGGCAACGCCGAGCCGCTGCATACGCCGGGAATTGTGGTGTATCAGCACGGGCCGTTCg                                                                 >  1:5618075/1‑71 (MQ=255)
     aTCGAAACGCTGGGCAACGCCGAGCCGCTGCATACGCCGGGAATTGTGGTGTATCAGCACGGGCCGTTCGc                                                                >  1:5418683/1‑71 (MQ=255)
          aaCGCTGGGCAACGCCGAGCCGCTGCATACGCCGGGAATTGTGGGGTATCAGCACGGGCCGTTCGc                                                                <  1:5356699/66‑1 (MQ=255)
           aCGCTGGGCAACGCCGAGCCGCTGCATACGCCGGGAATTGTGGTGTATCAGCACGGGCCGTTCGCCTgggg                                                          >  1:4287228/1‑71 (MQ=255)
                 ggCAACGCCGAGCCGCTGCATACGCCGGGAATTGTGGTGTATCAGCACGGGCCGTTCGCCTGGGGGaaa                                                      <  1:431422/69‑1 (MQ=255)
                     aCGCCGAGCCGCTGCATACGCCGGGAATTGTGGTGTATCAGCACGGGCCGTTCGCCTGGGGGAAAGATGCt                                                >  1:1056150/1‑71 (MQ=255)
                      cGCCGAGCCGCTGCATACGCCGGGAATTGTGGTGTATCAGCACGGGCCGTTCGCCTGGGGGAAAGATGCt                                                >  1:412421/1‑70 (MQ=255)
                       gCCGAGCCGCTGCATACGCCGGGAATTGTGGTGTATCAGCACGGGCCGTTCGCCTGGGGGAAAGATGCTCa                                              >  1:2834798/1‑71 (MQ=255)
                       gCCGAGCCGCTGCATACGCCGGGAATTGTGGTGTATCAGCACGGGCCGTTCGCCTGGGGGAAAGATGCTCa                                              >  1:3175748/1‑71 (MQ=255)
                       gCCGAGCCGCTGCATACGCCGGCAATTGTGGTGTATCAGCACGGGCCGTTGGCCTGGGGGAAAGATGCTCa                                              >  1:2461073/1‑71 (MQ=255)
                        ccGAGCCGCTGCATACGCCGGGAATTGTGGTGTATCAGCACGGGCCGTTCGCCTGGGGGAAAGATGCTCAc                                             <  1:943868/71‑1 (MQ=255)
                        ccGAGCCGCTGCATACGCCGGGAATTGTGGTGTATCAGCACGGGCCGTTCGCCTGGGGGAAAGATGCTCAc                                             <  1:157381/71‑1 (MQ=255)
                            gCCGCTGCATACGCCGGGAATTGTGGTGTATCAGCACGGGCCGTTCGCCTGGGGGAAAGATGCTc                                               >  1:5437778/1‑65 (MQ=255)
                            gCCGCTGCATACGCCGGGAATTGTGGTGTATCAGCACGGGCCGTTCGCCTGGGGGAAAGATGCTCACGATg                                         >  1:2041557/1‑71 (MQ=255)
                                  gCATACGCCGGGAATTGTGGTGTATCAGCACGGGCCGTTCGCCTgggg                                                          <  1:956092/48‑1 (MQ=255)
                                  gCATACGCCGGGAATTGTGGTGTATCAGCACGGGCCGTTCGCCTGGGGGaaa                                                      >  1:2585054/1‑52 (MQ=255)
                                    aTACGCCGGGAATTGTGGTGTATCAGCACGGGCCGTTCGCCTGGGGGAAAGATGCTCACGATGCGGTGCAt                                 <  1:2226036/71‑1 (MQ=255)
                                    aTACGCCGGGAATTGTGGTGTATCAGCACGGGCCGTTCGCCTGGGGGAAAGATGCTCACGATGCGGTGCAt                                 <  1:3846571/71‑1 (MQ=255)
                                    aTACGCCGGGAATTGTGGTGTATCAGCACGGGCCGTTCGCCTGGGGGAAAGATGCTCACGATGCGGTGCAt                                 <  1:3823442/71‑1 (MQ=255)
                                     tACGCCGGGAATTGTGGTGTATCAGCACGGGCCGTTCGCCTGGGGGAAAGATGCTCACGATGCGGTGCAt                                 >  1:875985/1‑70 (MQ=255)
                                         ccGGGAATTGTGGTGTATCAGCACGGGCCGTTCGCCt                                                              >  1:807407/1‑37 (MQ=255)
                                         ccGGGAATTGTGGTGTATCAGCACGGGCCGTTCGCCt                                                              >  1:3344338/1‑37 (MQ=255)
                                             gAATTGTGGTGTATCAGCACGGGCCGTTCGCCTGGGGGAAAGATGCTCACGATg                                         >  1:5521392/1‑54 (MQ=255)
                                             gAATTGTGGTGTATCAGCACGGGCCGTTCGCCTGGGGGAAAGATGCTCACGATg                                         >  1:22326/1‑54 (MQ=255)
                                              aaTTGTGGTGTATCAGCACGGGCCGTTCGCCTGGGGGaaa                                                      <  1:4718489/40‑1 (MQ=255)
                                                ttGTGGTGTATCAGCACGGGCCGTTCGCCTGGGGGAAAGATGCTCACGATGCGGTGCATAACGCggtggt                      >  1:923546/1‑70 (MQ=255)
                                                ttGTGGTGTATCAGCACGGGCCGTTCGCCTGGGGGAAAGATGCTCACGATGCGGTGCATAACGCggtggt                      >  1:2306317/1‑70 (MQ=255)
                                                  gtggtgTATCAGCACGGGCCGTTCGCCTGGGGGAAAg                                                     >  1:3587105/1‑37 (MQ=255)
                                                  gtggtgTATCAGCACGGGCCGTTCGCCTGGGGGAAAg                                                     >  1:2071874/1‑37 (MQ=255)
                                                    ggtgTATCAGCACGGGCCGTTCGCCTGGGGGAAAGATGCTCACGATGCGGTGCATAACGCGGTGGTGATg                  <  1:2067678/70‑1 (MQ=255)
                                                     gtgtATCAGCACGGGCCGTTCGCCTGGGGGAAAGATGCTCACGATGCGGTGCATAACGCGGTGGTGATgg                 <  1:2211007/70‑1 (MQ=255)
                                                      tgtATCAGCACGGGCCGTTCGCCTGGGGGAAAGATGCTCAc                                             <  1:1700472/41‑1 (MQ=255)
                                                        tATCAGCACGGGCCGTTCGCCTGGGGGAAAGATGCTCa                                              <  1:3323955/38‑1 (MQ=255)
                                                        tATCAGCACGGGCCGTTCGCCTGGGGGAAAGATGCTCa                                              <  1:1127238/38‑1 (MQ=255)
                                                        tATCAGCACGGGCCGTTCGCCTGGGGGAAAGATGCTCACGATGc                                        <  1:2000538/44‑1 (MQ=255)
                                                        tATCAGCACGGGCCGTTCGCCTGGGGGAAAGATGCTCACGATGCGGTGCATAACGCggtgg                       >  1:4805118/1‑61 (MQ=255)
                                                        tATCAGCACGGGCCGTTCGCCTGGGGGAAAGATGCTCACGATGCGGTGCATAACGCggtgg                       >  1:3860845/1‑61 (MQ=255)
                                                        tATCAGCACGGGCCGTTCGCCTGGGGGAAAGATGCTCACGATGCGGTGCATAACGCggtgg                       >  1:2581399/1‑61 (MQ=255)
                                                        tATCAGCACGGGCCGTTCGCCTGGGGGAAAGATGCTCACGATGCGGTGCATAACGCGGTGGTGATGgaaga             <  1:4604790/71‑1 (MQ=255)
                                                        tATCAGCACGGGCCGTTCGCCTGGGGGAAAGATGCTCACGATGCGGTGCATAACGCGGTGGTGATGgaag              >  1:4285458/1‑70 (MQ=255)
                                                          tCAGCACGGGCCGTTCGCCTGGGGGAAAGATGCTCACGATGCGGTGCATAAcgcg                           >  1:1051697/1‑55 (MQ=255)
                                                          tCAGCACGGGCCGTTCGCCTGGGGGAAAGATGCTCACGATGCGGTGCATAAcgcg                           >  1:1819863/1‑55 (MQ=255)
                                                          tCAGCACGGGCCGTTCGCCTGGGGGAAAGATGCTCACGATGCGGTGCATAACGCGGTGGTGATGgaagaa            <  1:4547312/70‑1 (MQ=255)
                                                            cccACGGCCCGTTTGCATGGGGCAAAAATGCCGAAGATGCGGTGCATAACGCCATCGTGCTGgaa               >  1:848685/3‑65 (MQ=255)
                                                            cccACGGCCCGTTTGCATGGGGCAAAAATGCCGAAGATGCGGTGCATAACGCCATCGTGCTGGAAGAGGtc         >  1:4479915/3‑70 (MQ=255)
                                                            cccACGGCCCGTTTGCATGGGGCAAAAATGCCGAAGATGCGGTGCATAACGCCATCGTGCTGGAAGAGGtc         >  1:2887968/3‑70 (MQ=255)
                                                            cccACGGCCCGTTTGCATGGGGCAAAAATGCCGAAGATGCGGTGCATAACGCCATCGTGCTGGAAGAGGt          >  1:1172157/3‑70 (MQ=255)
                                                            aGCACGGGCCGTTCGCCTGGGGGAAAGATGCTCACGATGCGGTGCATAACGCGGTGGTGATGGAAGAAGt          <  1:4557762/70‑1 (MQ=255)
                                                            aGCACGGGCCGTTCGCCTGGGGGAAAGATGCTCACGATGCGGTGCATAACGCGGTGGTGATGGAAGAAGTg         >  1:957472/1‑71 (MQ=255)
                                                            aGCACGGGCCGTTCGCCTGGGGGAAAGATGCTCACGATGCGGTGCATAACGCGGTGGTGATGGAAGAAGTg         <  1:3622580/71‑1 (MQ=255)
                                                            aGCACGGGCCGTTCGCCTGGGGGAAAGATGCTCACGATGCGGTGCATAACGCGGTGGTGATGGAAGAAGTg         >  1:374889/1‑71 (MQ=255)
                                                            aGCACGGGCCGTTCGCCTGGGGGAAAGATGCTCACGATGCGGTGCATAACGCGGTGGTGATGGAAGAAGTg         <  1:4289954/71‑1 (MQ=255)
                                                             ccACGGCCCGTTTGCATGGGGCAAAAATGCCGAAGATGCGGTGCATAACGCCATCGTGCTGGAAGAGGtc         >  1:3668662/2‑69 (MQ=255)
                                                              cACGGCCCGTTTGCATGGGGCAAAAATGCCGAAGATGCGGTGCATAACGCCATCGTGCTGGAAGAGGTCGc       <  1:1034563/71‑1 (MQ=255)
                                                               aCGGGCCGTTCGCCTGGGGGAAAGATGCTCACGATGCGGTGCATAACGCGGTGGTGATGGAAGAAGt          >  1:3283436/1‑67 (MQ=255)
                                                               aCGGGCCGTTCGCCTGGGGGAAAGATGCTCACGATGCGGTGCATAACGCGGTGGTGATGGAAGAAGt          >  1:5679899/1‑67 (MQ=255)
                                                               aCGGCCCGTTTGCATGGGGCAAAAATGCCGAAGATGCGGTGCATAACGCCATCGTGCTGGAAGAGGTCGct      >  1:648511/1‑70 (MQ=255)
                                                               aCGGCCCGTTTGCATGGGGCAAAAATGCCGAAGATGCGGTGCATAACGCCATCGTGCTGGAAGAGGTCGct      <  1:744477/71‑2 (MQ=255)
                                                               aCGGCCCGTTCGCATGGGGTAAAAACGCCGCCGATGCCGTGCATAAcgccg                          <  1:2246496/51‑3 (MQ=255)
                                                                cGGGCCGTTCGCCTGGGGGAAAGATGCTCACGATGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCg      >  1:3421720/1‑70 (MQ=255)
                                                                cGGGCCGTTCGCCTGGGGGAAAGATGCTCACGATGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGa     <  1:5533142/71‑1 (MQ=255)
                                                                cGGGCCGTTCGCCTGGGGGAAAGATGCTCACGATGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGa     <  1:394521/71‑1 (MQ=255)
                                                                 gggCCGTTCGCCTGGGGGAAAGATGCTCACGATGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGa     >  1:80711/1‑70 (MQ=255)
                                                                   gCCGTTCGCCTGGGGGAAAGATGCTCACGATGCGGTGCATAACGCGGTGGTGAt                   >  1:4939775/1‑54 (MQ=255)
                                                                   gCCGTTCGCCTGGGGGAAAGATGCTCACGATGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGaaaa  >  1:209393/1‑71 (MQ=255)
                                                                   gCCGTTCGCCTGGGGGAAAGATGCTCACGATGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGaaaa  <  1:4340449/71‑1 (MQ=255)
                                                                   |                                                                      
TGATTATCGAAACGCTGGGCAACGCCGAGCCGCTGCATACGCCGGGAATTGTGGTGTATCAGCACGGGCCGTTCGCCTGGGGGAAAGATGCTCACGATGCGGTGCATAACGCGGTGGTGATGGAAGAAGTGGCGAAAA  >  Exported/2872178‑2872315

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: