Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,884,356 G→T 15.3% A84A (GCC→GCA cpdB ← 2':3'‑cyclic‑nucleotide 2'‑phosphodiesterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,884,3560GT15.3% 168.8 / 16.6 78A84A (GCC→GCAcpdB2':3'‑cyclic‑nucleotide 2'‑phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base G (43/23);  new base T (10/2);  total (53/25)
Fisher's exact test for biased strand distribution p-value = 3.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.27e-01

TAATGCCTTATAGACCGGGTGAATATCACCTGCTTTTAATCCTTTCGCCGACATGTAATCGGCCAGCGGACTCCCCTGAATCAAATCGCCGTTATCAACCAGTACGCTGTTTTTCACTTCATTGCGGG  >  minE/2884296‑2884423
                                                            |                                                                   
tAATGCCTTATAGACCGGGTGAATATCACCTGCTTTTAATCCTTTCGCCGACATGTAATCGGCCAGCg                                                              <  1:2685373/68‑1 (MQ=255)
    gCCTTATAGACCGGGTGAATATCACCTGCTTTTAATCCTTTCGCCGACATGTAATCGGCCAGCGGACTc                                                         <  1:2112700/69‑1 (MQ=255)
      cTTATAGACCGGGTGAATATCACCTGCTTTTAATCCTTTCGCCGACATGTAATCGGCCAGCGGACTccc                                                       >  1:1176027/1‑69 (MQ=255)
      cTTATAGACCGGGTGAATATCACCTGCTTTTAATCCTTTCGCCGACATGTAATCGGCCAGCGGACTccc                                                       >  1:2243807/1‑69 (MQ=255)
          taGACCGGGTGAATATCACCTGCTTTTAATCCTTTCGCCGACATGTAATCGGCCAGCGGACTCCCCTGa                                                   >  1:847711/1‑69 (MQ=255)
          taGACCGGGTGAATATCACCTGCTTTTAATCCTTTCGCCGACATGTAATCGGCCAGCGGACTCCCCTGa                                                   >  1:1384846/1‑69 (MQ=255)
          taGACCGGGTGAATATCACCTGCTTTTAATCCTTTCGCCGACATGTAATCGGCCAGCGGACTCCCCTGa                                                   >  1:841724/1‑69 (MQ=255)
                    gAATATCACCTGCTTTTAATCCTTTCGCCGACATGTAATCGGCCAGCGGACTCCCCTGAATCAAATCGc                                         >  1:2159266/1‑69 (MQ=255)
                    gAATATCACCTGCTTTTAATCCTTTCGCCGACATGTAATCGGCCAGCGGACTCCCCTGAATCAAATCGc                                         >  1:2301570/1‑69 (MQ=255)
                    gAATATCACCTGCTTTTAATCCTTTCGCCGACATGTAATCGGCCAGCGGACTCCCCTGAATCAAATCGc                                         >  1:3024391/1‑69 (MQ=255)
                    gAATATCACCTGCTTTTAATCCTTTCGCCGACATGTAATCGGCCAGCGGACTCCCCTGAATCAAATCGc                                         >  1:3042548/1‑69 (MQ=255)
                    gAATATCACCTGCTTTTAATCCTTTCGCCGACATGTAATCGGCCAGCGGACTCCCCTGAATCAAATCGc                                         >  1:3236401/1‑69 (MQ=255)
                    gAATATCACCTGCTTTTAATCCTTTCGCCGACATGTAATCGGCCAGCGGACTCCCCTGAATCAAATCGc                                         >  1:1643628/1‑69 (MQ=255)
                    gAATATCACCTGCTTTTAATCCTTTCGCCGACATGTAATCGGCCAGCGGACTCCCCTGAATCAAATCGc                                         >  1:1403544/1‑69 (MQ=255)
                       tatCACCTGCTTTTAATCCTTTCGCCGACATGTAATCGGCCAGCGGACTCCCCTGAATCAAATCGcc                                        >  1:355273/1‑67 (MQ=255)
                       tatCACCTGCTTTTAATCCTTTCGCCGACATGTAATCGGCCAGCGGACTCCCCTGAATCAAATCGCCGt                                      >  1:98400/1‑69 (MQ=255)
                       tatCACCTGCTTTTAATCCTTTCGCCGACATGTAATCGGCCAGCGGACTCCCCTGAATCAAATCGCCGt                                      >  1:2130345/1‑69 (MQ=255)
                         tCACCTGCTTTTAATCCTTTCGCCGACATGTAATCGGCCAGCGGACTcc                                                        <  1:2442246/49‑1 (MQ=255)
                         tCACCTGCTTTTAATCCTTTCGCCGACATGTAATCGGCCAGCGGACTCCCCTGAATCAAATCGCCGtt                                     >  1:2915674/1‑68 (MQ=255)
                         tCACCTGCTTTTAATCCTTTCGCCGACATGTAATCGGCCAGCGGACTCCCCTGAATCAAATCGCCGtt                                     >  1:61532/1‑68 (MQ=255)
                         tCACCTGCTTTTAATCCTTTCGCCGACATGTAATCGGCCAGCGGACTCCCCTGAATCAAATCGCCGtt                                     >  1:2334072/1‑68 (MQ=255)
                         tCACCTGCTTTTAATCCTTTCGCCGACATGTAATCGGCCAGCGGACTCCCCTGAATCAAATCGCCGtt                                     >  1:232790/1‑68 (MQ=255)
                         tCACCTGCTTTTAATCCTTTCGCCGACATGTAATCGGCCAGCGGACTCCCCTGAATCAAATCGCCGtt                                     >  1:1119224/1‑68 (MQ=255)
                         tCACCTGCTTTTAATCCTTTCGCCGACATGTAATCGGCCAGCGGACTCCCCTGAATCAAATCGCCGtt                                     >  1:987538/1‑68 (MQ=255)
                         tCACCTGCTTTTAATCCTTTCGCCGACATGTAATCGGCCAGCGGACTCCCCTGAATCAAATCGCCGtt                                     >  1:1354241/1‑68 (MQ=255)
                         tCACCTGCTTTTAATCCTTTCGCCGACATGTAATCGGCCAGCGGACTCCCCTGAATCAAATCGCCGtt                                     >  1:1964757/1‑68 (MQ=255)
                         tCACCTGCTTTTAATCCTTACGCCGACATGTAATCGGCCAGCGGACTCCCCTGAATCAAATCGCCGtt                                     >  1:2263049/1‑68 (MQ=255)
                           aCCTGCTTTTAATCCTTTCGCCGACATGTAATCTGCCAGCGGACTCCCCTGAATCAAATCGCCGTTATc                                  <  1:961699/69‑1 (MQ=255)
                           aCCTGCTTTTAATCCTTTCGCCGACATGTAATCTGCCAGCGGACTCCCCTGAATCAAATCGCCGTTATc                                  <  1:1501248/69‑1 (MQ=255)
                           aCCTGCTTTTAATCCTTTCGCCGACATGTAATCGGCCAGCGGACTCCCCTGAATCAAATCGCCGTTAt                                   >  1:2230597/1‑68 (MQ=255)
                           aCCTGCTTTTAATCCTTTCGCCGACATGTAATCGGCCAGCGGACTCCCCTGAATCAAATCGCCGTTATc                                  >  1:1287190/1‑69 (MQ=255)
                           aCCTGCTTTTAATCCTTTCGCCGACATGTAATCGGCCAGCGGACTCCCCTGAATCAAATCGCCGTTATc                                  >  1:1422257/1‑69 (MQ=255)
                           aCCTGCTTTTAATCCTTTCGCCGACATGTAATCGGCCAGCGGACTCCCCTGAATCAAATCGCCGTTATc                                  >  1:543909/1‑69 (MQ=255)
                           aCCTGCTTTTAATCCTTTCGCCGACATGTAATCGGCCAGCGGACTCCCCTGAATCAAATCGCCGTTATc                                  >  1:3151958/1‑69 (MQ=255)
                           aCCTGCTTTTAATCCTTTCGCCGACATGTAATCGGCCAGCGGACTCCCCTGAATCAAATCGCCGTTATc                                  >  1:2987382/1‑69 (MQ=255)
                           aCCTGCTTTTAATCCTTTCGCCGACATGTAATCGGCCAGCGGACTCCCCTGAATCAAATCGCCGTTATc                                  >  1:1957822/1‑69 (MQ=255)
                           aCCTGCTTTTAATCCTTTCGCCGACATGTAATCGGCCAGCGGACTCCCCTGAATCAAATCGCCGTTATc                                  >  1:2645874/1‑69 (MQ=255)
                           aCCTACTTTTAATCCTTTCGCCGACATGTAATCGGCCAGCGGACTCCCCTGAATCAAATCGCCGTTATc                                  >  1:544760/1‑69 (MQ=255)
                                  tttAATCCTTTCGCCGACATGTAATCGGCCAGCGGACTCCCCTGAATCAAATCGCCGTTATCAACCAg                            <  1:581265/68‑1 (MQ=255)
                                        ccTTTCGCCGACATGTAATCGGCCAGCGGACTCCCCTGAATCATATCGCCGTTATCAACCAGTACGCTg                     >  1:2975078/1‑69 (MQ=255)
                                        ccTTTCGCCGACATGTAATCGGCCAGCGGACTCCCCTGAATCAAATCGCCGTTATCAACCAGTACTCTg                     <  1:2022916/69‑1 (MQ=255)
                                        ccTTTCGCCGACATGTAATCGGCCAGCGGACTCCCCTGAATCAAATCGCCGTTATCAACCAGTACGCTg                     <  1:1659723/69‑1 (MQ=255)
                                        ccTTTCGCCGACATGTAATCGGCCAGCGGACTCCCCTGAATCAAATCGCCGTTATCAACCAGTACGCTg                     <  1:3125188/69‑1 (MQ=255)
                                        ccTTTCGCCGACATGTAATCGGCCAGCGGACTCCCCTGAATCAAATCGCCGTTATCAACCAGTACGCTg                     <  1:3217750/69‑1 (MQ=255)
                                        ccTTTCGCCGACATGTAATCGGCCAGCGGACTCCCCTGAATCAAATCGCCGTTATCAACCAGTACGCTg                     <  1:341824/69‑1 (MQ=255)
                                        ccTTTCGCCGACATGTAATCGGCCAGCGGACTCCCCTGAATCAAATCGCCGTTATCAACCAGTACGCTg                     <  1:2824107/69‑1 (MQ=255)
                                        ccTTTCGCCGACATGTAATCGGCCAGCGGACTCCCCTGAATCAAATCGCCGTTATCAACCAGTACGCTg                     >  1:494917/1‑69 (MQ=255)
                                        ccTTTCGCCGACATGTAATCGGCCAGCGGACTCCCCTGAATCAAATCGCCGTTATCAACCAGTACGCTg                     <  1:770607/69‑1 (MQ=255)
                                        ccTTTCGCCGACATGTAATCGGCCAGCGGACTCCCCTGAATCAAATCGCCGTTATCAACCAGTACGCTg                     <  1:2598102/69‑1 (MQ=255)
                                        ccTTTCGCCGACATGTAATCGGCCAGCGGACTCCCCTGAATCAAATCGCCGTTATCAACCAGTACGCTg                     <  1:2594164/69‑1 (MQ=255)
                                        ccTTTCGCCGACATGTAATCGGCCAGCGGACTCCCCTGAATCAAATCGCCGTTATCAACCAGTACGCTg                     <  1:1315335/69‑1 (MQ=255)
                                         cTTTCGCCGAGATGTAATCGGCCAGCTGACTCCCCTGAATCAAATCGCCGTTATCAACCAGTACGCTg                     <  1:1649224/68‑1 (MQ=255)
                                            tCGCCGACATGTAATCGGCCAGCGGACTCCCCTGAATCAAATCGCCGTTATCAACCAGTACGCTg                     <  1:2679470/65‑1 (MQ=255)
                                              gCCGACATGTAATCGGCCAGCGGACTCCCCTGAATCAAATCGCCGTTATCAACCAGTACGCTGttttt                <  1:761907/68‑1 (MQ=255)
                                              gCCGACATGTAATCGGCCAGCGGACTCCCCTGAATCAAATCGCCGTTATCAACCAGTACGCTGttttt                <  1:850291/68‑1 (MQ=255)
                                              gCCGACATGTAATCGGCCAGCGGACTCCCCTGAATCAAATCGCCGTTATCAACCAGTACGCTGttttt                <  1:2138751/68‑1 (MQ=255)
                                              gCCGACATGTAATCGGCCAGCGGACTCCCCTGAATCAAATCGCCGTTATCAACCAGTACGCTGttttt                <  1:2554882/68‑1 (MQ=255)
                                              gCCGACATGTAATCGGCCAGCGGACTCCCCTGAATCAAATCGCCGTTATCAACCAGTACGCTGttttt                <  1:2628010/68‑1 (MQ=255)
                                              gCCGACATGTAATCGGCCAGCGGACTCCCCTGAATCAAATCGCCGTTATCAACCAGTACGCTGttttt                <  1:195005/68‑1 (MQ=255)
                                              gCCGACATGTAATCGGCCAGCGGACTCCCCTGAATCAAATCGCCGTTATCAACCAGTACGCTGttttt                <  1:2661350/68‑1 (MQ=255)
                                                 gACATGTAATCTGCCAGCGGACTCCCCTGAATCAAATCGCCGTTATCGCCCAGTACGc                       >  1:663483/1‑58 (MQ=255)
                                                 gACATGTAATCTGCCAGCGGACTCCCCTGAATCAAATCGCCGTTATCGACCAGTACg                        >  1:1689073/1‑57 (MQ=255)
                                                 gACATGTAATCTGCCAGCGGACTCCCCTGAATCAAATCGCCGTTATCGACCAGTACGc                       >  1:1976656/1‑58 (MQ=255)
                                                 gACATGTAATCTGCCAGCGGACTCCCCTGAATCAAATCGCCGTTATCGACCAGTACGc                       >  1:2276994/1‑58 (MQ=255)
                                                 gACATGTAATCTGCCAGCGGACTCCCCTGAATCAAATCGCCGTTATCGACCAGTACGc                       >  1:2283024/1‑58 (MQ=255)
                                                 gACATGTAATCTGCCAGCGGACTCCCCTGAATCAAATCGCCGTTATCGACCAGTACGc                       >  1:2092194/1‑58 (MQ=255)
                                                 gACATGTAATCTGCCAGCGGACTCCCCTGAATCAAATCGCCGTTATCGACCAGTACGc                       >  1:2562450/1‑58 (MQ=255)
                                                 gACATGTAATCTGCCAGCGGACTCCCCTGAATCAAATCGCCGTTATCGACCAGTACGc                       >  1:379908/1‑58 (MQ=255)
                                                 gACATGTAATCTGCCAGCGGACTCCCCTGAATCAAATCGCCGTTATCGACCAGTACGc                       >  1:2395914/1‑58 (MQ=255)
                                                 gACATGTAATCTGCCAGCGGACTCCCCTGAATCAAATCGCCGTTATCGACCAGTACGc                       >  1:3172558/1‑58 (MQ=255)
                                                       tAATCGGCCAGCGGACTGCCCTGAATCAAATCGCCGTTATCAACCAGTACGCTGTTTTTCACTTCATTg      >  1:1567059/1‑69 (MQ=255)
                                                       tAATCGGCCAGCGGACTCCCCTGAATCAAATCGCCGTTATCAACCAGTACGCTGTTTTTCACTTCATTg      >  1:2317335/1‑69 (MQ=255)
                                                       tAATCGGCCAGCGGACTCCCCTGAATCAAATCGCCGTTATCAACCAGTACGCTGTTTTTCACTTCATTg      >  1:2699592/1‑69 (MQ=255)
                                                         aTCGGCCAGCGGACTCCCCTGAATCAAATCGCCGTTATCAACCAGTACGCTGTTTTTCACTTCATTGCg    >  1:965076/1‑69 (MQ=255)
                                                         aTCGGCCAGCGGACTCCCCTGAATCAAATCGCCGTTATCAACCAGTACGCTGTTTTTCACTTCATTGCg    >  1:1960814/1‑69 (MQ=255)
                                                         aTCGGCCAGCGGACTCCCCTGAATCAAATCGCCGTTATCAACCAGTACGCTGTTTTTCACTTCATTGCg    >  1:2406006/1‑69 (MQ=255)
                                                         aTCGGCCAGCGGACTCCCCTGAATCAAATCGCCGTTATCAACCAGTACGCTGTTTTTCACTTCATTGCg    >  1:2685159/1‑69 (MQ=255)
                                                         aTCGGCCAGCGGACTCCCCTGAATCAAATCGCCGTTATCAACCAGTACGCTGTTTTTCACTTCATTGCg    >  1:2404071/1‑69 (MQ=255)
                                                            ggCCAGCGGACTCCCCTGAATCAAATCGCCGTTATCAACCAGTACg                        >  1:2991739/1‑46 (MQ=255)
                                                            ggCCAGCGGACTCCCCTGAATCAAATCGCCGTTATCAACCAGTACGCTGTTTTTCACTTCATTGCggg  >  1:276829/1‑68 (MQ=255)
                                                            |                                                                   
TAATGCCTTATAGACCGGGTGAATATCACCTGCTTTTAATCCTTTCGCCGACATGTAATCGGCCAGCGGACTCCCCTGAATCAAATCGCCGTTATCAACCAGTACGCTGTTTTTCACTTCATTGCGGG  >  minE/2884296‑2884423

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: