Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,533,437 A→G Y139H (TAC→CAC)  yddH ← flavin reductase‑like protein YddH

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,533,4370AG95.2% 63.2 / ‑4.0 21Y139H (TAC→CAC) yddHflavin reductase‑like protein YddH
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (10/10);  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.95e-01

CAAACTGCCAGCGGCCTTCGACAAATACCCGTGCGTCTGCCGCTGCTGATACTACTTCGCCAAACAGCGTGTCGTATTCTTCTTGCGCAGAAGTCGCAGGTAGCAATCGACACTCCATCCACGCCAGACATTTTTCCTCGACCAGAGGC  >  NC_000913/1533362‑1533510
                                                                           |                                                                         
caAACTGCCAGCGGCCTTCGACAAATACCCGTGCGTCTGCCGCTGCTGATACTACTTCGCCAAACAGCGTGTCGTGttctt                                                                      >  2:287418/1‑81 (MQ=255)
    cTGCCAGCGGCCTTCGACAAATACCCGTGCGTCTGCCGCTGCTGATACTACTTCGCCAAACAGCGTGTCGTGTTCTTCTTg                                                                  <  1:287393/81‑1 (MQ=255)
        cAGCGGCCTTCGACAAATACCCGTGCGTCTGCCGCTGCTGATACTACTTCGCCAAACAGCGTGTCGTGTTCTTCTTGCGCa                                                              <  1:287394/81‑1 (MQ=255)
                   gACAAATACCCGTGCGTCTGCCGCTGCTGATACTACTTCGCCAAACAGCGTGTCGTGTTCTTCTTGCGCAGAAGTCGCAgg                                                   >  1:287428/1‑81 (MQ=255)
                        aTACCCGTGCGTCTGCCGCTGCTGATACTACTTCGCCAAACAGCGTGTCGTGTTCTTCTTGCGCAGAAGTCGCAGGTAGCa                                              <  2:287395/81‑1 (MQ=255)
                          aCCCGTGCGTCTGCCGCTGCTGATACTAATTCGCCAAACAGCGTGTCGTGTTCTTCTTGCGCAGAAGTCGCAGGTAGCAAt                                            >  2:287419/1‑81 (MQ=255)
                            ccGTGCGTCTGCCGCTGCTGATACTACTTCGCCAAACAGCGTGTCGTGTTCTTCTTGCGCAGAAGTCGCAGGTAGCAATCg                                          <  1:287396/81‑1 (MQ=255)
                                     tGCCGCTGCTGATACTACTTCGCCAAACAGCGTGTCGTGTTCTTCTTGCGCAGAAGTCGCAGGTAGCAATCGACACtccat                                 >  2:287426/1‑81 (MQ=255)
                                       ccGCTGCTGATACTACTTCGCCAAACAGCGTGTCGTGTTCTTCTTGCGCAGAAGTCGCAGGTAGCAATCGACACtccatcc                               <  2:287397/81‑1 (MQ=255)
                                             ctgATACTACTTCGCCAAACAGCGTGTCGTGTTCTTCTTGCGCAGAAGTCGCAGGTAGCAATCGACACTCCATCCACGc                           >  2:287420/1‑79 (MQ=255)
                                                aTACTACTTCGCCAAACAGCGTGTCGTGTTCTTCTTGCGCAGAAGTCGCAGGTAGCAATCGACACTCCATCCACGCCAGa                       >  2:287422/1‑80 (MQ=255)
                                                   ctactTCGCCAAACAGCGTGTCGTGTTCTTCTTGCGCAGAAGTCGCAGGTAGCAATCGACACTCCATCCACGCCAGACAtt                   <  1:287399/81‑1 (MQ=255)
                                                   ctactTCGCCAAACAGCGTGTCGTGTTCTTCTTGCGCAGAAGTCGCAGGTAGCAATCGACACTCCATCCACGCCAGACAtt                   <  1:287398/81‑1 (MQ=255)
                                                          gCCAAACAGCGTGTCGTGTTCTTCTTGCGCAGAAGTCGCAGGTAGCAATCGACACTCCATCCACGCCAGACATTTTTCCTc            <  1:287400/81‑1 (MQ=255)
                                                           ccAAACAGCGTGTCGTGTTCTTCTTGCGCAGAAGTCGCAGGTAGCAATCGACACTCCATCCACGCCAGACATTTTTCCTc            >  1:287423/1‑80 (MQ=255)
                                                           ccAAACAGCGTGTCGTGTTCTTCTTGCGCAGAAGTCGCAGGTAGCAATCGACACTCCATCCACGCCAGACATTTTTCCTc            >  2:287427/1‑80 (MQ=255)
                                                             aaaCAGCGTGTCGTATTCTTCTTGCGCAGAAGTCGCAGGTAGCAATCGACACTCCATCCACGCCAGACATTTTTCCTCGAc         >  2:287425/1‑81 (MQ=255)
                                                              aaCAGCGTGTCGTGTTCTTCTTGCGCAGAAGTCGCAGGTAGCAATCGACACTCCATCCACGCCAGACATTTTTCCTCGAc         >  1:287434/1‑80 (MQ=255)
                                                              aaCAGCGTGTCGTGTTCTTCTTGCGCAGAAGTCGCAGGTAGCAATCGACACTCCATCCACGCCAGACATTTTTCCTCGAc         >  2:287432/1‑80 (MQ=255)
                                                                cAGCGTGTCGTGTTCTTCTTGCGCAGAAGTCGCAGGTAGCAATCGACACTCCATCCACGCCAGACATTTTTCCTCGACCag      <  2:287401/81‑1 (MQ=255)
                                                                   cgtgtcgtgtTCTTCTTGCGCAGAAGTCGCAGGTAGCAATCGACACTCCATCCACGCCAGACATTTTTCCTCGACCAGAgg   >  2:287424/1‑81 (MQ=255)
                                                                    gtgtCGTGTTCTTCTTGCGCAGAAGTCGCAGGTAGCAATCGACACTCCATCCACGCCAGACATTTTTCCTCGACCAGAGGc  <  1:287402/81‑1 (MQ=255)
                                                                           |                                                                         
CAAACTGCCAGCGGCCTTCGACAAATACCCGTGCGTCTGCCGCTGCTGATACTACTTCGCCAAACAGCGTGTCGTATTCTTCTTGCGCAGAAGTCGCAGGTAGCAATCGACACTCCATCCACGCCAGACATTTTTCCTCGACCAGAGGC  >  NC_000913/1533362‑1533510

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