Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,132,787 A→C I204S (ATC→AGC)  wcaA ← putative colanic acid biosynthesis glycosyl transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,132,7870AC97.1% 95.5 / NA 34I204S (ATC→AGC) wcaAputative colanic acid biosynthesis glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base C (15/18);  total (16/18)
Fisher's exact test for biased strand distribution p-value = 4.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGATCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGGTT  >  NC_000913/2132647‑2132926
                                                                                                                                            |                                                                                                                                           
ggATCTGGTAGAGGGTAAACAACTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTAATTGCTGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGATCACGAGATGATATGCATATCACGATGATAGATAAGCAGCAta                                                                                                                                           >  2:383245/1‑142 (MQ=255)
      gggagagGGTAAACAGCGGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTc                                                                                                                                     <  1:213409/140‑1 (MQ=255)
            gggTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGGTATGCAGGCTCTGCGTCGCctct                                                                                                                               >  1:14495/1‑143 (MQ=255)
                gaaCCGCTGGTATTTTTTGCTGGCGCGGTCGAATTTGTCGTGGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTa                                                                                                                           <  2:53589/142‑1 (MQ=255)
                     gCTGACATTTTGTGCTGGCGCGGGCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTc                                                                                                                      <  2:276791/143‑1 (MQ=255)
                      cTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCCTTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTcc                                                                                                                     >  1:112383/1‑143 (MQ=255)
                                   cTGGCGCGGTCGAATTTGTCTTTGTACTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCCTCTACTTTCCACGGTTCGCCGtg                                                                                                        >  2:676327/1‑142 (MQ=255)
                                    tGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCATAc                                                                                                       >  2:245455/1‑143 (MQ=255)
                                    tGGCGCGGGCGAATTTGTCTTCGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGGGATTGATATGCAGGCTCTGCGTCGCCGCGTCTACTTGCTCCGGTTCGCCGTAc                                                                                                       >  2:471914/1‑143 (MQ=255)
                                        gcgGGCAAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTcca                                                                                                   <  1:440218/143‑1 (MQ=255)
                                         cgGTCGCAGTTGTCTTTGTGCTCGCGGTAAAAGTGGAAGTACCGGGAGAATTTGTTCGGTGACGAGGTGATCTGCATGTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTccac                                                                                                  <  1:562676/143‑1 (MQ=255)
                                               aaTTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCAt                                                                                            <  1:241787/143‑1 (MQ=255)
                                               aaTTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCAt                                                                                            <  1:247195/143‑1 (MQ=255)
                                                aTTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATc                                                                                           >  2:444439/1‑143 (MQ=255)
                                                      cTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGTTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAgg                                                                                     <  2:519542/143‑1 (MQ=255)
                                                        ttGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATCCTTTCGGTGACGAGGTGATCTGCATTTCCCCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGaa                                                                                   <  1:684053/143‑1 (MQ=255)
                                                        tggtgCTTGGGGTAAACGTGGAAGTACCCGGAGGATGTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGaa                                                                                   <  2:452716/141‑1 (MQ=255)
                                                          gtgCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTCCTTTCCACGGTTCGCCGTACTCCACCACCATCCGCACGCAGa                                                                                 >  2:267849/1‑143 (MQ=255)
                                                              ttGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGAGCTGCATTTCACCGTGAGTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTAGCCCTACTACACCACCAGCGGCAGGAAGATGTc                                                                             >  1:536170/1‑143 (MQ=255)
                                                               gcgggTAAAAGTGGACGTGCCCGGAGAAGTTGTTCGGTGACGAGGTGATCTGCAGTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCGCTTCTACTGTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCg                                                                            <  2:472983/140‑1 (MQ=255)
                                                                         gTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCtgtg                                                                  <  1:207364/143‑1 (MQ=255)
                                                                                tACCCGGAGAAGTTGTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCttt                                                           <  1:400697/143‑1 (MQ=255)
                                                                                    cGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCg                                                                                                             <  2:671371/89‑1 (MQ=255)
                                                                                    cGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCg                                                                                                             >  1:671371/1‑89 (MQ=255)
                                                                                       agaATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACCGTTCGCCGTACTCCACCACCATCCGCAGGAAGAGGTCGTAATCCTGTGCCGCTTTCAGTTCg                                                    >  1:366418/1‑143 (MQ=255)
                                                                                       agaATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACAGTTCGCCGTACTCCACCACCATCCGCAGGCAGATGTCGGAATCCTGTGCCGCTTTCAGTTCg                                                    >  1:365490/1‑143 (MQ=255)
                                                                                         aaTTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGt                                                  <  1:336910/143‑1 (MQ=255)
                                                                                            ttttcCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATc                                               <  1:670728/143‑1 (MQ=255)
                                                                                                 cGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGaaca                                          >  1:436952/1‑143 (MQ=255)
                                                                                                           gTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGGAACCTTGTGCCGCTTGCAGTTCGGTATCGAACAAACAcgctcg                                >  2:716551/1‑141 (MQ=255)
                                                                                                                 tGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACg                          <  1:444439/143‑1 (MQ=255)
                                                                                                                           tgtgagtgaTATGCAGGCGCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGGTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAg                <  2:469213/142‑1 (MQ=255)
                                                                                                                                     atGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCt      >  1:326629/1‑143 (MQ=255)
                                                                                                                                         aGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGGGTTCGGTTCGGTATCGAACACACACTCGCTGAACCGCCATGCCCCGGCAAAGACGCGGtt  >  2:253103/1‑143 (MQ=255)
                                                                                                                                            |                                                                                                                                           
GGATCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGGTT  >  NC_000913/2132647‑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: