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.1% 115.9 / ‑5.9 42I204S (ATC→AGC) wcaAputative glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base C (21/18);  minor base T (1/1);  total (23/19)
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                                                                                                                                          <  1:94711/143‑1 (MQ=255)
     tGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTGGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGGCGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGt                                                                                                                                     >  2:224989/1‑143 (MQ=255)
      ggTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGGGGTGGTCTGCATTTCACCGTGATTGATATGCAGGCTCTGCggc                                                                                                                                    >  2:633581/1‑141 (MQ=255)
           agGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGGTCTGCGTCGccgc                                                                                                                               >  1:691253/1‑141 (MQ=255)
              gTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCtcttc                                                                                                                            >  2:25231/1‑143 (MQ=255)
              gTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACACGGAGAATTTTTTCGGTGACGAGGGGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCtcttc                                                                                                                            >  1:655491/1‑143 (MQ=255)
              gTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCGTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCtcttc                                                                                                                            <  2:433087/143‑1 (MQ=255)
                     gttgactggtgtggCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTc                                                                                                                     <  1:77914/130‑1 (MQ=255)
                     gCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGTTCTGCATGTGACCTTGATTTATATGCAGGTTCTGCTACGTCTTTCCTACTTTc                                                                                                                     >  1:324544/1‑143 (MQ=255)
                       tGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGACTTTTTTCGGTGACGAGGTGATCTGCCTTTCACCGTGCTTGCTATGCAGGCTCTGCGTCGCCTCTTCTACTTTccg                                                                                                                   >  2:290436/1‑142 (MQ=255)
                        gagagtttgtGCCGGCGCGGTCGAATTTGTCTTTGAGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCAc                                                                                                                  <  2:285786/138‑1 (MQ=255)
                                                  ttGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCCCCATCCg                                                                                        >  2:236603/1‑143 (MQ=255)
                                                  ttGTCTGTGGGCTTGCGGTAAAAGTGGAATTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTACGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCg                                                                                        <  2:93803/143‑1 (MQ=255)
                                                       tttGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCATCTACAACCTCTTTTTCGCCGTCCTCCACCACCATCCGCCggc                                                                                   >  1:151364/1‑142 (MQ=255)
                                                          gtgCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGa                                                                                >  2:314365/1‑143 (MQ=255)
                                                             cTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCAGCACCATCCGCAGGAAGATGt                                                                             >  2:313744/1‑143 (MQ=255)
                                                              ttGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTc                                                                            >  2:218720/1‑143 (MQ=255)
                                                                  gctAAGAGTGGAGGTCCCCGGAGAATTTTTTCTGTGCCGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTaa                                                                        <  1:528979/141‑1 (MQ=255)
                                                                   gTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGCTGTCGTAAt                                                                       >  1:477840/1‑143 (MQ=255)
                                                                      aaaGTGGCAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGCGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCt                                                                    <  2:591178/143‑1 (MQ=255)
                                                                       aaGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCtg                                                                   >  2:413890/1‑143 (MQ=255)
                                                                        aGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCtgt                                                                  <  2:618175/143‑1 (MQ=255)
                                                                         gTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCtgtg                                                                 <  1:25231/143‑1 (MQ=255)
                                                                          tGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTCCTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGCAGATGTCGTAACCCTGTGc                                                                >  1:436939/1‑143 (MQ=255)
                                                                             aaGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGc                                                             <  2:257210/143‑1 (MQ=255)
                                                                              aGTACCCGGAGAATCTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTAGCCACCTCGACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCt                                                            <  2:655491/143‑1 (MQ=255)
                                                                               gTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCtt                                                           <  2:729243/143‑1 (MQ=255)
                                                                                      gagaTTTTTTTCTGTGATGCGGTGATCTGCATTTTGCCGTGATTGATATGCAGGTAATGCGCCGCTCTTTTTACTGTTCTTGGTTCGCCGTACTTCACATTCATTTGCACGCGTATTTCGTAATCCTGTGCCGCTTTCAGTTc                                                    <  1:359000/143‑1 (MQ=255)
                                                                                      gagaATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTCCAGTTc                                                    <  1:559920/143‑1 (MQ=255)
                                                                                       agaATGTTTTCGGTGACGAGGTGATCTGCATGTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACTGTTCGCCGGACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCg                                                   <  2:214348/143‑1 (MQ=255)
                                                                                                   gTGACGATGTGATCTGCATTTCACCGTGATTCATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGaacaaa                                       <  2:366008/143‑1 (MQ=255)
                                                                                                   gTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGaacaaa                                       <  2:510894/143‑1 (MQ=255)
                                                                                                        gAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCGGTACTCCACCACCATCCGCAGGAAGATGTCGTAATACTGTGCCGCTTTCAGTTCGGTATCGAGCAAACAgtc                                  >  2:241200/1‑143 (MQ=255)
                                                                                                         aGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCCGTATCGAACAACCActct                                 >  1:580205/1‑143 (MQ=255)
                                                                                                              aTCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCACGAAGCTGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTCGaa                            >  1:467242/1‑143 (MQ=255)
                                                                                                                 tGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCAAGGGTCGCCGTGCTCCAACCCCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACg                         >  2:615800/1‑143 (MQ=255)
                                                                                                                  gCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGACGTACTCCGCCACCATCCGCAGGACGATGTCGTAATCCTGTGCCGCTTTCAGTTCGCTATCGAACAACCACTCTTTGAACCGc                        >  1:654376/1‑143 (MQ=255)
                                                                                                                      ttCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATg                    >  2:11190/1‑143 (MQ=255)
                                                                                                                      ttCACCGTGATTGATAGGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATg                    <  2:144871/143‑1 (MQ=255)
                                                                                                                      ttCACCGTGATTGAGATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCAGCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATg                    <  2:145893/143‑1 (MQ=255)
                                                                                                                                      tccAGGCTCTGCGTCGCCTCTTCTACTTTCGCCGGTTCGCCGTACTCCACCCCCATCCGCAGCAAGATGTCGTACCCCTTTGCCGACTCCAGTTCGGCAGCGGAGAAACACTGCTTGCAACGCCATGCCCAGGTAGACACCTg    >  1:235362/3‑143 (MQ=255)
                                                                                                                                        cAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGAGGTTGTCCTACACCATCCGCAGGAAGATGTCGTCCTCCTGTGCCGATTTCCGTTAGGTATCGAACAAACAATCTTGGAGACGCGATGCTCAGGTAAAGACCAggg  >  1: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: