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.9% 138.0 / NA 48I204S (ATC→AGC) wcaAputative glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base C (23/24);  minor base T (1/0);  total (24/24)
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.99e-01

TCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGGTT  >  NC_000913/2132650‑2132926
                                                                                                                                         |                                                                                                                                           
tCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCCGCATTTCTCCGTGATTGATATGCAGGCTCTGc                                                                                                                                        >  2:137394/1‑143 (MQ=255)
    ggagacGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCg                                                                                                                                    <  1:835946/141‑1 (MQ=255)
    gTAGAGGGTAAACAGCTGATATTTTTTTTTTTATTTTTCGAATTTGTCTTTGTGCTTGCGTTTAAAGTGGTAGTACCCGGAGAATTTTTTCGGTGACGATGTGATCTGCGTTTCACCGTGATTGATATGCTGGCTTTGCGTCg                                                                                                                                    >  1:529245/1‑143 (MQ=255)
     tAGAGGGTCAACAGCTGATATTCTTTGCTGGCGCGGTCGAATTTGCCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGc                                                                                                                                   <  2:670896/143‑1 (MQ=255)
        agGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCctc                                                                                                                                >  1:605011/1‑143 (MQ=255)
         gggTAAACAGCTGATATTTTTTGCTGGCGCGGTAGAATTTGTCTATGTGCTTGCGGTAAAAGTGGAAGTACCCGAAGAATTTTTTCGGTGACGAGGTGATCTGCAATTCACCGTGATTGATATGAAGACTCTGCGTCgaatag                                                                                                                               >  1:98095/1‑138 (MQ=255)
                cAGCTGATATTGTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACtt                                                                                                                        <  2:31595/143‑1 (MQ=255)
                    tGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCa                                                                                                                    <  2:873011/143‑1 (MQ=255)
                    tGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCCCCGCGGCGCCTCTTGTGCTGTCCa                                                                                                                    >  2:900562/1‑143 (MQ=255)
                      atacgttttCCCGGCCCGGTCGCATCTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGACTTCTTTCGGGGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACg                                                                                                                  <  2:870883/138‑1 (MQ=255)
                         ttttttGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGtt                                                                                                               <  2:333657/143‑1 (MQ=255)
                         ttttctGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGtt                                                                                                               <  2:340535/143‑1 (MQ=255)
                         agtgttGCTGGCGCGGTCGAATTTGTATTTGTGCTTCCGGTAAAAGTGTACGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGtt                                                                                                               <  1:620385/139‑1 (MQ=255)
                            ttggcgggCGCGGTCGAATTTGTCTTTGTGCTTGCGGGAAACGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGc                                                                                                            <  1:971234/137‑1 (MQ=255)
                                cTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTa                                                                                                        <  2:605011/143‑1 (MQ=255)
                                   gcgcGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGCTTCGCCgcccgc                                                                                                     >  1:120174/1‑138 (MQ=255)
                                       ggTTGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTAACCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTccacc                                                                                                 <  2:475842/143‑1 (MQ=255)
                                          cgtgtttGGCTTTGTGGTTGCGGTAAAAGTGGGGGTCCCCGGAGAATTTTTTCGGGGACGCGGTGATCTGCATTTCACCGTGATTGTTATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTccaccacc                                                                                              <  1:556853/139‑1 (MQ=255)
                                               ttGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCAGCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTCCTCCACCACCATCCg                                                                                         >  1:338811/1‑143 (MQ=255)
                                               ttGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCg                                                                                         >  2:766256/1‑143 (MQ=255)
                                               ttGTCTTTGTCCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCg                                                                                         <  2:536065/143‑1 (MQ=255)
                                                tGTCTGTGTTGTTGCGGTAAAGGTCCTAGTATCCGGAGAATTTTTTCGGTGACGAGGTTATCTGCATTTCACCCTGCACGTTATTCATGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGc                                                                                        <  1:404731/143‑1 (MQ=255)
                                                     ttGTGCTTGCGGTAAAAGTGGAAGCACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGaa                                                                                   <  2:943944/143‑1 (MQ=255)
                                                         gCTTGCGGTAAAAGTGGTAGTACCCGTAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATg                                                                               <  1:515025/143‑1 (MQ=255)
                                                             gCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTTTTTTTTTTTTTAATCTGAATTTGACCGTGATTTTTATGCTGGCTCTGCGTCGCCTCTCCTACTTTCTTCCGTTCTCCGTTCTCCACTACCATCTGCCGGAAGATGTctt                                                                           >  1:409603/1‑141 (MQ=255)
                                                                gTAAAAGTGGAAGTACCCGGAGAATTTTTTTTTTGTTTATCTAATCTGCATCTCACTGCGATTGTTTTGCTTGTTTTTCGTCGCCTCTTCTTCTTTCTTCCGTTTCTCGTACTCCTGCACCACCCGCAGGatcttttcttttt                                                                        >  1:771382/1‑131 (MQ=255)
                                                                       tGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGc                                                                 >  2:992945/1‑143 (MQ=255)
                                                                             tACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCttt                                                           <  1:845448/143‑1 (MQ=255)
                                                                                ccGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCCGGCAGATGTCCTAATCCCGTGCCGCTTTccg                                                        >  2:831873/1‑141 (MQ=255)
                                                                                caggAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTACCACGATACGCCGTACTGCACCACCTTCACCAAGCACATGTCGTAAGCCTCTGACgacatcag                                                        >  1:716078/3‑136 (MQ=255)
                                                                                        ttttttCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACGTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGACGAAGACTTATTCCTACGCCGCTTTCGGTGCGGtgg                                                >  1:8497/1‑141 (MQ=255)
                                                                                        ttttttCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTACTTCGCCTATTAAAATATCAACAATTCGACGTACTCCACCACCATCCGCAGGAAGATGTCATAATCCTGTGCCGCTTTCAGTTCGGTAt                                                >  2:438069/1‑143 (MQ=255)
                                                                                         ttttgCGGTGGCGACGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATc                                               <  2:220640/143‑1 (MQ=255)
                                                                                             tCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCACCTCTTCTACTTTCCACGGTTCGCCGTACTACACCACCATCCGCAGGAAGATGCCGTAAGCCTGTGCCGCTTTCAGTTCGGTGTCggga                                           >  2:683755/1‑140 (MQ=255)
                                                                                              cGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGAGTCGaaca                                          >  2:223728/1‑143 (MQ=255)
                                                                                                gTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGCTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGaacaaa                                        >  2:767972/1‑143 (MQ=255)
                                                                                                  gACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTGCACCACCATCCGCAGGACGATGTGATAATCCTGTCCCGCTTTCAGTTCTCGAGCGaacaaaca                                      >  1:1014853/1‑143 (MQ=255)
                                                                                                        gTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCCGTTCGGCATCGaacacacccccctt                                >  1:751921/1‑136 (MQ=255)
                                                                                                        gTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCttt                                <  2:362659/143‑1 (MQ=255)
                                                                                                           ttCTGCATTTCACCGTGGTTGCCATGCAGGCTCTGCGTCGCCTCTTCTCCTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGaa                             <  1:53907/142‑1 (MQ=255)
                                                                                                            tCTTCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGaaa                            <  2:658315/143‑1 (MQ=255)
                                                                                                               gCATTTCCCGGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGc                         <  2:120174/143‑1 (MQ=255)
                                                                                                                  tttCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGCAACGCCAt                      >  2:447476/1‑143 (MQ=255)
                                                                                                                   ttCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTGGCAACGCCATg                     >  2:623311/1‑143 (MQ=255)
                                                                                                                             ttgatATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTACCGAACAAACACTCTTTGAAACGCCATGCCCAGGTaaa           <  2:80310/143‑1 (MQ=255)
                                                                                                                              tgatATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAACCGCCATGCCCAGGTAAAg          >  1:288924/1‑143 (MQ=255)
                                                                                                                                 tatGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGAcc       <  2:529245/143‑1 (MQ=255)
                                                                                                                                      aGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGACGCTGTCGTACTCCTGTGCCCCTTTCCGCTCGGTATCGACCACACACTCTTTGAACCCCCATACCCCAGTCAAGACATCGtt  >  1:289781/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: ‑