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.3% 82.5 / NA 27L50L (CTC→CTAykfHuncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base T (16/10);  minor base A (1/0);  total (17/10)
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:824497/1‑68 (MQ=255)
gTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGAc                                                            >  1:796880/1‑68 (MQ=255)
gTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGAGTTCGGCGTATACCGCTAGGAc                                                            >  2:971455/1‑68 (MQ=255)
 tGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACt                                                           <  1:984032/68‑1 (MQ=255)
         atcaGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACAGCTAGGACTCCGGCGat                                                   >  2:370935/1‑68 (MQ=255)
           caggcaggTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGatat                                                 <  2:571087/68‑1 (MQ=255)
               caggTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCgg                                             >  1:1040037/1‑68 (MQ=255)
               caggTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAAGACTCCGGCGATATCAgg                                             >  2:413680/1‑68 (MQ=37)
                    ataGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGt                                        >  2:238428/1‑68 (MQ=255)
                        ggCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCt                                    >  2:91435/1‑68 (MQ=255)
                         gCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCAAGGACTCCGGCGATATCCGGAACGTCCCtc                                   >  1:479221/1‑67 (MQ=255)
                           ccATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGatat                                                 <  1:893572/52‑1 (MQ=255)
                           ccATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGatat                                                 >  2:893572/1‑52 (MQ=255)
                              ttGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGgt                              <  2:159694/68‑1 (MQ=255)
                              ttGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGgt                              <  1:596375/68‑1 (MQ=255)
                              ttGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACGCCGGCGATATCCGGAACGTCCCTGCCGgt                              <  1:824497/68‑1 (MQ=255)
                               tGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACGCCGGCGATAACCGGATCGTCCCTGCCGgtg                             <  2:672103/68‑1 (MQ=21)
                                    tCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATg                        >  1:616282/1‑68 (MQ=255)
                                     ccTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGa                       >  1:147491/1‑68 (MQ=255)
                                          cggcgTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGAc                  >  2:131660/1‑68 (MQ=255)
                                                    cGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGc        <  1:468303/68‑1 (MQ=255)
                                                     gTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGcc       <  1:642862/68‑1 (MQ=255)
                                                         aCCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTg   >  1:492317/1‑68 (MQ=255)
                                                         aCCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTg   >  1:361225/1‑68 (MQ=255)
                                                          ccGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTAc             >  1:85301/1‑57 (MQ=255)
                                                          ccGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTAc             <  2:85301/57‑1 (MQ=255)
                                                          ccGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTgg  >  1:77078/1‑68 (MQ=255)
                                                              |                                                               
GTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCGAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGG  >  NC_000913/264036‑264161

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