Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,132,787 A→C I204S (ATC→AGC)  wcaA ← putative glycosyl transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,132,7870AC95.8% 64.1 / NA 25I204S (ATC→AGC) wcaAputative glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base C (11/12);  total (13/12)
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.96e-01

ATCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGGTT  >  NC_000913/2132649‑2132926
                                                                                                                                          |                                                                                                                                           
aTCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTg                                                                                                                                         >  2:531318/1‑143 (MQ=255)
 tCTGGTAGAGGGTAACCAGCTGATTTTTTTTGCTGGCGCGGTCTAATTTGTCTTTGTGCTTGCGGTAAAAGTTGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGc                                                                                                                                        <  2:383295/143‑1 (MQ=255)
     gTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATCTGACCGTGATTGATATGCAGGCTCTGCGTCg                                                                                                                                    >  1:131659/1‑143 (MQ=255)
         agGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTGCCCGGAGAGTTTTTTCGGTGAAGAGGTGATCTGCATTTTACTGGGATTGATAGGCAGGGTATGCGTTGCctc                                                                                                                                >  1:544221/1‑143 (MQ=255)
          gggTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCctct                                                                                                                               <  1:548893/143‑1 (MQ=255)
                     tGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCAACGTGATTGATATGCAGGCTCTGCGTAGGCTCTTCTTCtttacc                                                                                                                    >  1:80554/1‑140 (MQ=255)
                      gATATTTTTTGCTGGCGCGGTCGCATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGGGTCGCCTTTTCTACTTTCCAc                                                                                                                   >  2:623033/1‑143 (MQ=255)
                              ttGCTGGCGCGGTCGACGTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATCTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGcc                                                                                                           <  1:440633/143‑1 (MQ=255)
                                 cTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTAGCATGAGAATATTTGCGGTGACGAGGTGAGCAGCAGTTCACCGTGATTCAAACACAGGATCTGCGTCGCCTCTTCTACTTTCGACGGCTCCCCGTa                                                                                                        >  2:670065/1‑143 (MQ=255)
                                  tGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATACAGGCTCTGCGTCGCCTCTTATCCTTTCCACGGTTAGCCGTAc                                                                                                       >  2:176722/1‑143 (MQ=255)
                                         gTCGACGTTGTCTTTGTGCTTGCGGTAGAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGTTCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTccacca                                                                                                <  1:182425/143‑1 (MQ=255)
                                                   tCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAg                                                                                      <  1:687978/143‑1 (MQ=255)
                                                         tgCTTGCGGTAAAAGTGCAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGAt                                                                                <  1:166962/143‑1 (MQ=255)
                                                              gCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCCGGACGCTGCCGt                                                                           >  1:64888/1‑143 (MQ=255)
                                                                                cccGGAGACTTTTTTCGTTGACGAGGTGACCTGCATCTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCCTCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCa                                                         <  2:63497/143‑1 (MQ=255)
                                                                                         ttttttCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACACCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTAt                                                <  1:377447/143‑1 (MQ=255)
                                                                                                ggTGACGCGGTGCTCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGaacaa                                         <  1:637889/143‑1 (MQ=255)
                                                                                                  tGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGaaccaac                                       >  1:315420/1‑143 (MQ=255)
                                                                                                   gACGGGGTGAGCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGaacaaaca                                      <  1:13280/143‑1 (MQ=255)
                                                                                                            aTCTGCATTTCCCCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTCTGCCGCTTTCAGTGCGGTATCGAACCAACCCTCTTTGaa                             >  2:327456/1‑143 (MQ=255)
                                                                                                              cTGCATTTCACCGTGATTGATATGCAGGCTCTCCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTCTGCCGCTTTCAGTTCCGTATCGAACAAACACTCtgtgcgac                           >  2:380983/1‑139 (MQ=255)
                                                                                                               tGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGCTGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATGGAACAAACACTCTTTGAAACg                          >  1:665525/1‑143 (MQ=255)
                                                                                                                       aCCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACGAACACTCTTTGAAACGCAATGccc                  >  2:94463/1‑143 (MQ=255)
                                                                                                                                 atatGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACACCCGGCCCCAATAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGAc        <  1:531318/143‑1 (MQ=255)
                                                                                                                                       aGGCTCTTCGTCGCCTCTGCTGCGTTCCACGGTTCGCCGTACGCCACCGCCTTCCGCTGGAAGATGTCGTGATCCTGTTCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCATGTAAAGACCTGGtt  <  2:133786/143‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                           
ATCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGGTT  >  NC_000913/2132649‑2132926

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