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