Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,919,571 T→C 100% K42R (AAA→AGA)  argI ← ornithine carbamoyltransferase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,919,5710TC77.8% 118.8 / 29.5 63K42R (AAA→AGA) argIornithine carbamoyltransferase 1
Reads supporting (aligned to +/- strand):  ref base T (0/14);  new base C (49/0);  total (49/14)
Fisher's exact test for biased strand distribution p-value = 2.67e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

ATGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCTTCGCGGCTAACTGCAGCAGGC  >  minE/2919544‑2919631
                           |                                                            
aTGAGCGCGTTGTTTTTGCCAGTGAGTCTGGCTTCTTCttt                                                 >  1:360664/1‑41 (MQ=37)
aTGAGCGCGATGTTTTTGCCAGTGAGTCTGGcttcttctt                                                  >  1:1557901/1‑40 (MQ=38)
aTGAGCGCGATGTTTTTGCCAGTGAGTCTGGcttcttctt                                                  >  1:66640/1‑40 (MQ=38)
aTGAGCGCGATGTTTTTGCCAGTGAGTCTGGcttcttctt                                                  >  1:613931/1‑40 (MQ=38)
aTGAGCGCGATGTTTTTGCCAGTGAGTCTGGcttcttctt                                                  >  1:116548/1‑40 (MQ=38)
aTGAGCGCGATGTTTTTGCCAGTGAGTCTGGcttcttctt                                                  >  1:527450/1‑40 (MQ=38)
aTGAGCGCGATGTTTTTGCCAGTGAGTCTGGcttcttctt                                                  >  1:1212901/1‑40 (MQ=38)
aTGAGCGCGATGTTTTTGCCAGTGAGTCTGGcttcttctt                                                  >  1:53665/1‑40 (MQ=38)
aTGAGCGCGATGTTTTTGCCAGTGAGTCTGGCTTCTTCttt                                                 >  1:2210391/1‑41 (MQ=38)
aTGAGCGCGATGTTTTTGCCAGTGAGTCTGGCTTCTTCttt                                                 >  1:2296235/1‑41 (MQ=38)
aTGAGCGCGATGTTTTTGCCAGTGAGTCTGGCTTCTTCttt                                                 >  1:2258778/1‑41 (MQ=38)
aTGAGCGCGATGTTTTTGCCAGTGAGTCTGGCTTCTTCttt                                                 >  1:2256123/1‑41 (MQ=38)
aTGAGCGCGATGTTTTTGCCAGTGAGTCTGGCTTCTTCttt                                                 >  1:2232217/1‑41 (MQ=38)
aTGAGCGCGATGTTTTTGCCAGTGAGTCTGGCTTCTTCttt                                                 >  1:454496/1‑41 (MQ=38)
aTGAGCGCGATGTTTTTGCCAGTGAGTCTGGCTTCTTCttt                                                 >  1:2117717/1‑41 (MQ=38)
aTGAGCGCGATGTTTTTGCCAGTGAGTCTGGCTTCTTCttt                                                 >  1:2111080/1‑41 (MQ=38)
aTGAGCGCGATGTTTTTGCCAGTGAGTCTGGCTTCTTCttt                                                 >  1:2080597/1‑41 (MQ=38)
aTGAGCGCGATGTTTTTGCCAGTGAGTCTGGCTTCTTCttt                                                 >  1:1013570/1‑41 (MQ=38)
aTGAGCGCGATGTTTTTGCCAGTGAGTCTGGCTTCTTCttt                                                 >  1:202405/1‑41 (MQ=38)
aTGAGCGCGATGTTTTTGCCAGTGAGTCTGGCTTCTTCttt                                                 >  1:284280/1‑41 (MQ=38)
aTGAGCGCGATGTTTTTGCCAGTGAGTCTGGCTTCTTCttt                                                 >  1:388616/1‑41 (MQ=38)
aTGAGCGCGATGTTTTTGCCAGTGAGTCTGGCTTCTTCttt                                                 >  1:2025158/1‑41 (MQ=38)
aTGAGCGCGATGTTTTTGCCAGTGAGTCTGGCTTCTTCttt                                                 >  1:509899/1‑41 (MQ=38)
aTGAGCGCGATGTTTTTGCCAGTGAGTCTGGCTTCTTCttt                                                 >  1:53164/1‑41 (MQ=38)
aTGAGCGCGATGTTTTTGCCAGTGAGTCTGGCTTCTTCttt                                                 >  1:565806/1‑41 (MQ=38)
aTGAGCGCGATGTTTTTGCCAGTGAGTCTGGCTTCTTCttt                                                 >  1:607752/1‑41 (MQ=38)
aTGAGCGCGATGTTTTTGCCAGTGAGTCTGGCTTCTTCttt                                                 >  1:850098/1‑41 (MQ=38)
aTGAGCGCGATGTTTTTGCCAGTGAGTCTGGCTTCTTCttt                                                 >  1:962167/1‑41 (MQ=38)
aTGAGCGCGATGTTTTTGCCAGTGAGTCTGGCTTCTTCttt                                                 >  1:992119/1‑41 (MQ=38)
aTGAGCGCGATGTTTTTGCCAGTGAGTCTGGCTTCTTCttt                                                 >  1:992155/1‑41 (MQ=38)
aTGAGCGCGATGTTTTTGCCAGTGAGTCTGGCTTCTTCttt                                                 >  1:1128819/1‑41 (MQ=38)
aTGAGCGCGATGTTTTTGCCAGTGAGTCTGGCTTCTTCttt                                                 >  1:1473353/1‑41 (MQ=38)
aTGAGCGCGATGTTTTTGCCAGTGAGTCTGGCTTCTTCttt                                                 >  1:139021/1‑41 (MQ=38)
aTGAGCGCGATGTTTTTGCCAGTGAGTCTGGCTTCTTCttt                                                 >  1:152758/1‑41 (MQ=38)
aTGAGCGCGATGTTTTTGCCAGTGAGTCTGGCTTCTTCttt                                                 >  1:1344682/1‑41 (MQ=38)
aTGAGCGCGATGTTTTTGCCAGTGAGTCTGGCTTCTTCttt                                                 >  1:1240503/1‑41 (MQ=38)
aTGAGCGCGATGTTTTTGCCAGTGAGTCTGGCTTCTTCttt                                                 >  1:1173056/1‑41 (MQ=38)
aTGAGCGCGATGTTTTTGCCAGTGAGTCTGGCTTCTTCttt                                                 >  1:1722757/1‑41 (MQ=38)
aTGAGCGCGATGTTTTTGCCAGTGAGTCTGGCTTCTTCttt                                                 >  1:1146203/1‑41 (MQ=38)
aTGAGCGCGATGTTTTTGCCAGTGAGTCTGGCTTCTTCttt                                                 >  1:1805614/1‑41 (MQ=38)
aTGAGCGCGATGTTTTTGCCAGTGAGTCTGGCTTCTTCttt                                                 >  1:1822033/1‑41 (MQ=38)
aTGAGCGCGATGTTTTTGCCAGTGAGTCTGGCTTCTTCttt                                                 >  1:1037778/1‑41 (MQ=38)
aTGAGCGCGATGTTTTTGCCAGTGAGTCTGGCTTCTTCttt                                                 >  1:1827288/1‑41 (MQ=38)
aTGAGCGCGATGTTTTTGCCAGTGAGTCTGGCTTCTTCttt                                                 >  1:187952/1‑41 (MQ=38)
aTGAGCGCGATGTTTTTGCCAGTGAGTCTGGCTTCTTCttt                                                 >  1:19255/1‑41 (MQ=38)
aTGAGCGCGATGTTTTTGCCAGTGAGTCTGGCTTCTTCttt                                                 >  1:1967904/1‑41 (MQ=38)
aTGAGCGCGATGTTTTTGCCAGTGAGTCTGGCTTATTCttt                                                 >  1:2132281/1‑41 (MQ=37)
aTGAGCGCGATGTTTTTGCCAGTCAGTCTGGCTTCTTCttt                                                 >  1:377754/1‑41 (MQ=37)
aTGAGCGCGATGTTTTTGCCAGTCAGTCTGGCTTCTTCttt                                                 >  1:380389/1‑41 (MQ=37)
                    aGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCTTCGCGGCTAACTGCAGCAgg   <  1:1638785/67‑1 (MQ=255)
                    aGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCTTCGCGGCTAACTGCAGCAgg   <  1:104568/67‑1 (MQ=255)
                     gTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCTTCGCGGCTAACTGCAGCAGGc  <  1:1988430/67‑1 (MQ=255)
                     gTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCTTCGCGGCTAACTGCAGCAGGc  <  1:1313996/67‑1 (MQ=255)
                     gTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCTTCGCGGCTAACTGCAGCAGGc  <  1:2315487/67‑1 (MQ=255)
                     gTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCTTCGCGGCTAACTGCAGCAGGc  <  1:1223348/67‑1 (MQ=255)
                     gTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCTTCGCGGCTAACTGCAGCAGGc  <  1:1207609/67‑1 (MQ=255)
                     gTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCTTCGCGGCTAACTGCAGCAGGc  <  1:16787/67‑1 (MQ=255)
                     gTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCTTCGCGGCTAACTGCAGCAGGc  <  1:1745997/67‑1 (MQ=255)
                     gTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCTTCGCGGCTAACTGCAGCAGGc  <  1:2234761/67‑1 (MQ=255)
                     gTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCTTCGCGGCTAACTGCAGCAGGc  <  1:915302/67‑1 (MQ=255)
                     gTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCTTCGCGGCTAACTGCAGCAGGc  <  1:1038779/67‑1 (MQ=255)
                     gTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCTTCGCGGCTAACTGCAGCAGGc  <  1:2121788/67‑1 (MQ=255)
                     gTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCTTCGCGGCTAACTGCAGCAGGc  <  1:1027492/67‑1 (MQ=255)
                           |                                                            
ATGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCTTCGCGGCTAACTGCAGCAGGC  >  minE/2919544‑2919631

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