Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,454,025 T→A S458S (TCA→TCT gspA ← Type II secretion system protein GspA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,454,0250TA100.0% 56.9 / NA 17S458S (TCA→TCTgspAType II secretion system protein GspA
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base A (10/7);  total (10/7)

CCGGGGAAATATTTGCTTCATCCCGATATAAGTAGGCACTCTCCCCTGCCACCTGCCAGAGATGTACTAACGTTGAAAAGCCAACAACACCGTCAGTTTTCAAGTGGTGGCTTTTTTGAAATTGCTTTATTTTTTCTACCAACAAGGGACG  >  NC_000913/3453952‑3454102
                                                                         |                                                                             
ccGGGGAAATATTTGCTTCATCCCGATATAAGTAGGCACTCTCCCCTGCCACCTGCCAGAGATGTACTAACGTAGAAAAGc                                                                        >  1:610137/1‑81 (MQ=255)
ccGGGGAAATATTTGCTTCATCCCGATATAAGTAGGCACTCTCCCCTGCCACCTGCCAGAGATGTACTAACGTAGAAAAGc                                                                        >  2:610130/1‑81 (MQ=255)
     gAAATATTTGCTTCATCCCGATATAAGTAGGCACTCTCCCCTGCCACCTGCCAGAGATGTACTAACGTAGAAAAGCcaac                                                                    >  1:610126/1‑80 (MQ=255)
           tttGCTTCATCCCGATATAAGTAGGCACTCTCCCCTGCCACCTGCCAGAGATGTACTAACGTAGAAAAGCCAACAACAcc                                                              <  1:610112/80‑1 (MQ=255)
                          tataAGTAGGCACTCTCCCCTGCCACCTGCCAGAGATGTACTAACGTAGAAAAGCCAACAACACCGTCAGTTTTCAAgtgg                                              >  1:610143/1‑81 (MQ=255)
                              aGTAGGCACTCTCCCCTGCCACCTGCCAGAGATGTACTAACGTAGAAAAGCCAACAACACCGTCAGTTTTCAAGTGGTGGc                                          >  2:610136/1‑81 (MQ=255)
                              aGTAGGCACTCTCCCCTGCCACCTGCCAGAGATGTACTAACGTAGAAAAGCCAACAACACCGTCAGTTTTCAAGTGGTGGc                                          <  2:610113/81‑1 (MQ=255)
                                    cACTCTCCCCTGCCACCTGCCAGAGATGTACTAACGTAGAAAAGCCAACAACACCGTCAGTTTTCAAGTGGTGGCtttttt                                    >  2:610135/1‑81 (MQ=255)
                                        ctcCCCTGCCACCTGCCAGAGATGTACTAACGTAGAAAAGCCAACAACACCGTCAGTTTTCAAGTGGTGGCTTTTTTGaaa                                >  2:610139/1‑81 (MQ=255)
                                                ccacctgccaGAGATGTACTAACGTAGAAAAGCCAACAACACCGTCAGTTTTCAAGTGGTGGCTTTTTTGAAATTGCtt                          <  2:610114/79‑1 (MQ=255)
                                                ccacctgccaGAGATGTACTAACGTAGAAAAGCCAACAACACCGTCAGTTTTCAAGTGGTGGCTTTTTTGAAATTGCTTTa                        <  1:610115/81‑1 (MQ=255)
                                                    ctgccaGAGATGTACTAACGTAGAAAAGCCAACAACACCGTCAGTTTTCAAGTGGTGGCTTTTTTGAAATTGCTTTAtttt                    <  2:610116/81‑1 (MQ=255)
                                                           agaTGTACTAACGTAGAAAAGCCAACAACACCGTCAGTTTTCAAGTGGTGGCTTTTTTGAAATTGCTTTATTTTTTCTAc              >  2:610142/1‑80 (MQ=255)
                                                                  cTAACGTAGAAAAGCCAACAACACCGTCAGTTTTCAAGTGGTGGCTTTTTTGAAATTGCTTTATTTTTTCTACCAACAAgg      >  1:610147/1‑81 (MQ=255)
                                                                     aCGTAGAAAAGCCAACAACACCGTCAGTTTTCAAGTGGTGGCTTTTTTGAAATTGCTTTATTTTTTCTACCAACAAGGGAc   <  1:610118/81‑1 (MQ=255)
                                                                     aCGTAGAAAAGCCAACAACACCGTCAGTTTTCAAGTGGTGGCTTTTTTGAAATTGCTTTATTTTTTCTACCAACAAGGGAc   <  1:610117/81‑1 (MQ=255)
                                                                      cGTAGAAAAGCCAACAACACCGTCAGTTTTCAAGTGGTGGCTTTTTTGAAATTGCTTTATTTTTTCTACCAACAAGGGACg  >  2:610154/1‑81 (MQ=255)
                                                                         |                                                                             
CCGGGGAAATATTTGCTTCATCCCGATATAAGTAGGCACTCTCCCCTGCCACCTGCCAGAGATGTACTAACGTTGAAAAGCCAACAACACCGTCAGTTTTCAAGTGGTGGCTTTTTTGAAATTGCTTTATTTTTTCTACCAACAAGGGACG  >  NC_000913/3453952‑3454102

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