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,7870AC95.0% 95.1 / NA 40I204S (ATC→AGC) wcaAputative glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base A (2/0);  new base C (18/20);  total (20/20)
Fisher's exact test for biased strand distribution p-value = 4.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.10e-01

TAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCT  >  NC_000913/2132655‑2132922
                                                                                                                                    |                                                                                                                                       
tAGCGGGTAAACAGCTGATATTGTGTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGTAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGc                                                                                                                               <  1:294926/143‑1 (MQ=255)
tAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGCGATCTGCACTTCATCCAGCCTGACATGCAGCCTCTGCGACGc                                                                                                                               >  1:603717/1‑143 (MQ=255)
     ggTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAGGTGGAAGTACCCGGAGCCTTTTTTCCGTGCCGTGGTGATCTACCTTACACAGAAACTGACATACAAGATCtacataaacactt                                                                                                                          >  2:5937/1‑131 (MQ=255)
            aGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGCTTGATGTGCAGGCTCTGCGTCGCCTCTTCTACttt                                                                                                                   >  2:506801/1‑143 (MQ=255)
            aGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGAAAGAACTGCGTCGCCACTTATACttg                                                                                                                   >  2:422252/1‑142 (MQ=255)
                 atatTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATACGCCGGCTCTGCGTCGCCTCGTCTCCTTTCCACg                                                                                                              >  2:634554/1‑143 (MQ=255)
                      ttttGCTGGCGTGGTTGTATTTGTCTTTGTGCTTGTGGTAAAAGTGGAAGTACCCGGAGAAGGTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCg                                                                                                         <  2:345141/143‑1 (MQ=255)
                       tttGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGACAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTACGTGGCCTCTTCTACTATCGACggaccga                                                                                                        >  1:393374/1‑138 (MQ=255)
                            tGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTAc                                                                                                   <  2:260008/143‑1 (MQ=255)
                              gcgcTGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTTATCTGCATTTCACCGTGATTGATATGCAGGCTCTGTGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTc                                                                                                 <  2:582421/143‑1 (MQ=255)
                                   gTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTAc                                                                                                   >  2:612339/1‑136 (MQ=255)
                                   gTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTAc                                                                                                   <  1:612339/136‑1 (MQ=255)
                                            ttCTTTGTGCTTGCGGTAGAAGCGGAAGTACCCGGAGGAGCTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCa                                                                                   <  2:43355/142‑1 (MQ=255)
                                              cTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACCGTTCGCCCGACTccaccacccacccccgg                                                                                 >  1:554816/1‑134 (MQ=255)
                                               tttGTGCTTGCGGTAAAAGTGGAAGTACCCGTAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGa                                                                                >  1:263991/1‑143 (MQ=255)
                                                  gtgCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCACGGTCCTCTCCCACCATCCGCCggcaca                                                                             >  1:111962/1‑139 (MQ=255)
                                                      ttGCGGTAAAATTGTAAGTACTCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTc                                                                         <  2:503417/143‑1 (MQ=255)
                                                      ttGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTc                                                                         >  2:248305/1‑143 (MQ=255)
                                                      ttGCGGAAAAAGAAGAAGGAACCGGAAAATAGTTTCGGAGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTc                                                                         <  2:157798/143‑1 (MQ=255)
                                                        gCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCCTTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACAATCCGCAGGAAGCTGTCGt                                                                       >  2:545392/1‑143 (MQ=255)
                                                           gTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTCCTCCACCACCCTCCGCACGACGATGTCGGAAt                                                                    >  2:732121/1‑143 (MQ=255)
                                                                aGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTGCTCCACCACCATCCGCAGGCAGATGAAGTAACCCcgt                                                               >  1:650640/1‑143 (MQ=255)
                                                                aGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCtgt                                                               >  2:310175/1‑143 (MQ=255)
                                                                    taaGTACCCGGAGAGTTTTTTCGGTGACTAGGTGATCTTCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTGCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCg                                                           <  1:670228/142‑1 (MQ=255)
                                                                      cgttcccGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCGTCTACTTTCCACGGTTCGCCGTACTCCACCACCTTCCGCAGGAAGATGTCGTAATCCTGTGCCGCt                                                         <  1:491273/139‑1 (MQ=255)
                                                                                gaATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCAGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACAATCCGCAGGAAGCTGTCGTAATCATGTACCGcatcaagctcga                                               >  1:171011/1‑132 (MQ=255)
                                                                                    tttttCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATc                                           <  1:603678/143‑1 (MQ=255)
                                                                                     ttttCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCGCGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCg                                          <  1:310175/143‑1 (MQ=255)
                                                                                     ttttCGGAGACGAAGTGATCTGCATCTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCg                                          <  2:442070/143‑1 (MQ=255)
                                                                                                  ggTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTCCGTCGCCTCTTCTTCTGTCGACGGTTCGCCGTACTCCACCGCCATCCCCTGGAAGATGTCGTATTCCTTTGCCGCTTTCAGTTCGGTATCGAACAAACACTCtt                             <  1:346220/143‑1 (MQ=255)
                                                                                                      aTCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTNATCCTGTGCCGCTTTCAGTTCGGCATCTCACAAacacacaaagaa                         >  1:383161/1‑135 (MQ=255)
                                                                                                        cTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCCCCATCCGCAGCAAGATGTCGTAATCCTCTGCCGCGTTTGGTTCGGTATCGAGCAAGCACTCGTTGagca                       >  1:146130/1‑140 (MQ=255)
                                                                                                         tGCATTTCACCGTGATTGAGATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACg                      <  1:288206/143‑1 (MQ=255)
                                                                                                           gttctcgtcGTGATTGATATGCAGGCTCTGCGGCTCATCTTCTACTGTCCGCGGTTCGCCGTACTCCACCGCCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGcc                    <  1:570959/135‑1 (MQ=255)
                                                                                                             tttCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTCCTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTCTCAGTGCGGCATCGAACCAACACTCCTTGAGACGCCAt                  >  1:118856/1‑143 (MQ=255)
                                                                                                                   cgcgattggtgcgCAGGCTCTGCGTCTCCTCCTTTACTTTCTACGGTTCGCCGTGCTCCACCACCATCCGTAGGAATATGTCGTAATCCTGTGCCGCTTTCCGTTCGGTATCGAACAAACCCTCTTTGAAACGCCATGCCCAg            <  1:175502/131‑1 (MQ=255)
                                                                                                                   cGTGATTGATTTGCATGCTTCGCGTCTCCTCTTCTACTTTCCACGGTTCGTCGTGCTCCACCACGATCCGCACGAAGATTCCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAg            <  1:308173/143‑1 (MQ=255)
                                                                                                                        ttgacATGCAGGCTCTTCGTCTCTTCTTATACTGTCCACGTTTCGCCTTACTCCACCACTATCCTCATGAAGATTTCGTAATCCTGTGCCGCTGTCAGTTCGGTATTGAACAAACACTTTTTGAAACGCCATGCCCAGGTaaa       <  1:731206/143‑1 (MQ=255)
                                                                                                                           atatGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGCCTATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTAACGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGAc    >  2:225431/1‑143 (MQ=255)
                                                                                                                             attcAGGCTTTGCGTCGCCTCTTCTCCTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCt  <  1:248305/140‑1 (MQ=255)
                                                                                                                                    |                                                                                                                                       
TAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCT  >  NC_000913/2132655‑2132922

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: