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% 82.3 / NA 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 G (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 = 9.98e-01

TCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTG  >  NC_000913/2132650‑2132923
                                                                                                                                         |                                                                                                                                        
tCTGGTAGAGGGTATACAGCTGATATTTTGTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGCAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGc                                                                                                                                     <  2:496935/143‑1 (MQ=255)
    ggagagGCCAAGCAGCTGAGATTTTGTGGTGGCCCGGTCGAATTTGTCGGTGTGCTGGCGGTAAAAGTGGAGGTACCCGGAGAATGTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCg                                                                                                                                 <  2:247336/141‑1 (MQ=255)
      cgagGGTAGACGGCTGATATCGCTTGCGGGCGCGGTCGAATGTGGCTGTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGcc                                                                                                                               <  1:179521/142‑1 (MQ=255)
      agagGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTCTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGcc                                                                                                                               >  1:599058/1‑143 (MQ=255)
       gagGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCct                                                                                                                              >  1:348188/1‑143 (MQ=255)
         gggTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGCTATGCAGGCTCTGCGTCGCctct                                                                                                                            >  2:623112/1‑143 (MQ=255)
             aaaCAGCTGATAGCTGTCGCGGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTa                                                                                                                        <  2:581274/143‑1 (MQ=255)
                      atatTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACg                                                                                                               <  1:249908/143‑1 (MQ=255)
                          tttttGCTGGCGCGGTCGAATTTGGCTATGCGCTTGCGGTAAAAGTGGTAGTACCCGGCGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTc                                                                                                           <  1:76818/143‑1 (MQ=255)
                                           gAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCAGTAACGAGGTGCTTTGCATCTCACCGTGATTCATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGCTTGGCCGTACTccaccacca                                                                                          >  2:379081/1‑143 (MQ=255)
                                                 gTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGAACTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCATCTACTTTCCACGGTGCGCCGGACTCCCAGACCATCCGCa                                                                                    >  2:500377/1‑143 (MQ=255)
                                                          cTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGt                                                                           >  1:70946/1‑143 (MQ=255)
                                                           ttGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCCTACTCCACCACCCGCCCCAGGAAGAGGTc                                                                          >  2:418709/1‑143 (MQ=255)
                                                           gtGCGGGAAAGGTGGACGTGCCCGGAGAGTTTTTTCGGTGACGGGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTGTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTc                                                                          <  2:342887/142‑1 (MQ=255)
                                                                     gggggAAGTGCCCGGAGAATTTTTTCGGTGACGGGGTGATCTGCATTTCACCGTGATTGATGTGCAGGCGCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCtgt                                                                <  2:517516/140‑1 (MQ=255)
                                                                     aGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTGACCGTAATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCCCCATACTACACCACCATACTCACAAAAATCTCGTAAACCtgt                                                                >  2:531999/1‑143 (MQ=255)
                                                                      gTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCtgtg                                                               >  1:186907/1‑143 (MQ=255)
                                                                       tGGAAGTACCCGGAGAATTTTTTCGGTGACGCGGTGATCTGCATTCCACCGTGATTGATATGCAGGCTCTGCGTCGCCGCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCCGGAAGCTGTCGTCATCCTGTGc                                                              >  2:476046/1‑143 (MQ=255)
                                                                        gggagggcccGGAGAACGTTTGCGGGGACGAGGTGATCTGCAGCCCCCCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACGCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGcc                                                             <  1:60300/136‑1 (MQ=255)
                                                                               cgcGGAGAATTTGTTCTGTGACGAGGTGATCTGGATTTCACCGTGATTGATATGCAGGCTCTGCCTCGCCTTTTCCGCTTTCCACGGTTCGCCGTACTCCCCCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCa                                                      <  2:563465/141‑1 (MQ=255)
                                                                                       aTTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCCGGGTCTGAAACCCCTCTTCGACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTa                                              <  1:623112/143‑1 (MQ=255)
                                                                                        ttttttCGGTGACGAGGTGATCTGCATTTCACCGTGATTGCTATCCCGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCCTCCGCAGGAAGGTGTCCTAATCCTGTGCCGCTCTCCGTTCGGTAt                                             >  1:39240/1‑143 (MQ=255)
                                                                                          ttttCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCg                                           <  1:344402/143‑1 (MQ=255)
                                                                                           tttCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCCTCCGCAGGAAGATGTTGTAATCCTGTGCCGCTTTCAGTTCGGTATCGa                                          <  1:173507/143‑1 (MQ=255)
                                                                                              cGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCCGGCAGATGTCCTAAGCCTGTGCAGCTTTCAGTTCCGTATcagaca                                       >  2:457055/1‑138 (MQ=255)
                                                                                                gTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGaacaac                                     >  1:597974/1‑142 (MQ=255)
                                                                                                        gTGATCTCCATTGCACCGTGCTCGGGCTGCAGGCTCTGCTTCGCCCCTTCTACTCTCCACCGTTCTCCTTCCGCCTCCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCttt                             <  2:554735/143‑1 (MQ=255)
                                                                                                                        cttgcttgatACGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAg             <  1:315927/138‑1 (MQ=255)
                                                                                                                                   tGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTACTCCTGTGCCGCTTTCAGCTCGGTATCGAACCAACACTCTTTGCAACGCCATGCCCACGTAAAGACATg  >  2:428714/1‑143 (MQ=255)
                                                                                                                                   tGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTg  <  2:70946/143‑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 < 22 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑