Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 2,294,105 T→G 19.1% intergenic (+132/‑143) gldA → / → yijF glycerol dehydrogenase, NAD/conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported2,294,1050TG19.1% 87.0 / 18.1 42intergenic (+132/‑143)gldA/yijFglycerol dehydrogenase, NAD/conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (22/12);  new base G (4/4);  total (26/16)
Fisher's exact test for biased strand distribution p-value = 4.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

ATCAGGCCTACGCAATCTCTGCAATATATTGAATTTGCGTGCTTTTGTAGGCCGGATAAGGCGTTCACTCCGCATCCGGCACTCTGTGCCAGTAATCTGAAACTCCCACGTTGTTCAGGAATTTTTTATCCGCTTCT  >  Exported/2294037‑2294173
                                                                    |                                                                    
aTCAGGCCTACGCAATCTCTGCAATATATTGAATTTGCGTGCTTTTGTAGGCCGGATAAGGCGTTCACTcc                                                                    >  1:1376865/1‑71 (MQ=255)
  cAGGCCTACGCAATCTCTGCAATATATTGAATTTGCGTGCTTTTGTAGGCCGGATAAGGCGTTCACTCCg                                                                   >  1:4720792/1‑70 (MQ=255)
   aGGCCTACGCAATCTCTGCAATATATTGAATTTGCGTGCTTTTGTAGGCCGGATAAGGCGTTCACTCCGCa                                                                 <  1:3097563/71‑1 (MQ=255)
     gCCTACGCCATCTCTGCAATATATTGAATTTGCTTGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATc                                                               <  1:5565589/71‑1 (MQ=14)
     gCCTACGCCATCTCTGCAATATATTGAATTTGCGTGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATc                                                               <  1:350415/71‑1 (MQ=14)
       cTACGCCATCTCTGCAATATATTGAATTTGCGTGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCg                                                             >  1:3390699/1‑71 (MQ=11)
          cGCCATCTCTGCAATATATTGAATTTGCGTGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCgg                                                            <  1:2120277/69‑1 (MQ=11)
          cGCCATCTCTGCAATATATTGAATTTGCGTGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCgg                                                            <  1:1258126/69‑1 (MQ=11)
          cGCCATCTCTGCAATATATTGAATTTGCGTGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCa                                                          >  1:3472742/1‑71 (MQ=11)
          cGCCATCTCTGCAATATATTGAATTTGCGTGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCa                                                          >  1:2113870/1‑71 (MQ=11)
          cGCCATCTCTGCAATATATTGAATTTGCGTGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCa                                                          >  1:4381388/1‑71 (MQ=11)
          cGCAATCTCTGCAATATAATGAATTTGCGTGCTTTTGTAGGCCGGATAAGGCGTTCACTCCGCTTCCgg                                                            <  1:4109549/69‑1 (MQ=255)
                ctctGCAATATATTGAATTTGCGTGCTTTTGTAAGCCGGATAAGGCGTTCACTCCGCATCCGGCACTCtgt                                                    >  1:2420119/1‑71 (MQ=255)
                 tctGCAATATATTGAATTTGCGTGCTTTTGTAGGCCGGATAAGGCGTTCACTCCGCATCCGGCa                                                          >  1:4102370/1‑64 (MQ=40)
                  ctGCAATATATTGAATTTGCGTGCTTTTTTAGGCCGGATAAGGCGTTCACTCCGCATCCGGCACTCTGTGc                                                  >  1:5678167/1‑71 (MQ=255)
                   tGCAATATATTGAATTTGCGTGCTTTTGTAGGCCGGATAAGGCGTTCACTCCGCATCCGGCACTCTGTGc                                                  >  1:4100520/1‑70 (MQ=255)
                      aaTATATTGAATTTGCGTGCTTTTGTAGGCCGGATAAGGCGTTCACTCCGCATCCGGCActct                                                      >  1:3623469/1‑63 (MQ=255)
                       atatatTGAATTTGCGTGCTTTTGTAGGCCGGATAAGGCGTTCACTCCGCATCCGGCACTCTGTGCCAGt                                              <  1:1771631/70‑1 (MQ=255)
                       atatatTGAATTTGCGTGCTTTTGTAGGCCGGATAAGGCGTTCACTCCGCATCCGGCACTCTGTGCCAGt                                              >  1:1044988/1‑70 (MQ=255)
                       atatatTGAATTTGCGTGCTTTTGTAGGCCGGATAAGGCGTTCACTCCGCATCCGGCACTCTGTGCCAGTa                                             <  1:2666049/71‑1 (MQ=255)
                       atatatTGAATTTGCGTGCTTTTGTAGGCCGGATAAGGCGTTCACTCCGCATCCGGCACTCTGTGCCAGTa                                             <  1:3342011/71‑1 (MQ=255)
                            ttGAATTTGCGTGCTTTTGTAGGCCGGATAAGGCGTTCACTCCGCATCCGGCACTCTGTGCCAGTAATCt                                         >  1:540453/1‑70 (MQ=255)
                              gAATTTGCGTGCTTTTGTAGGCCTGATAAGGCGTTCACTCCGCATCCGGCACTCTGTGCCAGTAATCTGa                                       <  1:770221/70‑1 (MQ=255)
                              gAATTTGCGTGCTTTTGTAGGCCGGATAAGGCGTTCACTCCGCATCCGGCACTCTGTGCCAGTAATCTGa                                       <  1:2242911/70‑1 (MQ=255)
                                aTTTGCGTGCTTTTGTAGGCCGGATAAGGCGTTCACTCCGCATCCGGCACTCTGTGCCAGTAATCt                                         >  1:3457347/1‑66 (MQ=255)
                                aTTTGCGTGCTTTTGTAGGCCGGATAAGGCGTTCACTCCGCATCCGGCACTCTGTGCCAGTAATCt                                         >  1:3459319/1‑66 (MQ=255)
                                 tttGCGTGCTTTTGTAGGCCGGATAAGGCGTTCACTc                                                                     <  1:150415/37‑1 (MQ=37)
                                  ttGCGTGCTTTTGTAGGCCGGATAAGGCGTTCACTCCGCATCCGGCACTCTGTGCCAGTAAt                                           >  1:815692/1‑62 (MQ=255)
                                  ttGCGTGCTTTTGTAGGCCGGATAAGGCGTTCACTCCGCATCCGGCACTCTGTGCCAGTAAt                                           >  1:2507610/1‑62 (MQ=255)
                                   tGCGTGCTTTTGTAGGCCGGATAAGGCGTTCACTCCGCATcc                                                              >  1:3887064/1‑42 (MQ=38)
                                             tGTAGGCCGGATAAGGCGTTCACTCCGCATCCGGCACTCTGTGCCAGTAATCTGAAACTCCCACgttgtt                        >  1:543986/1‑70 (MQ=255)
                                              gTAGGCCGGATAAGGCGTTCACTCCGCATCCGGCActc                                                       <  1:1144228/38‑1 (MQ=40)
                                                   ccGGATAAGGCGTTCACTCCGCATCCGGCACTCTGTGCCAGTAATCTGAAACTCCCACGTTGTTc                       >  1:1958157/1‑65 (MQ=255)
                                                    cGGATAAGGCGTTCACTCCGCATCCGGCACTCTGTGCCAGTAATCTGAAACTCCCACGTTGTTCAGGAAtt                >  1:5140115/1‑71 (MQ=255)
                                                     ggATAAGGCGTTCACTCCGCATCCGGCACTCTGTGCCAGTAATCTGAAACTCCCACGTTGTTCAGGAAttt               <  1:3460149/71‑1 (MQ=255)
                                                     ggATAAGGCGTTCACTCCGCATCCGGCACTCTGTGCCAGTAATCTGAAACTCCCACGTTGTTCAGGAAttt               <  1:2508091/71‑1 (MQ=255)
                                                        tAAGGCGTTCACTCCGCATCCGGCACTCTGTGCCAGTAATCTGAAACTCCCACGTTGTTCAGGAAttttt             >  1:2000649/1‑70 (MQ=255)
                                                        tAAGGCGTTCACTCCGCATCCGGCACTCTGTGCCAGTAATCTGAAACTCCCACGTTGTTCAGGAAttttt             >  1:1721810/1‑70 (MQ=255)
                                                                 cACTCCGCATCCGGCACTCTGTGCCAGTAATCTGAAACTCCCACGTTGTTCAg                     >  1:1550965/1‑53 (MQ=255)
                                                                 cACTCCGCATCCGGCACTCTGTGCCAGTAATCTGAAACTCCCACGTTGTTCAGGAAtttttt            >  1:2396311/1‑62 (MQ=255)
                                                                 cACTCCGCATCCGGCACTCTGTGCCAGTAATCTGAAACTCCCACGATGTTCAGGAAttttt             >  1:3135170/1‑61 (MQ=255)
                                                                  aCTCCGCATCCGGCACTCTGTGCCAGTAATCTGAAACTCCCACGTTGTTCAGGAATTTTTTATCCGCTTCt  <  1:3776280/71‑1 (MQ=255)
                                                                    |                                                                    
ATCAGGCCTACGCAATCTCTGCAATATATTGAATTTGCGTGCTTTTGTAGGCCGGATAAGGCGTTCACTCCGCATCCGGCACTCTGTGCCAGTAATCTGAAACTCCCACGTTGTTCAGGAATTTTTTATCCGCTTCT  >  Exported/2294037‑2294173

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