Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 2,336,250 T→C 5.9% G212G (GGT→GGC fdoH → formate dehydrogenase‑O, Fe‑S subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported2,336,2500TC 5.9% 203.9 / 3.0 68G212G (GGT→GGCfdoHformate dehydrogenase‑O, Fe‑S subunit
Reads supporting (aligned to +/- strand):  ref base T (23/41);  new base C (2/2);  total (25/43)
Fisher's exact test for biased strand distribution p-value = 6.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.71e-01

AGAGTCGATGAAAACGCTGGCGAGCGAGCGCGTTGCTGAGCTGAAAACCCGCGGTTACGACAATGCGGGTCTGTACGATCCGGCAGGCGTCGGTGGTACACACGTCATGTACGTGCTGCACCATGCTGACAAGCCAAAT  >  Exported/2336181‑2336319
                                                                     |                                                                     
agagTCGATGAAAACGCTGGCGAGCGAGCGCGTTGCTGAGCTGAAAACCCGCGGTTACGACAATGCGGGTc                                                                      <  1:2723309/71‑1 (MQ=255)
agagTCGATGAAAACGCTGGCGAGCGAGCGCGTTGCTGAGCTGAAAACCCGCGGTTACGACAATGCGGGTc                                                                      <  1:1563383/71‑1 (MQ=255)
   gTCGATGAAAACGCTGGCGAGCGAGCGCGTTGCTGAGCTGAAAACCCGCGGTTACGACAATGCGGGTCTGt                                                                   <  1:3052345/71‑1 (MQ=255)
   gTCGATGAAAACGCTGGCGAGCGAGCGCGTTGCTGAGCTGAAAACCCGCGGTTACGACAATGCGGGTCTGt                                                                   <  1:2390730/71‑1 (MQ=255)
   gTCGATGAAAACGCTGGCGAGCGAGCGCGTTGCTGAGCTGAAAACCCGCGGTTACGACAATGCGGGTCTGt                                                                   <  1:1381232/71‑1 (MQ=255)
   gTCGATGAAAACGCTGGCGAGCGAGCGCGTTGCTGAGCTGAAAACCCGCGGTTACGACAATGCGGGTCTGt                                                                   <  1:2323148/71‑1 (MQ=255)
   gTCGATGAAAACGCTGGCGAGCGAGCGCGTTGCTGAGCTGAAAACCCGCGGTTACGACAATGCGGGTCTGt                                                                   <  1:2105256/71‑1 (MQ=255)
      gATGAAAACGCTGGCGAGCGAGCGCGTTGCTGAGCTGAAAACCCGCGGTTACGACAATGCGGGTCTGTa                                                                  >  1:236035/1‑69 (MQ=255)
      gATGAAAACGCTGGCGAGCGAGCGCGTTGCTGAGCTGAAAACCCGCGGTTACGACAATGCGGGTCTGTAc                                                                 >  1:2362534/1‑70 (MQ=255)
      gATGAAAACGCTGGCGAGCGAGCGCGTTGCTGAGCTGAAAACCCGCGGTTACGACAATGCGGGTCTGTAc                                                                 >  1:824140/1‑70 (MQ=255)
        tGAAAACGCTGGCGAGCGAGCGCGTTGCTGAGCTGAAAACCCGCGGTTACGACAATGCGGGTCTGTACGa                                                               <  1:1940860/70‑1 (MQ=255)
        tGAAAACGCTGGCGAGCGAGCGCGTTGCTGAGCTGAAAACCCGCGGTTACGACAATGCGGGTCTGTACGAt                                                              >  1:1734742/1‑71 (MQ=255)
          aaaaCGCTGGCGAGCGAGCGCGTTGCTGAGCTGAAAACCCGCGGTTACGACAATGCGGGTCTGTACGATcc                                                            <  1:453209/71‑1 (MQ=255)
          aaaaCGCTGGCGAGCGAGCGCGTTGCTGAGCTGAAAACCCGCGGTTACGACAATGCGGGTCTGTACGATcc                                                            <  1:1836182/71‑1 (MQ=255)
              cGCTGGCGAGCGAGCGCGTTGCTGAGCTGAAAACCCGCGGTTACGACAATGCGGGTCTGTACGATCCGGCa                                                        <  1:466198/71‑1 (MQ=255)
                  ggCGAGCGAGCGCGTTGCTGAGCTGAAAACCCGCGGTTACGACAATGCGGGTCTGTACGATCCGGCAGGCg                                                    >  1:3371406/1‑71 (MQ=255)
                   gcgagcgaGCGCGTTGCTGAGCTGAAAACCCGCGGTTACGACAATGCGGGTCTGTACGATCCGGCAGGCGt                                                   <  1:2292254/71‑1 (MQ=255)
                    cgagcgagCGCGTTGCTGAGCTGAAAACCCGCGGTTACGACAATGCGGGTCTGTACGATCCGGCAGGCGTc                                                  <  1:400909/71‑1 (MQ=255)
                    cgagcgagCGCGTTGCTGAGCTGAAAACCCGCGGTTACGACAATGCGGGTCTGTACGATCCGGCAGGCGTc                                                  <  1:3335136/71‑1 (MQ=255)
                      agcgagcgCGTTGCTGAGCTGTAAACCCGCGGTTACGACAATGCGGGTc                                                                      <  1:1727840/49‑1 (MQ=255)
                          agcgCGTTGCTGAGCTGAAAACCCGCGGTTACGACAATGCGGGTCTGTACGATCCGGCAGGCGTCggtgg                                             <  1:420473/70‑1 (MQ=255)
                          agcgCGTTGCTGAGCTGAAAACCCGCGGTTACGACAATGCGGGTCTGTACGATCCGGCAGGCGTCgg                                                <  1:795188/67‑1 (MQ=255)
                          agcgCGTTGCTGAGCTGAAAACCCGCGGTTACGACAATGCGGGTCTGTACGATCCGGCAGGCGTCgg                                                <  1:496638/67‑1 (MQ=255)
                           gcgcGTTGCTGAGCTGAAAACCCGCGGTTACGACAATGCGGGTCTGTACGATCCGGCAGGCGTCGGTGGTa                                           >  1:2849245/1‑71 (MQ=255)
                                ttGCTGAGCTGAAAACCCGCGGTTACGACAATGCGGGTCTGTACGATCCGGCAGGCGTCGGTGGTa                                           <  1:3036386/66‑1 (MQ=255)
                                 tGCTGAGCTGAAAACCCGCGGTTACGACAATGCGGGTCTGTACGATCCGGCa                                                        >  1:1760008/1‑52 (MQ=255)
                                  gctgagctgaAAACCCGCGGTTACGACAATGCGGGTCTGTACGATCCGGCAGGCGTCGGTGGTACACACg                                     <  1:754223/70‑1 (MQ=255)
                                  gctgagctgaAAACCCGCGGTTACGACAATGCGGGTCTGTACGATCCGGCAGGCGTCGGTGGTACACACg                                     <  1:1162454/70‑1 (MQ=255)
                                  gctgagctgaAAACCCGCGGTTACGACAATGCGGGTCTGTACGATCCGGCAGGCGTCGGTGGTACACACg                                     <  1:1127972/70‑1 (MQ=255)
                                  gctgagctgaAAACCCGCGGTTACGACAATGCGGGTCTGTACGATCCGGCAGGCGTCGGTGGTACACACGt                                    >  1:2402768/1‑71 (MQ=255)
                                     gagctgaAAACCCGCGGTTACGACAATGCGGGTCTGTACGATCCGGCAGGCGTCGGTGGTa                                           <  1:329613/61‑1 (MQ=255)
                                      agctgaAAACCCGCGGTTACGACAATGCGGGTCTGTACGATCCGGCAGGCGTCGGTGGTACACACGTCa                                  <  1:874254/69‑1 (MQ=255)
                                      agctgaAAACCCGCGGTTACGACAATGCGGGTCTGTACGATCCGGCAGGCGTCGGTGGTACACACGTCAt                                 <  1:1038687/70‑1 (MQ=255)
                                      agctgaAAACCCGCGGTTACGACAATGCGGGTCTGTACGATCCGGCAGGCGTCGGTGGTACACACGTCATg                                >  1:2259923/1‑71 (MQ=255)
                                       gctgaAAACCCGCGGTTACGACAATGCGGGTCTGTACGATCCGGCAGGCGTCGGTGGTACACACGTCATg                                >  1:1785591/1‑70 (MQ=255)
                                       gctgaAAACCCGCGGTTACGACAATGCGGGTCTGTACGATCCGGCAGGCGTCGGTGGTACACACGTCATGt                               <  1:607436/71‑1 (MQ=255)
                                       gctgaAAACCCGCGGTTACGACAATGCGGGTCTGTACGATCCGGCAGGCGGCGGTGGTACACACGTCATGt                               <  1:1684724/71‑1 (MQ=255)
                                        ctgaAAACCCGCGGTTACGACAATGCGGGTCTGTACGATCCGGCAGGCGTCGGTGGTACACACGTCATGTa                              >  1:2765900/1‑71 (MQ=255)
                                        ctgaAAACCCGCGGTTACGACAATGCGGGTCTGTACGATCCGGCAGGCGTCGGTGGTACACACGTCATGTa                              >  1:3256528/1‑71 (MQ=255)
                                        ctgaAAACCCGCGGTTACGACAATGCGGGTCTGTACGATCCGGCAGGCGTCGGTGGTACACACGTCATGTa                              >  1:84511/1‑71 (MQ=255)
                                           aaaaCCCGCGGTTACGACAATGCGGGTCTGTACGATCCGGCAGGCGTCGGTGGTACACACGTCATGTACg                            <  1:3298997/70‑1 (MQ=255)
                                                  gcgGTTACGAGAACGCCGGCCTGTACGACCCCGAAGGGGTGGGCGGCACCCACGTGATGTACGTGCTGCAc                    <  1:130527/71‑1 (MQ=2)
                                                   cgGTTACGAGAACGCCGGCCTGTACGACCCCGAAGGGGTGGGCGGCACCCACGTGATGTACGTGCTGCAc                    <  1:489394/70‑1 (MQ=2)
                                                          gACAATGCGGGTCTGTACGATCCGGCAGGCGTCGGTGGTACACACGTCATGTACGtgc                         <  1:1227424/58‑1 (MQ=255)
                                                           aCAATGCGGGTCTGTACGATCCGGCAGGCGTCGGTGGTACACACGTCATGTACGTGCTGCACCATGCTGAc           <  1:3214027/71‑1 (MQ=255)
                                                            cAATGCGGGTCTGTACGATCCGGCAGGCGTCGGTGGTACACACGTCATGTACGtgctg                       >  1:475401/1‑58 (MQ=255)
                                                            cAAAGCGGGTCTGTACGATCCGGCAGGCGTCGGTGGTACACACGTCATGTACGtgctg                       >  1:280079/1‑58 (MQ=255)
                                                             aaTGCGGGTCTGTACGATCCGGCAGGCGTCGGTGGTACACACGTCATGTACGTGCTGCACCATGCTGAc           <  1:3378435/69‑1 (MQ=255)
                                                                gCGGGTCTGTACGATCCGGCAGGCGTCGGTGGTaca                                         >  1:1864720/1‑36 (MQ=255)
                                                                gCGGGTCTGTACGATCCGGCAGGCGTCGGTGGTaca                                         >  1:357662/1‑36 (MQ=255)
                                                                gCGGGTCTGTACGATCCGGCAGGCGTCGGTGGTac                                          >  1:2816471/1‑35 (MQ=255)
                                                                gCGGGTCTGTACGATCCGGCAGGCGTCGGTGGTACACACGTCATGTACGTGCTGCACCATGCTGACAAGc       >  1:547527/1‑70 (MQ=255)
                                                                gCGGGTCTGTACGATCCGGCAGGCGTCGGTGGTACACACGTCATGTACGTGCTGCACCATGCTGACAAGc       >  1:68255/1‑70 (MQ=255)
                                                                gCGGGTCTGTACGATCCGGCAGGCGTCGGTGGTACACACGTCATGTACGTGCTGCACCATGCTGACAAGc       >  1:1687797/1‑70 (MQ=255)
                                                                gCGGGTCTGTACGATCAGGCAGGCGTCGGTGGTaca                                         >  1:916027/1‑36 (MQ=255)
                                                                gCGGGCCTGTACGATCCGGCAGGCGTCGGTGGTACACACGTCATGTACGTGCTGCACCATGCTGACAAGc       >  1:2272288/1‑70 (MQ=255)
                                                                gCGGGCCTGTACGATCCGGCAGGCGTCGGTGGTACACACGTCATGTACGTGCTGCACCATGCTGACAAGc       >  1:2268237/1‑70 (MQ=255)
                                                                  gggTCTGTACGATCCGGCAGGCGTCGGTGGTACACACGTCATGTACGTGCTGCACCATGCTGACAAGCCaa    <  1:487339/71‑1 (MQ=255)
                                                                  gggTCTGTACGATCCGGCAGGCGTCGGTGGTACACACGTCATGTACGTGCTGCACCATGCTGACAAGCCa     <  1:2774485/70‑1 (MQ=255)
                                                                   ggTCTGTACGATCCGGCAGGCGTCGGTGGTACACACGTCATGt                               >  1:3418311/1‑43 (MQ=255)
                                                                   ggTCTGTACGATCCGGCAGGCGTCGGTGGTACACACGTCATGTACGTGCTGCACCATGCTGACAAGCCaaa   <  1:918192/71‑1 (MQ=255)
                                                                   ggTCTGTACGATCCGGCAGGCGTCGGTGGTACACACGTCATGTACGTGCTGCACCATGCTGACAAGCCaaa   <  1:865126/71‑1 (MQ=255)
                                                                   ggTCTGTACGATCCGGCAGGCGTCGGTGGTACACACGTCATGTACGTGCTGCACCATGCTGACAAGCCaaa   <  1:835361/71‑1 (MQ=255)
                                                                   ggTCTGTACGATCCGGCAGGCGTCGGTGGTACACACGTCATGTACGTGCTGCACCATGCTGACAAGCCaaa   <  1:1919486/71‑1 (MQ=255)
                                                                   ggTCTGTACGATCCGGCAGGCGTCGGTGGTACACACGTCATGTACGTGCTGCACCATGCTGACAAGCCaa    <  1:1142144/70‑1 (MQ=255)
                                                                    gTCTGTACGATCCGGCAGGGGTCGGTGGTACACACGTCATGTACGTGCTGCACCATGCTGACAAGCCAAAt  <  1:1888776/71‑1 (MQ=255)
                                                                    gTCTGTACGATCCGGCAGGCGTCGGTGGTACACACGTCATGTACGTGCTGCACCATGCTGACAAGCCAAAt  <  1:2718619/71‑1 (MQ=255)
                                                                    gTCTGTACGATCCGGCAGGCGTCGGTGGTACACACGTCATGTACGTGCTGCACCATGCTGACAAGCCAAAt  <  1:3288822/71‑1 (MQ=255)
                                                                     |                                                                     
AGAGTCGATGAAAACGCTGGCGAGCGAGCGCGTTGCTGAGCTGAAAACCCGCGGTTACGACAATGCGGGTCTGTACGATCCGGCAGGCGTCGGTGGTACACACGTCATGTACGTGCTGCACCATGCTGACAAGCCAAAT  >  Exported/2336181‑2336319

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: