Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,135,755 C→T 14.3% E192E (GAG→GAA pnp ← polynucleotide phosphorylase/polyadenylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,135,7550CT14.3% 96.4 / 4.9 55E192E (GAG→GAApnppolynucleotide phosphorylase/polyadenylase
Reads supporting (aligned to +/- strand):  ref base C (24/23);  new base T (4/4);  total (28/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.72e-01

CAGATTCAACCATCAGTACAGCGGCTTCAGTACCGGCAACAACCAGATCCAGTTTGCTCTCTTTCAGCTCGTCCTGAGTCGGGTTCAGTACGTACTGGTCATTGATGTAACCTACGCGGGCAGCACCA  >  minE/2135688‑2135815
                                                                   |                                                            
cAGATTCAACCATCAGTACAGCGGCTTCAGTACCGGCAACAACCAGATCCAGTTTGCTCTCTTTCAGCt                                                             >  1:1056027/1‑69 (MQ=255)
cAGATTCAACCATCAGTACAGCGGCTTCAGTACCGGCAACAACCAGATCCAGTTTGCTCTCTTTCAGCt                                                             >  1:606833/1‑69 (MQ=255)
cAGATTCAACCATCAGTACAGCGGCTTCAGTACCGGCAACAACCAGATCCAGTTTGCTCTCTTTCAGCt                                                             >  1:2439493/1‑69 (MQ=255)
cAGATTCAACCATCAGTACAGCGGCTTCAGTACCGGCAACAACCAGATCCAGTTTGCTCTCTTTCAGCt                                                             >  1:3149275/1‑69 (MQ=255)
cAGATTCAACCATCAGTACAGCGGCTTCAGTACCGGCAACAACCAGATCCAGTTTGCTCTCTTTCAGCt                                                             >  1:3323665/1‑69 (MQ=255)
     tCAACCATCAGTACAGCGGCTTCAGTACCGGCAACAACCAGATCCAGTTTGCTCTCTTTCAGCTCGTc                                                         >  1:3280371/1‑68 (MQ=255)
     tCAACCATCAGTACAGCGGCTTCAGTACCGGCAACAACCAGATCCAGTTTGCTCTCTTTCAGCTCGTc                                                         >  1:2723163/1‑68 (MQ=255)
     tCAACCATCAGTACAGCGGCTTCAGTACCGGCAACAACCAGATCCAGTTTGCTCTCTTTCAGCTCGTc                                                         >  1:1025391/1‑68 (MQ=255)
     tCAACCATCAGTACAGCGGCTTCAGTACCGGCAACAACCAGATCCAGTTTGCTCTCTTTCAGCTCGTc                                                         >  1:2089714/1‑68 (MQ=255)
     tCAACCATCAGTACAGCGGCTTCAGTACCGGCAACAACCAGATCCAGTTTGCTCTCTTTCAGCTCGTc                                                         >  1:3085583/1‑68 (MQ=255)
     tCAACCATCAGTACAGCGGCTTCAGTACCGGCAACAACCAGATCCAGTTTGCTCTCTTTCAGCTCGTc                                                         >  1:1927721/1‑68 (MQ=255)
     tCAACCATCAGTACAGCGGCTTCAGTACCGGCAACAACCAGATCCAGTTTGCTCTCTTTCAGCTCGTc                                                         >  1:1683121/1‑68 (MQ=255)
     tCAACCATCAGTACAGCGGCTTCAGTACCGGCAACAACCAGATCCAGTTTGCTCTCTTTCAGCTCGTc                                                         >  1:1601358/1‑68 (MQ=255)
     tCAACCATCAGTACAGCGGCTTCAGTACCGGCAACAACCAGATCCAGTTTGCTCTCTTTCAGCTCGTc                                                         >  1:1461641/1‑68 (MQ=255)
     tCAACCATCAGTACAGCGGCTTCAGTACCGGCAACAACCAGATCCAGTTTGCTCTCTTTCAGCTCGTc                                                         >  1:156297/1‑68 (MQ=255)
     tCAACCATCAGTACAGCGGCTTCAGTACCGGCAACAACCAGATCCAGTTTGCTCTCTTTCAGCTCGTc                                                         >  1:1200139/1‑68 (MQ=255)
     tCAACCATCAGTACAGCGGCTTCAGTACCGGCAACAACCAGATCCAGTTTGCTCTCTTTCAGCTCGTc                                                         >  1:584514/1‑68 (MQ=255)
                tACAGCGGCTTCAGTACCGGCAACAACCAGATCCAGTTTGCTCTCTTTCAGCTCGTCCTGAGTCggg                                               >  1:1402851/1‑67 (MQ=255)
                tACAGCGGCTTCAGTACCGGCAACAACCAGATCCAGTTTGCTCTCTTTCAGCTCGTCCTGAGTCggg                                               >  1:2489553/1‑67 (MQ=255)
                tACAGCGGCTTCAGTACCGGCAACAACCAGATCCAGATTGCTCTCTTTCAGCTCGTCCTGAGTCgg                                                >  1:2985779/1‑66 (MQ=255)
                 aCAGCGGCTTCAGTACCGGCAACAACCAGATCCAGTTTGCTCTCTTTCAGCTCGTCCTGAGTCGGGTTc                                            <  1:2602134/69‑1 (MQ=255)
                 aCAGCGGCTTCAGTACCGGCAACAACCAGATCCAGTTTGCTCTCTTTCAGCTCGTCCTGAGTCGGGTTc                                            <  1:2510639/69‑1 (MQ=255)
                 aCAGCGGCTTCAGTACCGGCAACAACCAGATCCAGTTTGCTCTCTTTCAGCTCGTCCTGAGTCGGGTTc                                            <  1:2133808/69‑1 (MQ=255)
                 aCAGCGGCTTCAGTACCGGCAACAACCAGATCCAGTTTGCTCTCTTTCAGCTCGTCCTGAGTCGGGTTc                                            <  1:509243/69‑1 (MQ=255)
                 aCAGCGGCTTCAGTACCGGCAACAACCAGATCCAGTTTGCTCTCTTTCAGCTCGTCCTGAGTCGGGTTc                                            <  1:353479/69‑1 (MQ=255)
                 aCAGCGGCTTCAGTACCGGCAACAACCAGATCCAGTTTGCTCTCTTTCAGCTCGTCCTGAGTCGGGTTc                                            <  1:1380444/69‑1 (MQ=255)
                  cAGCGGCTTCAGTACCGGCAACAACCAGATCCAGTTTGCTCTCTTTCAGCTCGTCCTGAGTCGGGtt                                             >  1:609463/1‑67 (MQ=255)
                  cAGCGGCTTCAGTACCGGCAACAACCAGATCCAGTTTGCTCTCTTTCAGCTCGTCCTGAGTCGGGTTc                                            >  1:2253520/1‑68 (MQ=255)
                  cAGCGGCTTCAGTACCGGCAACAACCAGATCCAGTTTGCTCTCTTTCAGCTCGTCCTGAGTCGGGTTc                                            >  1:2996831/1‑68 (MQ=255)
                    gcagCTTCTGTACCGGCAACAACCAGGTCCAGTTTGCTCTCTTTCAGTTCGTCCTGAGTCGGGTTCAg                                          <  1:2853806/65‑1 (MQ=255)
                    gcagCTTCTGTACCGGCAACAACCAGGTCCAGTTTGCTCTCTTTCAGTTCGTCCTGAGTCGGGTTCAg                                          <  1:1610732/65‑1 (MQ=255)
                     cagCTTCTGTACCCGCAACAACCAGGTCCAGTTTGCTCTCTTTCAGTTCGTCCTGAGTCGGGTTCAg                                          <  1:1189672/65‑1 (MQ=255)
                      agCTTCTGTACCGGCAACAACCAGGTCCAGTTTGCTCTCTTTCAGTTCGTCCTGAGTCg                                                 >  1:247428/2‑59 (MQ=255)
                      agCTTCTGTACCGGCAACAACCAGGTCCAGTTTGCTCTCTTTCAGTTCGTCCTGAGTCGGGTTCAgta                                        >  1:2486931/2‑68 (MQ=255)
                      agCTTCTGTACCGGCAACAACCAGGTCCAGTTTGCTCTCTTTCAGTTCGTCCTGAGTCGGGTTCAgta                                        >  1:639637/2‑68 (MQ=255)
                      agCTTCTGTACCGGCAACAACCAGGTCCAGTTTGCTCTCTTTCAGTTCGTCCTGAGTCGGGTTCAgta                                        >  1:2452657/2‑68 (MQ=255)
                       gCTTCAGTACCGGCAACAACCAGATCCAGTTTGCTCTCTTTCAGCTCGTCCTGAGTCGGGTTCAgtacg                                      <  1:998097/69‑1 (MQ=255)
                       gCTTCAGTACCGGCAACAACCAGATCCAGTTTGCTCTCTTTCAGCTCGTCCTGAGTCGGGTTCAgtacg                                      <  1:1070261/69‑1 (MQ=255)
                        cTTCAGTACCGGCAACAACCAGATCCAGTTTGCTCTCTTTCAGCTCGTCCTGAGTCGGGTTCAgtacgt                                     <  1:1347906/69‑1 (MQ=255)
                        cTTCAGTACCGGCAACAACCAGATCCAGTTTGCTCTCTTTCAGCTCGTCCTGAGTCGGGTTCAgtacgt                                     <  1:978676/69‑1 (MQ=255)
                                       caacCAGATCCAGTTTGCTCTCTTTCAGCTCGTCCTGAGTCGGGTTCAgtacgt                                     <  1:535915/54‑1 (MQ=255)
                                        aacCAGATCCAGTTTGCTCTCTTTCAGCTCGTCCTGAGTCGGGTTCAGTACGTACTGGTCATTGATGTa                     <  1:3001081/69‑1 (MQ=255)
                                        aacCAGATCCAGTTTGCTCTCTTTCAGCTCGTCCTGAGTCGGGTTCAGTACGTACTGGTCATTGATGTa                     <  1:1655022/69‑1 (MQ=255)
                                        aacCAGATCCAGTTTGCTCTCTTTCAGCTCGTCCTGAGTCGGGTTCAGTACGTACTGGTCATTGATGTa                     <  1:1625837/69‑1 (MQ=255)
                                        aacCAGATCCAGTTTGCTCTCTTTCAGCTCGTCCTGAGTCGGGTTCAGTACGTACTGGTCATTGATGTa                     <  1:3381566/69‑1 (MQ=255)
                                        aacCAGATCCAGTTTGCTCTCTTTCAGCTCGTCCTGAGTCGGGTTCAGTACGTACTGGTCATTGATGTa                     <  1:1284892/69‑1 (MQ=255)
                                        aacCAGATCCAGTTTGCTCTCTTTCAGCTCGTCCTGAGTCGGGTTCAGTACGTACTGGTCATTGATGTa                     <  1:680677/69‑1 (MQ=255)
                                        aacCAGATCCAGTTTGCTCTCTTTCAGCTCGTCCTGAGTCGGGTTCAGTACGTACTGGTCATTGATGTa                     <  1:2905116/69‑1 (MQ=255)
                                             ggtCCAGTTTGCTCTCTTTCAGTTCGTCCTGAGTCGGGTTCAg                                          <  1:1483105/41‑1 (MQ=39)
                                                   gTTTGCTCTCTTTCAGCTCGTCCTGAGTCGGGTTCAGTACGTACTGGTCATTGATGTAACCTAcgcg            >  1:3256661/1‑67 (MQ=255)
                                                          ctctTTCAGCTCGTCCTGAGTCGGGTTCAGTACGTACTGGTCATTGATGTAACCTACGCGGGCAGCAcc   <  1:1970165/69‑1 (MQ=255)
                                                             tttCAGCTCGTCCTGAGTCGGGTTCAGTACGTACTGGTCATTGATGTAACCTACGCGGGCAGCACCa  <  1:2370092/67‑1 (MQ=255)
                                                             tttCAGCTCGTCCTGAGTCGGGTTCAGTACGTACTGGTCATTGATGTAACCTACGCGGGCAGCACCa  <  1:1090434/67‑1 (MQ=255)
                                                             tttCAGCTCGTCCTGAGTCGGGTTCAGTACGTACTGGTCATTGATGTAACCTACGCGGGCAGCACCa  <  1:124337/67‑1 (MQ=255)
                                                             tttCAGCTCGTCCTGAGTCGGGTTCAGTACGTACTGGTCATTGATGTAACCTACGCGGGCAGCACCa  <  1:2533018/67‑1 (MQ=255)
                                                                   |                                                            
CAGATTCAACCATCAGTACAGCGGCTTCAGTACCGGCAACAACCAGATCCAGTTTGCTCTCTTTCAGCTCGTCCTGAGTCGGGTTCAGTACGTACTGGTCATTGATGTAACCTACGCGGGCAGCACCA  >  minE/2135688‑2135815

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