Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,132,104 A→G V151A (GTG→GCG)  wcaB ← putative colanic acid biosynthesis acyl transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,132,1040AG93.5% 89.5 / ‑1.8 32V151A (GTG→GCG) wcaBputative colanic acid biosynthesis acyl transferase
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (11/18);  total (12/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.95e-01

GACGCCAGCCCCTGTTGCAGCGCACGCTGGTGGAGATCTAACGCCACACCTGCTGCCCCGCCTTCCGCCAGTCGCACATTAAATTGCAAAATATTCATTTAATTACCTTCACTCGCGCTTTTTCTCCCACCACCAGCGCGTTGTCCGGGACAGAGTCGAGCACCACGCTGCCCGCGCCCACGGTGACGTTGTTACCAAGCGTGATATCACCAAGAATAATGACGTTGGCACCGAGTTCGACGCCGTTGCCAATGTGTGGACATGCC  >  NC_000913/2131973‑2132238
                                                                                                                                   |                                                                                                                                      
gACGCCAGCCCCTGTTTCAGCGCACGCTGGTCGGGTTCTAACGCCACACCTGCTGCCCCGCCTTCCGCCAGTCGCACATTAAATTGCAAAATATTCATTTAATTACCTTCACTCGCGCTTTTTCTCCCACCGCCAGcgcg                                                                                                                                <  2:344910/140‑1 (MQ=255)
  cGCCAGCCCCTGTTGCATCGCATGCTGGTGGTGATGTAATGCCACACCCGCTGCCCCGCCTTCCGCCAGTCGCACATTAAATTGCAAAATATTCATTTAATTACCTTCACTCGCGCTTTTTCTCCCACCGCCAGCGCGtt                                                                                                                              <  2:331686/140‑1 (MQ=255)
     cAGCCCCTGTTGCAGCGCACGCTGGTGGAGATCTAACGCCACACCTGCTGCCCCGCCTTCCGCCAGTCGCACATTAAATTGCAAAATATTCATTTAATTACCTTCACTCGCGCTTTTTCTCCCACCGCCAGCGCGTTGTc                                                                                                                           >  2:416302/1‑140 (MQ=255)
     cAGCCCCTGTTGCAGCGCACGCTGGTGGAGATCTAACGCCACACCTGCTGCCCCGCCTTCCGCCAGTCGCACATTAAATTGCAAAATATTCATTTAATTACCTTCACTCGCGCTTTTTCTCCCACCGCCAGCGCGTTGTc                                                                                                                           <  1:636405/140‑1 (MQ=255)
     cAGCCCCTGTTGCAGCGCACGCTGGTGGAGATCTAACGCCACACCTGCTGCCCCGCCTTCCGCCAGTCGCACATTAAATTGCAAAATATTCATTTAATTACCCTCACTCGTGCTGTTGCTCCCATCGCCAGCGCGTTGTc                                                                                                                           >  1:344910/1‑140 (MQ=255)
            tGTTGCAGCGCACGCTGGTGGATATCTAACGCCACCCCTTCTGCCCCGCCTTCCGCCAGTCGCACATTAAATTGCAAAATATTCATTTAATTACCTTCACTCGCGCTTTTTCTCCCACCGCCAGCGCGTTGTCCGGGACa                                                                                                                    <  2:278427/140‑1 (MQ=255)
                gCAGCGCACGCTGGTGGAGATCTAACGCCACACCTGCTGCCCCGCCTTCCGCCAGTCGCACATTAAATTGCAAAATATTCATTTAATTACCTTCACTCGCGCATTTTCTCTCACGGCCAGCGCGGTGTCCGGGACAGAGt                                                                                                                >  1:233109/1‑140 (MQ=255)
                 cAGCGCACGCTGGTGGAGATCTAACGCCACACCTGCTGCCCCGCCTTCCGCCAGTCGCACATTAAATTGCAAAATATTCATTTAATTACCTTCACTCGCGCTTTTTCTCCCACCGCCAGCGCGTTGTCCGGGACAGAGTc                                                                                                               <  1:669806/140‑1 (MQ=255)
                 cAGCGCACGCTGGTGGAGATCTAACGCCACACCTGCTGCCCCGCCTTCCGCCAGTCGCACATTAAATTGCAAAATATTCATTTAATTACCTTCACTCGCGCTTTTTCTCCCACCGCCAGCGCGTTGTCCGGGACAGAGTc                                                                                                               <  1:619391/140‑1 (MQ=255)
                                gagaTCTAACGCCACACCTGCTGCCCCGCCTTCCGCCAGTCGCACATTAAATTGCAAAATATTCATTTAATTACCTTCACTCGCGCTTTTTCTCCCACCGCCAGCGCGTTGTCCGGGACAGAGTCGAGCACCACGCTGcc                                                                                                <  1:416302/140‑1 (MQ=255)
                                    tCTAACGCCACACCTGCCGCCCGGCCTTCCGTCAGGCGCGCATTAAATTGCAAAATATTCATTTAATTACCTTCACTCGCGCTTTTTCTCCCACCGCCAGCGCGTTGTCCGGGACAGAGTCGAGCACCACGCTGCCcgcg                                                                                            <  2:233109/140‑1 (MQ=255)
                                     cTAACGCCACACCTGCTGCCCCGCCTTCCGCCAGTCGCACATTAAATTGCAAAATATTCATTTAATTACCTTCACTCGCGCTTTTTCTCCCACCGCCAGCGCTTTGTCCGCGACTGAATCGAGCGCCACGCTGCCCgcgc                                                                                           >  1:210597/1‑140 (MQ=255)
                                       aaCGCCACACCTGCTGCCCCGCCTTCCGCCAGTCGCACATTAAATTGCAAAATATTCATTTAATTACCTTCACTCGCGCTTTTTCTCCCACCGCCAGCGCGTTGTCCGGGACAGAGTCGAGCACCACGCTGCCCGCGccc                                                                                         >  1:37287/1‑140 (MQ=255)
                                          gCCACACCTGCTGCCCCGCCTTCCGCCAGTCGCACATTAAATTGCAAAATATTCATTTAATTACCTTCACTCGCGCTTTTTCTCCCACCACCAGCGCGTTGTCCGGGACAGAGTCGAGCACCACGCTGCCCGCGCCCACg                                                                                      >  1:394229/1‑140 (MQ=255)
                                              cacCTGCTGCCCCGCCTTCCGCCAGTCGCACATTAAATTGCAAAATATTCATTTAATTACCTTCACTCGCGCTTTTTCTCCCACCGCCAGCGCGTTGTCCGGGACAGAGTCGAGTACCACGCTGACCGCGCCCACGGTGa                                                                                  >  1:355193/1‑140 (MQ=255)
                                              cacCTGCTGCCCCGCCTTCCGCCAGTCGCACATTAAATTGCAAAATATTCATTTAATTACCTTCACTCGCGCTTTTTCTCCCACCGCCAGCGCGTTGTCCGGGACAGAGTCGAGCACCACGCTGCCCGCGCCCACGGTGa                                                                                  <  1:659729/140‑1 (MQ=255)
                                              cacCTGCTGCCCCGCCGTCCGCCAGTCGCACGTTAAATTGCAAAATATTCATTTAATTACCTTCACTCGCGCTTTTTCTCCCACCGCCAGCGCGTTGTCCGGGACAGAGTCGAGCACCACGCTGCCCGCGCCCACGGTGa                                                                                  <  2:328168/140‑1 (MQ=255)
                                                ccTGCTGCCCCGCCTTCCGCCAGTCGCACATTAAATTGCAAAATATTCATTTAATTACCTTCACTCGCGCTTTTTCTCCCACCGCCAGCGCGTTGTCCGGGACAGAGTCGAGCACCACGCTGCCCGCGCCCACGGTGACg                                                                                <  2:677380/140‑1 (MQ=255)
                                                  tgctgcCCCGCCTTCGGCCAGTTGCCCATTAAATTGCAAAATATTCATTTAATTACCTTCACTCGCGCTTTTTCTCCCACCGCCAGCGCGTTGTCCGGGACAGAGTCGAGCACCACGCTGCCCGCGCCCACGGTGACgtt                                                                              <  1:320951/140‑1 (MQ=255)
                                                    ctgcCCCGCCTTCCGCCAGTCGCACATTAAATTGCAAAATATTCATTTAATTACCTTCACTCGCGCTTTTTCTCCCACCGCCAGCGCGTTGTCCGGGACAGAATCGAGCAACACCCTCCCCGCGCCCACGGTGACgttgt                                                                            >  1:300088/1‑140 (MQ=255)
                                                           gCCTTCCGCCAGTCGCACATTAAATTGCAAAATATTCATTTAATTACCTTCACTCGCGCTTTTTCTCCCACCGCCAGCGCGTTGTCCGGGACAGAGTCGAGCACCACGCTGCCCGCGCCCACGGTGACGTTGTTACCAAg                                                                     >  2:68917/1‑140 (MQ=255)
                                                               tCCGCCAGTCGCACATTCAATTCCAAAATATTCATCTAACTACCTTCACTCGCGCTTTTTCTCCCACCGCCAGCGCGTTGTCCGGGACAGAGTCGAGCACCACGCTGCCCGCGCCCACGGTGACGTTGTTACCAAGCGTg                                                                 <  2:210597/140‑1 (MQ=255)
                                                                      gTCGCACATTAAATTGCAAAATATTCATTTAATTACCTTCACTCGCGCTTTTTCTCCCACCGCCAGCGCGTTGTCCGGGACAGAGTCGAGCACCACGCTGCCCGCGCCCACGGTGACGTTGTTACCAAGCGTGATATCAc                                                          >  1:411672/1‑140 (MQ=255)
                                                                       tCGCACATTAAATTGCAAAATATTCATTTAATTACCTTCACTCGCGCTTTTTCTCCCACCGCCAGCGCGTTGTCCGGGACAGAGTCGAGCACCACGCTGCCCGCGCCCACGGTGACGTTGTTACCAAGCGTGATATCAcc                                                         <  2:411672/140‑1 (MQ=255)
                                                                       tCGCACATTAAATTGCAAAATATTCATTTAATTACCTTCACTCGCGCTTTTTCTCCCACCGCCAGCGCGTTGTCCGGGACAGAGTCGAGCACCACGCTGCACGCCCCCACGGTGACGATGTTACCAAGCGTGATATCAcc                                                         >  2:246320/1‑140 (MQ=255)
                                                                                aaaTTGCAAAATATTCAGTTAATTACCTTCAGCCGTGCTTTTTCTCCCACCGCCAGCGCGTTGTCCGGGACAGAGTCGAGCACCACGCTGCCCGCGCCCACGGTGACGTTGTTACCAAGCGTGATATCACCAAGaataat                                                <  2:295586/140‑1 (MQ=255)
                                                                                       aaaaTATTCATTTAACTACCTTTACCCGCTCTTTTTCCCCCACCGCCAGCGCGTTGTCCGGGACAGAGTCGAGCACCACGCTGCCCGCGCCCACGGTGACGTTGTTACCAAGCGTGATATCACCAAGAATAATGACGTTg                                         <  1:215503/140‑1 (MQ=255)
                                                                                          atatTCATTTAATTACCTTCACTCGCGCTTTTTCTCCCACCACCAGCGCGTTGTCCGGGACAGAGTCGAGCACCACGCTGCCCGCGCCCACGGTGACGTTGTTACCAAGCGTGATATCACCAAGAATAATGACGTTGGCa                                      <  2:394229/140‑1 (MQ=255)
                                                                                                     aTTACCTTCACTCGCGCTTTTTCTCCCACCGCCAGCGCGTTGTCCGGGACAGAGTCGAGCACCACGCTGCCCGCGCCCACGGTGACGTTGTTACCAAGCGTGATATCACCAAGAATAATGACGTTGGCACCGAGTTCGAc                           >  1:60127/1‑140 (MQ=255)
                                                                                                            tCACTCGCGCTTTTTCTCCCACCGCCAGCGCGTTGTCCGGGACAGAGTCGAGCACCACGCTGCCCGCGCCCACGGTGACGTTGTTACCAAGCGTGATATCACAAAGAATAATGAAGTCGGCACCTAGTTCGACGCCGTTg                    >  1:249433/1‑140 (MQ=255)
                                                                                                              aCTCGCGCTTTTTCTCCCACCGCCAGCGCGTTGTCCGGGACAGAGTCGAGCACCACGCTGCCCGCGCCCACGGTGACGTTGTTACCAAGCGTGATATCACCAAGAATAATGACGTTGGCACCGAGTTCGACGCCGTTGcc                  <  2:37287/140‑1 (MQ=255)
                                                                                                               cTCGCGCTTTTTCTCCCACCGCCAGCGCGTTGTCCGGGACAGAGTCGAGCACCACGCTGCCCGCGCCCACGGTGACGTTGTTACCAAGCGTGATATCACCAAGAATAATGACGTTGGCACCGAGTTCGACGCCGTTGCCa                 <  1:146915/140‑1 (MQ=255)
                                                                                                                           ctcCCACCTCCAGCGCTCTGTCCGGGACAGACTCGAGCACCACGCTGCCCGCGCCCACGGTGACGTTGTTACCAAGCGTGATATCACCAAGAATAATGACGTTGGCACCGAGTTCGACGCCGTTGCCAATGTGTGGACAt     <  2:355193/140‑1 (MQ=255)
                                                                                                                              ccaccgccaGCGCGTTGTCCGGCACAGAGTCGAGCGCCAGGCTGCCCGCGCCCACGGTGACGTTGTTACCAAGCGTGATATCACCAAGAATAATGACGTTGGCACCGAGTTCGACGCCGTTGCCAATGTGTGGACATGcc  <  2:249433/140‑1 (MQ=255)
                                                                                                                                   |                                                                                                                                      
GACGCCAGCCCCTGTTGCAGCGCACGCTGGTGGAGATCTAACGCCACACCTGCTGCCCCGCCTTCCGCCAGTCGCACATTAAATTGCAAAATATTCATTTAATTACCTTCACTCGCGCTTTTTCTCCCACCACCAGCGCGTTGTCCGGGACAGAGTCGAGCACCACGCTGCCCGCGCCCACGGTGACGTTGTTACCAAGCGTGATATCACCAAGAATAATGACGTTGGCACCGAGTTCGACGCCGTTGCCAATGTGTGGACATGCC  >  NC_000913/2131973‑2132238

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: