Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 774,341 C→T 68.2% intergenic (‑135/‑37) ymfC ← / → icd 23S rRNA pseudouridine synthase/isocitrate dehydrogenase, specific for NADP+

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE774,3410CT68.2% 36.1 / 21.3 47intergenic (‑135/‑37)ymfC/icd23S rRNA pseudouridine synthase/isocitrate dehydrogenase, specific for NADP+
Reads supporting (aligned to +/- strand):  ref base C (8/7);  new base T (17/15);  total (25/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAAAAATTGCTCTAAAGCATCCGTATCGCAGGACGCAAACGCATATGCAACGTGGTGGCAGACGAGCAAACCAGTAGCGCTCGAAGGAGAGGTGAATGGAAAGTAAAGTAGTTG  >  minE/774283‑774396
                                                          |                                                       
aaaaaaTTGCTCTAAAGCATCCGTATCGCAGGACGCAAACGCATATGCAACGTGGTGGTAGACGAGCaa                                               <  1:2084544/69‑1 (MQ=255)
aaaaaaTTGCTCTAAAGCATCCGTATCGCAGGACGCAAACGCATATGCAACGTGGTGGTAGACGAGCaa                                               <  1:2832753/69‑1 (MQ=255)
aaaaaaTTGCTCTAAAGCATCCGTATCGCAGGACGCAAACGCATATGCAACGTGGTGGTAGACGAGCaa                                               <  1:2957401/69‑1 (MQ=255)
aaaaaaTTGCTCTAAAGCATCCGTATCGCAGGACGCAAACGCATATGCAACGTGGTGGTAGACGAGCaa                                               <  1:1272238/69‑1 (MQ=255)
aaaaaaTTGCTCTAAAGCATCCGTATCGCAGGACGCAAACGCATATGCAACGTGGTGGTAGACGAGCaa                                               <  1:2196238/69‑1 (MQ=255)
aaaaaaTTGCTCTAAAGCATCCGTATCGCAGGACGCAAACGCATATGCAACGTGGTGGTAGACGAGCaa                                               <  1:2147241/69‑1 (MQ=255)
aaaaaaTTGCTCTAAAGCATCCGTATCGCAGGACGCAAACGCATATGCAACGTGGTGGTAGACGAGCaa                                               <  1:1885019/69‑1 (MQ=255)
aaaaaaTTGCTCTAAAGCATCCGTATCGCAGGACGCAAACGCATATGCAACGTGGTGGTAGACGAGCaa                                               <  1:3296927/69‑1 (MQ=255)
aaaaaaTTGCTCTAAAGCATCCGTATCGCAGGACGCAAACGCATATGCAACGTGGTGGTAGACGAGCaa                                               <  1:1502208/69‑1 (MQ=255)
aaaaaaTTGCTCTAAAGCATCCGTATCGCAGGACGCAAACGCATATGCAACGTGGTGGTAGACGAGCaa                                               <  1:978819/69‑1 (MQ=255)
aaaaaaTTGCTCTAAAGCATCCGTATCGCAGGACGCAAACGCATATGCAACGTGGTGGTAGACGAGCaa                                               <  1:903885/69‑1 (MQ=255)
aaaaaaTTGCTCTAAAGCATCCGTATCGCAGGACGCAAACGCATATGCAACGTGGTGGTAGACGAGCaa                                               <  1:3386129/69‑1 (MQ=255)
aaaaaaTTGCTCTAAAGCATCCGTATCGCAGGACGCAAACGCATATGCAACGTGGTGGTAGACGAGCaa                                               <  1:838381/69‑1 (MQ=255)
aaaaaaTTGCTCTAAAGCATCCGTATCGCAGGACGCAAACGCATATGCAACGTGGTGGTAGACGAGCaa                                               <  1:389100/69‑1 (MQ=255)
                  aTCCGTATCGCAGGACGCAAACGCATATGCAACGTGGTGGCAGACGAGCAAACCAGTAGTGCTCGAAgg                             <  1:2442944/69‑1 (MQ=255)
                      gTATCGCAGGACGCAAACGCATATGCAACGTGGTGGTAGACGAGCaa                                               <  1:920547/47‑1 (MQ=255)
                         tCGCAGGACGCAAACGCATATGCAACGTGGTGGCAGACGAGCAAACCAGTAGCGCTCGAAGGAGAGGTg                      <  1:3043225/69‑1 (MQ=255)
                         tCGCAGGACGCAAACGCATATGCAACGTGGTGGCAGACGAGCAAACCAGTAGCGCTCGAAGGAGAGGTg                      <  1:2777183/69‑1 (MQ=255)
                                 cGCAAACGCATATGCAACGTGGTGGCAGACGAGCAAACCAGTAGCGCTCGAAGGAGAGGTGAATGGaa               >  1:2057/1‑68 (MQ=255)
                                 cGCAAACGCATATGCAACGTGGTGGCAGACGAGCAAACCAGTAGCGCTCGAAGGAGAGGTGAATGGaa               >  1:342819/1‑68 (MQ=255)
                                    aaaCGCATATGCAACGTGGTGGCAGACGAGCAAACCAGTAGCGCTCGAAGGAGAGGTg                      <  1:2267949/58‑1 (MQ=255)
                                           tatGCAACGTGGTGGCAGACGAGCAAACCAGTAGCGCTCGAAGGAGAGGTGAATGGAAAGTAAagtagt    <  1:3064395/69‑1 (MQ=255)
                                           tatGCAACGTGGTGGCAGACGAGCAAACCAGTAGCGCTCGAAGGAGAGGTGAATGGAAAGTAAagtagt    <  1:249219/69‑1 (MQ=255)
                                            atGCAACGTGGTGGTAGACGAGCAAACCAGTAGCGCTCGAAGGAGAGGTGAATGGaaagta           >  1:2647771/1‑61 (MQ=255)
                                            atGCAACGTGGTGGTAGACGAGCAAACCAGTAGCGCTCGAAGGAGAGGTGAATGGaaagta           >  1:3009353/1‑61 (MQ=255)
                                            atGCAACGTGGTGGTAGACGAGCAAACCAGTAGCGCTCGAAGGAGAGGTGAATGGaaagta           >  1:3241251/1‑61 (MQ=255)
                                            atGCAACGTGGTGGTAGACGAGCAAACCAGTAGCGCTCGAAGGAGAGGTGAATGGaaagta           >  1:2998203/1‑61 (MQ=255)
                                            atGCAACGTGGTGGTAGACGAGCAAACCAGTAGCGCTCGAAGGAGAGGTGAATGGaaagta           >  1:2990388/1‑61 (MQ=255)
                                            atGCAACGTGGTGGTAGACGAGCAAACCAGTAGCGCTCGAAGGAGAGGTGAATGGaaagta           >  1:2264219/1‑61 (MQ=255)
                                            atGCAACGTGGTGGTAGACGAGCAAACCAGTAGCGCTCGAAGGAGAGGTGAATGGaaagta           >  1:2170340/1‑61 (MQ=255)
                                            atGCAACGTGGTGGTAGACGAGCAAACCAGTAGCGCTCGAAGGAGAGGTGAATGGaaagta           >  1:2120528/1‑61 (MQ=255)
                                            atGCAACGTGGTGGTAGACGAGCAAACCAGTAGCGCTCGAAGGAGAGGTGAATGGaaagta           >  1:1282939/1‑61 (MQ=255)
                                            atGCAACGTGGTGGTAGACGAGCAAACCAGTAGCGCTCGAAGGAGAGGTGAATGGaaagta           >  1:1306010/1‑61 (MQ=255)
                                            atGCAACGTGGTGGTAGACGAGCAAACCAGTAGCGCTCGAAGGAGAGGTGAATGGaaagta           >  1:1754180/1‑61 (MQ=255)
                                            atGCAACGTGGTGGTAGACGAGCAAACCAGTAGCGCTCGAAGGAGAGGTGAATGGaaagta           >  1:1739204/1‑61 (MQ=255)
                                            atGCAACGTGGTGGTAGACGAGCAAACCAGTAGCGCTCGAAGGAGAGGTGAATGGaaagta           >  1:1358412/1‑61 (MQ=255)
                                            atGCAACGTGGTGGTAGACGAGCAAACCAGTAGCGCTCGAAGGAGAGGTGAATGGaaagta           >  1:1706255/1‑61 (MQ=255)
                                            atGCAACGTGGTGGTAGACGAGCAAACCAGTAGCGCTCGAAGGAGAGGTGAATGGaaagta           >  1:1601073/1‑61 (MQ=255)
                                            atGCAACGTGGTGGTAGACGAGCAAACCAGTAGCGCTCGAAGGAGAGGTGAATGGaaagta           >  1:1576951/1‑61 (MQ=255)
                                            atGCAACGTGGTGGTAGACGAGCAAACCAGTAGCGCTCGAAGGAGAGGTGAATGGaaag             >  1:2594977/1‑59 (MQ=255)
                                            atGCAACGTGGTGGCAGACGAGCAAACCAGTAGCGCTCGAAGGAGAGGTGAATGGAAAGTAAagtagt    <  1:2794227/68‑1 (MQ=255)
                                             tGCAACGTGTTGGCAGACGAGCAAACCAGTAGCGCTCGAAGGAGAGGTGAATGGAAAGTAAAGTAgttg  >  1:2723372/1‑69 (MQ=255)
                                             tGCAACGTGGTGGCAGACGAGCAAACCAGTAGCGCTCGAAGGAGAGGTGAATGGAAAGTAAAGTAgttg  >  1:501949/1‑69 (MQ=255)
                                             tGCAACGTGGTGGCAGACGAGCAAACCAGTAGCGCTCGAAGGAGAGGTGAATGGAAAGTAAAGTAgttg  >  1:3418106/1‑69 (MQ=255)
                                             tGCAACGTGGTGGCAGACGAGCAAACCAGTAGCGCTCGAAGGAGAGGTGAATGGAAAGTAAAGTAgttg  >  1:2002188/1‑69 (MQ=255)
                                             tGCAACGTGGTGGCAGACGAGCAAACCAGTAGCGCTCGAAGGAGAGGTGAATGGAAAGTAAAGTAgttg  >  1:976929/1‑69 (MQ=255)
                                             tGCAACGTGGTGGCAGACGAGCAAACCAGTAGCGCTCGAAGGAGAGGTGAATGGAAAGTAAAGTAgttg  >  1:2072995/1‑69 (MQ=255)
                                                          cAGACGAGCAAACCAGTAGCGCTCGAAGGAGAGGTg                      <  1:641079/36‑1 (MQ=255)
                                                          |                                                       
AAAAAATTGCTCTAAAGCATCCGTATCGCAGGACGCAAACGCATATGCAACGTGGTGGCAGACGAGCAAACCAGTAGCGCTCGAAGGAGAGGTGAATGGAAAGTAAAGTAGTTG  >  minE/774283‑774396

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: