Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,919,561 A→G 46.7% G45G (GGT→GGC argI ← ornithine carbamoyltransferase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,919,5610AG46.7% 2.3 / 50.2 60G45G (GGT→GGCargIornithine carbamoyltransferase 1
Reads supporting (aligned to +/- strand):  ref base A (19/13);  new base G (17/11);  total (36/24)
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.90e-01

TGCGGCAACTTCGAAAGAGCATCGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCTTCGCGGCTAACTGC  >  minE/2919495‑2919624
                                                                  |                                                               
tGCGGCAACTTCGAAAGAGCATCGGGTACGAGTCGAGTCTTTTTCGAAGCTGAGCGCGATGTTTTTAcc                                                               >  1:319881/1‑69 (MQ=255)
tGCGGCAACTTCGAAAGAGCATCGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTAcc                                                               >  1:1033731/1‑69 (MQ=255)
tGCGGCAACTTCGAAAGAGCATCGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTAcc                                                               >  1:723839/1‑69 (MQ=255)
tGCGGCAACTTCGAAAGAGCATCGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTAcc                                                               >  1:2490669/1‑69 (MQ=255)
tGCGGCAACTTCGAAAGAGCATCGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTAcc                                                               >  1:1991696/1‑69 (MQ=255)
tGCGGCAACTTCGAAAGAGCATCGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTAcc                                                               >  1:1959147/1‑69 (MQ=255)
tGCGGCAACTTCGAAAGAGCATCGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTAcc                                                               >  1:3065189/1‑69 (MQ=255)
tGCGGCAACTTCGAAAGAGCATCGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTAcc                                                               >  1:1655505/1‑69 (MQ=255)
tGCGGCAACTTCGAAAGAGCATCGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTAc                                                                <  1:1116561/68‑1 (MQ=255)
tGCGGCAACTTCGAAAGAGCATCGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTAc                                                                >  1:1190938/1‑68 (MQ=255)
tGCGGCAACTTCGAAAGAGCATCGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTAc                                                                >  1:496386/1‑68 (MQ=255)
tGCGGCAACTTCGAAAGAGCATCGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTAc                                                                <  1:1841693/68‑1 (MQ=255)
           cGAAAGAGCATCGGGTACGGGTCGGGTCTTTTTCGAAGATGAGCGCGATGTTTTTGCCAGTGAGTCTgg                                                    >  1:2183428/1‑69 (MQ=255)
           cGAAAGAGCATCGGGTACGGGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTGCTAGTGAGTCTgg                                                    >  1:1053918/1‑69 (MQ=255)
           cGAAAGAGCATCGGGTACGGGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTGCCAGTGAGTCTgg                                                    >  1:1666106/1‑69 (MQ=255)
           cGAAAGAGCATCGGGTACGGGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTGCCAGTGAGTCTgg                                                    >  1:3086861/1‑69 (MQ=255)
           cGAAAGAGCATCGGGTACGGGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTGCCAGTGAGTCTgg                                                    >  1:1890075/1‑69 (MQ=255)
           cGAAAGAGCATCGGGTACGGGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTGCCAGTGAGTCTgg                                                    >  1:2729582/1‑69 (MQ=255)
           cGAAAGAGCATCGGGTACGGGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTGCCAGTGAGTCTgg                                                    >  1:1296098/1‑69 (MQ=255)
           cGAAAGAGCATCGGGTACGGGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTGCCAGTGAGTCTgg                                                    >  1:602478/1‑69 (MQ=255)
           cGAAAGAGCATCGGGTACGGGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTGCCAGTGAGTCTgg                                                    >  1:715691/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:2271943/1‑69 (MQ=255)
                           aCGGGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTGCCAGTGAGTCTGGCTTCTTCTTTACCGCt                                    >  1:403703/1‑69 (MQ=255)
                           aCGGGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTGCCAGTGAGTCTGGCTTCTTCTTTACCGCt                                    >  1:3374042/1‑69 (MQ=255)
                           aCGGGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTGCCAGTGAGTCTGGCTTCTTCTTTACCGCt                                    >  1:650730/1‑69 (MQ=255)
                           aCGGGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTGCCAGTGAGTCTGGCTTCTTCTTTACCGCt                                    >  1:1213226/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:194625/1‑68 (MQ=255)
                                         tttCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCt                       >  1:2569552/1‑68 (MQ=255)
                                         tttCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCt                       >  1:2510445/1‑68 (MQ=255)
                                         tttCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCt                       >  1:1166408/1‑68 (MQ=255)
                                         tttCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCt                       >  1:437937/1‑68 (MQ=255)
                                         tttCGAAGATGAGCGCGATGTTTTTACCAGTCAGTTTGGCTTCTTCTTTACCGCTTTTCTAATCTGc                        >  1:1449409/1‑67 (MQ=255)
                                           tCGAAGATGAGCGCGATGTTTTTGCCAGTGAGTCTGGCTTCTTCTTTACCGCTTTTCTTATCGGCtt                      <  1:783885/67‑1 (MQ=255)
                                           tCGAAGATGAGCGCGATGTTTTTGCCAGTGAGTCTGGCTTCTTCTTTACCGCTTTTCTTATCGGCtt                      <  1:661562/67‑1 (MQ=255)
                                           tCGAAGATGAGCGCGATGTTTTTGCCAGTGAGTCTGGCTTCTTCTTTACCGCTTTTCTTATCGGCtt                      <  1:1277570/67‑1 (MQ=255)
                                           tCGAAGATGAGCGCGATGTTTTTGCCAGTGAGTCTGGCTTCTTCTTTACCGCTTTTCTTATCGGCtt                      <  1:2482426/67‑1 (MQ=255)
                                           tCGAAGATGAGCGCGATGTTTTTGCCAGTGAGTCTGGCTTCTTCTTTACCGCTTTTCTTATCGGCtt                      <  1:499666/67‑1 (MQ=255)
                                           tCGAAGATGAGCGCGATGTTTTTGCCAGTGAGTCTGGCTTCTTCTTTACCGCTTTTCTTATCGGCtt                      <  1:1375494/67‑1 (MQ=255)
                                           tCGAAGATGAGCGCGATGTTTTTGCCAGTGAGTCTGGCTTCTTCTTTACCGCTTTTCTTATCGGCtt                      <  1:177684/67‑1 (MQ=255)
                                           tCGAAGATGAGCGCGATGTTTTTGCCAGTGAGTCTGGCTTCTTCTTTACCGCTTTTCTTATCGGCtt                      <  1:193973/67‑1 (MQ=255)
                                           tCGAAGATGAGCGCGATGTTTTTGCCAGTGAGTCTGGCTTCTTCTTTACCGCTTTTCTTATCGGCtt                      <  1:2717360/67‑1 (MQ=255)
                                           tCGAAGATGAGCGCGATGTTTTTGCCAGTGAGTCTGGCTTCTTCTTTACCGCTTTTCTTATCGGCtt                      <  1:2653980/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)
                                                                  |                                                               
TGCGGCAACTTCGAAAGAGCATCGGGTACGAGTCGAGTCTTTTTCGAAGATGAGCGCGATGTTTTTACCAGTGAGTTTGGCTTCTTCTTTACCGCTTTTCTTATCGGCTTTCAGCTTCGCGGCTAACTGC  >  minE/2919495‑2919624

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: