Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 982,677 C→G 15.6% A41A (GCC→GCG ydfA → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported982,6770CG15.6% 152.9 / 19.2 64A41A (GCC→GCGydfAhypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (20/34);  new base G (4/6);  total (24/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.65e-01

CAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGGAAAGAAACTCAGGTGAGTGATATGGATTTCGACACAATCATGGAAAAGGCTTACGAA  >  Exported/982607‑982747
                                                                      |                                                                      
cAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGcc                                                                        >  1:578358/1‑71 (MQ=255)
cAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGcc                                                                        >  1:4718072/1‑71 (MQ=255)
cAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGcc                                                                        >  1:5012760/1‑71 (MQ=255)
  aaCCAGGATGGAACGTTCACTGCCATGACGTATACCAAAAGCAAAACATTTAAAACCGAAACTGGTGcgcg                                                                      <  1:4637559/71‑1 (MQ=255)
  aaCCAGGATGGAACGTTCACTGCCATGACGTATACCAAAAGCAAAACATTTAAAACCGAAACTGGTGcgcg                                                                      <  1:5188924/71‑1 (MQ=255)
  aaCCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCg                                                                      >  1:4722547/1‑71 (MQ=255)
      aGGATGGAACGTTCACTGCCATGACGTATACCAAAAGCAAAACATTTAAAACCGAAACTGGTGCGCGCCga                                                                  <  1:3936469/71‑2 (MQ=255)
      aGGATGGAACGTTCACTGCCATGACGTATACCAAAAGCAAAACATTTAAAACCGAAACTGGTGCGCGCCga                                                                  <  1:3247907/71‑2 (MQ=255)
         aCGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCgctgg                                                               <  1:3627447/71‑1 (MQ=255)
             aaCGTTCACTGCCATGACGTATACCAAAAGCAAAACATTTAAAACCGAAACTGGTGCGCGCCGatggttgg                                                           <  1:3552788/70‑1 (MQ=255)
               cGTTCACTGCCATGACGTATACCAAAAGCAAAACATTTAAAACCGAAACTGGTGCGCGCCgatggttg                                                            >  1:5157372/1‑68 (MQ=255)
               cGTTCACTGCCATGACGTATACCAAAAGCAAAACATTTAAAACCGAAACTGGTGCGCGCCgatggttg                                                            >  1:3223200/1‑68 (MQ=255)
               cGTTCACTGCCATGACGTATACCAAAAGCAAAACATTTAAAACCGAAAATGGTGCGCGCCgatggttg                                                            >  1:3871693/1‑68 (MQ=255)
               cGTTCACCGCAATGACGTATACCAATAGCAAAACGTTTAAAACCGAAAATGGTGcc                                                                        <  1:8980/56‑1 (MQ=255)
                   cACTGCCATGACGTATACCAAAAGCAAAACATTTAAAACCGAAACTGGTGCGCGCCGATGGTTGgag                                                         >  1:5012133/1‑66 (MQ=255)
                   cACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGgaaaga                                                      <  1:4847846/70‑1 (MQ=255)
                   cACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGgaaaga                                                      <  1:2573547/70‑1 (MQ=255)
                   cACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGgaaaga                                                      <  1:4180364/70‑1 (MQ=255)
                      ctcAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGGAAAGAAACt                                                  >  1:2670830/3‑71 (MQ=255)
                          aTGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGgaaa                                                        <  1:5006571/61‑1 (MQ=255)
                          aTGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGgaaa                                                        <  1:1546652/61‑1 (MQ=255)
                           tGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGGAAAGAAACTCAGgt                                             <  1:3351675/71‑1 (MQ=255)
                           tGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGGAAAGAAACTCAGgt                                             <  1:2549657/71‑1 (MQ=255)
                                  taCCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCGTc                                                                    <  1:3813525/41‑1 (MQ=255)
                                  taCCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCGTc                                                                    <  1:5434580/41‑1 (MQ=255)
                                  taCCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCGTc                                                                    <  1:5663943/41‑1 (MQ=255)
                                  taCCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGGAAAGAAACTCAg                                               <  1:1387796/62‑1 (MQ=255)
                                    ccAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGGAAAGAAACTCAGgtgagtga                                       <  1:3495523/68‑1 (MQ=255)
                                    ccAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGGAAAGAAACTCAGGTGAGTGatat                                    >  1:2952438/1‑71 (MQ=255)
                                      aaaaGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGGAAAGAAACTCAGGTGAGTGATATg                                   >  1:1078776/1‑70 (MQ=255)
                                       aaaGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGGAAAGAAACTCAGGTGAGTGATATGGa                                 >  1:3235177/1‑71 (MQ=255)
                                        aaGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGgaa                                                         <  1:4209019/46‑1 (MQ=255)
                                        aaGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGgaa                                                         <  1:1216691/46‑1 (MQ=255)
                                        aaGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGGAAAGAAACTCAGGTGAGTGATATGGAt                                >  1:3058826/1‑71 (MQ=255)
                                        aaGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGGAAAGAAACTCAGGTGAGTGATATGGAt                                >  1:4786817/1‑71 (MQ=255)
                                          gCAAAACGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGGAAAGAAACTCAGgtgag                                          >  1:824947/1‑59 (MQ=255)
                                          gCAAAACGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGGAAAGAAACTCAGgtgag                                          >  1:2541855/1‑59 (MQ=255)
                                          gCAAAACGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGGAAAGAAACTCAGgtgag                                          >  1:4030641/1‑59 (MQ=255)
                                          gCAAAACATTTAAAACCGAAACTGGTGCGCGCCgatgg                                                               <  1:904566/38‑1 (MQ=16)
                                          gCAAAACATTTAAAACCGAAACTGGTGCGCGCCgatgg                                                               <  1:138582/38‑1 (MQ=16)
                                          gCAAAACATTTAAAACCGAAACTGGTGCGCGCCgatgg                                                               <  1:5535686/38‑1 (MQ=16)
                                           cAAAACGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGGAAAGAAACTCAGGTGAGTGATATGGATTTc                             >  1:3982701/1‑71 (MQ=255)
                                            aaaaCGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGGAAAGAAACTCAGGTGAGTGatat                                    <  1:1819954/63‑1 (MQ=255)
                                            aaaaCGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGGAAAGAAACTCAGGTGAGTGatat                                    <  1:1281384/63‑1 (MQ=255)
                                             aaaCGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGGAAAGAAACTc                                                 <  1:4594842/49‑1 (MQ=255)
                                             aaaCGTTTAAAACCGAAAATGGTGCCCGGCgctggctg                                                            <  1:5362432/38‑1 (MQ=255)
                                              aaCGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGgaaagaaa                                                    >  1:1557849/1‑45 (MQ=255)
                                              aaCGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGgaaaga                                                      <  1:4373710/43‑1 (MQ=255)
                                               aCGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGga                                                          >  1:848033/1‑38 (MQ=255)
                                                cGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGGAAAGAAACTCAGGTGAGTGATATGGATTTc                             <  1:2853523/66‑1 (MQ=255)
                                                 gTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGGAAAGAAACTCAGGTGAGTGATATGGATTTCGacac                        >  1:1360133/1‑70 (MQ=255)
                                                  tttAAAACCGAAAATGGTGCCCGTCGCTGTCTGga                                                          <  1:1823361/35‑1 (MQ=255)
                                                  tttAAAACCGAAAATGGTGCCCGTCGCTGGCTGga                                                          <  1:52607/35‑1 (MQ=255)
                                                  tttAAAACCGAAAATGGTGCCCGTCGCTGGCTGga                                                          <  1:3863878/35‑1 (MQ=255)
                                                  tttAAAACCGAAAATGGTGCCCGTCGCTGGCTGga                                                          <  1:4712045/35‑1 (MQ=255)
                                                  tttAAAACCGAAAATGGTGCCCGTCGCTGGCTGga                                                          <  1:462759/35‑1 (MQ=255)
                                                   ttAAAACCGAAAATGGTGCCCGTCGCTGGCTGGAAAGAAACTCAGGTGAGTGATATGGATTTCGACACAAt                     >  1:231596/1‑71 (MQ=255)
                                                    tAAAACCGAAAATGGTGCCCGTCGCTGGCTGgaaagaaa                                                    <  1:2010744/39‑1 (MQ=255)
                                                     aaaaCCGAAAATGGTGCCCGTCGCTGGCTGGAAAGAAACTCAgg                                              >  1:2804462/1‑44 (MQ=255)
                                                      aaaCCGAAAATGGTGCCCGTCGCTGGCTGGAAAGAAACTCAGGTGAGTGATATGGATTTCGACACAATCa                   <  1:2700682/70‑1 (MQ=255)
                                                        aCCGAAAATGGTGCCCGTCGCTGGCTGGAAAGAAACTCAgg                                              >  1:4197133/1‑41 (MQ=255)
                                                           gAAAATGGTGCCCGTCGCTGGCTGGAAAGAAACTc                                                 >  1:3943956/1‑35 (MQ=255)
                                                           gAAAATGGTGCCCGTCGCTGGCTGGAAAGAAACTCAg                                               >  1:3515714/1‑37 (MQ=255)
                                                           gAAAATGGTGCCCGTCGCTGGCTGGAAAGAAACTCAGGTGAGTGATATGGATTTCGACACAATCATGGaaa             <  1:1496753/71‑1 (MQ=255)
                                                           gAAAATGGTGCCCGTCGCTGGCTGGAAAGAAACTCAGGTGAGTGATATGGATTTCGACACAATCATGGa               <  1:4864487/69‑1 (MQ=255)
                                                            aaaaTGGTGCCCGTCGCTGGCTGGAAAGAAACTCAGgt                                             <  1:3779239/38‑1 (MQ=255)
                                                            aaaaTGGTGCCCGTCGCTGGCTGGAAAGAAACTCAGgt                                             <  1:3168484/38‑1 (MQ=255)
                                                              aaTGGTGCCCGTCGCTGGCTGGAAAGAAACTCAGGTGAGTGata                                     >  1:2644295/1‑44 (MQ=255)
                                                               aTGGTGCCCGTCGCTGGCTGGAAAGAAACTCAGgtgag                                          <  1:2108554/38‑1 (MQ=255)
                                                               aTGGTGCCCGTCGCTGGCTGGAAAGAAACTCAGgtgag                                          <  1:4265683/38‑1 (MQ=255)
                                                                     cccGTCGCTGGCTGGAAAGAAACTCAGgtgagtga                                       <  1:2217944/35‑1 (MQ=255)
                                                                      ccGTCGCTGGCTGGAAAGAAACTCAGGTGAGTGATATGGATTTCGACACAATCATGGAAAAGGCTTACgaa  <  1:2910793/71‑1 (MQ=255)
                                                                      |                                                                      
CAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGGAAAGAAACTCAGGTGAGTGATATGGATTTCGACACAATCATGGAAAAGGCTTACGAA  >  Exported/982607‑982747

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: