Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,132,787 A→C 100% 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,7870AC98.9% 76.0 / ‑6.7 27I204S (ATC→AGC) wcaAputative colanic acid biosynthesis glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base C (13/13);  total (13/14)
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.89e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GCGGATCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGG  >  NC_000913/2132645‑2132924
                                                                                                                                               |                                                                                                                                         
gCGGATCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTAGAATTTGTCTTTGTGCTTGCGGTACAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGCGAGTGACATGCGGGa                                                                                                                                           >  2:129710/1‑143 (MQ=255)
          tAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTTCCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGc                                                                                                                                 <  1:391937/143‑1 (MQ=255)
                   aaCAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTAc                                                                                                                        >  1:338050/1‑143 (MQ=255)
                        cTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTGCTTTcc                                                                                                                   >  1:289569/1‑143 (MQ=255)
                             atTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTGCCCGGAGAATTTTTTCGGTGACGGGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGt                                                                                                              <  2:315359/143‑1 (MQ=255)
                              ttttttGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCGCCGTGATTGATATGCAGGCTCTGCGTCGCCTCGTCTACTTTCCACGGtt                                                                                                             >  1:309727/1‑143 (MQ=255)
                                   tGCTGGCGCGGTCGAATGTGTCTGTGGGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCg                                                                                                        <  2:201106/143‑1 (MQ=255)
                                      tGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGATAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTAc                                                                                                     <  1:28028/143‑1 (MQ=255)
                                         cgcgGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGTCGCTGCGTCGCCTCTTCTACTTCCTACGGTTAGCCGTACTcc                                                                                                  >  2:572249/1‑143 (MQ=255)
                                                gAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGCCGAGGTGATCTGCATTTCACCGTGATTGTTATGCAGGCGCGGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTccaccacca                                                                                           >  2:70803/1‑143 (MQ=255)
                                                 aaTTTGTCTTTGTGCTTGCGGTATAAGTTTGAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTCCTCCACCCCCAt                                                                                          >  2:546440/1‑143 (MQ=255)
                                                 aaTTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTACCCGGTTCGCGGTACTCCAGCACCAt                                                                                          >  2:342066/1‑143 (MQ=255)
                                                   tttGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGCGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATcc                                                                                        <  1:90222/143‑1 (MQ=255)
                                                          ttGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGaa                                                                                 >  1:307503/1‑143 (MQ=255)
                                                                                    tgcGGAGAATTTTTTCGAGGACGAGGCGACCAGCAGGTCACCGCGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTGCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCa                                                       <  1:228048/141‑1 (MQ=255)
                                                                                     actgagaATTTTTTCGTCTACGAGTTCATCTGCATTTCACCTTTATTGATATGCAGGCTCTGCTTCGCTTCTTCTCCTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAg                                                      <  2:81080/140‑1 (MQ=255)
                                                                                       gcagaATTTTTTCGGTGACGGGGAGATCGCCATGTCACCGAGGTCAGTATCCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGtt                                                    <  1:84460/141‑1 (MQ=255)
                                                                                            atctgttCGGTGACGAGGTGAGCTGCCTTTCCCCGGGATTGCGAGCCCGGCCCGGCGCCGCCTCCTCCACTTTCCACGGTTCGCCGTACTCCACCTCCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTa                                               <  1:23103/138‑1 (MQ=255)
                                                                                             ttttttCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTACACttagatct                                              >  2:284769/1‑137 (MQ=255)
                                                                                              ttgttgggTGACGAGGTGATCTGCCTTTCCCCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATc                                             <  2:309727/137‑1 (MQ=255)
                                                                                                 ttCGGTGACGAGGTGCTCTGCTTCGATTGTTGGTGGATATGGAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGaa                                          <  1:84857/143‑1 (MQ=255)
                                                                                                              tGATCTCCCTTTCACCGTGATTGCTATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTg                             <  2:561051/143‑1 (MQ=255)
                                                                                                                    gCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCACCCGCCGGAAGATGTCGTAATCCTGTGCCCCGTTCGGTTCGCTCTCGAACAAACACTCATAGAAAccc                       >  1:35014/1‑141 (MQ=255)
                                                                                                                     cATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCACGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGcc                      >  1:312201/1‑143 (MQ=255)
                                                                                                                        ttCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCAACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATg                   >  2:23958/1‑143 (MQ=255)
                                                                                                                                 gattgatATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTGATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACCCTCTTTGAAACGCCATGCCCAggcc           >  2:199442/1‑141 (MQ=255)
                                                                                                                                   ttgatGTGCAGGCTTTGTGTCGCCTCTTCTACTGTCCAGGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTaaa         <  2:253620/143‑1 (MQ=255)
                                                                                                                                          gCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTgg  <  1:104811/143‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                         
GCGGATCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGG  >  NC_000913/2132645‑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: