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% 112.4 / ‑5.8 41I204S (ATC→AGC) wcaAputative glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base C (21/17);  minor base T (1/1);  total (23/18)
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

GGATCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGGT  >  NC_000913/2132647‑2132925
                                                                                                                                            |                                                                                                                                          
ggATCTGGTGGAGGGTAAACGGCTGATATTTGTTGCTGGCGCGGTGGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACTCGGTGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTc                                                                                                                                          <  2:94711/143‑1 (MQ=255)
     tGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTGGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGGCGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGt                                                                                                                                     >  1:224989/1‑143 (MQ=255)
      ggTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGGGGTGGTCTGCATTTCACCGTGATTGATATGCAGGCTCTGCggc                                                                                                                                    >  1:633581/1‑141 (MQ=255)
           agGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGGTCTGCGTCGccgc                                                                                                                               >  2:691253/1‑141 (MQ=255)
              gTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCtcttc                                                                                                                            >  1:25231/1‑143 (MQ=255)
              gTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACACGGAGAATTTTTTCGGTGACGAGGGGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCtcttc                                                                                                                            >  2:655491/1‑143 (MQ=255)
              gTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCGTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCtcttc                                                                                                                            <  1:433087/143‑1 (MQ=255)
                     gCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGTTCTGCATGTGACCTTGATTTATATGCAGGTTCTGCTACGTCTTTCCTACTTTc                                                                                                                     >  2:324544/1‑143 (MQ=255)
                       tGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGACTTTTTTCGGTGACGAGGTGATCTGCCTTTCACCGTGCTTGCTATGCAGGCTCTGCGTCGCCTCTTCTACTTTccg                                                                                                                   >  1:290436/1‑142 (MQ=255)
                        gagagtttgtGCCGGCGCGGTCGAATTTGTCTTTGAGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCAc                                                                                                                  <  1:285786/138‑1 (MQ=255)
                                                  ttGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCCCCATCCg                                                                                        >  1:236603/1‑143 (MQ=255)
                                                  ttGTCTGTGGGCTTGCGGTAAAAGTGGAATTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTACGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCg                                                                                        <  1:93803/143‑1 (MQ=255)
                                                       tttGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCATCTACAACCTCTTTTTCGCCGTCCTCCACCACCATCCGCCggc                                                                                   >  2:151364/1‑142 (MQ=255)
                                                          gtgCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGa                                                                                >  1:314365/1‑143 (MQ=255)
                                                             cTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCAGCACCATCCGCAGGAAGATGt                                                                             >  1:313744/1‑143 (MQ=255)
                                                              ttGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTc                                                                            >  1:218720/1‑143 (MQ=255)
                                                                  gctAAGAGTGGAGGTCCCCGGAGAATTTTTTCTGTGCCGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTaa                                                                        <  2:528979/141‑1 (MQ=255)
                                                                   gTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGCTGTCGTAAt                                                                       >  2:477840/1‑143 (MQ=255)
                                                                      aaaGTGGCAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGCGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCt                                                                    <  1:591178/143‑1 (MQ=255)
                                                                       aaGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCtg                                                                   >  1:413890/1‑143 (MQ=255)
                                                                        aGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCtgt                                                                  <  1:618175/143‑1 (MQ=255)
                                                                         gTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCtgtg                                                                 <  2:25231/143‑1 (MQ=255)
                                                                          tGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTCCTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGCAGATGTCGTAACCCTGTGc                                                                >  2:436939/1‑143 (MQ=255)
                                                                             aaGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGc                                                             <  1:257210/143‑1 (MQ=255)
                                                                              aGTACCCGGAGAATCTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTAGCCACCTCGACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCt                                                            <  1:655491/143‑1 (MQ=255)
                                                                               gTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCtt                                                           <  1:729243/143‑1 (MQ=255)
                                                                                      gagaTTTTTTTCTGTGATGCGGTGATCTGCATTTTGCCGTGATTGATATGCAGGTAATGCGCCGCTCTTTTTACTGTTCTTGGTTCGCCGTACTTCACATTCATTTGCACGCGTATTTCGTAATCCTGTGCCGCTTTCAGTTc                                                    <  2:359000/143‑1 (MQ=255)
                                                                                      gagaATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTCCAGTTc                                                    <  2:559920/143‑1 (MQ=255)
                                                                                       agaATGTTTTCGGTGACGAGGTGATCTGCATGTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACTGTTCGCCGGACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCg                                                   <  1:214348/143‑1 (MQ=255)
                                                                                                   gTGACGATGTGATCTGCATTTCACCGTGATTCATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGaacaaa                                       <  1:366008/143‑1 (MQ=255)
                                                                                                   gTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGaacaaa                                       <  1:510894/143‑1 (MQ=255)
                                                                                                        gAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCGGTACTCCACCACCATCCGCAGGAAGATGTCGTAATACTGTGCCGCTTTCAGTTCGGTATCGAGCAAACAgtc                                  >  1:241200/1‑143 (MQ=255)
                                                                                                         aGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCCGTATCGAACAACCActct                                 >  2:580205/1‑143 (MQ=255)
                                                                                                              aTCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCACGAAGCTGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTCGaa                            >  2:467242/1‑143 (MQ=255)
                                                                                                                 tGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCAAGGGTCGCCGTGCTCCAACCCCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACg                         >  1:615800/1‑143 (MQ=255)
                                                                                                                  gCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGACGTACTCCGCCACCATCCGCAGGACGATGTCGTAATCCTGTGCCGCTTTCAGTTCGCTATCGAACAACCACTCTTTGAACCGc                        >  2:654376/1‑143 (MQ=255)
                                                                                                                      ttCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATg                    >  1:11190/1‑143 (MQ=255)
                                                                                                                      ttCACCGTGATTGATAGGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATg                    <  1:144871/143‑1 (MQ=255)
                                                                                                                      ttCACCGTGATTGAGATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCAGCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATg                    <  1:145893/143‑1 (MQ=255)
                                                                                                                                      tccAGGCTCTGCGTCGCCTCTTCTACTTTCGCCGGTTCGCCGTACTCCACCCCCATCCGCAGCAAGATGTCGTACCCCTTTGCCGACTCCAGTTCGGCAGCGGAGAAACACTGCTTGCAACGCCATGCCCAGGTAGACACCTg    >  2:235362/3‑143 (MQ=255)
                                                                                                                                        cAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGAGGTTGTCCTACACCATCCGCAGGAAGATGTCGTCCTCCTGTGCCGATTTCCGTTAGGTATCGAACAAACAATCTTGGAGACGCGATGCTCAGGTAAAGACCAggg  >  2:157332/1‑142 (MQ=255)
                                                                                                                                            |                                                                                                                                          
GGATCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGGT  >  NC_000913/2132647‑2132925

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: