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.0% 109.1 / NA 41I204S (ATC→AGC) wcaAputative glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base A (2/0);  new base C (18/20);  total (20/21)
Fisher's exact test for biased strand distribution p-value = 4.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

GGATCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGAC  >  NC_000913/2132647‑2132920
                                                                                                                                            |                                                                                                                                     
ggATCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTc                                                                                                                                     >  1:546235/1‑143 (MQ=255)
  gtcgggTAGAGGTTGAGCAGCTGAGAGTTTTGGCTGGCGGGGTGGAATTTGGCTTTGTGCTTGCTGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTg                                                                                                                                   <  1:387852/139‑1 (MQ=255)
  aTCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTg                                                                                                                                   >  2:541473/1‑143 (MQ=255)
       gTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCg                                                                                                                              >  1:467003/1‑143 (MQ=255)
        tggacggTAAGCTGCTGATATTTGTTGCTTGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGATGTATCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGc                                                                                                                             <  1:61410/138‑1 (MQ=255)
        tAGAGGGTAAACAGCTGATCTTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTCCCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGc                                                                                                                             <  2:487426/143‑1 (MQ=255)
           agGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCctc                                                                                                                          <  2:281337/143‑1 (MQ=255)
              gTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCtcttc                                                                                                                       <  2:546736/143‑1 (MQ=255)
                    aGCTGATATGTTTTGCTGGCGCGGTCGAATTCGTCTTTGTGCCGGCGGTAAAAGTGGAAGTACCCGGAGAATTTGTTCGGTGACGAGGATATCTCCATTTCACCGTGATTTATATGCAGGCTCTTCGTCGCCTCTTCTACttt                                                                                                                 <  2:205097/143‑1 (MQ=255)
                       tGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCa                                                                                                              >  2:188900/1‑143 (MQ=255)
                       tGATATTTTTTGCTGGCGCGGTCGAAGTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTAGCCGGAGAATTTTTTCGGTGACAAGGTGATCTGAAGTTCACCGTGAGTGATATGCAGGCTCTGCGTAGCCTCTTCTACttgaca                                                                                                              >  1:540546/1‑139 (MQ=255)
                          tatTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACgg                                                                                                           <  2:661535/143‑1 (MQ=255)
                               gttGCTGTCGCGGTCTAATTTGTCTTTGTGCTTGCGCTCAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGc                                                                                                      <  2:431199/142‑1 (MQ=255)
                                   cTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAAGTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGGGGGATCTGCGTCGCCTCTGCTACTTTCGACGGTACGCCATa                                                                                                  >  2:175271/1‑143 (MQ=255)
                                      gcgcGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCTGTGACGAGGTGATCTGCATTTCACCGTGAGTGATATGCAGGCTCTTCGTCGCCTCTTTGATTTTTCACGGTTCGTCGTACtt                                                                                               >  2:151971/1‑142 (MQ=255)
                                           gTCGATTTTGTCTTGGTGCTCGCGGTAAAAGTGGAAGGACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTAATTGATATGCAGGCTCTGCGTCGCCGCTTCTACTTTCCACGGTTCGCCGTACTccacca                                                                                          <  2:657019/143‑1 (MQ=255)
                                                gcttGTCTGTGTGCTTTCCGTAAAAGTGGAATTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATc                                                                                     <  2:434807/141‑1 (MQ=255)
                                                    gTCTTTGTGCTAGCGGTAAAAGTGGAAGCACCCGGAGAATTTTTTCGGTGATGAGGTGATCTGCATTTCACCGTGATTGGTATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCa                                                                                 <  1:957693/143‑1 (MQ=255)
                                                        ttGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGGTATGCAGGCTCTGCGTCGCCTCTCTCACTTTCCGCGGTGTGCCGACCCCCACCCCCATCAGCCGGaa                                                                             >  1:917938/1‑143 (MQ=255)
                                                         tgtgCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAg                                                                            <  2:201459/143‑1 (MQ=255)
                                                          gtgCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGATCTGAGTCGGCCCTTCCACTTTCCACGGTTCGCCGTACTCCCCCACCATCCGCAGGAAGa                                                                           >  2:947563/1‑143 (MQ=255)
                                                            gCTTGCGGAAAACGTGGCAGCACCCGGAGACTTTTTCCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCCCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATg                                                                         <  2:255161/143‑1 (MQ=255)
                                                              ttGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTACGCCGTACTCCACCACCATCACCCGGACGATGTc                                                                       >  1:535542/1‑143 (MQ=255)
                                                                   gTAAAAGTGGAAGTACCCGGAGACTTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGCCGCCTCTTCTACTTTCCACGGTACGCCGTCCTCCACCACCACCCGCAGCATGATGTCGTagt                                                                  >  1:535961/1‑141 (MQ=255)
                                                                   gTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGCAGATGCCGTact                                                                  >  1:881495/1‑141 (MQ=255)
                                                                    tAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCCGGCCGATGTCCTAATc                                                                 >  2:870585/1‑143 (MQ=255)
                                                                          tGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCCGGAAGATGTCGTAATCCTGTGc                                                           >  2:863076/1‑143 (MQ=255)
                                                                          tGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGCGGATGTCGTAATCCTGTGc                                                           >  2:517123/1‑143 (MQ=255)
                                                                                   ccGGAGAATTTTTTCGGTGACGCGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCCGGAACGCGTCGTAATCCCGTCCCGCTTTCAg                                                  >  1:444135/1‑143 (MQ=255)
                                                                                       agaATTTTTTCGGTGACGAGGTGAGCTGCATTTCACCGTGATTGATATGCAGGCGCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTGCg                                              >  2:499427/1‑143 (MQ=255)
                                                                                            tttttGTGTGACGAGGTGATCTGCATGTCACCGTGATGGAGATGCAGGCGCTGCGTCGCCGCGTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATc                                         <  2:734002/143‑1 (MQ=255)
                                                                                             ttttCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCg                                        <  2:395606/143‑1 (MQ=255)
                                                                                                   gTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGaacaaa                                  <  2:881055/143‑1 (MQ=255)
                                                                                                      aCGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTTGCGCGTCCCCCAGACTCACAAACTATCCCCTTGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAacac                               <  2:444135/143‑1 (MQ=255)
                                                                                                               gcTGCATTTCACCGTGAGTGTGATGCAGGCTCTGCGTCGCCTCTTCTCCTCTCCACGGGTCGGCGTACGCCACCACCATCCGCAGCAAGACGTCGTAAGCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGaaa                      <  1:333828/142‑1 (MQ=255)
                                                                                                                       tCACCGTGATTGATATGCAGGCTCTGCGTCGCCGCTCCTACTTTCCTCGGTTCGCCGTACTCCACCACCATACGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGc              <  1:881489/143‑1 (MQ=255)
                                                                                                                           cGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTTTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGACATACACAg          >  1:518178/1‑143 (MQ=255)
                                                                                                                           cGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAg          >  2:922606/1‑143 (MQ=255)
                                                                                                                            attattgttTTGCAGGGTCGGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTTCTCCACCACCATCCGCAGGTAGATGTCTTACTCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACGCTGTGAAACGCCATGCCCAgg         <  1:628406/140‑1 (MQ=255)
                                                                                                                                gtgatAAGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTCCTCCACCACCATCCGCAGGCAGATGTCGTAATCCTGTGCCGCTTTCAGTCCGGTGACGAACAAACACTCTTTGAACCGCCACGCCCAGGGaac     >  1:70914/2‑142 (MQ=255)
                                                                                                                                   atacGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGAc  <  2:724243/143‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                     
GGATCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGAC  >  NC_000913/2132647‑2132920

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: