Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 264,098 G→T L50L (CTC→CTA ykfH ← uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913264,0980GT95.2% 60.8 / NA 21L50L (CTC→CTAykfHuncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base T (15/5);  minor base A (1/0);  total (16/5)
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

GTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCGAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGG  >  NC_000913/264036‑264161
                                                              |                                                               
gTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGAc                                                            >  2:813028/1‑68 (MQ=255)
gTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGAc                                                            >  1:696886/1‑68 (MQ=255)
gTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGAc                                                            >  2:811202/1‑68 (MQ=255)
gTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGAc                                                            >  2:101073/1‑68 (MQ=255)
               caggTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCgg                                             <  1:379347/68‑1 (MQ=255)
               caggTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCgg                                             >  2:153038/1‑68 (MQ=255)
                    ataGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGt                                        >  2:844752/1‑68 (MQ=255)
                      aGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTcc                                      >  1:911061/1‑68 (MQ=255)
                      aGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTcc                                      <  2:696886/68‑1 (MQ=255)
                        ggCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCt                                    >  2:503829/1‑68 (MQ=255)
                        ggCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCt                                    >  2:355747/1‑68 (MQ=255)
                        ggCACATTGCGGTCCTGGCGGCGTTCGGCGTATACCGATAGGACTCCGGCGATATCCGGAAGGTccag                                    >  2:725564/1‑66 (MQ=21)
                         gaccATTGCGGTCCTGGCGGCGTTCGGCGTATACAGCAAGGACTCCGGCGATATCCGGAACGTCCCTg                                   >  2:201783/3‑68 (MQ=18)
                         gCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTg                                   >  1:461251/1‑68 (MQ=255)
                               tGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGgtg                             >  2:212871/1‑68 (MQ=255)
                                         gcggcgTTCGGCGTATTCCGCTAGGACTCCGGCGATATCCGGGAAGTACCTGCGGGTGTAAGGACAgg                   >  2:808978/1‑67 (MQ=18)
                                                   gCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTg         >  2:127800/1‑68 (MQ=255)
                                                      tataCCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCa      <  1:244742/68‑1 (MQ=255)
                                                      tataCCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCa      <  1:355747/68‑1 (MQ=255)
                                                       ataCCGCTAGGACGCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCAc     <  1:79960/68‑1 (MQ=255)
                                                          ccGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTgg  >  2:685218/1‑68 (MQ=255)
                                                              |                                                               
GTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCGAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGG  >  NC_000913/264036‑264161

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