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

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