Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,962,983 (T)5→4 coding (174/408 nt) ygdB ← DUF2509 domain‑containing protein YgdB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,962,9790T.90.5% 65.2 / 2.7 21coding (178/408 nt)ygdBDUF2509 domain‑containing protein YgdB
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base . (10/9);  total (11/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 = 9.14e-01

CAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAATGACTGACTTTCCATGCTCACGCGAGAGG  >  NC_000913/2962899‑2963053
                                                                                |                                                                          
cAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCAt                                                                            <  2:484260/81‑1 (MQ=255)
  gTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCAttt                                                                          <  1:484261/81‑1 (MQ=255)
  gTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCAttt                                                                          >  2:484302/1‑81 (MQ=255)
   tAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCAtttt                                                                         <  2:484262/81‑1 (MQ=255)
    aaaCGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCAttttc                                                                        >  1:484306/1‑80 (MQ=255)
       cGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTcccc                                                                     <  1:484263/81‑1 (MQ=255)
        gCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTcccca                                                                    >  1:484305/1‑80 (MQ=255)
          aaacaaacCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTCCCCACg                                                                 >  2:484299/1‑81 (MQ=255)
               aacCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTCCCCACGCCAgc                                                            >  1:484313/1‑81 (MQ=255)
                 ccTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTCCCCACGCCAgcgc                                                          >  2:484324/1‑81 (MQ=255)
                  cTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTCCCCACGCCAGCGcc                                                         <  2:484264/81‑1 (MQ=255)
                           ggTTTCAGCGTACTGCGAGCACTGAACTGCAGGCTGCGTCTGCCAGCAGTGCATTTTCCCCACGCCAGCGCCGACTGAAc                                                 <  2:484265/80‑1 (MQ=255)
                               ttaGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTCCCCACGCCAGCGCCGACTGAACGATg                                             <  2:484266/78‑1 (MQ=255)
                               tCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTCCCCACGCCAGCGCCGACTGAACGATgg                                            >  1:484307/1‑81 (MQ=255)
                                    gTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTgg                                       <  1:484267/81‑1 (MQ=255)
                                                        ctggctgCGTCTGCCAGCAGTGCATTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAATGACTGACttt                    <  1:484268/80‑1 (MQ=255)
                                                            ctgCGTCTGCCAGCAGTGCATTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAATGACTGACTTTCCATg               <  1:484269/81‑1 (MQ=255)
                                                                  cTGCCAGCAGTGCATTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAATGACTGACTTTCCATGCTCAcg         >  2:484308/1‑81 (MQ=255)
                                                                  cTGCCAGCAGTGCATTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAATGACTGACTTTCCATGCTCAcg         >  2:484316/1‑81 (MQ=255)
                                                                    gCCAGCAGTGCATTTTCCCCACGCCAGCGCCGACCGAACGATGGCCTGGCGGCCCAATGACTGACTTTCCATGCTCAcgc        >  1:484325/1‑80 (MQ=255)
                                                                     ccAGCAGTGCATTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAATGACTGACTTTCCATGCTCACGCga      <  1:484270/81‑1 (MQ=255)
                                                                      cagcagTGCATTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAATGACTGACTTTCCATGCTCAcgc        >  2:484309/1‑78 (MQ=255)
                                                                       agcagTGCATTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAATGACTGACTTTCCATGCTCACGCgaga    <  2:484271/81‑1 (MQ=255)
                                                                        gcagTGCATTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAATGACTGACTTTCCATGCTCACGCGagag   <  2:484272/81‑1 (MQ=255)
                                                                         cagTGCATTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAATGACTGACTTTCCATGCTCACGCGAGAgg  >  2:484318/1‑81 (MQ=255)
                                                                          agTGCATTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAATGACTGACTTTCCATGCTCACGCGAGAgg  >  1:484320/1‑80 (MQ=255)
                                                                                |                                                                          
CAGTAAACGCAAACAAACCTGGGCATCGGTTTCAGCGTACTGCGAGCACTGAACTGCTGGCTGCGTCTGCCAGCAGTGCATTTTTCCCCACGCCAGCGCCGACTGAACGATGGCCTGGCGGCGCAATGACTGACTTTCCATGCTCACGCGAGAGG  >  NC_000913/2962899‑2963053

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