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.7% 69.4 / ‑6.8 30I204S (ATC→AGC) wcaAputative colanic acid biosynthesis glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base C (14/15);  minor base T (0/1);  total (14/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTG  >  NC_000913/2132650‑2132923
                                                                                                                                         |                                                                                                                                        
tCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGc                                                                                                                                     <  1:644674/143‑1 (MQ=255)
    gTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGttg                                                                                                                                 >  1:122297/1‑141 (MQ=255)
           gTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCtcttc                                                                                                                          <  2:673217/143‑1 (MQ=255)
            tAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCtcttct                                                                                                                         >  1:502131/1‑143 (MQ=255)
                         ttttttGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGTGTCGCCTGTTCTACTTTCTACGGtt                                                                                                            >  1:199980/1‑143 (MQ=255)
                             ttgttgGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGcc                                                                                                        <  2:662555/143‑1 (MQ=255)
                                     gcgGTCGAATTTGTCTTTGGCATTGCGGGAAAAGTGGAAGGACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTcca                                                                                                <  2:148244/143‑1 (MQ=255)
                                                    tttGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCCggc                                                                                 >  2:507708/1‑142 (MQ=255)
                                                    tttGTCCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTCCTCCACCACCATCCGCAGGa                                                                                 >  2:533922/1‑143 (MQ=255)
                                                     ttGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACCAGGTAATAAATGTTTCGCCGTGATTCATATTCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGaa                                                                                <  1:793389/143‑1 (MQ=255)
                                                        tgCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGAt                                                                             <  2:502131/143‑1 (MQ=255)
                                                                gTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCACCGTACTCCACCACCATCCGCAGGAAGATGTCGTAAt                                                                     >  1:593082/1‑143 (MQ=255)
                                                                 tAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCCTCCGCCGTAAGACGTCGCAATc                                                                    >  2:412231/1‑143 (MQ=255)
                                                                     aGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCtgt                                                                >  1:722663/1‑143 (MQ=255)
                                                                       tGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGCAGCTGTCGTCACCCTGTGc                                                              >  2:736575/1‑143 (MQ=255)
                                                                            gTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCGTTTTCTCCTTTCCACGGTTCGCCGTACTCCACCACCGTCCGCAGGAAGATGTCGTAATCCTGTGACGCtt                                                         >  2:11751/1‑143 (MQ=255)
                                                                               cccGGAGGATTTTTTCGGTGACGAGGTGATCTGCAGTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCTCGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCa                                                      <  2:740720/143‑1 (MQ=255)
                                                                                  ggAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGCTGTTGTCTTTATCCTTTGCCTCTTTCAGtt                                                   >  1:455100/1‑143 (MQ=255)
                                                                                            tccGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGaa                                         <  1:575203/141‑1 (MQ=255)
                                                                                                 tGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCCGTCCGGTCTCGaacaaac                                    >  2:70828/1‑143 (MQ=255)
                                                                                                     gAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTCCTACTTAACACGGTCGGCGGTAGTCCACCACCATCCGCCGGAAGCTGTCGTACTCCTGTGCCGCCTTCAGTTCGGTATCGAACAAAAActc                                >  2:738068/1‑143 (MQ=255)
                                                                                                      aGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCACCTCTTCTACTTTCCACGGTTCCCCGTACTCCCCCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAGACAcgct                               >  1:439720/1‑143 (MQ=255)
                                                                                                          gATCTGCATTTCACCGTGATTGATATGCAGGCTTTGCTTCGGCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCCCCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGa                           <  2:562145/143‑1 (MQ=255)
                                                                                                               gCAGTTCACCGTGATTGAGATGCAGGTTCGGCGTCGCCGCTTCTACTGGCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGc                      <  1:67940/143‑1 (MQ=255)
                                                                                                                  tttCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACGTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCAt                   <  2:122297/143‑1 (MQ=255)
                                                                                                                   ttCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATg                  <  1:408046/143‑1 (MQ=255)
                                                                                                                        cGTGATTGATATGCAGGCTCTGCGTCGCCTGTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGATTTCAGTGCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAg             <  1:590602/143‑1 (MQ=255)
                                                                                                                                atatGCGGGCTCTGCGTCGCCCCTTCTACTGTGCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGAGGTCGTAATCCTGTGTCGCTTTTAGTTCGGTATCGAACAAACACTCTTGGAAACGCCATGCCCAGGTAAAGAc     <  2:481617/143‑1 (MQ=255)
                                                                                                                                   ggCAGGCTCTGCGTCGCCGCTTCTACTTTCCGCGGTTGTCCGTGCTCCACCACCATCCGCAGGAAGATGTCGGAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTg  <  2:237811/142‑1 (MQ=255)
                                                                                                                                   ggCAGGCTCTGCGTCGCCGCTGCTGCTTTCCACGGTTCGCCGTACTCCACCACCAGCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTGAGGTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTg  <  2:464877/142‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                        
TCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTG  >  NC_000913/2132650‑2132923

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: ‑