Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,244,132 A→C 100% A180A (GCT→GCG hemE ← uroporphyrinogen decarboxylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,244,1320AC80.6% 94.8 / 11.6 36A180A (GCT→GCGhemEuroporphyrinogen decarboxylase
Reads supporting (aligned to +/- strand):  ref base A (6/1);  new base C (0/29);  total (6/30)
Fisher's exact test for biased strand distribution p-value = 3.59e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.48e-01
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TGAGCACCGGCTTTAATCTGCGCATTCAGATACAAAGTGACGCTTTTCGCCAGTTTATCGAGCAGAGCGTGCAGCGCCTGCGGATCGGCATACATCATTTTTTTGATCACGGTGAAG  >  minE/2244067‑2244183
                                                                 |                                                   
tGAGCACCGGCTTTAATCTGCGCATTCAGATACAAAGTGACGCTTTTCGCCAGTTTATCGAGCagag                                                    <  1:1088013/67‑1 (MQ=255)
             tAATCTGCGCATTCAGATACAACGTGACGCGTTTCGCCAGTTTATCGAGTAGCGCGTGCa                                              <  1:805448/60‑1 (MQ=255)
             tAATCTGCGCATTCAGATACAAAGTGACGCTTTTCGCCAGTTTATCGAGTAGCGCGTGCa                                              <  1:2065273/60‑1 (MQ=255)
             tAATCTGCGCATTCAGATACAAAGTGACGCTTTTCGCCAGTTTATCGAGTAGCGCGTGCa                                              <  1:6425/60‑1 (MQ=255)
             tAATCTGCGCATTCAGATACAAAGTGACGCTTTTCGCCAGTTTATCGAGTAGCGCGTGCa                                              <  1:470870/60‑1 (MQ=255)
             tAATCTGCGCATTCAGATACAAAGTGACGCTTTTCGCCAGTTTATCGAGTAGCGCGTGCa                                              <  1:265784/60‑1 (MQ=255)
             tAATCTGCGCATTCAGATACAAAGTGACGCTTTTCGCCAGTTTATCGAGTAGCGCGTGCa                                              <  1:250742/60‑1 (MQ=255)
             tAATCTGCGCATTCAGATACAAAGTGACGCTTTTCGCCAGTTTATCGAGTAGCGCGTGCa                                              <  1:645838/60‑1 (MQ=255)
             tAATCTGCGCATTCAGATACAAAGTGACGCTTTTCGCCAGTTTATCGAGTAGCGCGTGCa                                              <  1:2164805/60‑1 (MQ=255)
             tAATCTGCGCATTCAGATACAAAGTGACGCTTTTCGCCAGTTTATCGAGTAGCGCGTGCa                                              <  1:2160741/60‑1 (MQ=255)
             tAATCTGCGCATTCAGATACAAAGTGACGCTTTTCGCCAGTTTATCGAGTAGCGCGTGCa                                              <  1:2130334/60‑1 (MQ=255)
             tAATCTGCGCATTCAGATACAAAGTGACGCTTTTCGCCAGTTTATCGAGTAGCGCGTGCa                                              <  1:2107875/60‑1 (MQ=255)
             tAATCTGCGCATTCAGATACAAAGTGACGCTTTTCGCCAGTTTATCGAGTAGCGCGTGCa                                              <  1:2101727/60‑1 (MQ=255)
             tAATCTGCGCATTCAGATACAAAGTGACGCTTTTCGCCAGTTTATCGAGTAGCGCGTGCa                                              <  1:2092161/60‑1 (MQ=255)
             tAATCTGCGCATTCAGATACAAAGTGACGCTTTTCGCCAGTTTATCGAGTAGCGCGTGCa                                              <  1:2073282/60‑1 (MQ=255)
             tAATCTGCGCATTCAGATACAAAGTGACGCTTTTCGCCAGTTTATCGAGTAGCGCGTGCa                                              <  1:2070598/60‑1 (MQ=255)
             tAATCTGCGCATTCAGATACAAAGTGACGCTTTTCGCCAGTTTATCGAGTAGCGCGTGCa                                              <  1:2043029/60‑1 (MQ=255)
             tAATCTGCGCATTCAGATACAAAGTGACGCTTTTCGCCAGTTTATCGAGTAGCGCGTGCa                                              <  1:1987702/60‑1 (MQ=255)
             tAATCTGCGCATTCAGATACAAAGTGACGCTTTTCGCCAGTTTATCGAGTAGCGCGTGCa                                              <  1:1866633/60‑1 (MQ=255)
             tAATCTGCGCATTCAGATACAAAGTGACGCTTTTCGCCAGTTTATCGAGTAGCGCGTGCa                                              <  1:1852038/60‑1 (MQ=255)
             tAATCTGCGCATTCAGATACAAAGTGACGCTTTTCGCCAGTTTATCGAGTAGCGCGTGCa                                              <  1:181459/60‑1 (MQ=255)
             tAATCTGCGCATTCAGATACAAAGTGACGCTTTTCGCCAGTTTATCGAGTAGCGCGTGCa                                              <  1:653433/60‑1 (MQ=255)
             tAATCTGCGCATTCAGATACAAAGTGACGCTTTTCGCCAGTTTATCGAGTAGCGCGTGCa                                              <  1:1571927/60‑1 (MQ=255)
             tAATCTGCGCATTCAGATACAAAGTGACGCTTTTCGCCAGTTTATCGAGTAGCGCGTGCa                                              <  1:707252/60‑1 (MQ=255)
             tAATCTGCGCATTCAGATACAAAGTGACGCTTTTCGCCAGTTTATCGAGTAGCGCGTGCa                                              <  1:1414869/60‑1 (MQ=255)
             tAATCTGCGCATTCAGATACAAAGTGACGCTTTTCGCCAGTTTATCGAGTAGCGCGTGCa                                              <  1:1001581/60‑1 (MQ=255)
             tAATCTGCGCATTCAGATACAAAGTGACGCTTTTCGCCAGTTTATCGAGTAGCGCGGGCa                                              <  1:1663221/60‑1 (MQ=255)
             tAATCTGCGCATTCAGATACAAAGCGACGCTTTTCGCCAGTTTATCGAGTAGCGCGTGCa                                              <  1:1645659/60‑1 (MQ=255)
              aaTCTGCGCATTCAGATACAAAGTGACGCTTTTCGCCAGTTTATCGAGTAGCGCGTGCa                                              <  1:360350/59‑1 (MQ=255)
                 cTGCGCATTCAGATACAAAGTGACGCTTTTCGCCAGTTTATCGAGTAGCGCGTGCa                                              <  1:1072331/56‑1 (MQ=255)
                            ggtACAAAGTGACGCTTTTCGCCAGTTTATCGAGCAGAGCGTGCAGCGCCTGCGGATCgg                               >  1:52537/3‑60 (MQ=255)
                            gATACAAAGTGACGCTTTTCGCCAGTTTATCGAGCAGAGCGTGCAGCGCCTGCGGATCgg                               >  1:2215497/1‑60 (MQ=255)
                            gATACAAAGTGACGCTTTTCGCCAGTTTATCGAGCAGAGCGTGCAGCGCCTGCGGATCgg                               >  1:145673/1‑60 (MQ=255)
                            gATACAAAGTGACGCTTTTCGCCAGTTTATCGAGCAGAGCGTGCAGCGCCTGCGGATCgg                               >  1:1422637/1‑60 (MQ=255)
                            gATACAAAGTGACGCTTTTCGCCAGTTTATCGAGCAGAGCGTGCAGCGCCTGCGGATCgg                               >  1:903754/1‑60 (MQ=255)
                            gATACAAAGTGACGCTTTTCGCCAGTTTATCGAGCAGAGCGTGAAGCGCCTGCGGATCgg                               >  1:1681633/1‑60 (MQ=255)
                                                  cAGTTTATCGAGCAGAGCGTGCAGCGCCTGCGGATCGGCATACATCATTTTTTTGATCACGGTGAAg  <  1:624526/67‑1 (MQ=255)
                                                                 |                                                   
TGAGCACCGGCTTTAATCTGCGCATTCAGATACAAAGTGACGCTTTTCGCCAGTTTATCGAGCAGAGCGTGCAGCGCCTGCGGATCGGCATACATCATTTTTTTGATCACGGTGAAG  >  minE/2244067‑2244183

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