Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,919,571 T→C 54.8% 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,5710TC54.8% 8.2 / 40.4 53K42R (AAA→AGA) argIornithine carbamoyltransferase 1
Reads supporting (aligned to +/- strand):  ref base T (9/15);  new base C (17/12);  total (26/27)
Fisher's exact test for biased strand distribution p-value = 1.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.83e-01

CGAAAGAGCATCGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCTTCGCGGCTAACTGCAGCAGGCTGTTG  >  minE/2919506‑2919636
                                                                 |                                                                 
cGAAAGAGCATCGGGTACGGGTCGGGTCTTTTTCGAAGATGAGCGCGATGTTTTTGCCAGTGAGTCTgg                                                                >  1:2183428/1‑69 (MQ=255)
cGAAAGAGCATCGGGTACGGGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTGCTAGTGAGTCTgg                                                                >  1:1053918/1‑69 (MQ=255)
cGAAAGAGCATCGGGTACGGGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTGCCAGTGAGTCTgg                                                                >  1:2729582/1‑69 (MQ=255)
cGAAAGAGCATCGGGTACGGGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTGCCAGTGAGTCTgg                                                                >  1:3086861/1‑69 (MQ=255)
cGAAAGAGCATCGGGTACGGGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTGCCAGTGAGTCTgg                                                                >  1:602478/1‑69 (MQ=255)
cGAAAGAGCATCGGGTACGGGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTGCCAGTGAGTCTgg                                                                >  1:1890075/1‑69 (MQ=255)
cGAAAGAGCATCGGGTACGGGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTGCCAGTGAGTCTgg                                                                >  1:715691/1‑69 (MQ=255)
cGAAAGAGCATCGGGTACGGGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTGCCAGTGAGTCTgg                                                                >  1:1666106/1‑69 (MQ=255)
cGAAAGAGCATCGGGTACGGGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTGCCAGTGAGTCTgg                                                                >  1:1296098/1‑69 (MQ=255)
  aaaGAGCATCGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGc                                                               <  1:1767657/68‑1 (MQ=255)
  aaaGAGCATCGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGc                                                               <  1:1557469/68‑1 (MQ=255)
   aaGAGCATCGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGc                                                               <  1:3402101/67‑1 (MQ=255)
         aTCGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAGtt                                                                   >  1:1743898/1‑57 (MQ=255)
         aTCGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGcttcttctt                                                       >  1:2514961/1‑69 (MQ=255)
         aTCGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGcttcttctt                                                       >  1:89145/1‑69 (MQ=255)
             ggTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGc                                                               <  1:761766/57‑1 (MQ=255)
                aCGGGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTGCCAGTGAGTCTGGCTTCTTCTTTACCGCt                                                >  1:650730/1‑69 (MQ=255)
                aCGGGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTGCCAGTGAGTCTGGCTTCTTCTTTACCGCt                                                >  1:403703/1‑69 (MQ=255)
                aCGGGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTGCCAGTGAGTCTGGCTTCTTCTTTACCGCt                                                >  1:1213226/1‑69 (MQ=255)
                aCGGGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTGCCAGTGAGTCTGGCTTCTTCTTTACCGCt                                                >  1:2271943/1‑69 (MQ=255)
                aCGGGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTGCCAGTGAGTCTGGCTTCTTCTTTACCGCt                                                >  1:3374042/1‑69 (MQ=255)
                aCGGGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTGCCAGTGAGTCTGGCTTCTTCTTTACCGCt                                                >  1:1200243/1‑69 (MQ=255)
                aCGGGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTGCCAGTGAGTCTGGCTTCTTCTTTACCGCt                                                >  1:2909766/1‑69 (MQ=255)
                aCGGGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTGCCAGTGAGTCTGGCTTCTTCTTTACCGCt                                                >  1:280330/1‑69 (MQ=255)
                 cgagtcgagtCTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCtt                                               <  1:1815701/69‑1 (MQ=255)
                      cgagtcTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCtt                                          <  1:2038474/69‑1 (MQ=255)
                       gagtcTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCtt                                          <  1:841454/68‑1 (MQ=255)
                              tttCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCt                                   >  1:1166408/1‑68 (MQ=255)
                              tttCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCt                                   >  1:437937/1‑68 (MQ=255)
                              tttCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCt                                   >  1:2569552/1‑68 (MQ=255)
                              tttCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCt                                   >  1:2510445/1‑68 (MQ=255)
                              tttCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCt                                   >  1:194625/1‑68 (MQ=255)
                              tttCGAAGATGAGCGCGATGTTTTTACCAGTCAGTTTGGCTTCTTCTTTACCGCTTTTCTAATCTGc                                    >  1:1449409/1‑67 (MQ=255)
                                tCGAAGATGAGCGCGATGTTTTTGCCAGTGAGTCTGGCTTCTTCTTTACCGCTTTTCTTATCGGCtt                                  <  1:1277570/67‑1 (MQ=255)
                                tCGAAGATGAGCGCGATGTTTTTGCCAGTGAGTCTGGCTTCTTCTTTACCGCTTTTCTTATCGGCtt                                  <  1:1375494/67‑1 (MQ=255)
                                tCGAAGATGAGCGCGATGTTTTTGCCAGTGAGTCTGGCTTCTTCTTTACCGCTTTTCTTATCGGCtt                                  <  1:2653980/67‑1 (MQ=255)
                                tCGAAGATGAGCGCGATGTTTTTGCCAGTGAGTCTGGCTTCTTCTTTACCGCTTTTCTTATCGGCtt                                  <  1:783885/67‑1 (MQ=255)
                                tCGAAGATGAGCGCGATGTTTTTGCCAGTGAGTCTGGCTTCTTCTTTACCGCTTTTCTTATCGGCtt                                  <  1:177684/67‑1 (MQ=255)
                                tCGAAGATGAGCGCGATGTTTTTGCCAGTGAGTCTGGCTTCTTCTTTACCGCTTTTCTTATCGGCtt                                  <  1:661562/67‑1 (MQ=255)
                                tCGAAGATGAGCGCGATGTTTTTGCCAGTGAGTCTGGCTTCTTCTTTACCGCTTTTCTTATCGGCtt                                  <  1:2482426/67‑1 (MQ=255)
                                tCGAAGATGAGCGCGATGTTTTTGCCAGTGAGTCTGGCTTCTTCTTTACCGCTTTTCTTATCGGCtt                                  <  1:2717360/67‑1 (MQ=255)
                                tCGAAGATGAGCGCGATGTTTTTGCCAGTGAGTCTGGCTTCTTCTTTACCGCTTTTCTTATCGGCtt                                  <  1:193973/67‑1 (MQ=255)
                                tCGAAGATGAGCGCGATGTTTTTGCCAGTGAGTCTGGCTTCTTCTTTACCGCTTTTCTTATCGGCtt                                  <  1:499666/67‑1 (MQ=255)
                                            gcgATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCTTCGCGGCt                    <  1:3290114/69‑1 (MQ=255)
                                                  tttttACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCTTCGCGGCTAACTgc              <  1:3413708/69‑1 (MQ=255)
                                                  tttttACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCTTCGCGGCTAACTgc              <  1:3228571/69‑1 (MQ=255)
                                                  tttttACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCACCTTCGCGGCTAACTgc              <  1:2750634/69‑1 (MQ=255)
                                                   ttttGCGAGTGAGTCTTGCTTCTTCTTTACCGGTTTTCTTATCGGCtt                                  <  1:2108595/48‑1 (MQ=25)
                                                              aGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCTTCGCGGCTAACTGCAGCAGGCTGTTg  <  1:534741/69‑1 (MQ=255)
                                                              aGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCTTCGCGGCTAACTGCAGCAGGCTGTTg  >  1:2121637/1‑69 (MQ=255)
                                                              aGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCTTCGCGGCTAACTGCAGCAGGCTGTTg  <  1:3070605/69‑1 (MQ=255)
                                                              aGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCTTCGCGGCTAACTGCAGCAGGCTGTTg  <  1:3040356/69‑1 (MQ=255)
                                                              aGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCTTCGCGGCTAACTGCAGCAGGCTGTTg  <  1:2676655/69‑1 (MQ=255)
                                                              aGTCTGGCTTCTTCTTTACCGCTTTTCTTATCGGCtt                                  <  1:3358756/37‑1 (MQ=255)
                                                                 |                                                                 
CGAAAGAGCATCGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCTTCGCGGCTAACTGCAGCAGGCTGTTG  >  minE/2919506‑2919636

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: