Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,977,701 T→A 6.3% P59P (CCT→CCA wecF → TDP‑Fuc4NAc:lipidIIFuc4NAc transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,977,7010TA 6.3% 222.8 / 4.8 63P59P (CCT→CCAwecFTDP‑Fuc4NAc:lipidIIFuc4NAc transferase
Reads supporting (aligned to +/- strand):  ref base T (29/30);  new base A (2/2);  total (31/32)
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.71e-01

CTCACCATAACCGAACCGTTTTGCGGTTTTTCAATGACGCGCTGGCCGCGACGAGCGAGCACGCGCGCGAGTTTATGGTTGTTGGCAAGGACGACGGCTTAAGTGATAGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGGAAAAAATCGCTGGCGGAAGCGGTCATCGCGAAAGCAAAAGCTAACCGTCAGCAGCGTTTTTTCTTCCACGGTCAGTTCAATCCCACACTGTGGCTGGCTCTGCTGAGTGGTGGCATTAAGCCCAGCCAGTTTTTCTG  >  NC_000913/3977562‑3977841
                                                                                                                                           |                                                                                                                                            
cTCACCATAACCGAACCGTTTTGCGGTTTTTCAATGACGCGCTGGCCGCGACGAGCGAGCACGCGCGCGAGTTTATGGTTGTTGGCAAGGACGACGGCTTAAGTGATAGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGGaaaaa                                                                                                                                      >  1:90452/1‑148 (MQ=255)
 tCACCATAACCGAACCGTTTTGCGGTTTTTCAATGACGCGCTGGCCGCGACGAGCGAGCACGCGCGCGAGTTTATGGTTGTTGGCAAGGACGACGGCTTAAGTGATAGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGGaaaaaa                                                                                                                                     <  4:793716/148‑1 (MQ=255)
               ccgTTTTGCGGTTTTTCAATGACGCGCTGGCCGCGACGAGCGAGCACGCGCGCGAGTTTATGGTTGTTGGCAAGGACGACGGCTTAAGTGATAGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGGAAAAAATCGCTGGCGGAAGc                                                                                                                       <  4:373307/148‑1 (MQ=255)
                  ttttGCGGTTTTTCAATGACGCGCTGGCCGCGACGAGCGAGCACGCGCGCGAGTTTATGGTTGTTGGCAAGGACGACGGCTTAAGTGATAGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGGAAAAAATCGCTGGCGGAAGCGGt                                                                                                                    <  2:90452/148‑1 (MQ=255)
                                 aTGACGCGCTGGCCGCGACGAGCGAGCACGCGCGCGAGTTTATGGTTGTTGGCAAGGACGACGGCTTAAGTGATAGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGGAAAAAATCGc                                                                                                                                 >  2:369110/1‑120 (MQ=255)
                                 aTGACGCGCTGGCCGCGACGAGCGAGCACGCGCGCGAGTTTATGGTTGTTGGCAAGGACGACGGCTTAAGTGATAGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGGAAAAAATCGc                                                                                                                                 <  1:369110/120‑1 (MQ=255)
                                 aTGACGCGCTGGCCGCGACGAGCGAGCACGCGCGCGAGTTTATGGTTGTTGGCAAGGACGACGGCTTAAGTGATAGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCAGGGaaaaa                                                                                                                                      <  4:636846/115‑1 (MQ=255)
                                 aTGACGCGCTGGCCGCGACGAGCGAGCACGCGCGCGAGTTTATGGTTGTTGGCAAGGACGACGGCTTAAGTGATAGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCAGGGaaaaa                                                                                                                                      >  3:636846/1‑115 (MQ=255)
                                         cTGGCCGCGACGAGCGAGCACGCGCGCGAGTTTATGGTTGTTGGCAAGGACGACGGCTTAAGTGATAGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGGAAAAAATCGCTGGCGGAAGCGGTCATCGCGAAAGCAAAAGCTAAcc                                                                                             >  2:946457/1‑148 (MQ=255)
                                                           cACGCGCGCGAGTTTATGGTTGTTGGCAAGGACGACGGCTTAAGTGATAGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGGAAAAAATCGc                                                                                                                                 >  4:274900/1‑94 (MQ=255)
                                                           cACGCGCGCGAGTTTATGGTTGTTGGCAAGGACGACGGCTTAAGTGATAGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGGAAAAAATCGc                                                                                                                                 <  3:274900/94‑1 (MQ=255)
                                                           cACGCGCGCGAGTTTATGGTTGTTGGCAAGGACGACGGCTTAAGTGATAGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGGAAAAAATCGCTGGCGGAAGCGGTCATCGCGAAAGCAAAAGCTAAc                                                                                              <  2:254214/129‑1 (MQ=255)
                                                           cACGCGCGCGAGTTTATGGTTGTTGGCAAGGACGACGGCTTAAGTGATAGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGGAAAAAATCGCTGGCGGAAGCGGTCATCGCGAAAGCAAAAGCTAAc                                                                                              >  1:254214/1‑129 (MQ=255)
                                                             cgcgcgcgAGTTTATGGTTGTTGGCAAGGACGACGGCTTAAGTGATAGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGGAAAAAATCGCTGGCGGAAGCg                                                                                                                      <  3:1100601/103‑1 (MQ=255)
                                                             cgcgcgcgAGTTTATGGTTGTTGGCAAGGACGACGGCTTAAGTGATAGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGGAAAAAATCGCTGGCGGAAGCg                                                                                                                      >  4:1100601/1‑103 (MQ=255)
                                                                gcgcgAGTTTATGGTTGTTGGCAAGGACGACGGCTTAAGTGATAGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGGAAAAAATCGCTGGCGGAAGCGGTCATCGCGAAAGCAAAAGCTAACCGTCAGCAGCGTTTTTTCTTCCAc                                                                      <  1:946457/148‑1 (MQ=255)
                                                                            ggTTGTTGGCAAGGACGACGGCTTAAGTGATAGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGGAAGAAATCGCTGGCGGAAGCGGTCATCGCGAAAGCAAAAGCTAACCGTCAGCAGCGTTTTttctt                                                                          <  1:688743/132‑1 (MQ=255)
                                                                            ggTTGTTGGCAAGGACGACGGCTTAAGTGATAGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGGAAGAAATCGCTGGCGGAAGCGGTCATCGCGAAAGCAAAAGCTAACCGTCAGCAGCGTTTTttctt                                                                          >  2:688743/1‑132 (MQ=255)
                                                                             gttgttGGCAAGGACGACGGCTTAAGTGATAGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGGAAAAAATCGc                                                                                                                                 >  1:856426/1‑76 (MQ=255)
                                                                             gttgttGGCAAGGACGACGGCTTAAGTGATAGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGGAAAAAATCGc                                                                                                                                 <  2:856426/76‑1 (MQ=255)
                                                                                  tgGCAAGGACGACGGCTTAAGTGATAGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTg                                                                                                                                             >  1:1082675/1‑59 (MQ=255)
                                                                                  tgGCAAGGACGACGGCTTAAGTGATAGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTg                                                                                                                                             <  2:1082675/59‑1 (MQ=255)
                                                                                    gCAAGGACGACGGCTTAAGTGATAGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGGAAAAAATCGCTGGCGGAAGCGGTCATCGCGaa                                                                                                           <  1:1072489/91‑1 (MQ=255)
                                                                                    gCAAGGACGACGGCTTAAGTGATAGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGGAAAAAATCGCTGGCGGAAGCGGTCATCGCGaa                                                                                                           >  2:1072489/1‑91 (MQ=255)
                                                                                     cAAGGACGACGGCTTAAGTGATAGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGGAAAAAATCGCTGGCGGAAGCGGTCATCGCGAAAGCAAAAGCTAACCGTcagcag                                                                                     >  3:116862/1‑112 (MQ=255)
                                                                                     cAAGGACGACGGCTTAAGTGATAGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGGAAAAAATCGCTGGCGGAAGCGGTCATCGCGAAAGCAAAAGCTAACCGTcagcag                                                                                     <  4:116862/112‑1 (MQ=255)
                                                                                     cAAGGACGACGGCTTAAGTGATAGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGGAAAAAATCGCTGGCGGAAGCGGTCATCGCGAAAGCAAAAGCTAACCGTCAGCAGCGTTTTTTCTTCCACGGTCAGTTCa                                                            >  1:1032988/1‑137 (MQ=255)
                                                                                     cAAGGACGACAGCTTAAGTGATAGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGGAAAAAATCGCTGGCGGAAGCGGTCATCGCGAAAGCAAAAGCTAACCGTCAGCAGCGTTTTTTCTTCCACGGTCAGTTCa                                                            <  2:1032988/137‑1 (MQ=255)
                                                                                          acgacgGCTTAAGTGATAGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGGAAAAAATCGCTGGCGGAAGCGGTCATCGCGAAAGCAAAAGCTAACCGTCAGCAGCGTTTTTTCTTCCACGGTCAGTTCAATc                                                         <  4:259624/135‑1 (MQ=255)
                                                                                          acgacgGCTTAAGTGATAGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGGAAAAAATCGCTGGCGGAAGCGGTCATCGCGAAAGCAAAAGCTAACCGTCAGCAGCGTTTTTTCTTCCACGGTCAGTTCAATc                                                         >  3:259624/1‑135 (MQ=255)
                                                                                             acgGCTTAAGTGATAGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGGAAAAAATCGCTGGCGGaa                                                                                                                         >  3:624796/1‑68 (MQ=255)
                                                                                             acgGCTTAAGTGATAGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGGAAAAAATCGCTGGCGGaa                                                                                                                         <  4:624796/68‑1 (MQ=255)
                                                                                                   tAAGTGATAGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGGAAAAAATCGCTGGCGGAAGCGGTCATCGCGAAAGCaaa                                                                                                     <  2:20090/82‑1 (MQ=255)
                                                                                                   tAAGTGATAGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGGAAAAAATCGCTGGCGGAAGCGGTCATCGCGAAAGCaaa                                                                                                     >  1:20090/1‑82 (MQ=255)
                                                                                                   tAAGTGATAGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGGAAAAAATCGCTGGCGGAAGCGGTCATCGCGAAAGCAAAAGCTaa                                                                                               <  4:970149/88‑1 (MQ=255)
                                                                                                   tAAGTGATAGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGGAAAAAATCGCTGGCGGAAGCGGTCATCGCGAAAGCAAAAGCTaa                                                                                               >  3:970149/1‑88 (MQ=255)
                                                                                                       tGATAGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGGAAAAAATCGCTGGCGGAAGCg                                                                                                                      >  2:784407/1‑61 (MQ=255)
                                                                                                       tGATAGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGGAAAAAATCGCTGGCGGAAGCg                                                                                                                      <  1:784407/61‑1 (MQ=255)
                                                                                                       tGATAGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGGAAAAAATCGCTGGCGGAAGCGGTCATCGCGAAAGCAAAAGCTAACCGTCAGCAGCGTTTTTTCTTCCACGGTCAGTTCAATCCCACACTGTGGCTGGCTCTGCTGa                                 <  2:168740/146‑1 (MQ=255)
                                                                                                       tGATAGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGGAAAAAATCGCTGGCGGAAGCGGTCATCGCGAAAGCAAAAGCTAACCGTCAGCAGCGTTTTTTCTTCCACGGTCAGTTCAATCCCACACTGTGGCTGGCTCTGCTGa                                 >  1:168740/1‑146 (MQ=255)
                                                                                                       tGATAGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGGAAAAAATCGCTGGCGGAAGCGGTCATCGCGAAAGCAAAAGCTAACCGTCAGCAGCGTTTTTTCTTCCACGGTCAGTTCAATCCCACACTGTGGCTGGCTCTGCTGAgt                               >  3:925864/1‑148 (MQ=255)
                                                                                                        gATAGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGGAAAAAATCGCTGGCGGAAGCGGTCATCGCGAAAGCaaa                                                                                                     >  1:623383/1‑77 (MQ=255)
                                                                                                        gATAGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGGAAAAAATCGCTGGCGGAAGCGGTCATCGCGAAAGCaaa                                                                                                     <  2:623383/77‑1 (MQ=255)
                                                                                                          tAGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTg                                                                                                                                             <  1:116758/35‑1 (MQ=255)
                                                                                                          tAGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTg                                                                                                                                             <  1:256161/35‑1 (MQ=255)
                                                                                                          tAGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTg                                                                                                                                             >  2:256161/1‑35 (MQ=255)
                                                                                                          tAGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTg                                                                                                                                             >  2:116758/1‑35 (MQ=255)
                                                                                                            gCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGGAAAAAATCg                                                                                                                                  <  1:753856/44‑1 (MQ=255)
                                                                                                            gCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGGAAAAAATCg                                                                                                                                  >  2:753856/1‑44 (MQ=255)
                                                                                                             cTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGGAAAAAATCGCTGGCGGAAGCGGTCATCGCGa                                                                                                            <  1:195029/65‑1 (MQ=255)
                                                                                                             cTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGGAAAAAATCGCTGGCGGAAGCGGTCATCGCGa                                                                                                            >  2:195029/1‑65 (MQ=255)
                                                                                                                    gcgcTTTCTGTGCAATTTTTCCCTGGGAAAAAATCGCTGGCGGAAGCGGTCATCGCGAAAGCAAAAGCTTACCGTCAGCAGCGTTTTTTCTTCCACGGTCAGTTCAATCCCACACTGTGGCTGGCTctg                                     <  3:662696/129‑1 (MQ=255)
                                                                                                                    gcgcTTTCTGTGCAATTTTTCCCTGGGAAAAAATCGCTGGCGGAAGCGGTCATCGCGAAAGCAAAAGCTAACCGTCAGCAGCGTTTTTTCTTCCACGGTCAGTTCAATCCCACACTGTGGCTGGCTctg                                     <  3:999726/129‑1 (MQ=255)
                                                                                                                    gcgcTTTCTGTGCAATTTTTCCCTGGGAAAAAATCGCTGGCGGAAGCGGTCATCGCGAAAGCAAAAGCTAACCGTCAGCAGCGTTTTTTCTTCCACGGTCAGTTCAATCCCACACTGTGGCTGGCTctg                                     >  4:999726/1‑129 (MQ=255)
                                                                                                                    gcgcTTTCTGTGCAATTTTTCCCTGGGAAAAAATCGCTGACGGAAGCGGTCATCGCGAAAGCAAAAGCTAACCGTCAGCAGCGTTTTTTCTTCCACGGTCAGTTCAATCCCACACTGTGGCTGGCTctg                                     >  4:662696/1‑129 (MQ=255)
                                                                                                                       cTTTCTGTGCAATTTTTCCCTGGGAAAAAATCGCTGGCGGAAGCGGTCATCGCGAAAGCAAAAGCTAACCGTCAGCAGCGTTTTTTCTTCCACGGTCAGTTCAATCCCACACTGTGGCTGGCTCTGCTGAGTGGTGGCATTAAGCcca               >  3:581577/1‑148 (MQ=255)
                                                                                                                          tCTGTGCAATTTTTCCCTGGGAAAAAATCGCTGGCGGAAGCGGTCATCGCGAAAGCAAAAGCTAACCGTCAGCAGCGTTTTTTCTTCCACGGTCAGTTCAATCCCACACTGTGGCTGGCTCTGCTGAGTGGTGGCATTAAGCccagcc            <  4:581577/148‑1 (MQ=255)
                                                                                                                             gtgCAATTTTTCCCTGGGAAAAAATCGCTGGCGGAAGCGGTCATCGCGAAAGCAAAAGCTAACCGTCAGCAGCGTTTTTTCTTCCACGGTCAGTTCAATCCCACACTGTGGCTGGCTCTGCTGAGTGGTGGCATTAAGCCCAGCCAGt         >  2:350238/1‑148 (MQ=255)
                                                                                                                               gCAATTTTTCCCTGGGAAAAAATCGCTGGCGGAAGCGGTCATCGCGaa                                                                                                           >  2:891054/1‑48 (MQ=255)
                                                                                                                               gCAATTTTTCCCTGGGAAAAAATCGCTGGCGGAAGCGGTCATCGCGaa                                                                                                           <  1:891054/48‑1 (MQ=255)
                                                                                                                               cgattttttCCCAGGGAAAAATCGCTGGCGGAAGCGGTCATCGCGAAAGCAAAAGCTAACCGTCAGCAGCGTTTTTTCTTCCACGGTCAGTTCAATCCCACACTGTGGCTGGCTCTGCTGAGTGGTGGCAt                       >  3:939586/5‑131 (MQ=255)
                                                                                                                               cgattttttCCCAGGGAAAAATCGCTGGCGGAAGCGGTCATCGCGAAAGCAAAAGCTAACCGTCAGCAGCGTTTTTTCTTCCACGGTCAGTTCAATCCCACACTGTGGCTGGCTCTGCTGAGTGGTGGCAt                       <  4:939586/127‑1 (MQ=255)
                                                                                                                                    ttttCCCTGGGAAAAAATCGCTGGCGGAAGCGGTCATCGCGAAAGCAAAAGCTAACCGTCAGCAGCGTTTTTTCTTCCACGGTCAGTTCAATCCCACACTGTGGCTGGCTCTGCTGAGTGGTGGCATTAAGCCCAGCCAGTTTTTCTg  <  1:350238/148‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                            
CTCACCATAACCGAACCGTTTTGCGGTTTTTCAATGACGCGCTGGCCGCGACGAGCGAGCACGCGCGCGAGTTTATGGTTGTTGGCAAGGACGACGGCTTAAGTGATAGCTGTCCGGCGCTTTCTGTGCAATTTTTCCCTGGGAAAAAATCGCTGGCGGAAGCGGTCATCGCGAAAGCAAAAGCTAACCGTCAGCAGCGTTTTTTCTTCCACGGTCAGTTCAATCCCACACTGTGGCTGGCTCTGCTGAGTGGTGGCATTAAGCCCAGCCAGTTTTTCTG  >  NC_000913/3977562‑3977841

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