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,7870AC96.5% 75.9 / NA 29I204S (ATC→AGC) wcaAputative colanic acid biosynthesis glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base C (12/16);  minor base T (0/1);  total (12/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

CAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGG  >  NC_000913/2132666‑2132924
                                                                                                                         |                                                                                                                                         
cAGCTGATATTTTTTGCTTGTGCAGTCGATTTTGTCTTTTTGTTTGTGGTGGACGAGCAAATACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACtt                                                                                                                      <  1:385364/143‑1 (MQ=255)
  gCTGATATTGTTTGCTGGTGCGGTCGAATTTGTCTTTGTGCTTGCGGGAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCCGCGTCGCCTCTTCTACTTTc                                                                                                                    <  1:120097/143‑1 (MQ=255)
                cgggggcgGTCGAGTTTGTCGTTGGGCTTGCGGTAAAAGTGGAAGTAGCCGGAGAATTTTTTCGGTGACGCGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTa                                                                                                      <  1:817796/138‑1 (MQ=255)
                cTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGCTCGCCggc                                                                                                      >  2:753788/1‑141 (MQ=255)
                       ggTCGAATTTGTCTGTGTGCTTGCGGTGAGAGTGGATGTGCCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCGCCGTGATTGATATGCAGGCTCTGCGTCGTCTCTTCTACTTTCCACGGTTCGCCGTACTccacc                                                                                               <  2:318535/143‑1 (MQ=255)
                            aaTTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAAGCTCTGCATCACCTCAAATACTTTCCACGGTTCACCGTACTCCACCACCAt                                                                                          >  1:224920/1‑143 (MQ=255)
                             tgttGTCTTGGTGCCTGCGGGACGAGTGGACGTCCCCGGAGAATTTTTTCGGTGACGATGTGATCTGCCTTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATc                                                                                         <  2:595797/141‑1 (MQ=255)
                                 gTCTTTGTGCTTGCGGTAAAAGTGGAAGTATCCGGAGATTATTGGCATAAACGAGGTGATCTGCATATCACCGTGATTGATATGCAGGCTCTGCGTCGCCTTTTCTATTTTCCACGGTTCGCCGTACTCCACCACCATCCGCa                                                                                     <  1:113841/143‑1 (MQ=255)
                                 gTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTAACCCACCCCCATCCGCa                                                                                     >  2:808694/1‑143 (MQ=255)
                                                gTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGTTTCGCCGTACTCCACCACCATCCGACAGAAGATGTCGTAAt                                                                      >  1:606263/1‑143 (MQ=255)
                                                      gtgtatgTATCCGGAGATTTTTTTCGGTGACGTGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCtgtg                                                                <  1:695743/137‑1 (MQ=255)
                                                       tGGGAGCACCCGGAGAATTTTTGCGGTGACGAGTGGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGc                                                               <  2:774954/143‑1 (MQ=255)
                                                                ccGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTGGTAATCCTGTGCCGCTGTCAg                                                      >  2:424012/1‑143 (MQ=255)
                                                                   gtgtgttttttCTGTGTTGATGTGATCTACATTTCACAGTGATTGATATGCACGCTATGAGTCGTCGCTACTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCACGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTc                                                   <  1:429054/138‑1 (MQ=255)
                                                                              cGGTGACGCGGTGATCTGCATTTCTCCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGaaca                                        <  2:603107/143‑1 (MQ=255)
                                                                                  gACGGGGTGATCTGCATGTCACCGTGATTGAGATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGGAATCCTGTGCCGCTTTCAGTTCGATACCGaactgaaa                                    <  1:808694/143‑6 (MQ=255)
                                                                                  gACGAGGTTATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTCCTCCACCACCATACGCAGGAAAATGTCGTAATCCTGTGACGCGTTCAGTTCCGTAGCGaacaaaca                                    >  2:18329/1‑143 (MQ=255)
                                                                                      aGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCCTGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGAATCGAGCAAACActct                                >  2:481261/1‑143 (MQ=255)
                                                                                       ggTGCTCTGCATTTCACCGTGATTGATGTGGAGGCTCTGCGTCGCCTCTTCTACGTTCCACGGTTCGCCGTACTCCACCCCCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCtt                               <  2:536529/143‑1 (MQ=255)
                                                                                        gTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCttt                              >  1:449209/1‑143 (MQ=255)
                                                                                        gTGATCTGCATTTCACCGTGATTGATATGAAGACTCTGTGTCGCCTCTTCTACTTTCCACAGTTCGCCGTACTCCACCACCATCCCCAGGAAGATGTCGTGATCCTGTGGCGGTTTCAGTTCGGTAGCGACGAAACActctct                              >  2:297685/1‑141 (MQ=255)
                                                                                           aTCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCAGCACCATCCGCAGGAGGATGTCGTACTCCTGTGCCTCTGTCTGTTCGGTATCGAACAAACACCCTTTGaa                           >  1:127024/1‑143 (MQ=255)
                                                                                            gcTGCATTTCACCGTGATTGATATTCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGaaa                          <  2:353445/142‑1 (MQ=255)
                                                                                                       gcGTGATTGATATGCAGGCTCTGCGTCGTCTCTTCTACTTTCCCCGGTTCGCCGGACTCCACCACCAGCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTTAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCa               <  2:608684/142‑1 (MQ=255)
                                                                                                            attgatATGCAGTTTCTGCGTTGCTTCTTCTACTTTTCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGGAATCCTGTGCCGCTTTCAGTTCGGTATCTAACAAACACTCTTTGAAACGCCATGCCCAGGTaa          <  1:84659/143‑1 (MQ=255)
                                                                                                             ttgatATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCACGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGACCAAACACTCTCTGAAACGCCATGCCCAGGTaaa         >  2:37080/1‑143 (MQ=255)
                                                                                                              tgagacgCAGGCTCTGCGTCGCCTCTTCTACTGTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTTTGCCGCTTGCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAg        <  2:178382/137‑1 (MQ=255)
                                                                                                                  atGCAGGCTCTGGGGCTCCTCTTCTACTGTCCACGGTTCGTCGTACTCCACCACCAGCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCt    <  2:164970/143‑1 (MQ=255)
                                                                                                                    gCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCTCCATCCGCAGTCAGATGTCGTAATACTATGCTGCTTTCAGTTCGGTATCGAAAAAACACTCTTTGAAAAGCCATGCCCAGGTAAAGATCTgg  >  1:82025/1‑143 (MQ=255)
                                                                                                                         |                                                                                                                                         
CAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGG  >  NC_000913/2132666‑2132924

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: