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,0980GT100.0% 104.2 / NA 31L50L (CTC→CTAykfHuncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (19/12);  total (19/12)

GGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCGAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGT  >  NC_000913/264034‑264162
                                                                |                                                                
gggTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGGCCTGGCGGCGTTCGGCGTATACCGCTAgg                                                               <  2:83726/68‑1 (MQ=255)
  gTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAggcc                                                             >  2:170651/1‑66 (MQ=255)
       gCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCgg                                                        >  1:780009/1‑68 (MQ=255)
            tcaGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGata                                                   >  2:1388401/1‑68 (MQ=255)
            tcaGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGata                                                   >  1:642541/1‑68 (MQ=255)
             caggcaggTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACCCCGGCGatat                                                  <  2:814621/68‑1 (MQ=255)
                 caggTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCgg                                              >  2:247171/1‑68 (MQ=255)
                 caggTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCgg                                              >  1:589952/1‑68 (MQ=255)
                 caggTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCgg                                              >  1:83353/1‑68 (MQ=255)
                     tataGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACg                                          <  1:1141630/68‑1 (MQ=255)
                      ataGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGt                                         <  2:902683/68‑1 (MQ=255)
                      ataGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGt                                         <  2:83352/68‑1 (MQ=255)
                       taGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGata                                                   >  2:1260585/1‑57 (MQ=255)
                       taGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGata                                                   <  1:1260589/57‑1 (MQ=255)
                          ggCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCt                                     >  2:701415/1‑68 (MQ=255)
                          ggCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCt                                     >  2:473547/1‑68 (MQ=255)
                           gCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTg                                    >  2:838323/1‑68 (MQ=255)
                           gCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTg                                    >  1:1171306/1‑68 (MQ=255)
                           gCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTg                                    >  1:456470/1‑68 (MQ=255)
                           gCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTg                                    >  2:1305975/1‑68 (MQ=255)
                           gCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTg                                    >  2:1271167/1‑68 (MQ=255)
                           gCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTg                                    >  1:1082222/1‑68 (MQ=255)
                                cttcGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGgt                               <  2:1237985/65‑1 (MQ=255)
                                  gCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGgtgt                             >  2:116965/1‑68 (MQ=255)
                                      tCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCTGAACTTCCCTGCCGGTGTAATg                         <  1:541285/68‑1 (MQ=18)
                                      tCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATg                         <  2:830259/68‑1 (MQ=255)
                                      tCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATg                         >  1:622178/1‑68 (MQ=255)
                                       ccTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGa                        <  2:428138/68‑1 (MQ=255)
                                                gTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTa               <  2:714037/68‑1 (MQ=255)
                                                     gCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTg          <  2:622174/68‑1 (MQ=255)
                                                           aCCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTg    >  1:1358611/1‑68 (MQ=255)
                                                             cgctAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGt  <  1:110355/68‑1 (MQ=255)
                                                             cgctAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGt  <  2:325085/68‑1 (MQ=255)
                                                             cgctAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGt  <  1:237585/68‑1 (MQ=255)
                                                             cgctAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGt  <  2:1125615/68‑1 (MQ=255)
                                                                |                                                                
GGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCGAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGT  >  NC_000913/264034‑264162

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: