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,7870AC94.9% 106.6 / NA 39I204S (ATC→AGC) wcaAputative glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base C (14/23);  minor base G (1/1);  total (15/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 = 1.00e+00

TAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGGT  >  NC_000913/2132655‑2132925
                                                                                                                                    |                                                                                                                                          
tggagGGTAAACAGCTGATAGTGTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGc                                                                                                                                  <  1:672151/141‑1 (MQ=255)
  gagGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCCGGCTCTGCGTCGCct                                                                                                                                >  2:402881/1‑143 (MQ=255)
  gagGGGAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGAGCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCct                                                                                                                                <  2:424153/143‑1 (MQ=255)
        taacatcgtctatTCTTTGTTGGCGCGGTCTATGTTGTTTTTGTGCTTGCGGTAAGGGTGGAAGTACCCGGAGCATTGTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTa                                                                                                                          <  1:49916/133‑1 (MQ=255)
        caaCGGCTGATATTTTTTGCTCGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTa                                                                                                                          <  2:331797/142‑1 (MQ=255)
            aGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACttt                                                                                                                      >  2:250002/1‑143 (MQ=255)
               tGATCTTTTTTCCTGTCGCGGTCGCTTTTGTCTTTGTGCTTGCGGTAATAGTGGACGTACCCGGAGAATTTTTTCGGTGACGGGGTGATCTGCATTTCCCCGTGATTGTTATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCa                                                                                                                   <  1:61770/143‑1 (MQ=255)
                  tatTGTTTGCGTGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACgg                                                                                                                <  1:709395/143‑1 (MQ=255)
                       tgtGCTGGCGCGGTCGAGTGTGTCTTTGAGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGc                                                                                                           <  1:306606/141‑1 (MQ=255)
                          gCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTGTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGt                                                                                                        >  1:589085/1‑143 (MQ=255)
                            tGGCGCGGTCGAATTGGTCTTTGTGCTTGCGGTAAAAGTGGACGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTAc                                                                                                      <  1:264127/143‑1 (MQ=255)
                                     cGAACTGGTCTTTGGGCTTGCGGTAAAAGTGGAGGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTccaccacc                                                                                             <  1:295696/143‑1 (MQ=255)
                                        aTTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATGTCACCGTGATTGATATGCAGGGTCTGGGTCCCCTATTCTACTTAGCAAGGTTTGCGGTACTCCACAGCCATc                                                                                          >  1:432115/1‑143 (MQ=255)
                                           tGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGc                                                                                       <  2:209769/143‑1 (MQ=255)
                                                ttGTGCTTGCGGTAAAAGTGGAAGCACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGaa                                                                                  <  2:567201/143‑1 (MQ=255)
                                                      ggcaggTAAAGGTGGAAGTAGCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTc                                                                            <  1:626753/139‑1 (MQ=255)
                                                           gTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACCTTCCACGGTTCGCCGTACTCCCCCACCCTCCGCAGGCACACGTCGTaac                                                                       >  1:290439/1‑142 (MQ=255)
                                                             aaaaGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCCTCCGCAGGAAGATGTCGTAATcc                                                                     >  2:286542/1‑143 (MQ=255)
                                                                aGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCtgt                                                                  >  1:5631/1‑143 (MQ=255)
                                                                 gTGGAAGTACCCGGGTAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCTGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCtgtg                                                                 <  2:210075/143‑1 (MQ=255)
                                                                 gTGGAAGTACCCGGAGAATTTTTTCGATAACGAGGTGAACAGCATAGCACCGTGAATGATATGCAGGCTCTGCGTCGCCTGTTCTACTTTCCACGGTGCGCCGGACACGACCCCCATCCGCCGGGAGACGTCGTAGTcatctg                                                                 >  1:295311/1‑138 (MQ=255)
                                                                      aGTACCCGGAGAATTTTATCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCt                                                            <  2:650131/143‑1 (MQ=255)
                                                                              gagaATTTTGTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTGCTACTTTCCACAGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTc                                                    <  2:571964/143‑1 (MQ=255)
                                                                               agcagtgtttCGGTGACGGGGTGATCTGCATTTCACCGTGAGTGATCTGCAGGCTCCGCTCCGCCTCTTTTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCCTAATCCTGTGCCGCTTTCAGTTCg                                                   <  2:534885/136‑1 (MQ=255)
                                                                                 aaTTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGGCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGGTCTCAGCTCgga                                                 >  1:431684/1‑142 (MQ=255)
                                                                                    tttttCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGAAGGCGCGGCGTCGGCTCTTCTACTTTCCGCGGTTGGTCGTACTCTAGCACCATCCGCAGGAAGGTGTCGTAATGCGGTGCGCGTTACAGTGCGGtgta                                              >  1:273203/1‑140 (MQ=255)
                                                                                     ttttCGGTGACCAGGTGATCTGCATTCCACCGTGCTCGATATGCAGGCTCTGCGTCGCCTCTTCTCCTTTCCACGCTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCg                                             <  1:563542/143‑1 (MQ=255)
                                                                                         cGGTGACGGGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGaaca                                         <  1:121774/143‑1 (MQ=255)
                                                                                          ggTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGaccaa                                        >  2:597088/1‑143 (MQ=255)
                                                                                           gTGCCGAGGTGATCTGCTTTTCACCGGGATTGATAGGCCGGCTCTGCGTCGCCGCTTCGACTTTCCACGGTCCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGaacaaa                                       <  1:269995/143‑1 (MQ=255)
                                                                                             taCGAGCTTATCTGTATTTCAACGTGATGGAGATGCAGGGGCAGCGTCGCCTCTTCTACTGTCCACGGTGCGGCGTACTCCACCACCATCCGCAGGACGATGTCGTAATCCTGTTCCGCTTTCAGTTCGGTATCGaacaaaca                                     <  2:97961/142‑1 (MQ=255)
                                                                                             gACGCGGTGTTCTGCATTTCACCGTGATTGATGTGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGaacaaaca                                     <  1:196607/143‑1 (MQ=255)
                                                                                              ccGAGGTGCTCTGCATTTCGCCGTGGTTGATATGTAGGCTCTGCGTCGCCTCTTCTACGTTCCGCGGTGCGCCGTAATCCACCTCCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAacac                                    <  1:10247/142‑1 (MQ=255)
                                                                                                  ggTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATGCTGTGCCGCTTTCAGTGCGGCATCGAACAAACACTGtt                                >  2:716255/1‑143 (MQ=255)
                                                                                                         tGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACg                         >  2:496872/1‑143 (MQ=255)
                                                                                                                   cctgattgatATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACGCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAg               <  1:321991/141‑1 (MQ=255)
                                                                                                                    gcggctgaggtGCAGGCTCCGCGTCGCCGCTTCTACGTCCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAgg              <  1:662599/133‑1 (MQ=255)
                                                                                                                            tatGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGAcc      >  1:415314/1‑143 (MQ=255)
                                                                                                                                cAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGCAACGCCCCGCCCAGGTAAAGCCCCGGt  >  1:75157/1‑143 (MQ=255)
                                                                                                                                    |                                                                                                                                          
TAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGGT  >  NC_000913/2132655‑2132925

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: