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.6% 102.6 / NA 37I204S (ATC→AGC) wcaAputative glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base C (14/21);  minor base G (1/1);  total (15/22)
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                                                                                                                                  <  2:672151/141‑1 (MQ=255)
  gagGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCCGGCTCTGCGTCGCct                                                                                                                                >  1:402881/1‑143 (MQ=255)
  gagGGGAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGAGCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCct                                                                                                                                <  1:424153/143‑1 (MQ=255)
        caaCGGCTGATATTTTTTGCTCGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTa                                                                                                                          <  1:331797/142‑1 (MQ=255)
            aGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACttt                                                                                                                      >  1:250002/1‑143 (MQ=255)
               tGATCTTTTTTCCTGTCGCGGTCGCTTTTGTCTTTGTGCTTGCGGTAATAGTGGACGTACCCGGAGAATTTTTTCGGTGACGGGGTGATCTGCATTTCCCCGTGATTGTTATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCa                                                                                                                   <  2:61770/143‑1 (MQ=255)
                  tatTGTTTGCGTGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACgg                                                                                                                <  2:709395/143‑1 (MQ=255)
                       tgtGCTGGCGCGGTCGAGTGTGTCTTTGAGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGc                                                                                                           <  2:306606/141‑1 (MQ=255)
                          gCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTGTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGt                                                                                                        >  2:589085/1‑143 (MQ=255)
                            tGGCGCGGTCGAATTGGTCTTTGTGCTTGCGGTAAAAGTGGACGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTAc                                                                                                      <  2:264127/143‑1 (MQ=255)
                                     cGAACTGGTCTTTGGGCTTGCGGTAAAAGTGGAGGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTccaccacc                                                                                             <  2:295696/143‑1 (MQ=255)
                                        aTTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATGTCACCGTGATTGATATGCAGGGTCTGGGTCCCCTATTCTACTTAGCAAGGTTTGCGGTACTCCACAGCCATc                                                                                          >  2:432115/1‑143 (MQ=255)
                                           tGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGc                                                                                       <  1:209769/143‑1 (MQ=255)
                                                ttGTGCTTGCGGTAAAAGTGGAAGCACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGaa                                                                                  <  1:567201/143‑1 (MQ=255)
                                                      ggcaggTAAAGGTGGAAGTAGCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTc                                                                            <  2:626753/139‑1 (MQ=255)
                                                           gTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACCTTCCACGGTTCGCCGTACTCCCCCACCCTCCGCAGGCACACGTCGTaac                                                                       >  2:290439/1‑142 (MQ=255)
                                                             aaaaGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCCTCCGCAGGAAGATGTCGTAATcc                                                                     >  1:286542/1‑143 (MQ=255)
                                                                aGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCtgt                                                                  >  2:5631/1‑143 (MQ=255)
                                                                 gTGGAAGTACCCGGGTAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCTGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCtgtg                                                                 <  1:210075/143‑1 (MQ=255)
                                                                 gTGGAAGTACCCGGAGAATTTTTTCGATAACGAGGTGAACAGCATAGCACCGTGAATGATATGCAGGCTCTGCGTCGCCTGTTCTACTTTCCACGGTGCGCCGGACACGACCCCCATCCGCCGGGAGACGTCGTAGTcatctg                                                                 >  2:295311/1‑138 (MQ=255)
                                                                      aGTACCCGGAGAATTTTATCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCt                                                            <  1:650131/143‑1 (MQ=255)
                                                                              gagaATTTTGTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTGCTACTTTCCACAGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTc                                                    <  1:571964/143‑1 (MQ=255)
                                                                               agcagtgtttCGGTGACGGGGTGATCTGCATTTCACCGTGAGTGATCTGCAGGCTCCGCTCCGCCTCTTTTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCCTAATCCTGTGCCGCTTTCAGTTCg                                                   <  1:534885/136‑1 (MQ=255)
                                                                                 aaTTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGGCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGGTCTCAGCTCgga                                                 >  2:431684/1‑142 (MQ=255)
                                                                                    tttttCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGAAGGCGCGGCGTCGGCTCTTCTACTTTCCGCGGTTGGTCGTACTCTAGCACCATCCGCAGGAAGGTGTCGTAATGCGGTGCGCGTTACAGTGCGGtgta                                              >  2:273203/1‑140 (MQ=255)
                                                                                     ttttCGGTGACCAGGTGATCTGCATTCCACCGTGCTCGATATGCAGGCTCTGCGTCGCCTCTTCTCCTTTCCACGCTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCg                                             <  2:563542/143‑1 (MQ=255)
                                                                                         cGGTGACGGGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGaaca                                         <  2:121774/143‑1 (MQ=255)
                                                                                          ggTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGaccaa                                        >  1:597088/1‑143 (MQ=255)
                                                                                           gTGCCGAGGTGATCTGCTTTTCACCGGGATTGATAGGCCGGCTCTGCGTCGCCGCTTCGACTTTCCACGGTCCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGaacaaa                                       <  2:269995/143‑1 (MQ=255)
                                                                                             taCGAGCTTATCTGTATTTCAACGTGATGGAGATGCAGGGGCAGCGTCGCCTCTTCTACTGTCCACGGTGCGGCGTACTCCACCACCATCCGCAGGACGATGTCGTAATCCTGTTCCGCTTTCAGTTCGGTATCGaacaaaca                                     <  1:97961/142‑1 (MQ=255)
                                                                                             gACGCGGTGTTCTGCATTTCACCGTGATTGATGTGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGaacaaaca                                     <  2:196607/143‑1 (MQ=255)
                                                                                              ccGAGGTGCTCTGCATTTCGCCGTGGTTGATATGTAGGCTCTGCGTCGCCTCTTCTACGTTCCGCGGTGCGCCGTAATCCACCTCCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAacac                                    <  2:10247/142‑1 (MQ=255)
                                                                                                  ggTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATGCTGTGCCGCTTTCAGTGCGGCATCGAACAAACACTGtt                                >  1:716255/1‑143 (MQ=255)
                                                                                                         tGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACg                         >  1:496872/1‑143 (MQ=255)
                                                                                                                   cctgattgatATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACGCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAg               <  2:321991/141‑1 (MQ=255)
                                                                                                                            tatGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGAcc      >  2:415314/1‑143 (MQ=255)
                                                                                                                                cAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGCAACGCCCCGCCCAGGTAAAGCCCCGGt  >  2: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: