Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 2,409,844 +C coding (102/933 nt) yfdC → predicted inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6062,409,8401.C100.0% 56.8 / NA 17L33P (CTC→CCC) yfdCpredicted inner membrane protein
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base C (8/9);  total (8/9)

AACACAGCGACGAAATTGAAGTTGAGAGCGAAGAAAAAGAGCGCGGCAAAAAAATAGAAATAGATGAAGACCGACTCCCCTCCCGGGCGATGGCAATTCATGAGCATATCCGCCAGGATGGTGAAAAAGAGCTGGAACGCGACGCA  >  REL606/2409765‑2409910
                                                                             |                                                                      
aaCACAGCGACGAAATTGAAGTTGAGAGCGAAGAAAAAGAGCGCGGCAAAAAAATAGAAATAGATGAAGACCGAct                                                                         >  2:345033/1‑76 (MQ=255)
       cgacgaAATTGAAGTTGAGCGCGAAGAAAAAGAGCGCGGCAAAAAAATAGAAATAGATGAAGACCGActccccctc                                                                 >  1:184539/1‑73 (MQ=255)
       cgacgaAATTGAAGTTGAGAGCGAAGAAAAAGAGCGCGGCAAAAAAATAGAAATAGATGAAGACCGActccccctc                                                                 >  1:184544/1‑73 (MQ=255)
          cgaAATTGAAGTTGAGAGCGAAGAAAAAGAGCGCGGCAAAAAAATAGAAATAGATGAAGACCGACCCCCCCTCCCg                                                               <  1:1152615/76‑1 (MQ=255)
          cgaAATTGAAGTTGAGAGCGAAGAAAAAGAGCGCGGCAAAAAAATAGAAATAGATGAAGACCGACTCCCCCTCCCg                                                               <  1:1152619/76‑1 (MQ=255)
                 gAAGTTGAGAGCGAAGAAAAAGAGCGCGGCAAAAAAATAGAAATAGATGAAGACCGACTCCCCCTCCCGGGCGATg                                                        >  2:518771/1‑76 (MQ=255)
                 gAAGTTGAGAGCGAAGAAAAAGAGCGCGGCAAAAAAATAGAAATAGATGAAGACCGACTCCCCCTCCCGGGCGATg                                                        <  2:1197872/76‑1 (MQ=255)
                   aGTTGAGAGCGAAGAAAAAGAGCGCGGCAAAAAAATAGAAATAGATGAAGACCGACTCCCCCTCCCGGGCGATGGc                                                      >  2:989455/1‑76 (MQ=255)
                            cGAAGAAAAAGAGCGCGGCAAAAAAATAGAAATAGATGAAGACCGACTCCCCCTCCCGGGCGATGGCAATTCATGa                                             <  1:486017/76‑1 (MQ=255)
                                 aaaaaGAGCGCGGCAAAAAAATAGAAATAGATGAAGACCGACTCCCCCTCCCGGGCGATGGCAATTCATGAGCata                                        >  1:1218839/1‑76 (MQ=255)
                                    aaGAGCGCGGCAAAAAAATAGAAATAGATGAAGACCGACTCCCCCTCCCGGGCGATGGCAATTCATGAGCATATcc                                     <  1:1119030/76‑1 (MQ=255)
                                            ggCAAAAAAATAGAAATAGATGAAGACCGACTCCCCCTCCCGGGCGATGGCAATTCATGAGCATATCCGCCAGGAt                             <  2:778886/76‑1 (MQ=255)
                                                aaaaaaTAGAAATAGATGAAGACCGACTCCCCCTCCCGGGCGATGGCAATTCATGAGCATATCCGCCAGGATGGTg                         >  2:1205065/1‑76 (MQ=255)
                                                   aaatagaaatagATGAAGACCGACTCCCCCTCCCGGGCGATGGCAATTCATGAGCATATCCGCCAGGATGGTGaaa                      <  1:364219/76‑1 (MQ=255)
                                                    aatagaaatagaTGAAGACCGACTCCCCCTCCCGGGCGATGGCAATTCATGAGCATATCCGCCAGGATGGTGaaaa                     <  1:745360/76‑1 (MQ=255)
                                                      tagaaatagaTGAAGACCGACTCCCCCTCCCGGGCGATGGCAATTCATGAGCATATCCGCCAGGATGGTGAAAAag                   <  2:1152983/76‑1 (MQ=255)
                                                              gaTGAAGACCAACTCCCCCTCCCGGGCGATGGCAATTCATGAGCATATCCGCCAGGATGGTGAAAAAGAGCTGGaa           <  2:926703/76‑1 (MQ=255)
                                                                   aGACCGACTCCCCCTCCCGGGCGATGGCAATTCATGAGCATATCCGCCAGGATGGTGAAAAAGAGCTGGAACGCGa      >  1:991118/1‑76 (MQ=255)
                                                                    gACCGACTCCCCCTCCCGGGCGATGGCAATTCATGAGCATATCCGCCAGGATGGTGAAAAAGAGCTGGAACGCGAc     >  1:181270/1‑76 (MQ=255)
                                                                     aCCGACTCCCCCTCCCGGGCGATGGCAATTCATGAGCATATCCGCCAGGATGGTGAAAAAGAGCTGGAACGCGACg    <  1:1211926/76‑1 (MQ=255)
                                                                      ccGACTCCCCCTCCCGGGCGATGGCAATTCATGAGCATATCCGCCAGGATGGTGAAAAAGAGCTGGAACGCGACGc   >  1:15582/1‑76 (MQ=255)
                                                                       cGACTCCCCCTCCCGGGCGATGGCAATTCATGAGCATATCCGCCAGGATGGTGAAAAAGAGCTGGAACGCGACGCa  >  2:1028483/1‑76 (MQ=255)
                                                                             |                                                                      
AACACAGCGACGAAATTGAAGTTGAGAGCGAAGAAAAAGAGCGCGGCAAAAAAATAGAAATAGATGAAGACCGACTCCCCTCCCGGGCGATGGCAATTCATGAGCATATCCGCCAGGATGGTGAAAAAGAGCTGGAACGCGACGCA  >  REL606/2409765‑2409910

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: