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.3% 96.5 / NA 37I204S (ATC→AGC) wcaAputative glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base C (18/18);  total (18/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.76e-01

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