Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,476,830 T→C Y563H (TAT→CAT)  ydbD → DUF2773 domain‑containing protein YdbD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,476,8300TC88.2% 39.7 / ‑0.7 17Y563H (TAT→CAT) ydbDDUF2773 domain‑containing protein YdbD
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base C (7/8);  total (8/9)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.78e-01

GCACGTCGATGAAACAAAAAATTTGGCAACGTGCTTTATCGCATCGGAAAAATAATCCGAAAGCAGATTCAGATGCATATGAAACCAGCGCTGATATGATTTTGTCTGAGCTGATTAGCCACGGCGAAGTTGATGACCAGATGCTAT  >  NC_000913/1476753‑1476899
                                                                             |                                                                     
gCACGTCGATGAAACAAAAAATTTGGCAACGTGCTTTATCGCATCGGAAAAATAATCCGAAAGCAGATTCAGATGCACATg                                                                    <  1:279324/81‑1 (MQ=255)
gCACGTCGATGAAACAAAAAATTTGGCAACGTGCTTTATCGCATCGGAAAAATAATCCGAAAGCAGATTCAGATGCACATg                                                                    <  1:279323/81‑1 (MQ=255)
  aCGTCGATGAAACAAAAAATTTGGCAACGTGCTTTATCGCATCGGAAAAATAATCCGAAAGCAGATTCAGATGCATATGaa                                                                  <  1:279325/81‑1 (MQ=255)
  aCGTCGATGAAACAAAAAATTTGGCAACGTGCTTTATCGCATCGGAAAAATAATCCGAAAGCAGATTCAGATGCACATGa                                                                   >  2:279349/1‑80 (MQ=255)
     tCGATGAAACAAAAAATTTGGCAACGTGCTTTATCGCATCGGAAAAATAATCCGAAAGCAGATTCAGATGCATATGAAAc                                                                >  2:279342/1‑80 (MQ=255)
     tCGATGAAACAAAAAATTTGGCAACGTGCTTTATCGCATCGGAAAAATAATCCGAAAGCAGATTCAGATGCACATGAAAc                                                                >  2:279340/1‑80 (MQ=255)
        aTGAAACAAAAAATTTGGCAACGTGCTTTATCGCATCGGAAAAATAATCCGAAAGCAGATTCAGATGCACATGAAACCAgc                                                            <  1:279326/81‑1 (MQ=255)
           aaaCAAAAAATTTGGCAACGTGCTTTATCGCATCGGAAAAATAATCCGAAAGCAGATTCAGATGCACATGAAACCAgcgc                                                          <  1:279327/80‑1 (MQ=255)
           aaaCAAAAAATTTGGCAACGTGCTTTATCGCATCGGAAAAATAATCCGAAAGCAGATTCAGATGCACATGAAACCAGCGCt                                                         >  2:279347/1‑81 (MQ=255)
                        ggCAACGTGCTTTATCGCATCGGAAAAATAATCCGAAAGCAGATTCAGATGCACATGAAACCAGCGCTGATATGATTTTg                                             <  2:279328/80‑1 (MQ=255)
                              gTGCTTTATCGCATCGGAAAAATAATCCGAAAGCAGATTCAGATGCACATGAAACCAGCGCTGATATGATTTTGTCTGAGc                                      >  2:279354/1‑81 (MQ=255)
                                       cGCATCGGAAAAATAATCCGAAAGCAGATTCAGATGCACATGAAACCAGCGCTGATATGATTTTGTCTGAGCTGATTAGc                              >  1:279348/1‑80 (MQ=255)
                                         cATCGGAAAAATAATCCGAAAGCAGATTCAGATGCACATGAAACCAGCGCTGATATGATTTTGTCTGAGCTGATTAGCCAc                           <  2:279329/81‑1 (MQ=255)
                                                      atCCGAAAGCAGATTCAGATGCACATGAAACCAGCGCTGATATGATTTTGTCTGAGCTGATTAGCCACGGCGAAGTtgat               <  2:279330/80‑1 (MQ=255)
                                                          gAAAGCAGATTCAGATGCACATGAAACCAGCGCTGATATGATTTTGTCTGAGCTGATTAGCCACGGCGAAGTTGATGAc            >  2:279351/1‑79 (MQ=255)
                                                          gAAAGCAGATTCAGATGCACATGAAACCAGCGCTGATATGATTTTGTCTGAGCTGATTAGCCACGGCGAAGTTGATGACCa          <  2:279331/81‑1 (MQ=255)
                                                                  aTTCAGATGCACATGAAACCAGCGCTGATATGATTTTGTCTGAGCTGATTAGCCACGGCGAAGTTGATGACCAGATGCTAt  >  1:279358/1‑81 (MQ=255)
                                                                             |                                                                     
GCACGTCGATGAAACAAAAAATTTGGCAACGTGCTTTATCGCATCGGAAAAATAATCCGAAAGCAGATTCAGATGCATATGAAACCAGCGCTGATATGATTTTGTCTGAGCTGATTAGCCACGGCGAAGTTGATGACCAGATGCTAT  >  NC_000913/1476753‑1476899

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