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.5% 108.1 / NA 40I204S (ATC→AGC) wcaAputative glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base C (19/20);  minor base T (0/1);  total (19/21)
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.68e-01

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