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

TCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTG  >  NC_000913/2132650‑2132923
                                                                                                                                         |                                                                                                                                        
tCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGGCGAAAGTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGc                                                                                                                                     <  1:498920/143‑1 (MQ=255)
   ggTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGGGGTCGAATTTGTCGTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTGGCCGTGATTGATATGCAGGCTCTGCGTc                                                                                                                                  >  2:149491/1‑143 (MQ=255)
                  gCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTc                                                                                                                   >  1:650892/1‑143 (MQ=255)
                  gCTGATATCTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTc                                                                                                                   <  1:269184/143‑1 (MQ=255)
                    tGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCa                                                                                                                 >  1:274820/1‑143 (MQ=255)
                      agtttttttGCTGGCTCGGTCGAATTTGTCTTTGTGCTTGCGTTACGCGTGTCAGTACCCGGCGACTTTTTGCGGTGAGGGGGTGACCTGCATTTCACCGTTATTGTTATGCAGGTTCTGCGTCGCCTACTCTCCTGTCCACg                                                                                                               <  2:226749/140‑1 (MQ=255)
                         ttttttGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATCTGCAGGCTCTGCGTCGACCCTTCTACTTTCCACgtgt                                                                                                            >  1:232851/1‑140 (MQ=255)
                                cTGTCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTa                                                                                                     <  1:588241/143‑1 (MQ=255)
                                cTGGCGCGGTCGGAGTTGTCTTTGTGCTGGCGGTAGAGGTGGAAGTGCCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGGTGTGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTa                                                                                                     <  2:95912/143‑1 (MQ=255)
                                     gcgGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTcca                                                                                                >  1:136093/1‑143 (MQ=255)
                                         tCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCGCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTccaccac                                                                                            >  2:601741/1‑143 (MQ=255)
                                          cGACTTTGTGTGTTTGCTTGCGGTAAAAGTGTAAGTCCCCGTTGAATTCTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTccaccacc                                                                                           <  2:428560/143‑1 (MQ=255)
                                              tttGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATcc                                                                                       >  2:776480/1‑143 (MQ=255)
                                                 gTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCa                                                                                    >  1:220049/1‑143 (MQ=255)
                                                     ttGTGCTTGCTTTATAAGTGGATGTTCCCGGAGTATTTTTTCGGTGTCGAGGTGATCTGCATTGCACCGGGATAGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGaa                                                                                <  1:530041/143‑1 (MQ=255)
                                                     ttGTGCTTGCGGTATATTTGGTATTTTCCGGAGTATTTTTTCGGTTTCGTGGTGATTTGCATTTCACCATGATTGATATGCAGGCTCTGCGGCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGaa                                                                                <  1:326323/143‑1 (MQ=255)
                                                          cTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCAGTACTCCAACACCATCCACAGGAAGATGt                                                                           >  2:865176/1‑143 (MQ=255)
                                                                      gTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCtgtg                                                               >  1:262792/1‑143 (MQ=255)
                                                                      gTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCtgtg                                                               >  1:710279/1‑143 (MQ=255)
                                                                      gTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCtgtg                                                               >  1:717292/1‑143 (MQ=255)
                                                                                    agaCTTTCTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCg                                                 <  1:273152/143‑1 (MQ=255)
                                                                                            ttCGGTGACGAGGTGATCTGCATTTCACCGTGAATGATATGCAGGCTCTGCGTCGCCTCTTCTACCTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGaa                                         <  1:715836/143‑1 (MQ=255)
                                                                                                    cgcggTCATCTGCATTTCACCGTGATTGATCTGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACAct                                 <  1:315935/140‑1 (MQ=255)
                                                                                                    ccaGGTGATCTGCATTTCCCCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACAct                                 <  1:256952/141‑1 (MQ=255)
                                                                                                      aGGTGAGCGGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACActct                               <  1:624718/143‑1 (MQ=255)
                                                                                                       ggTGATCTGCATTTCACCGTGTTTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCtt                              <  2:686216/143‑1 (MQ=255)
                                                                                                        gTGATCTGCCTTGCACCGTGATTGATGTGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCttt                             <  2:449965/143‑1 (MQ=255)
                                                                                                              tGCATTTCACCGTGATTGATATGCAGGCGCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGAGGGCGTAATGCTGTGCCGCCGTCAGTGCGGAATAGAACAAGCACTCTTTGcaaag                       >  2:850131/1‑141 (MQ=255)
                                                                                                                  gtgcACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTCTCCACGGTTCGCCGTGCTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCAt                   <  1:785611/140‑1 (MQ=255)
                                                                                                                      aCCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGccc               >  1:843539/1‑143 (MQ=255)
                                                                                                                        cgggattgatATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACCCCACCACCATCGGCAGGAAGAGGTCGTAATCACGTGCCGCTCTCAGTTCGATATCGACCGAACACTCTTTGAAACGCTATGCCAAg             >  2:316735/4‑143 (MQ=255)
                                                                                                                        cGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCac             >  2:596857/1‑142 (MQ=255)
                                                                                                                        cGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAg             >  2:486888/1‑143 (MQ=255)
                                                                                                                             ttgatATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTCCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTaaa        <  2:523797/143‑1 (MQ=255)
                                                                                                                              tgatATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCTCAGGAAGCTGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAACCCCCATGCCCCGGTAAAg       >  1:482095/1‑143 (MQ=255)
                                                                                                                              tgatATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACACTTTGAAACGCCATGCCCAGGTAAAg       >  2:597740/1‑143 (MQ=255)
                                                                                                                              tgacATGCGGGCTCTGCGTCGCCTATTCGAGTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAg       <  2:693784/143‑1 (MQ=255)
                                                                                                                                atatGCAGGCTCTGCGTCGCCACTTCTACTTGCCACGGGTCGCCGTACTCCACCGCCGTCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGAc     <  1:326112/143‑1 (MQ=255)
                                                                                                                                   tGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTg  <  1:601741/143‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                        
TCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTG  >  NC_000913/2132650‑2132923

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: