Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,255,335 A→G S206S (AGT→AGC ycgV ← putative autotransporter adhesin YcgV

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,255,3350AG96.4% 95.4 / ‑6.6 28S206S (AGT→AGCycgVputative autotransporter adhesin YcgV
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (14/13);  total (15/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATCGACAACAGAGTTTTTCTGTACGTTTATCCCCATGGCGCTGTAGCCCTGAACATCGATTGTCAGGTCTGTAGCCGTAAAACGCGCAGCACCATTGGCGTTATTGCCGTTGAGACCACCGATATAAACACCTGTACTTCCATCGGTCGTGATTTTACTTCCGCTTCCAAGATCGACACTGGTGCCATAGTCATTGATGGTTAGGCCTATACCGTTCGAGTTTTCAATGGTGAATTGAGTGGCTGTTAGTGTTGAGCTATGCCCAATAATGA  >  NC_000913/1255199‑1255471
                                                                                                                                        |                                                                                                                                        
gaTCGACAACAGAGTTTTTCTGTACGTTTATCCCCATGGCGCTGTAGCCCTGAACATCGATTGTCAGGTCTGTAGCCGTAAAACGCGCAGCACCATTGGCGTTATTGCCGTTGAGACCACCGATATAAACACCTGTGCt                                                                                                                                        <  2:316953/139‑1 (MQ=255)
gaTCGACAACAGAGTTTTTCTGTACGTTTATCCCCATGGCGCTGTAGCCCTGAACATCGATTGTCAGGTCTGTAGCCGTAAAACGCGCAGCACCATTGGCGTTATTGCCGTTGAGACCACCGATATAAACACCTGTGCt                                                                                                                                        >  1:316953/1‑139 (MQ=255)
    gACAACAGAGTTTTTCTGTACGTTTATCCCCATGGCGCTGTAGCCCTGAACATCGATTGTCAGGTCTGTAGCCGTAAAACGCGCAGCACCATTGGCGTTATTGCCGTTGAGACCACCGATATAAACACCTGTGCTTCCAt                                                                                                                                   >  1:308510/1‑140 (MQ=255)
    gACAACAGAGTTTTTCTGTACGTTTATCCCCATGGCGCTGTAGCCCTGAACATCGATTGTCAGGTCTGTAGCCGTAAAACGCGCAGCACCATTGGCGTTATTGCCGTTGAGACCACCGATATAAACACCTGTACTTCCAt                                                                                                                                   >  1:268716/1‑140 (MQ=255)
                            tATCCCCATGGCGCTGTAGCCCTGAACATCGATTGTCAGGTCTGTAGCCGTAAAACGCGCAGCACCATTGGCGTTATTGCCGTTGAGACCACCGATATAAACACCTGTGCTTCCATCGGTCGTGATTTTACTTCCGCTTc                                                                                                           >  1:170417/1‑140 (MQ=255)
                                  cATGGCGCTGTAGCCCTGAACATCGATTGTCAGGTCTGTAGCCGTAAAACGCGCAGCACCATTGGCGTTATTGCCGTTGAGACCACCGATATAAACACCTGTGCTTCCATCGGTCGTGATTTTACTTCCGCTTCCAAGAt                                                                                                     >  1:260917/1‑140 (MQ=255)
                                    tGGCGCTGTAGCCCTGAACATCGATTGTCAGGTCTGTAGCCGTAAAACGCGCAGCACCATTGGCGTTATTGCCGTTGAGACCACCGATATAAACACCTGTGCTTCCATCGGTCGTGATTTTACTTCCGCTTCCAAGATCg                                                                                                   >  1:342704/1‑140 (MQ=255)
                                    tGGCGCTGTAGCCCTGAACATCGATTGTCAGGTCTATAGCCGTAAAACGCGCAGCACCATTGGCGTTATTGCCGTTGAGACCACCGATATAAACACCTGTGCTTCCATCGGTCGTGATTTTACTTCCGCTTCCAAGATCg                                                                                                   >  1:54370/1‑140 (MQ=255)
                                     ggCGCTGTAGCCCTGAACATCGATTGTCAGGTCTGTAGCCGTAAAACGCGCAGCACCATTGGCGTTATTGCCGTTGAGACCACCGATATAAACACCTGTGCTTCCATCGGTCGTGATTTTACTTcc                                                                                                                <  1:389562/126‑1 (MQ=255)
                                     ggCGCTGTAGCCCTGAACATCGATTGTCAGGTCTGTAGCCGTAAAACGCGCAGCACCATTGGCGTTATTGCCGTTGAGACCACCGATATAAACACCTGTGCTTCCATCGGTCGTGATTTTACTTcc                                                                                                                >  2:389562/1‑126 (MQ=255)
                                           gTAGCCCTGAACATCGATTGTCAGGTCTGTAGCCGTAAAACGCGCAGCACCATTGGCGTTATTGCCGTTGAGACCACCGATATAAACACCTGTGCTTCCATCGGTCGTGATTTTACTTCCGCTTCCAAGATCGACACTgg                                                                                            <  2:308510/140‑1 (MQ=255)
                                                ccTGAACATCGATTGTCAGGTCTGTAGCCGTAAAACGCGCAGCACCATTGGCGTTATTGCCGTTGAGACCACCGATATAAACACCTGTGCTTCCATCGGTCGTGATTTTACTTCCGCTTCCAAGATCGACACTGGTGCCa                                                                                       <  2:342704/140‑1 (MQ=255)
                                                ccTGAACATCGATTGTCAGGTCTATAGCCGTAAAACGCGCAGCACCATTGGCGTTATTGCCGTTGAGACCACCGATATAAACACCTGTGCTTCCATCGGTCGTGATTTTACTTCCGCTTCCAAGATCGACACTGGTGCCa                                                                                       <  2:54370/140‑1 (MQ=255)
                                                    aaCATCGATTGTCAGGTCTGTAGCCGTAAAACGCGCAGCACCATTGGCGTTATTGCCGTTGAGACCACCGATATAAACACCTGTGCTTCCATCGGTCGTGATTTTACTTCCGCTTCCAAGATCGACACTGGTGCCATAGt                                                                                   >  2:381016/1‑140 (MQ=255)
                                                         cGATTGTCAGGTCTGTAGCCGTAAAACGCGCAGCACCATTGGCGTTATTGCCGTTGAGACCACCGATATAAACACCTGTGCTTCCATCGGTCGTGATTTTACTTCCGCTTCCAAGATCGACACTGGTGCCATAGTCATTg                                                                              <  1:407583/140‑1 (MQ=255)
                                                              gTCAGGTCTGTAGCCGTAAAACGCGCAGCACCATTGGCGTTATTGCCGTTGAGACCACCGATATAAACACCTGTGCTTCCATCGGTCGTGATTTTACTTCCGCTTCCAAGATCGACACTGGTGCCATAGTCATTGATGGt                                                                         <  1:411349/140‑1 (MQ=255)
                                                              gTCAGGTCTGTAGCCGTAAAACGCGCAGCACCATTGGCGTTATTGCCGTTGAGACCACCGATATAAACACCTGTGCTTCCATCGGTCGTGATTTTACTTCCGCTTCCAAGATCGACACTGGTGCCATAGTCATTGATGGt                                                                         <  1:471900/140‑1 (MQ=255)
                                                              gTCAGGTCTGTAGCCGTAAAACGCGCAGCACCATTGGCGTTATTGCCGTTGAGACCACCGATATAAACACCTGTGCTTCCATCGGTCGTGATTTTACTTCCGCTTCCAAGATCGACACTGGTGCCATAGTCATTGATGGt                                                                         <  2:260917/140‑1 (MQ=255)
                                                                                     cgcAGCACCATTGGCGTTATTGCCGTTGAGACCACCGATATAAACACCTGTGCTTCCATCGGTCGTGATTTTACTTCCGCTTCCAAGATCGACACTGGTGCTATAGTCATTGATGGTTAGGCCTATACCGTTCGAGtttt                                                  >  2:531476/1‑140 (MQ=255)
                                                                                             cATTGGCGTTATTGCCGTTGAGACCACCGATATAAACACCTGTGCTTCCATCGGTCGTGATTTTACTTCCGCTTCCAAGATCGACACTGGTGCCATAGTCATTGATGGTTAGGCCTATACCGTTCGAGTTTTCAATGGTg                                          <  1:276593/140‑1 (MQ=255)
                                                                                                            cGTTGAGACCACCGATATAAACACCTGTGCTTCCATCGGTCGTGATTTTACTTCCGCTTCCAAGATCGACACTGGTGCCATAGTCATTGATGGTTAGGCCTATACCGTTCGAGTTTTCAATGGTGAATTGAGTGGCTGtt                           >  1:444970/1‑140 (MQ=255)
                                                                                                                     caccGATATAAACACCTGTGCTTCCATCGGTCGTGATTTTACTTCCGCTTCCAAGATCGACACTGGTGCCATAGTCATTGATGGTTAGGCCTATACCGTTCGAGTTTTCAATGGTGAATTGAGTGGCTGTTAGTGTTGAg                  >  1:432652/1‑140 (MQ=255)
                                                                                                                       ccGATATAAACACCTGTGCTTCCATCGGTCGTGATTTTACTTCCGCTTCCAAGATCGACACTGGTGCCATAGTCATTGATGGTTAGGCCTATACCGTTCGAGTTTTCAATGGTGAATTGAGTGGCTGTTAGTGTTGAGCt                <  2:260500/140‑1 (MQ=255)
                                                                                                                         gATATAAACACCTGTGCTTCCATCGGTCGTGATTTTACTTCCGCTTCCAAGATCGACACTGGTGCCATAGTCATTGATGGTTAGGCCTATACCGTTCGAGTTTTCAATGGTGAATTGAGTGGCTGTTAGTGTTGAGCTAt              <  2:170417/140‑1 (MQ=255)
                                                                                                                                  acCTGTGCTTCCATCGGTCGTGATTTTACTTCCGCTTCCAAGATCGACACTGGTGCCATAGTCATTGATGGTTAGGCCTATACCGTTCGAGTTTTCAATGGTGAATTGAGTGGCTGTTAGTGTTGAGCTATGccc          >  1:173372/1‑135 (MQ=255)
                                                                                                                                  acCTGTGCTTCCATCGGTCGTGATTTTACTTCCGCTTCCAAGATCGACACTGGTGCCATAGTCATTGATGGTTAGGCCTATACCGTTCGAGTTTTCAATGGTGAATTGAGTGGCTGTTAGTGTTGAGCTATGccc          <  2:173372/135‑1 (MQ=255)
                                                                                                                                   ccTGTGCTTCCATCGGTCGTGATTTTACTTCCGCTTCCAAGATCGACACTGGTGCCATAGTCATTGATGGTTAGGCCTATACCGTTCGAGTTTTCAATGGTGAATTGAGTGGCTGTTAGTGTTGAGCTATGCCCaataat    >  1:433110/1‑140 (MQ=255)
                                                                                                                                    cTGTGCTTCCATCGGTCGTGATTTTACTTCCGCTTCCAAGATCGACACTGGTGCCATAGTCATTGATGGTTAGGCCTATACCGTTCGAGTTTTCAATGGTGAATTGAGTGGCTGTTAGTGTTGAGCTATGCCCAATAatg   >  1:203236/1‑140 (MQ=255)
                                                                                                                                     tgtgCTTCCATCGGTCGTGATTTTACTTCCGCTTCCAAGATCGACACTGGTGCCATAGTCATTGATGGTTAGGCCTATACCGTTCGAGTTTTCAATGGTGAATTGAGTGGCTGTTAGTGTTGAGCTATGCCCAATAatga  >  1:284362/1‑140 (MQ=255)
                                                                                                                                        |                                                                                                                                        
GATCGACAACAGAGTTTTTCTGTACGTTTATCCCCATGGCGCTGTAGCCCTGAACATCGATTGTCAGGTCTGTAGCCGTAAAACGCGCAGCACCATTGGCGTTATTGCCGTTGAGACCACCGATATAAACACCTGTACTTCCATCGGTCGTGATTTTACTTCCGCTTCCAAGATCGACACTGGTGCCATAGTCATTGATGGTTAGGCCTATACCGTTCGAGTTTTCAATGGTGAATTGAGTGGCTGTTAGTGTTGAGCTATGCCCAATAATGA  >  NC_000913/1255199‑1255471

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