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,0980GT96.0% 74.4 / NA 25L50L (CTC→CTAykfHuncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base T (9/15);  minor base A (1/0);  total (10/15)
Fisher's exact test for biased strand distribution p-value = 4.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCGAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCC  >  NC_000913/264034‑264156
                                                                |                                                          
gggTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAgg                                                         <  2:715774/68‑1 (MQ=255)
  gTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGAc                                                       >  1:402483/1‑68 (MQ=255)
  gTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGAc                                                       >  1:631241/1‑68 (MQ=255)
   tGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGGTCGGCGTATACCGCTAGGACt                                                      <  2:552603/68‑1 (MQ=255)
    gATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCAAGGCCTc                                                     >  2:743832/1‑68 (MQ=255)
       gCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCgg                                                  >  1:52966/1‑68 (MQ=255)
       gCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCgg                                                  >  2:277313/1‑68 (MQ=255)
             caggcaggTATAGGGCCCTTTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGatat                                            <  1:1029/68‑1 (MQ=255)
             caggcaggTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGatat                                            <  2:34372/68‑1 (MQ=255)
                   ggTATAGGGCCAATTGCGGTCCAGGCGGCGTACGGCGTATAACGCTAGGACTCCGGCGATAACCGGaa                                      >  1:804967/1‑68 (MQ=255)
                      ataGGGCCCATTGCGGTCCTGGCTGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGt                                   <  2:268380/68‑1 (MQ=37)
                       taGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGTTATCCGGAACGTc                                  <  2:441314/68‑1 (MQ=255)
                          ggCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCt                               <  2:519055/68‑1 (MQ=255)
                           gCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTg                              >  1:768012/1‑68 (MQ=255)
                           gCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTg                              >  1:434108/1‑68 (MQ=255)
                           gCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTg                              >  2:506576/1‑68 (MQ=255)
                                ttGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGgt                         <  2:118477/68‑1 (MQ=255)
                                 tGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGgtg                        >  2:360242/1‑68 (MQ=255)
                                 tGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGgtg                        <  2:458867/68‑1 (MQ=255)
                                 ggCGGGCCTGGCGGCGTTCGTCGTATACCGCTAGGACGCCGGCGATATCCGGAACGTCCCTGCCgggg                        <  1:354006/67‑3 (MQ=255)
                                          ggcggcGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAg               <  2:768013/68‑1 (MQ=255)
                                            cggcgTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGAc             <  1:431734/68‑1 (MQ=255)
                                               cgTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCt          <  2:39164/68‑1 (MQ=255)
                                                  tCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTAcc       <  1:716112/68‑1 (MQ=255)
                                                       gTATACCGCTAGGACGCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGcc  <  2:705877/68‑1 (MQ=255)
                                                                |                                                          
GGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCGAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCC  >  NC_000913/264034‑264156

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