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.4% 113.4 / NA 39I204S (ATC→AGC) wcaAputative glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base C (16/22);  total (17/22)
Fisher's exact test for biased strand distribution p-value = 4.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGG  >  NC_000913/2132654‑2132924
                                                                                                                                     |                                                                                                                                         
gttgagGGTAAAGAGCGGATATTGTGTGGTGGGGCTGTCGAATGTGGCTTTGTGCTTGCGGTAAAAGTGGAAGTGCCCGGAGAATTTTGTCGGTGACGAGGTGATCTGCATTTGACCGTGATTGATATGCAGGCTCTGCGTCg                                                                                                                                  <  2:614405/140‑1 (MQ=255)
 tAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATGGATATGCAGGCTCTGCGTCGc                                                                                                                                 >  1:643432/1‑143 (MQ=255)
   gagGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGAAGAGGAGTGCTTCGTTTCACCGTGCTTGATATGCAGGCTCTGCGTCGCct                                                                                                                               >  1:676030/1‑143 (MQ=255)
      ggTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCGCTGCGTCGCCtctt                                                                                                                            >  2:529995/1‑143 (MQ=255)
       gTAAACCGCTGATATTTTGTGCTGGCGCGGTCGAATTTGTCTTTGCGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCtcttc                                                                                                                           <  1:72187/143‑1 (MQ=255)
             aGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGGACCCGGAGAGTTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCGGGCTCTGCGCCGCCTCCactactgc                                                                                                                     >  2:43710/1‑141 (MQ=255)
                 gATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGTTGATCTGCATTTCACCGTGATTGATATGCTGGATGTGCGGCGCGTGTTCTACTGTCCaa                                                                                                                 >  1:877197/1‑142 (MQ=255)
                     ttttttGCTGGCGCGGTCGAATTTGTATCTCAGCATTGGCTAAAATTTGAAGTACCCGGAGAATTTTTTCCGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTAATTTCCAAGGtt                                                                                                             >  1:640353/1‑143 (MQ=255)
                       tttgGCTGGCGCGGTCGAATTTGTCTTTGTGCGTGCGGTAAAAGTGGAAGTACCCGGAGAATGTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCg                                                                                                           <  2:687519/143‑1 (MQ=255)
                            cTGGCGCGGTCGAATTTGTCTCTGTGCATGCGGTAGAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTa                                                                                                      <  1:249677/143‑1 (MQ=255)
                             tGGCTCGGTCGCCGTGGTCTTCGTGGTTGCGGTTCAAGTGGAAGTACCCGGAGAATGTTTTCGGTGACGAGGTGATCTGCCTTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTGCGCCGTAc                                                                                                     <  2:24128/143‑1 (MQ=255)
                                     tCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTccaccac                                                                                             >  1:483333/1‑143 (MQ=255)
                                              tCTGTGTGCGTGCGGTAAAAGTGGAAGTACCCGTAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAg                                                                                    <  1:408768/143‑1 (MQ=255)
                                              gtctatagCTTGCGGTAAAAGTGGAAGTACCCGGTGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAg                                                                                    <  1:26178/136‑1 (MQ=255)
                                                tttGTGCTTACGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTGGCCTCTTCTACTTTACACGGTTCGCCGTACTCCACCACCATCCGCAgca                                                                                  >  2:69901/1‑141 (MQ=255)
                                                 ttGTGCTTGCGGTAAAAGGGGAGGTACCCGGAGAATGTGTGCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGaa                                                                                 <  1:338654/143‑1 (MQ=255)
                                                 ttGTGCTTACTGCAAAAGTTGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTAATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGaa                                                                                 <  1:764357/143‑1 (MQ=255)
                                                                  gTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCtgtg                                                                <  1:714714/143‑1 (MQ=255)
                                                                           cccGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCa                                                       <  1:902095/143‑1 (MQ=255)
                                                                            ccGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAg                                                      >  1:584722/1‑143 (MQ=255)
                                                                                cgaACTTTTTCGGCGCCGAGGTGATCTGCGTCTCACCGTGATTGATACGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCACCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCg                                                  <  1:691647/142‑1 (MQ=255)
                                                                                        ttCGGTGACGAGGTGATCTGCACTTCGCCGTGATTGATATGCAGGCGCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGaa                                          <  2:244595/143‑1 (MQ=255)
                                                                                            gTGACGAGGCGATCTGCATTTCACCCTGATTGATCTGCACGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTCCTCCACCCCCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGaacaaa                                      <  1:493091/143‑1 (MQ=255)
                                                                                             tGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCa                                                       >  1:438287/1‑125 (MQ=255)
                                                                                             tGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCa                                                       <  2:438287/125‑1 (MQ=255)
                                                                                                       aTCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGCATCGAACAAACACTCTTTGaa                           >  2:337536/1‑143 (MQ=255)
                                                                                                           gCATTTCACCGTGATTGATATGCAGGCTCTGCGTAGCCTCTTCTACTTTCCACGGATCGCCATACTCCACCACCATCCGCAGGAAGATGTAGTACACCTGTGCCGCTTTCAGTTCGGGCTGGAACGAACACTGTTTgcaccct                       >  1:206691/1‑137 (MQ=255)
                                                                                                              tttCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTGCTCTACCACCATCCGCAGGAAGCTGTCGTAATCCTGTGCCGCTTTCCGTTCGGTATCGAAAAAACCCTCTTTGAAACGCCAt                    >  2:155080/1‑143 (MQ=255)
                                                                                                              cttCACTGTGCTTGTCATGCAGTCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCTGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCAt                    <  2:737424/142‑1 (MQ=255)
                                                                                                                  tccGTGGTTGTTACGCGGGCTCTTCCTCGCCTCTTCTTCTTTCCACGTTTCGCCGTACTCCTCCACCATCCGCAGGAAGATGTCTGCTTTCTTTTCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGccc                <  2:230140/142‑1 (MQ=255)
                                                                                                                   ccGTTATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCa               <  2:851760/143‑1 (MQ=255)
                                                                                                                   ccGTCATTGATATGCAGGCTCTGCGTCGCCCCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCa               <  1:403427/143‑1 (MQ=255)
                                                                                                                    cGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCTTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAg              <  1:484424/143‑1 (MQ=255)
                                                                                                                    cGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTAATCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAg              <  1:834414/143‑1 (MQ=255)
                                                                                                                     gTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAgg             <  1:583783/143‑1 (MQ=255)
                                                                                                                             tatGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAACCGCCATGCCCAGGTACAGAcc     >  1:610685/1‑143 (MQ=255)
                                                                                                                              atGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCt    >  1:818520/1‑143 (MQ=255)
                                                                                                                               tGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTCCTCCTGTGCCGCTTTCCGTTCGGTATCGACCACACACTCTTTGAGACGCGATGCCCCGGTAAAGAcccg   >  2:663360/1‑141 (MQ=255)
                                                                                                                                gCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAGAGACCTgg  >  1:647083/1‑143 (MQ=255)
                                                                                                                                     |                                                                                                                                         
GTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGG  >  NC_000913/2132654‑2132924

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: