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,7870AC93.9% 88.0 / NA 33I204S (ATC→AGC) wcaAputative glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base A (2/0);  new base C (17/14);  total (19/14)
Fisher's exact test for biased strand distribution p-value = 4.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGATCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACC  >  NC_000913/2132647‑2132921
                                                                                                                                            |                                                                                                                                      
ggATCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGCCTTTTTTCGGTGACGAGGTCCTCTGCATATCAACGTGAATAATAAGCAAGAta                                                                                                                                      >  1:32841/1‑142 (MQ=255)
agATATGGTAGGGGGTAAACAGCTGATATGCGTTGCTGGCGCGGTCGAATTTGTCCTTGTGCTTGCTGTTAAAGTGGAAGTACCCGGACAATTTTTTCGGTGACGAGGAGATCGGCATTTCACCGTGATTGATATGCAGGCTc                                                                                                                                      <  1:296650/142‑1 (MQ=255)
 gATCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGCTCTGCATTTCACCGTGATTGATATGCAGGAGCt                                                                                                                                     >  1:658398/1‑143 (MQ=255)
      ggTAGAGGTTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTc                                                                                                                                <  1:84555/143‑1 (MQ=255)
      ggTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGTTCTGCATTCCACCGTGATCGATATGCAGGCCCTGCGTc                                                                                                                                >  2:150323/1‑143 (MQ=255)
       gTAGAGGGTAAACGGCGGATCGTGTTTGCTGGCGCGGTCGAATGTGTCTTTGTTCTTGCGGTAAAAGGGGAAGTTCTCGGTGAATTTTTTCGGGGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCg                                                                                                                               <  2:489867/143‑1 (MQ=255)
       gTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCg                                                                                                                               <  1:622407/143‑1 (MQ=255)
           agGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTTAACTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCctc                                                                                                                           >  1:570175/1‑143 (MQ=255)
             ggTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCtctt                                                                                                                         <  1:404192/143‑1 (MQ=255)
                    aGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTGTACttt                                                                                                                  >  2:546320/1‑143 (MQ=255)
                                     ggCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTATTCGGTGACGAGGTGATCTGCATTTCACCGTGATGGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGtccc                                                                                                 >  1:206625/1‑140 (MQ=255)
                                       cgcgGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGtcctcc                                                                                               >  2:278925/1‑143 (MQ=255)
                                           gTCTAATTTGTCTTTGTGCTTTCTGGAAAAGTGCAAGTACCCGGAGAATTTTTTCGGTGACGATGTGATCTGCCTTTCACTGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTccacca                                                                                           <  2:291393/143‑1 (MQ=255)
                                            tCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACGTTCCACGGTTCGCCGTACTccaccac                                                                                          >  1:367799/1‑143 (MQ=255)
                                              gAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGACGTACTccaccacca                                                                                        >  2:96506/1‑143 (MQ=255)
                                              gAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAAGCTCTGCGTCGCCTCTTCTACTGTCCCCGGTTCGCCGTACTccaccacga                                                                                        >  2:101264/1‑141 (MQ=255)
                                                   tGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGc                                                                                   >  1:233489/1‑143 (MQ=255)
                                                       tttGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCCCCATCCGCAGGa                                                                               >  1:211205/1‑143 (MQ=255)
                                                        ttGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGTTCTGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCCCAggcc                                                                              >  2:663633/1‑141 (MQ=255)
                                                                      aaaGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCt                                                                <  1:184276/143‑1 (MQ=255)
                                                                        aggggACGTCCCCGGAGAATTTTTGCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCtgt                                                              <  1:187253/140‑1 (MQ=255)
                                                                          tgctaGTACTCGGAGAATTTTTTCGGTGACGATGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCTACTGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGc                                                            <  2:109651/139‑1 (MQ=255)
                                                                               gTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCtt                                                       >  1:595965/1‑143 (MQ=255)
                                                                                  cccGGAGAATTTTTTCGGTGACGAGTTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTGCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCa                                                    <  1:481625/143‑1 (MQ=255)
                                                                                  cccGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCCCCTGCAGCTCTCGTAATCCCGcgccgcctcca                                                    >  2:719281/1‑138 (MQ=255)
                                                                                        gaATTTTTTCGGTGACGAGGTGATCTGCATTTCCCCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCgg                                              <  1:730447/143‑1 (MQ=255)
                                                                                            tttgtgggTGAGGAGGTGATCTGCATTTCACCTGGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATc                                          <  2:19948/137‑1 (MQ=255)
                                                                                                   gTGACGAGGTGATCTGCATTTCACCGTGATTGCTATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCTACCACCCTCCGCCGGACGATGTTGTGCTCCTTTCCCGCTTTCAGTTCGGTATCGaacaaa                                   <  2:54079/143‑1 (MQ=255)
                                                                                                   gTGACGAGGTGATCTGCATTTAACCGAGATTGCTATGCAGGCTCTGCGACGCCTCTTCTACTTTCAACGGTTCGCCGTACTGCACCACCATCCGCAGGAAGATGGCGAAATCCGATGCCGCTGTCAGTTCGGTATGGaacacc                                   >  2:381881/1‑141 (MQ=255)
                                                                                                     gagcgggTGATCTGCATTTCACGGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCGTGTGCCGCTTTCAGTTCGGTATCGaacaaaca                                 <  2:382548/138‑1 (MQ=255)
                                                                                                      aCGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTGGCGAACAAacac                                >  1:399636/1‑143 (MQ=255)
                                                                                                        gAGGTGATCTGCATTTCACAGTGATTGATATGCAGGCTCTGCGTCGCCTCATCTACTTTCCACGGTTCACCGTACTCCACCACCATCCGCAGGACGATGTCGTAATCCTGTGCGGCTTTCAATGCGGAAACGaacaaacggac                              >  2:249062/1‑139 (MQ=255)
                                                                                                                                    tatGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGCAAAGatc  >  1:193338/1‑141 (MQ=255)
                                                                                                                                            |                                                                                                                                      
GGATCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACC  >  NC_000913/2132647‑2132921

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: ‑