Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 982,680 T→C 15.1% R42R (CGT→CGC ydfA → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported982,6800TC15.1% 152.5 / 18.1 67R42R (CGT→CGCydfAhypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (20/36);  new base C (4/6);  total (24/43)
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.54e-01

AGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGGAAAGAAACTCAGGTGAGTGATATGGATTTCGACACAATCATGGAAAAGGCTTACGAAGAA  >  Exported/982613‑982750
                                                                   |                                                                      
aGGATGGAACGTTCACTGCCATGACGTATACCAAAAGCAAAACATTTAAAACCGAAACTGGTGCGCGCCga                                                                     <  1:3247907/71‑2 (MQ=255)
aGGATGGAACGTTCACTGCCATGACGTATACCAAAAGCAAAACATTTAAAACCGAAACTGGTGCGCGCCga                                                                     <  1:3936469/71‑2 (MQ=255)
   aCGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCgctgg                                                                  <  1:3627447/71‑1 (MQ=255)
       aaCGTTCACTGCCATGACGTATACCAAAAGCAAAACATTTAAAACCGAAACTGGTGCGCGCCGatggttgg                                                              <  1:3552788/70‑1 (MQ=255)
         cGTTCACTGCCATGACGTATACCAAAAGCAAAACATTTAAAACCGAAACTGGTGCGCGCCgatggttg                                                               >  1:3223200/1‑68 (MQ=255)
         cGTTCACTGCCATGACGTATACCAAAAGCAAAACATTTAAAACCGAAACTGGTGCGCGCCgatggttg                                                               >  1:5157372/1‑68 (MQ=255)
         cGTTCACTGCCATGACGTATACCAAAAGCAAAACATTTAAAACCGAAAATGGTGCGCGCCgatggttg                                                               >  1:3871693/1‑68 (MQ=255)
             cACTGCCATGACGTATACCAAAAGCAAAACATTTAAAACCGAAACTGGTGCGCGCCGATGGTTGgag                                                            >  1:5012133/1‑66 (MQ=255)
             cACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGgaaaga                                                         <  1:2573547/70‑1 (MQ=255)
             cACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGgaaaga                                                         <  1:4180364/70‑1 (MQ=255)
             cACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGgaaaga                                                         <  1:4847846/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:5434580/41‑1 (MQ=255)
                            taCCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCGTc                                                                       <  1:3813525/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:4786817/1‑71 (MQ=255)
                                  aaGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGGAAAGAAACTCAGGTGAGTGATATGGAt                                   >  1:3058826/1‑71 (MQ=255)
                                    gCAAAACGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGGAAAGAAACTCAGgtgag                                             >  1:2541855/1‑59 (MQ=255)
                                    gCAAAACGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGGAAAGAAACTCAGgtgag                                             >  1:4030641/1‑59 (MQ=255)
                                    gCAAAACGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGGAAAGAAACTCAGgtgag                                             >  1:824947/1‑59 (MQ=255)
                                    gCAAAACATTTAAAACCGAAACTGGTGCGCGCCgatgg                                                                  <  1:5535686/38‑1 (MQ=16)
                                    gCAAAACATTTAAAACCGAAACTGGTGCGCGCCgatgg                                                                  <  1:904566/38‑1 (MQ=16)
                                    gCAAAACATTTAAAACCGAAACTGGTGCGCGCCgatgg                                                                  <  1:138582/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:3863878/35‑1 (MQ=255)
                                            tttAAAACCGAAAATGGTGCCCGTCGCTGGCTGga                                                             <  1:4712045/35‑1 (MQ=255)
                                            tttAAAACCGAAAATGGTGCCCGTCGCTGGCTGga                                                             <  1:462759/35‑1 (MQ=255)
                                            tttAAAACCGAAAATGGTGCCCGTCGCTGGCTGga                                                             <  1:52607/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:4265683/38‑1 (MQ=255)
                                                         aTGGTGCCCGTCGCTGGCTGGAAAGAAACTCAGgtgag                                             <  1:2108554/38‑1 (MQ=255)
                                                               cccGTCGCTGGCTGGAAAGAAACTCAGgtgagtga                                          <  1:2217944/35‑1 (MQ=255)
                                                                ccGTCGCTGGCTGGAAAGAAACTCAGGTGAGTGATATGGATTTCGACACAATCATGGAAAAGGCTTACgaa     <  1:2910793/71‑1 (MQ=255)
                                                                 cGTCGCTGGCTGGAAAGAAACTCAGGTGAGTGATATGGATTTCGACACAATCATGGAAAAGGCTTACgaag    <  1:4749366/71‑1 (MQ=255)
                                                                 cGTCGCTGGCTGGAAAGAAACTCAGGTGAGTGATATGGATTTCGACACAATCATGGAAAAGGCTTACg       >  1:3124421/1‑68 (MQ=255)
                                                                   tCGCTGTCTGGAAAGAAACTCAGGTGAGTGATATGGATTTCGACACAATCATGGa                  <  1:2156247/55‑1 (MQ=255)
                                                                   tCGCTGGCTGGAAAGAAACTCAGGTGAGTGATATGGAttt                                 >  1:1215278/1‑40 (MQ=255)
                                                                   tCGCTGGCTGGAAAGAAACTCAGGTGAGTGATATGGAtt                                  <  1:3664461/39‑1 (MQ=255)
                                                                   tCGCTGGCTGGAAAGAAACTCAGGTGAGTGATATGGATTTCGACACAATCATGGa                  <  1:4582549/55‑1 (MQ=255)
                                                                   tCGCTGGCTGGAAAGAAACTCAGGTGAGTGATATGGATTTCGACACAATCATGGAAAAGGCTTACgaagaa  <  1:2944908/71‑1 (MQ=255)
                                                                   tCGCTGGCTGGAAAGAAACTCAGGTGAGTGATATGGATTTCGACACAATCATGGAAAAGGCTTACgaagaa  >  1:1168268/1‑71 (MQ=255)
                                                                   |                                                                      
AGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGGAAAGAAACTCAGGTGAGTGATATGGATTTCGACACAATCATGGAAAAGGCTTACGAAGAA  >  Exported/982613‑982750

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