Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,114,502 C→T *465* (TAG→TAA wcaM ← colanic acid biosynthesis protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,114,5020CT98.3% 189.0 / ‑6.6 59*465* (TAG→TAAwcaMcolanic acid biosynthesis protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (28/30);  minor base G (1/0);  total (29/30)
Fisher's exact test for biased strand distribution p-value = 4.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTCGATTTTTTTCGCCAGCATCAATTACCCTGAATTGATTACTGAATTACTTGTGATGTTACGCCGCTTCGTTGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATCTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATGGTGATTAATCCTGTCGATATCGACGGAGCTTTGTCCTTTTTCATTCACCGCATGAAC  >  NC_000913/2114369‑2114612
                                                                                                                                     |                                                                                                              
ttCGATTTTTTTCGCCAGCATCAATTACCCTGAATTGATTACTGAATTACTTGTGATGTTACGCCGCTTCGTTGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTc                                                                                                           >  1:264126/1‑139 (MQ=255)
  cGATTTTTTTCGCCAGCATCAATTACCCTGAATTGATTACTGAATTACTTGTGATGTTACGCCGCTTCGTTGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTccc                                                                                                         <  1:1306816/139‑1 (MQ=255)
   gATTTTTTTCGCCAGCATCAATTACCCTGAATTGATTACTGAATTACTTGTGATGTTACGCCGCTTCGTTGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCt                                                                                                        <  2:264126/139‑1 (MQ=255)
            cGCCAGCATCAATTACCCTGAATTGATTACTGAATTACTTGTGATGTTACGCCGCTTCGTTGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCgg                                                                                               <  1:1417505/139‑1 (MQ=255)
                           cccTGAATTGATTACTGAATTACTTGTGATGTTACGCCGCTTCGTTGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAAt                                                                                <  2:237485/139‑1 (MQ=255)
                            ccTGAATTGATTACTGAATTACTTGTGATGTTACGCCGCTTCGTTGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAAt                                                                                           <  2:1217743/127‑1 (MQ=255)
                            ccTGAATTGATTACTGAATTACTTGTGATGTTACGCCGCTTCGTTGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAAt                                                                                           >  1:1217743/1‑127 (MQ=255)
                               gAATTGATTACTGAATTACTTGTGATGTTACGCCGCTTCGTTGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAAtcttc                                                                            >  1:369350/1‑139 (MQ=255)
                                 aTTGATTACTGAATTACTTGTGATGTTACGCCGCTTCGTTGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTc                                                                                                       >  2:1175327/1‑110 (MQ=255)
                                 aTTGATTACTGAATTACTTGTGATGTTACGCCGCTTCGTTGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTc                                                                                                       <  1:1175327/110‑1 (MQ=255)
                                     aTTACTGAATTACTTGTGATGTTACGCCGCTTCGTTGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCt                                                                                         <  2:978737/120‑1 (MQ=255)
                                     aTTACTGAATTACTTGTGATGTTACGCCGCTTCGTTGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCt                                                                                         >  1:978737/1‑120 (MQ=255)
                                               tACTTGTGATGTTACGCCGCTTCGTTGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAAtctt                                                                             >  2:518815/1‑122 (MQ=255)
                                               tACTTGTGATGTTACGCCGCTTCGTTGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAAtctt                                                                             <  1:518815/122‑1 (MQ=255)
                                                 cTTGTGATGTTACGCCGCTTCGTTGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATg                                                           <  1:1445117/138‑1 (MQ=255)
                                                 cTTGTGATGTTACGCCGCTTCGTTGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATg                                                           >  2:1445117/1‑138 (MQ=255)
                                                    gtgATGTTACGCCGCTTCGTTGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCt                                                                                                        <  1:781762/90‑1 (MQ=255)
                                                    gtgATGTTACGCCGCTTCGTTGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCt                                                                                                        >  2:781762/1‑90 (MQ=255)
                                                      gATGTTACGCCGCTTCGTTGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATg                                                           >  2:624281/1‑133 (MQ=255)
                                                      gATGTTACGCCGCTTCGTTGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATg                                                           <  1:624281/133‑1 (MQ=255)
                                                                  cTTCGTTGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTTCTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATGGTGATTAATCCTGTCGat                                         >  2:206118/1‑139 (MQ=255)
                                                                  cTTCGTTGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTccc                                                                                                         <  2:783881/75‑1 (MQ=255)
                                                                  cTTCGTTGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTccc                                                                                                         >  1:783881/1‑75 (MQ=255)
                                                                  cTTCGTTGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTaa                                                                                       >  1:676935/1‑93 (MQ=255)
                                                                  cTTCGTTGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTaa                                                                                       <  2:676935/93‑1 (MQ=255)
                                                                  cTTCGTTGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTtctc                                                                          <  2:1298244/106‑1 (MQ=255)
                                                                  cTTCGTTGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTtctc                                                                          >  1:1298244/1‑106 (MQ=255)
                                                                       ttGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAAtcttc                                                                            >  2:1358334/1‑99 (MQ=255)
                                                                       ttGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAAtcttc                                                                            <  1:1358334/99‑1 (MQ=255)
                                                                       ttGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCcac                                                                       <  1:796166/104‑1 (MQ=255)
                                                                       ttGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCcac                                                                       >  2:796166/1‑104 (MQ=255)
                                                                               tGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATGGTGATTAATCCTGTCGATATCGACGGAGCtt                            >  2:1413780/1‑139 (MQ=255)
                                                                                    tAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAAt                                                                                           <  1:1344351/71‑1 (MQ=255)
                                                                                    tAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAAt                                                                                           >  2:1344351/1‑71 (MQ=255)
                                                                                         ttGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAata                                                             <  1:549501/96‑1 (MQ=255)
                                                                                         ttGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAata                                                             >  2:549501/1‑96 (MQ=255)
                                                                                                    tATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATGGTGATTAATCCTGTCGATATCGACGGAGCTTTGTCCTTTTTCATTCACCGCa       >  2:1126106/1‑139 (MQ=255)
                                                                                                      tGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCt                                                                                                        <  2:1063509/40‑1 (MQ=255)
                                                                                                      tGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCt                                                                                                        >  1:1063509/1‑40 (MQ=255)
                                                                                                      tGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATc                                                                                          >  1:1087110/1‑54 (MQ=255)
                                                                                                      tGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATc                                                                                          <  2:1087110/54‑1 (MQ=255)
                                                                                                      tGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATg                                                           >  2:1474817/1‑85 (MQ=255)
                                                                                                      tGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATg                                                           <  1:1474817/85‑1 (MQ=255)
                                                                                                       gACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATGGTGATTAATCCTGTCGATATCGACGGAGCTTTGTCCt                      >  2:906501/1‑121 (MQ=255)
                                                                                                       gACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATGGTGATTAATCCTGTCGATATCGACGGAGCTTTGTCCt                      <  1:906501/121‑1 (MQ=255)
                                                                                                         cTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATGGTGATTAATCCTGTCGATATCGACGGAGCTTTGTCCTTTTTCATTCACCGCATGAAc  <  2:369350/139‑1 (MQ=255)
                                                                                                         cTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATGGTGATTAATCCTGTCGATATCGACGGAGCTTTGTCCTTTTTCATTCACCGCATGAAc  <  1:1126106/139‑1 (MQ=255)
                                                                                                            taTTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATGGTGATTAATCCTGTCGATATCGACGGAGCTTTGTCCTTTttctttca            >  1:926721/2‑126 (MQ=255)
                                                                                                            cATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATGGTGATTAATCCTGTCGATATCGACGGAGCTTTGTCCTTTttcattca            <  2:926721/126‑1 (MQ=255)
                                                                                                                   ccAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATGGTGATTAATCCTGTCGatat                                       >  2:225194/1‑92 (MQ=255)
                                                                                                                   ccAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATGGTGATTAATCCTGTCGatat                                       <  1:225194/92‑1 (MQ=255)
                                                                                                                            gTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTa                                                                                  >  1:1115365/1‑40 (MQ=255)
                                                                                                                            gTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTa                                                                                  <  2:813907/40‑1 (MQ=255)
                                                                                                                            gTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTa                                                                                  >  1:813907/1‑40 (MQ=255)
                                                                                                                            gTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTa                                                                                  <  2:1115365/40‑1 (MQ=255)
                                                                                                                            gTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTtctc                                                                          <  2:1030016/48‑1 (MQ=255)
                                                                                                                            gTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACAt                                                                     <  2:105520/53‑1 (MQ=255)
                                                                                                                            gTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACAt                                                                     >  1:105520/1‑53 (MQ=255)
                                                                                                                            gTCGCAAATGTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTtctc                                                                          >  1:1030016/1‑48 (MQ=255)
                                                                                                                                     |                                                                                                              
TTCGATTTTTTTCGCCAGCATCAATTACCCTGAATTGATTACTGAATTACTTGTGATGTTACGCCGCTTCGTTGTGGATTGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATCTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATGGTGATTAATCCTGTCGATATCGACGGAGCTTTGTCCTTTTTCATTCACCGCATGAAC  >  NC_000913/2114369‑2114612

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