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,7870AC98.2% 157.7 / NA 54I204S (ATC→AGC) wcaAputative glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base C (28/25);  minor base T (0/1);  total (28/26)
Fisher's exact test for biased strand distribution p-value = 4.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGATCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGGTT  >  NC_000913/2132647‑2132926
                                                                                                                                            |                                                                                                                                           
ggATCTGGTAGAGGGGAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATGTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGTATAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTc                                                                                                                                           <  1:1039359/143‑1 (MQ=255)
 gATCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGctct                                                                                                                                          <  1:1045218/143‑1 (MQ=255)
   tCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGc                                                                                                                                        >  1:557612/1‑143 (MQ=255)
    cTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCg                                                                                                                                       >  1:528149/1‑143 (MQ=255)
        tAGAGGGTAGACAGCTGATATTTTGTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGGAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGc                                                                                                                                   <  2:418967/143‑1 (MQ=255)
             ggTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATATCACCGTGATTGATATGCAGGCTCTGCGTCGCCtctt                                                                                                                              >  1:758162/1‑143 (MQ=255)
              gTAGACAGCTGATATTTTTTGCTGGCGCGGTCGACTTTATCATTGTGCTTGCTGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCtcttc                                                                                                                             <  2:1063508/143‑1 (MQ=255)
               tAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCtcttct                                                                                                                            >  1:242263/1‑143 (MQ=255)
                aaaCAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTa                                                                                                                           <  1:863101/143‑1 (MQ=255)
                   cAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCCACtt                                                                                                                        >  1:611821/1‑143 (MQ=255)
                       tGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCGTCTACTTTCCa                                                                                                                    >  2:1087757/1‑143 (MQ=255)
                       tGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCGTCTACTTTCCa                                                                                                                    >  1:819121/1‑143 (MQ=255)
                       tGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGAGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCa                                                                                                                    >  1:755633/1‑143 (MQ=255)
                         atatTTTTTGCTGGCGCGGTCGAAGTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACg                                                                                                                  <  1:744605/143‑1 (MQ=255)
                         atatTTTTTGCTGGCGCGGTCGAACTTCTCTTTGTGCTTTCGGTAACAGTGGCAGTACCCGGAGAATTTTTTCGGTGACGGGGTGATCTGCATTTCACCGTGATTGCTATGCACGCCCTGCCGCGCCCCTCCTCCACTCCACg                                                                                                                  >  2:953347/1‑143 (MQ=255)
                           atTTTTTGCTGGCGCGGTCGAAATTGTCTTGGTGCTTGCGGTAAAAGTGGAACTACTAGGAGAATTTTTTGGTTGTCGAGGTGATCTGCATTTCACCGGGATTGCTATGCAGGCTCTGCGTCGCCTCTTCTGCTTTCCACGGt                                                                                                                >  2:778076/1‑143 (MQ=255)
                                    tGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTAc                                                                                                       >  1:419811/1‑143 (MQ=255)
                                      gcgcGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGGACtg                                                                                                     >  1:132557/1‑142 (MQ=255)
                                         cgGTCGAATTGGTCTGTGTGGCTGCTGTAAAAGGGGAAGTACCCTGATGATTTTGTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTccac                                                                                                  <  2:1046583/143‑1 (MQ=255)
                                             cGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTccaccacc                                                                                              <  1:668943/143‑1 (MQ=255)
                                              gagtttGTGTTCGTGGTTGGGGTAGGTGTGTTAGTGCCCGGAGAATCTTTTCGGTGAGGAGGTGATCTGCATTGCACCGTGATTGATATGCAGGCTCTCCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTccaccacca                                                                                             <  2:877722/140‑1 (MQ=255)
                                              gAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTccaccacca                                                                                             >  2:617571/1‑143 (MQ=255)
                                              gAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCACCGTACTccaccacca                                                                                             >  2:328267/1‑143 (MQ=255)
                                              gAATTTGTCCTGGTGCCTGCGGTAAAAGCGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTccaccacca                                                                                             <  1:66669/143‑1 (MQ=255)
                                                 tttGTCTTTGTGCTTGGTGTAAAGGTGGAAGTACCCGGAGAAGGGTTGCGGGGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCGTCTACTTTCCACGGTTCGCCGTACTCCACCACCATcc                                                                                          <  1:128288/143‑1 (MQ=255)
                                                   tGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCCTTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTAGCCGTACTGCACCACCATCCGc                                                                                        >  2:1011762/1‑143 (MQ=255)
                                                    gTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCa                                                                                       <  1:700858/143‑1 (MQ=255)
                                                    gTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCa                                                                                       <  1:161401/143‑1 (MQ=255)
                                                       tttGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACTGGTCGCCGTACTCCACCACCCTCCGCAGGa                                                                                    >  2:307370/1‑143 (MQ=255)
                                                        ttGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCCCGCTTCGCCGTACCGCACCCCCATcagcagcac                                                                                   >  2:849161/1‑140 (MQ=255)
                                                              ttGCGGTAACGGTGGAAGTCCCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTc                                                                             <  2:286194/143‑1 (MQ=255)
                                                                gCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGAAGTCGt                                                                           >  2:404284/1‑143 (MQ=255)
                                                                  ggTACAAGTGGAAGTCCCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTaa                                                                         <  1:934564/143‑1 (MQ=255)
                                                                  ggTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTaa                                                                         >  1:624129/1‑143 (MQ=255)
                                                                     aaaaGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATcc                                                                      >  2:591426/1‑143 (MQ=255)
                                                                        aGTGGACGCACCCGGAGAATTGTTTCGGTGACGAGGTGATCTGCACAACACCCTTATTGATATGCAGGCTCTGCGTCGCCTCTTCGACTTTCCACGGTTCGCCTTACTCCACCACCATCCGCAGTAAGATGTCGTAATCCtgt                                                                   <  1:351469/143‑1 (MQ=255)
                                                                             aaGTACCCGGAGAATTTTTTCGGTGAAGAAGTGATGTACATTTCACCGTGATTGATATGAAGGCTATACGGAGCCTCTTCTAATTTACCCGGTTCGCTGTTTCCGAGCAACATCCGCAAGGACAGGGCGTAATGCttggacgc                                                              >  2:377728/1‑136 (MQ=255)
                                                                              aGTCCCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCt                                                             <  1:778076/143‑1 (MQ=255)
                                                                                    cGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATTCAGGCTCTGCGTCGCCTTTTGTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCCGGAAGATGTCGTAATCCTGTGCCGCGTTCAGt                                                       >  1:83157/1‑143 (MQ=255)
                                                                                       agatctttttCCGTCACGAGGTCATCTGCTTTTCAGCCCGAACCATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCCCCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCg                                                    <  2:770601/138‑1 (MQ=255)
                                                                                        ggagtgtttCGGTGACGAGGTGATCTGCATTTCACCGTGGTTGATATGCAGGCTCTGCGTCGCCTCTTCGACTTTCCACGGTTCGCCGTACGCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCgg                                                   <  2:145315/137‑1 (MQ=255)
                                                                                            gttttCGGTGACGAGGTGATCTGCATTTCACCGTGACTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATc                                               <  2:927150/142‑1 (MQ=255)
                                                                                             ttttCGCTGACGAGGCGATCTGCAGTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCg                                              <  1:883743/143‑1 (MQ=255)
                                                                                                  ggTGACGACGTGAGCTGCACGCCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGaacaa                                         <  2:787390/143‑1 (MQ=255)
                                                                                                          ggTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAACATGTCGTAATCCTGTGCCGCTTTCAGTTCGGCATCGAACAAACCCTCtt                                 >  2:1064942/1‑143 (MQ=255)
                                                                                                                   cATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGcc                        <  1:621653/143‑1 (MQ=255)
                                                                                                                       tCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCCCCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTCTCGAACACACACTCTTTGCAACGGCATgg                    >  2:798807/1‑142 (MQ=255)
                                                                                                                           cGTGATTGATCTGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAg                <  2:755633/143‑1 (MQ=255)
                                                                                                                                ttgatATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCGGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTaaa           >  1:150536/1‑143 (MQ=255)
                                                                                                                                   tgtggCAGGTTCGGCGTCTCCTCGGCTACGTGCCACGGTTCGCCGTACTCCACCACCAGCCGCAGGAAGATGTTGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAACAACACTCTTTATCGCCATGCCCAGGTAAAGAc        <  2:390526/139‑1 (MQ=255)
                                                                                                                                   atatGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACGACGAACCGCAGGAGGATGTCGTAATCCTGTGCCGCTGTCAGTTCGGTATCGAGCAAACACTCTTTGAAACGCCGTGCCCAGGTAAAGAc        >  1:252438/1‑143 (MQ=255)
                                                                                                                                     atGCAGGCTCTCCATCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCTCCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGTCATGCCCAGGTAAAGACCt      <  2:60684/143‑1 (MQ=255)
                                                                                                                                        cAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCCGTTCGGTATCGAACAAACACTCTTGGAAACGCCCAGCCCCGGTACAGACCtcgg   >  1:528399/1‑140 (MQ=255)
                                                                                                                                         aGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCCGTTCGGTCTCGAACAAACACTCTTTTAAACGCCTTGCCCAGGTTACGCCATGGtt  >  1:844797/1‑143 (MQ=255)
                                                                                                                                         aGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGGtt  <  1:371999/143‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                           
GGATCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGGTT  >  NC_000913/2132647‑2132926

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: