Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,179,295:1 +A 100% coding (745/984 nt) ansA → cytoplasmic L‑asparaginase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,179,2951.A90.9% 62.5 / 1.4 22coding (745/984 nt)ansAcytoplasmic L‑asparaginase I
Reads supporting (aligned to +/- strand):  ref base . (0/2);  new base A (20/0);  total (20/2)
Fisher's exact test for biased strand distribution p-value = 4.33e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.
Rejected as polymorphism: Polymorphic indel expands or contracts a homopolymer stretch.

ACCGGTGAAAGCATTGATTCTGCGCTCCTATGGCGTGGGTAATGCGCCACAAAACAAAGCCTTCCTGCAGGAATTACAAGAAGCCAGCGATCGCGGTATTGTGGTGGTCAACCTGACACA  >  minE/1179246‑1179365
                                                  |                                                                      
aCCGGTGAAAGCATTGATTCTGCGCTCCTATGGCGTGGGTAATGCGCCACAAAACAAAGCCTTCCTGCAgg                                                   <  1:2772335/71‑1 (MQ=255)
aCCGGTGAAAGCATTGATTCTGCGCTCCTATGGCGTGGGTAATGCGCCACAAAACAAAGCCTTCCTGCAgg                                                   <  1:1454989/71‑1 (MQ=255)
         aGCATTGATTCTGCGCTCCTATGGCGTGGGTAATGCGCCACAAAAACAAAGCCTTCCTGCAGGAATTa                                              >  1:822623/1‑68 (MQ=255)
         aGCATTGATTCTGCGCTCCTATGGCGTGGGTAATGCGCCACAAAAACAAAGCCTTCCTGCAGGAATTAc                                             >  1:1828583/1‑69 (MQ=255)
         aGCATTGATTCTGCGCTCCTATGGCGTGGGTAATGCGCCACAAAAACAAAGCCTTCCTGCAGGAATTAc                                             >  1:2742641/1‑69 (MQ=255)
         aGCATTGATTCTGCGCTCCTATGGCGTGGGTAATGCGCCACAAAAACAAAGCCTTCCTGCAGGAATTAc                                             >  1:2784482/1‑69 (MQ=255)
         aGCATTGATTCTGCGCTCCTATGGCGTGGGTAATGCGCCACAAAAACAAAGCCTTCCTGCAGGAATTAc                                             >  1:2701672/1‑69 (MQ=255)
         aGCATTGATTCTGCGCTCCTATGGCGTGGGTAATGCGCCACAAAAACAAAGCCTTCCTGCAGGAATTAc                                             >  1:2538925/1‑69 (MQ=255)
         aGCATTGATTCTGCGCTCCTATGGCGTGGGTAATGCGCCACAAAAACAAAGCCTTCCTGCAGGAATTAc                                             >  1:2505701/1‑69 (MQ=255)
         aGCATTGATTCTGCGCTCCTATGGCGTGGGTAATGCGCCACAAAAACAAAGCCTTCCTGCAGGAATTAc                                             >  1:2470188/1‑69 (MQ=255)
         aGCATTGATTCTGCGCTCCTATGGCGTGGGTAATGCGCCACAAAAACAAAGCCTTCCTGCAGGAATTAc                                             >  1:2861959/1‑69 (MQ=255)
         aGCATTGATTCTGCGCTCCTATGGCGTGGGTAATGCGCCACAAAAACAAAGCCTTCCTGCAGGAATTAc                                             >  1:2160235/1‑69 (MQ=255)
         aGCATTGATTCTGCGCTCCTATGGCGTGGGTAATGCGCCACAAAAACAAAGCCTTCCTGCAGGAATTAc                                             >  1:342612/1‑69 (MQ=255)
         aGCATTGATTCTGCGCTCCTATGGCGTGGGTAATGCGCCACAAAAACAAAGCCTTCCTGCAGGAATTAc                                             >  1:449452/1‑69 (MQ=255)
         aGCATTGATTCTGCGCTCCTATGGCGTGGGTAATGCGCCACAAAAACAAAGCCTTCCTGCAGGAATTAc                                             >  1:726973/1‑69 (MQ=255)
         aGCATTGATTCTGCGCTCCTATGGCGTGGGTAATGCGCCACAAAAACAAAGCCTTCCTGCAGGAATTAc                                             >  1:1553341/1‑69 (MQ=255)
         aGCATTGATTCTGCGCTCCTATGGCGTGGGTAATGCGCCACAAAAACAAAGCCTTCCTGCAGGAATTAc                                             >  1:1455970/1‑69 (MQ=255)
         aGCATTGATTCTGCGCTCCTATGGCGTGGGTAATGCGCCACAAAAACAAAGCCTTCCTGCAGGAATTAc                                             >  1:762730/1‑69 (MQ=255)
         aGCATTGATTCTGCGCTCCTATGGCGTGGGTAATGCGCCACAAAAACAAAGCCTTCCTGCAGGAATTAc                                             >  1:1423533/1‑69 (MQ=255)
         aGCATTGATTCTGCGCTCCTATGGCGTGGGTAATGCGCCACAAAAACAAAGCCTTCCTGCAGGAATTAc                                             >  1:1133866/1‑69 (MQ=255)
         aGCATTGATTCTGCGCTCCTATGGCGTGGGTAATGCGCCACAAAAACAAAGCCTTCCTGCAGGAATTAc                                             >  1:1098102/1‑69 (MQ=255)
         aGCATTGATTCTGCGCTCCTATGGCGTGGGTAATGCGCCACAAAAACAAAGCCTTCCTGCAGGAATTAc                                             >  1:1091828/1‑69 (MQ=255)
                                                 caaaacaaaGCCTTCCTGCAGGAATTACAAGAAGCCAGCGATCGCGGTATTGTGGTGGTCAACCTGAcaca  <  1:2478063/71‑1 (MQ=255)
                                                 caaaacaaaGCCTTCCTGCAGGAATTACAAGAAGCCAGCGATCGCGGTATTGTGGTGGTCAACCTGAcaca  <  1:697600/71‑1 (MQ=255)
                                                 caaaacaaaGCCTTCCTGCAGGAATTACAAGAAGCCAGCGATCGCGGTATTGTGGTGGTCAACCTGAcaca  <  1:334184/71‑1 (MQ=255)
                                                 caaaacaaaGCCTTCCTGCAGGAATTACAAGAAGCCAGCGATCGCGGTATTGTGGTGGTCAACCTGAcaca  <  1:747397/71‑1 (MQ=255)
                                                 caaaacaaaGCCTTCCTGCAGGAATTACAAGAAGCCAGCGATCGCGGTATTGTGGTGGTCAACCTGAcaca  <  1:2824351/71‑1 (MQ=255)
                                                 caaaacaaaGCCTTCCTGCAGGAATTACAAGAAGCCAGCGATCGCGGTATTGTGGTGGTCAACCTGAcaca  <  1:2807804/71‑1 (MQ=255)
                                                 caaaacaaaGCCTTCCTGCAGGAATTACAAGAAGCCAGCGATCGCGGTATTGTGGTGGTCAACCTGAcaca  <  1:931007/71‑1 (MQ=255)
                                                 caaaacaaaGCCTTCCTGCAGGAATTACAAGAAGCCAGCGATCGCGGTATTGTGGTGGTCAACCTGAcaca  <  1:2710856/71‑1 (MQ=255)
                                                 caaaacaaaGCCTTCCTGCAGGAATTACAAGAAGCCAGCGATCGCGGTATTGTGGTGGTCAACCTGAcaca  <  1:2708267/71‑1 (MQ=255)
                                                 caaaacaaaGCCTTCCTGCAGGAATTACAAGAAGCCAGCGATCGCGGTATTGTGGTGGTCAACCTGAcaca  <  1:2571922/71‑1 (MQ=255)
                                                 caaaacaaaGCCTTCCTGCAGGAATTACAAGAAGCCAGCGATCGCGGTATTGTGGTGGTCAACCTGAcaca  <  1:1032248/71‑1 (MQ=255)
                                                 caaaacaaaGCCTTCCTGCAGGAATTACAAGAAGCCAGCGATCGCGGTATTGTGGTGGTCAACCTGAcaca  <  1:2396355/71‑1 (MQ=255)
                                                 caaaacaaaGCCTTCCTGCAGGAATTACAAGAAGCCAGCGATCGCGGTATTGTGGTGGTCAACCTGAcaca  <  1:2342944/71‑1 (MQ=255)
                                                 caaaacaaaGCCTTCCTGCAGGAATTACAAGAAGCCAGCGATCGCGGTATTGTGGTGGTCAACCTGAcaca  <  1:1863488/71‑1 (MQ=255)
                                                 caaaacaaaGCCTTCCTGCAGGAATTACAAGAAGCCAGCGATCGCGGTATTGTGGTGGTCAACCTGAcaca  <  1:184866/71‑1 (MQ=255)
                                                 caaaacaaaGCCTTCCTGCAGGAATTACAAGAAGCCAGCGATCGCGGTATTGTGGTGGTCAACCTGAcaca  <  1:1705894/71‑1 (MQ=255)
                                                 caaaacaaaGCCTTCCTGCAGGAATTACAAGAAGCCAGCGATCGCGGTATTGTGGTGGTCAACCTGAcaca  <  1:167547/71‑1 (MQ=255)
                                                 caaaacaaaGCCTTCCTGCAGGAATTACAAGAAGCCAGCGATCGCGGTATTGTGGTGGTCAACCTGAcaca  <  1:1472034/71‑1 (MQ=255)
                                                 caaaacaaaGCCTTCCTGCAGGAATTACAAGAAGCCAGCGATCGCGGTATTGTGGTGGTCAACCTGAcaca  <  1:1425993/71‑1 (MQ=255)
                                                 caaaacaaaGCCTTCCTGCAGGAATTACAAGAAGCCAGCGATCGCGGTATTGTGGTGGTCAACCTGAcaca  <  1:1320825/71‑1 (MQ=255)
                                                 caaaacaaaGCCTTCCTGCAGGAATTACAAGAAGCCAGCGATCGCGGTATTGTGGTGGTCAACCTGAcaca  <  1:1079109/71‑1 (MQ=255)
                                                  |                                                                      
ACCGGTGAAAGCATTGATTCTGCGCTCCTATGGCGTGGGTAATGCGCCACAAAACAAAGCCTTCCTGCAGGAATTACAAGAAGCCAGCGATCGCGGTATTGTGGTGGTCAACCTGACACA  >  minE/1179246‑1179365

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: