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% 98.5 / NA 38I204S (ATC→AGC) wcaAputative glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base C (19/18);  total (19/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 = 9.77e-01

TCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGGTT  >  NC_000913/2132650‑2132926
                                                                                                                                         |                                                                                                                                           
tCTGGGAGAGGGTATTCGGCTGAGGTTTTTGGGTGGCGCGGTCTAATTTGTGTTTGTGCTGGCGGTAAAAGTGGAGTTTCCCGTAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGc                                                                                                                                        <  1:258507/143‑1 (MQ=255)
      agagGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGCGATGCATGTCACGGTGATTGCTATGAAGGCTGTGCGTCGcc                                                                                                                                  >  1:196432/1‑143 (MQ=255)
            tAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTGGCCtcttcg                                                                                                                            >  2:169101/1‑142 (MQ=255)
                 aGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCGTCGTCTAGttt                                                                                                                       >  2:367327/1‑143 (MQ=255)
                 aGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACttt                                                                                                                       >  2:495299/1‑143 (MQ=255)
                    tcatagTTTTTGCTGGCGCGGTCGAATTTGTCGTTGTGCGTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCa                                                                                                                    <  1:751095/141‑1 (MQ=255)
                    tGATATTTTTTGACGGAGCGGTCGAAGGTGTTTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCa                                                                                                                    <  1:349016/143‑1 (MQ=255)
                          tttttGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGCTTGATATGGAGGCTCTGCGGCGCCTCTTCTACTTTCCACGGTTc                                                                                                              >  2:645002/1‑143 (MQ=255)
                          tttttGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTAGTTTCCACGGttt                                                                                                              >  2:892148/1‑142 (MQ=255)
                           tgttGCTCGCGCGCTCAAATTTGTCTTTGTGCTTGCGTTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTAATCTGCATTTCACCGTGATTGATATGCAGGCTGTGCGTCGCCTCTTCTACTTTCCACGGTTCg                                                                                                             <  2:951419/141‑1 (MQ=255)
                                cTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCCTCGCCTCTTCTACCTTCCACGGCTCACCGTa                                                                                                        >  1:955699/1‑143 (MQ=255)
                                  gggggggTCGAATTTGGCTTTGGGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGGGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACt                                                                                                      <  2:735768/138‑1 (MQ=255)
                                   gcgcGGTCGAATTTGTCTTAGTGCTTGAGGTCAAAGTGGAAGTACCCGGAGGATTTTTTCGGTGACGAGGTGATATGCATTTCACCGTGATTGCTATGGAAGCTCTACGACGCATCTTGTACTTTCCACGGTTCACCGTAcac                                                                                                     >  1:437366/1‑141 (MQ=255)
                                      cgGTCGAATTGTGCTTTGTGCTGGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGTTCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTccac                                                                                                  <  1:830795/143‑1 (MQ=255)
                                       cgTCAAAGTTGTCTTTGTGCTTGCGGTAAAATTGCAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTccacc                                                                                                 <  1:694804/142‑1 (MQ=255)
                                        gTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTccacca                                                                                                <  1:495299/143‑1 (MQ=255)
                                           gAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTGCGCAGGACTccaccacca                                                                                             >  1:20909/1‑143 (MQ=255)
                                                          cTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCCCGGGAAGgtgt                                                                              >  2:954611/1‑143 (MQ=255)
                                                            tGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCg                                                                            <  2:901461/143‑1 (MQ=255)
                                                               ggTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGCGCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCACCGTACTCCACCACCATCCgcaggcagccgccggaa                                                                         >  1:189645/1‑134 (MQ=255)
                                                                   aaaGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGCCGCCGCTTCTACTTTTCGCGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTCATCCt                                                                     >  2:849968/1‑143 (MQ=255)
                                                                                 cGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTTGGAATCCTGTGCCGCTTTCAGt                                                       >  2:122014/1‑143 (MQ=255)
                                                                                 cGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGt                                                       >  2:44992/1‑143 (MQ=255)
                                                                                  ggggaATTTTTGCGGTGACGAGGGGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGtt                                                      <  1:122014/140‑1 (MQ=255)
                                                                                  ggAGACTTTTTGCGGTGACGAGGTGATCTGCATGTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGtt                                                      <  1:504596/143‑1 (MQ=255)
                                                                                      aaTTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGt                                                  <  2:889775/143‑1 (MQ=255)
                                                                                          ttttCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGCCGGAATCCTGTGCCGCTTTCAGTTCGGCATCg                                              >  1:695457/1‑143 (MQ=255)
                                                                                             tctgggACGGGGTGGTCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCGCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGaac                                           <  1:331768/138‑1 (MQ=255)
                                                                                               ggTGACGTGGTGATCTGCAGTTCAGCGTGATGGTTATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCTTCCGCAGGAGGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGaacaa                                         <  1:742352/143‑1 (MQ=255)
                                                                                                gTGACGAGGTGATCTGCATTTCACCGTGATTGATAAGAAAAATGTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGaacaaa                                        <  2:792618/143‑1 (MQ=255)
                                                                                                   aCGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGCTATCGAACAAacac                                     >  1:895584/1‑143 (MQ=255)
                                                                                                      aGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACActct                                  <  2:817104/143‑1 (MQ=255)
                                                                                                        gTGACCTTCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCttt                                <  1:849968/143‑1 (MQ=255)
                                                                                                             cTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAAc                           >  1:217170/1‑143 (MQ=255)
                                                                                                                     cgccGTGATTGATCTGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGcc                   <  1:313169/141‑1 (MQ=255)
                                                                                                                             ttgatATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACTGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCCTAATCCTGTGTCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTaaa           <  1:257175/143‑1 (MQ=255)
                                                                                                                                atatGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGACACGCCCTGCCCAGGTAAAGAc        >  2:659645/1‑143 (MQ=255)
                                                                                                                                      aGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTCTCGCACAAACACTCTTTGAAACGCCATGCCCAGGTACTGACCTCGtt  >  2:469742/1‑143 (MQ=255)
                                                                                                                                         |                                                                                                                                           
TCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGGTT  >  NC_000913/2132650‑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: