Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,504,097 C→T R14W (CGG→TGG)  yjhB → putative sialic acid transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,504,0970CT96.0% 84.5 / ‑3.3 25coding (40/1218 nt)yjhBputative sialic acid transporter
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (18/6);  minor base . (0/1);  total (18/7)
Fisher's exact test for biased strand distribution p-value = 2.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGCAACAGCATGGTATAAACAAGTTAATCCACCACAACGGAAAGCTCTTTTTTCCGCATGGCTTGGATATGTATTTGATGGCTTTGATTTTATGATGATATTTTACATTCT  >  NC_000913/4504020‑4504170
                                                                             |                                                                         
aCTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGCAACAGCATGGTATAAACAAGTTAATCCACCACAATgg                                                                         <  2:905327/80‑1 (MQ=255)
aCTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGCAACAGCATGGTATAAACAAGTTAATCCACCACAATgg                                                                         >  2:905360/1‑80 (MQ=255)
     aaGGGCTATATTAGAATAACACAGGAAACAAATATGGCAACAGCATGGTATAAACAAGTTAATCCACCACAATGGAAAGct                                                                   >  2:905361/1‑81 (MQ=255)
     aaGGGCTATATTAGAATAACACAGGAAACAAATATGGCAACAGCATGGTATAAACAAGTTAATCCACCACAATGGAAAGc                                                                    <  2:905328/80‑1 (MQ=255)
        ggCTATATTAGAATAACACAGGAAACAAATATGGCAACAGCATGGTATAAACAAGTTAATCCACCACAATGGAAAGCTCtt                                                                <  2:905329/81‑1 (MQ=255)
              atTAGAATAACACAGGAAACAAATATGGCAACAGCATGGTATAAACAAGTTAATCCACCACAATGGAAAGCTCTTTTTTcc                                                          >  2:905357/1‑81 (MQ=255)
                  gAATAACACAGGAAACAAATATGGCAACAGCATGGTATAAACAAGTTAATCCACCACAATGGAAAGCTCTTTTTTCCGCAt                                                      >  2:905365/1‑81 (MQ=255)
                       acacAGGAAACAAATATGGCAACAGCATGGTATAAACAAGTTAATCCACCACAATGGAAAGCTCTTTTTTCCGCATGGCt                                                  <  1:905330/80‑1 (MQ=255)
                            ggAAACAAATATGGCAACAGCATGGTATAAACAAGTTAATCCACCACAATGGAAAGCTCTTTTTTCCGCATGGCTTGGata                                            >  2:905352/1‑81 (MQ=255)
                             gAAACAAATATGGCAACAGCATGGTATAAACAAGTTAATCCACCACAATGGAAAGCTCTTTTTTCCGCATGGCTTGGatat                                           <  2:905331/81‑1 (MQ=255)
                             gAAACAAATATGGCAACAGCATGGTATAAACAAGTTAATCCACCACAATGGAAAGCTCTTTTTTCCGCATGGCTTGGatat                                           >  2:905353/1‑81 (MQ=255)
                                       tGGCAACAGCATGGTATAAACAAGTTAATCCACCACAATGGAAAGCTCTTTTTTCCGCATGGCTTGGATATGTATTTGATg                                 >  2:905356/1‑81 (MQ=255)
                                         gCAACAGCATGGTATAAACAAGTTAATCCACCACAATGGAAAGCTCTTTTTTCCGCATGGCTTGGATATGTATTTGATGGc                               >  2:905351/1‑81 (MQ=255)
                                          cAACAGCATGGTATAAACAAGTTAATCCACCACAATGGAAAGCTCTTTTTTCCGCATGGCTTGGATATGTATTTGATGGCt                              >  2:905366/1‑81 (MQ=255)
                                              aGCATGGTATAAACAAGTTAATCCACCACAATGGAAAGCTCTTTTTTCCGCATGGCTTGGATATGTATTTGATGGCTTTGa                          >  1:905354/1‑81 (MQ=255)
                                                 aTGGTATAAACAAGTTAATCCACCACAATGGAAAGCTCTTTTTTCCGCATGGCTTGGATATGTATTTGATGGCTTTGAttt                       <  2:905332/81‑1 (MQ=255)
                                                   ggTATAAACAAGTTAATCCACCACAATGGAAAGCTCTTTTTTCCGCATGGCTTGGATATGTATTTGATGGCTTTGATTTTa                     >  1:905374/1‑81 (MQ=255)
                                                    gTATAAACAAGTTAATCCACCACAATGGAAAGCTCTTTTTTCCGCATGGCTTGGATATGTATTTGATGGCTTTGATTTTat                    >  2:905358/1‑81 (MQ=255)
                                                        aaaCAAGTTAATCCACCACAATGGAAAGCTCTTTTTTCCGCATGGCTTGGATATGTATTTGATGGCTTTGATTTTatgat                 >  2:905368/1‑80 (MQ=255)
                                                         aaCAAGTTAATCCACCACAATGGAAAGCTCTTTTTTCCGCATGGCTTGGATATGTATTTGATGGCTTTGATTTTAtgatga               >  1:905373/1‑81 (MQ=255)
                                                          aCAAGTTAATCCACCACAATGGAAAGCTCTTTTTTCCGCATGGCTTGGATATGTATTTGATGGCTTTGATTTTAtgatga               >  1:905363/1‑80 (MQ=255)
                                                                tAATCCACCACAATGGAAAGCTCTTTTTTCCGCATGGCTTGGATATGTATTTGATGGCTTTGATTTTATGATGATATTTTa        >  1:905383/1‑81 (MQ=255)
                                                                  aTCCACCACAATGGAAAGCTCTTTTTTCCGCATGGCTTGGATATGTATTTGATGGCTTTGATTTTATGATGATATTTTACa      >  2:905367/1‑81 (MQ=255)
                                                                     caccacAATGGAAAGCTCTTTTTTCCGCATGGCTTGGATATGTATTTGATGGCTTTGATTTTATGATGATATTTTACAttc   >  1:905372/1‑81 (MQ=255)
                                                                     caccacAATGAAAGCTCTTTTTTCCGCATGGCTTGGATATGTATTTGATGGCTTTGATTTTATGATGATATTTTACAttct  <  2:905333/81‑1 (MQ=255)
                                                                             |                                                                         
ACTTTAAGGGCTATATTAGAATAACACAGGAAACAAATATGGCAACAGCATGGTATAAACAAGTTAATCCACCACAACGGAAAGCTCTTTTTTCCGCATGGCTTGGATATGTATTTGATGGCTTTGATTTTATGATGATATTTTACATTCT  >  NC_000913/4504020‑4504170

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