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,7870AC95.6% 105.1 / NA 46I204S (ATC→AGC) wcaAputative colanic acid biosynthesis glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base A (2/0);  new base C (26/17);  total (29/17)
Fisher's exact test for biased strand distribution p-value = 5.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGATCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGGTT  >  NC_000913/2132647‑2132926
                                                                                                                                            |                                                                                                                                           
ggATCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCTGTAAAAGTGGAAGTACCCGGAGACTTTTTTCTGTGACGAGGTGTTCTGCATTTTACCGTAGCTGATATGCAGAATc                                                                                                                                           >  1:853386/1‑143 (MQ=255)
ggATCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGACTTTTTTCGGTGACGAGGTGTTCTGCCTTTCACCGTGTTTGTTCTGCAGGCTc                                                                                                                                           >  1:132156/1‑143 (MQ=255)
     tGGTAGAGGGTAACCAGGTGATATTGTTGTGTGGCGCGGTCGAATTTGTCCTTGTGCTTGCGGTAAAAGTGTAAGTACCCGGAGAATTGTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGt                                                                                                                                      <  1:268570/143‑1 (MQ=255)
      ggTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCgcc                                                                                                                                     >  1:143049/1‑141 (MQ=255)
      ggTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTACTTGCGGTAACAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCTTTTCACCGTGAGTGATATGCAGGCTCTGCGTc                                                                                                                                     >  1:608936/1‑143 (MQ=255)
        tAGAGGGTCACCGGCTGAGCGTTTTTGCTGGCGCGGTCGAATGTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGc                                                                                                                                   <  2:370259/143‑1 (MQ=255)
        tAGAGCGTCAACAGCTGCTATTTTTTGCTGGCGCGGTCGAATTTGTCCTTGTGCTTGCGGTAACAGTGGAAGTACCCGGAGAATTTTTTCGGTGCCGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGc                                                                                                                                   <  1:1541/143‑1 (MQ=255)
          gagGGTAAACAGCTGATATTTTTTTTTTTTTTTTTCGAATTTGTCTTTGTGCTTGCGTTCCATGTGGAAGTACCCGGAGAATTTTTTCGGTGCCGAGGTGCTCTGCATTTCACCGTGATTGATTTGCTGCTTCTGCGTCGCct                                                                                                                                 >  2:1099227/1‑143 (MQ=255)
             ggTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGCCGAGGTGATCTGCATTTCACCGGGATGGATATGCAGGCTCTGCGTCGCatctc                                                                                                                              >  1:534721/1‑142 (MQ=255)
             ggTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGAAGAGGTGATCTGCATTTCAACGTGATTGCTATGAAGGCTCTGCGTCGCcttct                                                                                                                              >  2:1019594/1‑140 (MQ=255)
               tAAACAGCTGATATTTTTTGCTGGCGCGGTCGAGTTTGTCTTTGTGCTTGCGGTAAAGGTGGAAGGACCCGGAGAATTTTTTCGGGGACGAGGTGATCTGCATTTCACCGTGGTGGATATGGAGGAGTTGCGTCGCctctggg                                                                                                                            >  1:20540/1‑140 (MQ=255)
                caaCAGCTGATCTTTGGTGCTGGCGCGGTCGAACTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTGCGGTGACGGGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTa                                                                                                                           <  2:688535/142‑1 (MQ=255)
                aaaCAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTa                                                                                                                           >  1:59577/1‑143 (MQ=255)
                  aCAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGTTAAAAGTGGAAGTACCCGGAGACTTTTTTCGGTGACGAGGGGATCTGCATTTCACCGTGATTGATCTGCAGGCTCGCCGTCGCCTCTTCTACt                                                                                                                         >  2:757559/1‑143 (MQ=255)
                     gCTGATCTTTTTTGCTGGCGCGGTCGATGTTGGCTTTGTGCTTGCGGTATAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTc                                                                                                                      <  2:1114577/143‑1 (MQ=255)
                     gCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGTGTCGCCTCTTCTACTTTc                                                                                                                      >  1:861595/1‑143 (MQ=255)
                               tttGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTGGCCTCTGCTACTTTCAACGGtgccg                                                                                                            >  2:30364/1‑139 (MQ=255)
                                   cTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATATGCATTTCACCGTCATTGATATGCAGGCTCTGCGTCGCCTCTACTACTTTCCACCGTTCGCCGTa                                                                                                        >  2:277692/1‑143 (MQ=255)
                                    tGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTTATTGATATGCAGGCTCTGCGTCCCCTCCTCTACTTTCCCCGGTTCGCCGTAc                                                                                                       >  2:420014/1‑143 (MQ=255)
                                      gcgcGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTc                                                                                                     >  1:342058/1‑143 (MQ=255)
                                       cgcgGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTGCGGTGACGAGGTGAGCTGCACTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTcc                                                                                                    <  1:209991/143‑1 (MQ=255)
                                                aTTTGTCTTTGTGCTTGCGGTAAAAGTAGAAGGAACCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATc                                                                                           >  1:768793/1‑143 (MQ=255)
                                                 tttGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACAACCATcc                                                                                          >  1:477912/1‑143 (MQ=255)
                                                     tCTTGGTGCTGGCGGTGAAAGTGGAAGTGCCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAg                                                                                      <  2:202360/143‑1 (MQ=255)
                                                             cTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCCGGAAGCTgc                                                                              >  2:552375/1‑142 (MQ=255)
                                                                   gTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAAt                                                                        >  1:210143/1‑143 (MQ=255)
                                                                      ataGTGGAAGTTCCCGGAGACGTTTTTCGCCGACGCTGTGATCTGCATTTCACCGTGATTGATATGCAAGCTCTGCGCCGTCTCTTCTACTTTCCACGGTTCACCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCt                                                                     <  2:280092/141‑1 (MQ=255)
                                                                             aaGTGCCCGGAGAATTGTTTCGGTGACGAGGGGATCTGCATCTCACCGGGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTGTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGc                                                              <  1:730390/143‑1 (MQ=255)
                                                                             aaGTACCCGGAGAATTTTTTCGGTGACGAGGTCATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGc                                                              <  1:831299/143‑1 (MQ=255)
                                                                              aGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCGCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGGAGCTGTCGTACTCCCGTGCCGcc                                                             >  2:604469/1‑142 (MQ=255)
                                                                                caCCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCttt                                                           <  1:964436/142‑1 (MQ=255)
                                                                                   ccGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAg                                                        >  1:860455/1‑143 (MQ=255)
                                                                                            tttttCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCAGCGGCAGCAAGACGTCGTAATCCTGTGCCGCTTTGAGTTCGGTATc                                               >  1:319703/1‑143 (MQ=255)
                                                                                             ttttCGGTGCCGGCGTGAGCTGCCTTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTCTCCACGGTTCGCCGTTCTCCACCACCAGCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCg                                              <  2:496778/143‑1 (MQ=255)
                                                                                                tCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGaac                                           <  1:307551/143‑1 (MQ=255)
                                                                                                ccGGGGACGAGGTCAGCTGCCTTTCACCGTGATTGGTATGCAGGCTCTGCGTCGCCTCTTCTACGTTCCGGGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGaac                                           <  2:608527/142‑1 (MQ=255)
                                                                                                   gTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGaacaaa                                        >  1:973023/1‑143 (MQ=255)
                                                                                                   gTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCCTTCAGTTCGGTATccacccaa                                        >  2:650990/1‑136 (MQ=255)
                                                                                                            tGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGAATCGAACAAACACTCTTTg                               >  2:1031464/1‑143 (MQ=255)
                                                                                                             gATCTGCATTTCACCGTGATTGATATGCAGGCCCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGa                              <  2:1052137/143‑1 (MQ=255)
                                                                                                                   cATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCAGCCGCAGGAAGGTGGCGGGAGCCTGTGCCGCGCTCAGTTCGGTATCGAACAAGCACTCTGTGAAACGcc                        >  1:206909/1‑143 (MQ=255)
                                                                                                                           cGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGGTCGGTATCGAACAAACACTCTTTGAAACGCCATCCCCAg                >  2:492232/1‑143 (MQ=255)
                                                                                                                                   atatGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCACCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTACCGCTTCCATTTCGGTATCCAACAAACACTCTTTGCAACACCACGCCCAGGTAAAGAc        >  2:463766/1‑143 (MQ=255)
                                                                                                                                     atGCAGGCTCTGCGTCGTCTCTTCTACGTTCCACGGTTCGCCGGACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCt      <  2:608936/143‑1 (MQ=255)
                                                                                                                                     atGCAGGCTCGGCGTCGCCTCGTCTGCTGTGCACGGTTCGCCGTACTCCGCCACCATCCGCAGTAAGCTGTCGTACAACTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCt      <  2:845679/143‑1 (MQ=255)
                                                                                                                                         aGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTCCTCCACCACCATCCGCAGGAAGATGCCGGAAGCCTGTGCCGCTTTCAGTCCGGTAGCGAACCAACACGCTTCGACACGCCACGACCACGACAAGCCCCGGtt  >  2:808791/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 < 22 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: