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

CTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGGTT  >  NC_000913/2132651‑2132926
                                                                                                                                        |                                                                                                                                            
cTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCg                                                                                                                                       >  2:245863/1‑143 (MQ=255)
  ggTAGAGGGTAAACGGCTGATAGTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTc                                                                                                                                     <  1:197511/143‑1 (MQ=255)
    tAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCCCCGTGATTGATATGCAGGCTCTGCGTTGc                                                                                                                                   >  1:171613/1‑143 (MQ=255)
       agGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCctc                                                                                                                                >  1:602983/1‑143 (MQ=255)
        gggTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCctct                                                                                                                               >  1:483515/1‑143 (MQ=255)
         ggTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACACGGAGAATTTTTTCGGTGACACAGTGATATACTTTTCACGGTGATTGATATGCAGGCTCTGCGTGGCCtctt                                                                                                                              >  2:335973/1‑143 (MQ=255)
          gTAGACAGCTGGTATTTTTTGCTGTGGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTGCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCtcttc                                                                                                                             <  2:577738/143‑1 (MQ=255)
               cAGCTGAGATCTTGTCCGGGCGCGGCCGAATCTGTAAATGTGCGTGCGTTATATGTGGAAGTACCCGGCGAATTTTGTCGGCGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACtt                                                                                                                        <  1:550549/143‑1 (MQ=255)
                aGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCGGACCGTGCGTCGCCTATTCTACttt                                                                                                                       >  2:510967/1‑143 (MQ=255)
                        ttttttGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGtt                                                                                                               >  2:155440/1‑143 (MQ=255)
                          tgttGCTGGCGAGGTCGAATTTGTCTTTGTGCTGGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCg                                                                                                             <  1:250760/141‑1 (MQ=255)
                                  gcgcGGTCTTATTTGTCTTTTTGTTTGCTGTAAAATTGGAAGTAGCCGGAGAATTTGTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTc                                                                                                     <  1:299685/143‑1 (MQ=255)
                                   cgcgGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCGACGGTTCGCCGTACTcc                                                                                                    >  2:879752/1‑143 (MQ=255)
                                   cgcgGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACGTTCCACGGTTCGCCGTACTcc                                                                                                    >  2:411032/1‑143 (MQ=255)
                                   cgcgGTCCACTTTCTCTTTGTGCTTGCGGTAAAGTGGGATGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCCACGTGATTTATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTcc                                                                                                    <  2:958444/143‑1 (MQ=255)
                                                   tttGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGa                                                                                    >  1:910436/1‑143 (MQ=255)
                                                                            ttcccGGAGTATTTTTTCGGTGACGGGGTGATCTGCATTTCACCTTGATTGATAGGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCTCCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCttt                                                           <  2:323786/141‑1 (MQ=255)
                                                                                 gtatacttttttCGGTGACGATGTCATCTGCATTTCTCCGTGTTTGATATGCACGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGtt                                                      <  1:418755/137‑1 (MQ=255)
                                                                                 ggAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCCCCACCATCCGCAGGAAGATGTCGTCATCCTGTGCCGCCTTCAttt                                                      >  1:542682/1‑143 (MQ=255)
                                                                                  gagaATGTTGTCTGTGTCGGGGGTGTCTGCATTTCACCCGGATTGTTATGCAGGCTCTGCGTCGCTGTTTCTACTTTCCACGGTTCTCCGTACTCCACCCCCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTc                                                     <  2:306098/143‑1 (MQ=255)
                                                                                         ttttCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCGGGCTCTGCTCCTCCGCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCg                                              <  1:447957/143‑1 (MQ=255)
                                                                                         ttgtCGGTGACGGGGTGATCTGCATTTCACCGTGATTGAGATGCAGGCTCTGCGTCGGCTCTTCTACTGGCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCg                                              <  2:228439/140‑1 (MQ=255)
                                                                                               gggatgAGGTGATCTGCACTTCACCGTGATTGAAATCCAGGCTCTGCGTCGCCTCGTCTACTTTCCACGGTTCGCCGAACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGaacaaa                                        <  1:415504/138‑1 (MQ=255)
                                                                                                     aGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACActct                                  <  2:602983/143‑1 (MQ=255)
                                                                                                      cgtgggcTGCTTTTCTCCGTGACTGCTATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCTCCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCtt                                 <  2:973533/137‑1 (MQ=255)
                                                                                                            cTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGCACTCCTGTGCCGCTTTCAGTTCGGTATCGAACAACCACCCTTTGaacc                           >  2:77611/1‑141 (MQ=255)
                                                                                                              gCATTTCACTGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTATGAAACGc                         >  1:773957/1‑143 (MQ=255)
                                                                                                               cATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGcc                        <  2:483515/143‑1 (MQ=255)
                                                                                                                aTTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCCTCCGCAGGAAGCTGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGCAACCCCa                       >  2:226934/1‑143 (MQ=255)
                                                                                                                       cGTGCTTGATATGCAGGCTCTGCGCCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAg                <  1:379214/143‑1 (MQ=255)
                                                                                                                          gattgatATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTa             <  1:155440/143‑1 (MQ=255)
                                                                                                                           attgatATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCtgtg                                                                  >  2:566954/1‑89 (MQ=255)
                                                                                                                           attgatATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCtgtg                                                                  <  1:566954/89‑1 (MQ=255)
                                                                                                                               atatGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCCCTTTCAGTTCGGTCTCGAACAACCACCCTTTGAAACGCCACGCCCCGGTACAGAc        >  1:163801/1‑143 (MQ=255)
                                                                                                                                   gCAGGCTCTGCGTCGCCTCCTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTgg    <  1:584950/143‑1 (MQ=255)
                                                                                                                                    cAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCACCGCAAGATGTCGTAATCCTGTGCCGCTTTCAGTACGCTATCGAACACACACTCTTCGCAACCCCATGCCCAGGTAAACACCTgac   >  2:532034/1‑141 (MQ=255)
                                                                                                                                     aGGCTCTGCGTCGCCTCTTCGACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGGtt  <  1:332038/143‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                            
CTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGGTT  >  NC_000913/2132651‑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: