Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,971,996 T→A 12.7% L246Q (CTG→CAG)  wecC → UDP‑N‑acetyl‑D‑mannosamine dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,971,9960TA12.7% 84.2 / 7.7 41L246Q (CTG→CAG) wecCUDP‑N‑acetyl‑D‑mannosamine dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (12/21);  new base A (3/2);  total (16/25)
Fisher's exact test for biased strand distribution p-value = 3.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.51e-01

GAGTGTGTCGTCACTAACTCGCGGACGGCGGAAATGTGTAAGCTCACCGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGTTTATCGTGGCACAGAACCCCCAGCAGGCGCGGCTTATCCGTACCGCGCGCGAAGTGAACGA  >  NC_000913/3971854‑3972137
                                                                                                                                              |                                                                                                                                             
gAGTGTGTCGTCACTAACGCGCGGACGGCGGAAATGTGTAAGCTCCCCGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCt                                                                                                                                               <  1:202104/143‑1 (MQ=255)
       tcgtcACTAACTCGCGGACGGCGGAAATGTGTAAGCTCACCGAAACACGCTTCGGCGATGTGGATATCGCTTTTGCTAATGAACTGTCGCCGCTATGACCCCGTCAGGGGATTGCCGGCTGGGCACTGGTTCGCCTGGCGact                                                                                                                                        >  2:264569/1‑141 (MQ=255)
             cTAACTCGCGGACGGCGGAAATGTGTAAGCTCACCGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACCGATTCGCCAGGCGAATcagttc                                                                                                                                  <  1:505668/143‑6 (MQ=255)
              tAACTCGCGGACGGCGGAAATGTGTAAGCTCACCGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCAcc                                                                                                                                 <  2:50850/143‑1 (MQ=255)
                      ggACGGCGGAAATGTGTAAGATCACCGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGACGCTGATTTGTGCCGATCAGGGGACTCAGACCTGAGAACTGATACGCCTGGCGCATCGCCACGCTcgcggc                                                                                                                         >  2:33701/1‑141 (MQ=255)
                          ggcggAAATGTGTAAGCTCACCGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAata                                                                                                                     <  1:191775/143‑1 (MQ=255)
                              gAAATGTGTACGCTCACCGAAAACAGCTGCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATAttct                                                                                                                 <  2:295005/143‑1 (MQ=255)
                               aattgtgTAAGCTCACCGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATAttctt                                                                                                                <  1:167766/140‑1 (MQ=255)
                               aaaggagtAAGCTCACCGAAAACAGCTTCCGCTATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGAGCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATAttctt                                                                                                                <  2:555966/137‑1 (MQ=255)
                               aaaTGTGTAAGCTCACCGAAAACAGCTTCCGCGATGTGAATATCGCTGTTGCTAATGAGTTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATAttctt                                                                                                                <  2:242350/143‑1 (MQ=255)
                                           tCACCGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGcctggccct                                                                                                    <  1:494245/143‑1 (MQ=255)
                                              ccGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGACCTGATCCGCAGGGCGGATCGTCATCCTTGCGTTAATATTCTTCAGcctagccgcggc                                                                                                 >  2:438870/1‑138 (MQ=255)
                                              ccGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTGACTATTCTTCAGcctggccctggc                                                                                                 >  2:610178/1‑143 (MQ=255)
                                                gAAACCAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAGTTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCGTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGt                                                                                               <  2:117372/143‑1 (MQ=255)
                                                 aaaaCAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGGTACTATGCTTCAGCCTGGCCCTGGCCTg                                                                                              >  2:100853/1‑143 (MQ=255)
                                                  aaaCGGCTTCCGCGATGTGAATATCGCTTTTGCGAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATGAACGTCTGGGAACTGATTCGCCGGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTgg                                                                                             <  2:417488/143‑1 (MQ=255)
                                                     caggctCCGCGTCGTGAAGATCGCGTGTGCTAATGAATTGTCGCTGATTTGTGCCGATCGGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTTAGTGCGGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCg                                                                                          <  2:47663/138‑1 (MQ=255)
                                                     cAGGTTCCGCGATGCGAATATCGCTTTGGCTAATGAGTTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGGCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCg                                                                                          <  1:656210/143‑1 (MQ=255)
                                                         ttCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGCGGGCGGTCa                                                                                      >  1:410766/1‑143 (MQ=255)
                                                         ttCCGCGAGGTGAATCTCGCTTTTGCTAATGAATTGTCGCTGCTTTGTGCCGATCAGGGGATTAACGTCTGGAAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCa                                                                                      <  2:128677/143‑1 (MQ=255)
                                                                  gtgAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCCGGCGTGGTCGGTCACTGAAATGc                                                                             >  1:518710/1‑143 (MQ=255)
                                                                   tgAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCt                                                                            <  2:579012/143‑1 (MQ=255)
                                                                      atatCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGtt                                                                         <  2:358535/143‑1 (MQ=255)
                                                                              tttGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTg                                                                 <  2:411412/143‑1 (MQ=255)
                                                                              cattcTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTg                                                                 <  1:146676/139‑1 (MQ=255)
                                                                                   tAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGACCTGGCGTGGGCGGTCACTGCATTGCTGTTAATCCGTGGTGTa                                                            >  2:197894/1‑143 (MQ=255)
                                                                                            gTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGTTTATCGTGGcac                                                   <  1:642795/143‑1 (MQ=255)
                                                                                                    tctgtgCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATCCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGATGTTCATCCGTGGTTTATCGTGGCACAGAAcccc                                           <  1:649644/141‑1 (MQ=255)
                                                                                                          ccGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGTTTATGGTGGCACAGAACACCcagcag                                     >  1:389713/1‑143 (MQ=255)
                                                                                                            gATCAGGGGATTAACGTCTGGGAACTGATTCGCCAGGCGAAt                                                                                                                                        >  2:281852/1‑42 (MQ=255)
                                                                                                            gATCAGGGGATTAACGTCTGGGAACTGATTCGCCAGGCGAAt                                                                                                                                        <  1:281852/42‑1 (MQ=255)
                                                                                                                aGGGGGTTAGCGTCTGGGAACTGATTCGCCTGGCGAATCGTCACGCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGTTTATCGTGGCACAGAACCCCCAGCAGGCGCgg                               <  2:642875/143‑1 (MQ=255)
                                                                                                                  gggATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATAATCTACAGCCTCGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGTTTATCGTGGCACAGAACCACCAGCAAGCGCGGCt                             >  1:95213/1‑143 (MQ=255)
                                                                                                                        aaCGTCTGGGAACTGAGGCGGCTGGCGAGTCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGGCACTGCATTGCTGTTGATCCGTGGTTTATCGTGGCACAGAACCCCCAGCAGGCGCGGCTTATCCg                       <  2:605632/143‑1 (MQ=255)
                                                                                                                           gTCTGGGAACTGATTCGCCTAGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGTTTATCGTGGCACAGAACCCCCAGCAGGCGCGGCTTATTCGTAc                    >  2:314360/1‑143 (MQ=255)
                                                                                                                           gTCCGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCGTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGTTTATCGTGGCACAGAACCCCCAGCAGGCGCGGCTTATCCGTAc                    <  1:264569/143‑1 (MQ=255)
                                                                                                                              tGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGTTTATCGTGGCACAGAACCCCCAGCAGGCGCGGCTTATCCGTACcgc                 >  1:70527/1‑143 (MQ=255)
                                                                                                                                 gAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGTTTATCGTGGCACAGAACCCCCAGCAGGCGCGGCTTATCCGTACcgcgcg              >  1:37618/1‑143 (MQ=255)
                                                                                                                                      gATTCGCCAGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGTTTATCGTGGCACAGAACCCCCAGCAGGCGCGGCTTATCCGGACCGCGCGCGAAg         >  2:622200/1‑143 (MQ=255)
                                                                                                                                      gATTCGCCAGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGTTTATCGTGGCACAGAACCCCCACCAGGCGCGGCTTATCCCTCCCCCGCGCGCAg         >  2:263616/1‑143 (MQ=255)
                                                                                                                                         gcGCCGGGCGAATCGTCAGCCGCGCGGTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGTTTATCGTGGCACAGAACCCCCAGCAGGCGCGGCTTATCCGTACCGCGCGCGAAGTGa      <  1:330405/142‑1 (MQ=255)
                                                                                                                                             cTGGCGAACCGTCACCCTCGCGTTAATATTCTCCAGCCCGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGTTTATCGTGGCACAGAACCCCCAGCAGGCGCGGCTTATCCGTACCGCGCGCGAAGTGAACGa  <  1:523708/143‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                             
GAGTGTGTCGTCACTAACTCGCGGACGGCGGAAATGTGTAAGCTCACCGAAAACAGCTTCCGCGATGTGAATATCGCTTTTGCTAATGAATTGTCGCTGATTTGTGCCGATCAGGGGATTAACGTCTGGGAACTGATTCGCCTGGCGAATCGTCACCCTCGCGTTAATATTCTTCAGCCTGGCCCTGGCGTGGGCGGTCACTGCATTGCTGTTGATCCGTGGTTTATCGTGGCACAGAACCCCCAGCAGGCGCGGCTTATCCGTACCGCGCGCGAAGTGAACGA  >  NC_000913/3971854‑3972137

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