Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,743,058 T→A 7.2% V109E (GTG→GAG)  yiaK → 2,3‑diketo‑L‑gulonate reductase, NADH‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,743,0580TA 7.2% 159.2 / 3.9 56V109E (GTG→GAG) yiaK2,3‑diketo‑L‑gulonate reductase, NADH‑dependent
Reads supporting (aligned to +/- strand):  ref base T (24/26);  new base A (2/2);  total (27/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.68e-01

GATATCATTCCTGATGCCCAACCCAAACGTATAACCAGCCTCGGCGCAATTGAACAGTGGGACGCCCAGCGTTCGATCGGTAACCTGACAGCGAAAAAGATGATGGATCGCGCCATTGAACTGGCTGCCGATCACGGTATTGGTCTGGTGGCACTACGTAATGCCAACCACTGGATGCGCGGCGGCAGCTACGGCTGGCAGGCGGCGGAAAAAGGCTATATTGGCATTTGCTGGACCAACTCCATCGCCGTAATGCCGCCGTGGGGCGCAAAAGAGTGTCGCATAGGCACTAACCC  >  NC_000913/3742910‑3743205
                                                                                                                                                    |                                                                                                                                                   
gATATCATTCCTGATGCCCAACCCAAACGTATAACCAGCCTCGGCGCATTTGAACAGTGGGACGCCCAGCGTTCGATCGGTAACCTGTCAGCGAAAAAGATGATGGATCGCGCCATTGAACTGGCTGCCGATCACGGTATTGGTCtggt                                                                                                                                                     <  5:202377/149‑1 (MQ=255)
gATATCATTCCTGATGACAAACCCAAACGTATAACAAGCCTCTGCGCAATTGAACAGTGGGACGCCCAGCTTTCGAACGGTAACCTGACAGCGAAAAAGATGATGGATCGCGCAATTGACATGGCTGCCTATCACTGTAATGGTCtggt                                                                                                                                                     >  4:114639/1‑149 (MQ=255)
      aTTCCTGATGCCCAACCCAAACGTATAACCAGCCTCGGCGCAATTGAACAGTGGGACGCCCAGCGTTCGATCGGTAACCTGACAGCGAAAAAGATGATGGATCGCGCCATTGAACTGGCTGCCGATCACGGTATTGGTCTGGTGGCAc                                                                                                                                                >  1:191622/1‑148 (MQ=255)
      aTTCCTGATGCCCAACCCAAACGTATAACCAGCCTCGGCGCAATTGAACAGTGGGACGCCCAGCGTTCGATCGGTAACCTGACAGCGAAAAAGATGATGGATCGCGCCATTGAACTGGCTGCCGATCACGGTATTGGTCTGGTGGCAc                                                                                                                                                <  5:181385/148‑1 (MQ=255)
         ccTGATGCCCAACCCAAACGTATAACCAGCCTCGGCGCAATTGAACAGTGGGACGCCCAGCGTTCGATCGGTAACCTGACAGCGAAAAAGATGATGGATCGCGCCATTGAACTGGCTGCCGATCACGGTATTGGTCTGGTGGCACTAc                                                                                                                                             <  1:150530/148‑1 (MQ=255)
         ccTGATGCCCAACCCAAACGTATAACCAGCCTCGGCGCAATTGAACAGTGGGACGCCCAGCGTTCGATCGGTAACCTGACAGCGAAAAAGATGATGGAACGCGCCATTGAACTGGCTGCCGATCACGGTATTGGTCTGGTGGCACTAc                                                                                                                                             <  5:141543/148‑1 (MQ=255)
                 ccaaaccaaaCGAATAACAAGCATATGCGCAATTCAACAGACCGACGCACAGCGATCGATCGGTAACCTGACAGCCAAAAACATGATGGATCGCGCCATAGACCTGGCTGCCGATCACGGTATTGGTCTGGTGGCACTACGTAATGcc                                                                                                                                     >  5:166305/1‑148 (MQ=255)
                  caacccaaACGTATAACCAGCCTCGGCGCAATTGAACAGTGGGACGCCCAGCGTTCGATCGGTAACCTGACAGCGAAAAAGATGATGGATCGCGCCATTGAACTGGCTGCCGATCACGGTATTGGTCTGGTGGCACTACGTAATGCCaa                                                                                                                                   >  4:225699/1‑149 (MQ=255)
                      ccaaACGTATAACCAGCCTCGGCGCAATTGAACAGTGGGACGCCCAGCGTTCGATCGGTAACCTGACAGCGAAAAAGATGATGGATCGCGCCATTGAACTGGCTGCCGATCACGGTATTGGTCTGGTGGCACTACGTAATGCCAACCAc                                                                                                                               >  3:199129/1‑149 (MQ=255)
                         aaCGTATAACCAGCCGCGGCGCAATTGAACAGTGGGACGCCCAGCGTTCGATCGGTAACCTGACAGCGAAAAAGATGATGGATCGCGCCATTGAACTGGCTGCCGATCACGGTATTGGTCTGGTGGCACTACGTAATGCCAACCACTgg                                                                                                                            <  8:211828/149‑1 (MQ=255)
                            gTATAACCCGCCTCGGCGCAATTGAACAGTGGGACGCCCAGCGTTCGATCGGTAACCTGACAGCGAAAAAGATGATGGATCGCGCCATTGAACTGGCTGCCGATCACGGTATTGGTCTGGTGGCACTACGTAATGCCAACCACTGGATg                                                                                                                         >  6:89934/1‑149 (MQ=255)
                            gTATAACCAGCCTCGGCGCAATTGAACAGTGGGACGCCCAGCGTTCGATCGGTAACCTGACAGCGAAAAAGATGATGGATCGCGCCATTGAACTGGCTGCCGATCACGGTATTGGTCTGGTGGCACTACGTAATGCCAACCACTGGATg                                                                                                                         >  6:149726/1‑149 (MQ=255)
                               taACCAGCCTCGGCGCAATTGAACAGTGGGACGCCCAGCGTTCGATCGGTAACCTGACAGCGAAAAAGATGATGGATCGCGCCATTGAACTGGCTGCCGATCACGGTATTGGTCTGGTGGCACTACGTAATGCCAACCACTGGATgcg                                                                                                                       <  6:159597/148‑1 (MQ=255)
                                aaCCAGCCTCGGCGCAATTGAACAGTGGGACGCCCAGCGTTCGATCGGTAACCTGACAGCGAAAAAGATGATGGATCGCGCCATTGAACTGGCTGCCGATCACGGTATTGGTCTGGTGGCACTACGTAATGCCAACCACTGGATGcgcg                                                                                                                     >  5:170555/1‑149 (MQ=255)
                                  ccAGCCTCGGCGCAATTGAACAGTGGGACGCCCAGCGTTCGATCGGTAACCTGACAGCGAAAAAGATGATGGATCGCGCCATTGAACTGGCTGCCGATCACGGTATTGGTCTGGTGGCACTACGTAATGCCAACCACTGGATGCgcggc                                                                                                                   <  4:43312/149‑1 (MQ=255)
                                    aGCCTCGGCGCAATTGAACAGTGGGACGCCCAGCGTTCGATCGGTAACCTGACAGCGAAAAAGATGATGGATCGCGCCATTGAACTGGCTGCCGATCACGGTATTGGTCTGGTGGCACTACGTAATGCCAACCACTGGATGCGcggcgg                                                                                                                 <  6:114932/149‑1 (MQ=255)
                                     gCCTCGGCGCAATTGAACAGTGGGACGCCCAGCGTTCGATCGGTAACCTGACAGCGAAAAAGATGATGGATCGCGCCATTGAACTGGCTGCCGATCACGGTATTGGTCTGGTGGCACTACGTAATGCCAACCACTGGATGCGcggcgg                                                                                                                 >  6:175686/1‑148 (MQ=255)
                                      ccTCGGCGCAATTGAACAGTGGGACGCCCAGCGTTCGATCGGTAACCTGACAGCGAAAAAGATGATGGATCGCGCCATTGAACTGGCTGCCGATCACGGTATTGGTCTGGTGGCACTACGTAATGCCAACCACTGGATGCGCGGCGGCa                                                                                                               <  2:99529/149‑1 (MQ=255)
                                               aaTTGAACAGTGGGACGCCCAGCGTTCGATCGGTAACCTGACAGCGAAAAAGATGATGGATCGCGCCATTGAACTGGCTGCCGATCACGGTATTGGTCTGGTGGCACTACGTAATGCCAACCACTGGATGCGCGGCGGCAGCTACGGCt                                                                                                      >  4:112543/1‑149 (MQ=255)
                                                            gACGCCCAGCGTTCGATCGGTAACCTGACAGCGAAAAAGATGATGGATCGCGCCATTGAACTGGCTGCCGATCACGGTATTGGTCTGGTGGCACTACGTAATGCCAACCACTGGATGCGCGGCGGCAGCTACGACTGGc                                                                                                   <  7:154569/139‑1 (MQ=255)
                                                            gACGCCCAGCGTTCGATCGGTAACCTGACAGCGAAAAAGATGATGGATCGCGCCATTGAACTGGCTGCCGATCACGGTATTGGTCTGGTGGCACTACGTAATGCCAACCACTGGATGCGCGGCGGCAGCTACGACTGGc                                                                                                   >  8:154569/1‑139 (MQ=255)
                                                                      gTTCGATCGGTAATTTGACAGCGTAATAGATGATGGATCGCGCTATTGAACTGGCTGCCGATCACGGTATTGGTCTGGTGGCACTTCGTAATGCCAACCTCTGGATGCGCGGCGGCAGTTTCGGCTGGCAGGCGGCGGaaata                                                                                     <  4:190228/143‑3 (MQ=255)
                                                                      gTTCGATCGGTAACCTGACAGCGAAAAAGATGATGGATCGCGCCATTGAACTGGCTGCCGATCACGGTATTGGTCTGGTGGCACTACGTAATGCCAACCACTGGATGCGCGGCGGCAGCTACGGCTGGCAGGCGGCGGaaaaa                                                                                     >  3:190228/1‑143 (MQ=255)
                                                                      gTTCGATCGGTAACCTGACAGCGAAAAAGATGATGGATCGCGCCATTGAACTGGCTGCCGATCACGGTATTGGTCTGGTGGCACTACGTAATGCCAACCACTGGATGCGCGGCGGCAGCTACGGCTGGCAGGCGGCGGAAAAAGGCta                                                                                >  6:87411/1‑148 (MQ=255)
                                                                          gATCGGTAACCTGACAGCGAAAAAGATGATGGATGGGGACATTGAACTGGCTGCCGATCACGGTATTTGTCTGGTGGCACTACGTAATGCCATCCACTGGATGCGCGGCGGCTGCTACGGGTGGCAGGCGGCGGAAAAATGCTATAtt                                                                            <  6:166305/148‑1 (MQ=255)
                                                                          gATCGGTAACCTGACAGCGAAAAAGATGATGGATCGCGCCATTGAACTGGCTGCCGATCACGGTATTGGTCTGGTGGCACTACGTAATGCCAACCACTGGATGCGCGGCGGCAGCTACGGCTGGCAGGCGGCGGAAAAAGGCTATATTg                                                                           >  1:102980/1‑149 (MQ=255)
                                                                          gATCGGTAACCTGACAGCGAAAAAGATGATGGATCGCGCCATTGAACTGGCTGCCGATCACGGTATTGGTCTGGTGGCACTACGTAATGCCAACCACTGGATGCGCGGCGGCAGCTACGGCTGGCAGGCGGCGGAAAAAGGCTATATTg                                                                           <  6:210353/149‑1 (MQ=255)
                                                                           agcGGTAACCTGACAGCGAAAACGATGATGGAGCGCGCCATTGCACGGGCTGCCGATCACGGTATTGGTCTGGGGGCACTACGTAAAGCCAACCACTGGATGCGCGGCGGCAGCTACGGCTGGCAGGCGGCGGAAAAAGGCTATATTgg                                                                          <  3:161393/147‑1 (MQ=255)
                                                                            tCGGTAACCTGACAGCGAAAAAGATGATGGATCGCGCCATTGAACTGGCTGCCGATCACGGTATTGGTCTGGAGGCACTACGTAATGCCAACCACTGGATGCGCGGCGGCAGCTACGGCTGGCAGGCGGCGGAAAAAGGCTATATTgg                                                                          >  7:34493/1‑148 (MQ=255)
                                                                            tCGGTAACCTGACAGCGAAAAAGATGATGGATCGCGCCATTGAACTGGCTGCCGATCACGGTATTGGTCTGGAGGCACTACGTAATGCCAACCACTGGATGCGCGGCGGCAGCTACGGCTGGCAGGCGGCGGAAAAAGGCTATATTgg                                                                          >  7:73320/1‑148 (MQ=255)
                                                                               gTAACCTGACAGCGAAAAAGATGATGGATCGCGCCATTGAACTGGCTGCCGATCACGGTATTGGTCTGGTGGCACTACGTAATGCCAACCACTGGATGCGCGGCGGCAGCTACGGCTGGCAggcgg                                                                                             >  8:219775/1‑126 (MQ=255)
                                                                                  aCCTGACAGCGAAAAAGATGATGGATCGCGCCATTGAACTGGCTGCCGATCACGGTATTGGTCTGGTGGCACTACGTAATGCCAACCACTGGATGCGCGGCGGCAGCTACGGCTGGCAGGCGGCGGAAAAAGGCTATATTGGCATTTGc                                                                   >  6:168125/1‑149 (MQ=255)
                                                                                    cTGACAGCGAAAAAGATGATGGATCGCGCCATTGAACTGGCTGCCGATCACGGTATTGGTCTGGTGGCACTACGTAATGCCAACCACTGGATGCGCGGCGGCAGCTACGGCTGGCAGGCGGCGGAAAAAGGCTATATTGGCATTTGCTg                                                                 <  5:168125/149‑1 (MQ=255)
                                                                                    cTGACAGCGAAAAAGATGATGGATAGAGCAATTGAACTGGCTGCCGATCACGGTATTGGTCTGGTGGCACTACGTAATGCCAACCACTGGATGCGCGGCGGCAGCTACGGCTGGCAGGCGGCGGAAAAAGGCTATATTTGCATTTGCTg                                                                 <  4:199129/149‑1 (MQ=255)
                                                                                     tGACAGCGAAAAAGATGATGGATCGCGCCATTGAACTGGCTGCCGATCACGGTATTGGTCTGGTGGCACTACGTAATGCCAACCACTGGATGCGCGGCGGCAGCTATGGCTGGCAGGCGGCGGAAAAAGGCTATATTGGCATTTGCTg                                                                 >  7:143948/1‑148 (MQ=255)
                                                                                       aCAGCGAAAAAGATGATGGATCGCGCCATTGAACTGGCTGCCGATCACGGTATTGGTCTGGTGGCACTACGTAATGCCAACCACTGGATGCGCGGCGGCAGCTACGGCTGGCAGGCGGCGGAAAAAGGCTATATTGGCATTTGCTGGa                                                               <  5:175686/148‑1 (MQ=255)
                                                                                       aCAGCGAAAAAGATGATGGATCGCGCCATTGAACTGGCTGCCGATCACGGTATTGGTCTGGTGGCACTACGTAATGCCAACCACTGGATGCGCGGCGGCAGCTACGGCTGGCAGGCGGCGGAAAAAGGCTATATTGGCATTTGCTGGa                                                               <  5:199703/148‑1 (MQ=255)
                                                                                       aCAGCGAAAAAGATGATGGATCGCGCCATTGAACTGGCTGCCGATCACGGTATTGGTCTGGAGGCACTACGTAATGCCAACCACTGGATGCGCGGCGGCAGCTACGGCTGGCAGGCGGCGGAAAAAGGCTATATTGGCATTTGCTgg                                                                <  8:73320/147‑1 (MQ=255)
                                                                                       aCAGCGAAAAAGATGATGGATCGCGCCATTGAACTGGCTGCCGATCACGGTATTGGTCTGGAGGCACTACGTAATGCCAACCACTGGATGCGCGGCGGCAGCTACGGCTGGCAGGCGGCGGAAAAAGGCTATATTGGCATTTGCTgg                                                                <  8:34493/147‑1 (MQ=255)
                                                                                            gAAAAAGATGATGGATCGCGCCATTGAACTGGCTGCCGATCACGGTATTGGTCTGGTGGCACTACGTAATGCCAACCACTGGATGCGCGGCGGCAGCTACGGCTGGCAGGCGGCGGAAAAAGGCTATATTGGCATTTGCTGGACCAACt                                                         >  1:184922/1‑149 (MQ=255)
                                                                                             aaaaaGATGATGGATCGCGCCATTGAACTGGCTGCCGATCACGGTATTGGTCTGGTGGCACTACGTAATGCCAACCACTGGATGCGCGGCGGCAGCTACGGCTGGCAGGCGGCGGAAAAAGGCTATATTGGCATTTGCTGGACCAACTc                                                        <  6:55830/149‑1 (MQ=255)
                                                                                              aaaaGATGATGGATCGCGCCATTGAACTGGCTGCCGATCACGGTATTGGTCTGGTGGCACTACGTAATGCCAACCACTGGATGCGCGGCGGCAGCTACGGCTGGCAGGCGGCGGAAAAAGGCTATATTGGCATTTGCTGGACCAACTcc                                                       <  6:170555/149‑1 (MQ=255)
                                                                                                      atgGATCGCGCCATTGAACTGGCTGCCGATCACGGTATTGGTCTGGTGGCACTACGTAATGCCAACCACTGGATGCGCGGCGGCAGCTACGGCTGGCAGGCGGCGGAAAAAGGCTATATTGGCATTTGCTGGACCAACTCCATCGCCg                                                >  1:108209/1‑148 (MQ=255)
                                                                                                        ggATCGCGCCATTGAACTGGCTGCCGATCACGGTATTGGTCTGGTGGCACTACGTAATGCCAACCACt                                                                                                                              >  6:102628/1‑68 (MQ=255)
                                                                                                        ggATCGCGCCATTGAACTGGCTGCCGATCACGGTATTGGTCTGGTGGCACTACGTAATGCCAACCACt                                                                                                                              <  5:102628/68‑1 (MQ=255)
                                                                                                          aTCGCGCCATTGAACTGGCTGCCGATCACGGTATTGGTCTGGCGGCACTACGTAATGCCAACCACTGGATGCGCGGCGGCAGCTACGGCTGGCAGGCGGCAGAAAAAGGCTATATTGGCATTTGCTGGACCAACTCCATCGCCGTAATg                                           >  2:202979/1‑149 (MQ=255)
                                                                                                                     gAACTGGCTGCCGATCACGGTATTTGTCTGGTGGCACTACGTAATGCCAACCACTGGATGCGCGGCGGCAGCTACGGCTGGCAGGCGGCGGAAAAAGGCTATATTGGCATTTGCTGGACCAACTCCATCGCCGTAATGCCGCCGTggg                                 <  3:114639/148‑1 (MQ=255)
                                                                                                                      aaCTGGCTGCCGATCACGGTATTGGTCTGGTGGCACTACGTAATGCCAACCACTGGATGCGCGGCGGCAGCTACGGCTGGCAGGCGGCGGAAAAAGGCTATATTGGCATTTGCTGGACCAACTCCATCGCCGGATTGCCGCCGTGGGgc                               >  6:168537/1‑149 (MQ=255)
                                                                                                                          ggCTGCCGATCACGGTATTGGTCTGGTGGCACTACGTAATGCCAACCACTGGATGCGCGGCGGCAGCTACGGCTGGCAGGCGGCGGAAAAAGGCTATATTGGc                                                                         <  1:163555/103‑1 (MQ=255)
                                                                                                                          ggCTGCCGATCACGGTATTGGTCTGGTGGCACTACGTAATGCCAACCACTGGATGCGCGGCGGCAGCTACGGCTGGCAGGCGGCGGAAAAAGGCTATATTGGc                                                                         >  2:163555/1‑103 (MQ=255)
                                                                                                                             tGCCGATCACGGTATTGGTCTGGTGGCACTACGTAATGCCAACCACTGGATGCGCGGCGGCAGCTACGGCTGGCAGGCGGCGGAAAAAGGCTATATTGGCATTTGCTGGACCAACTCCATCGCCGTAATGCCGCCGTGGGGCGCaa                           >  3:64354/1‑146 (MQ=255)
                                                                                                                                   tCACGGTATTGGTCTGGTGGCACTACGTAATGCCAACCACTGGATGCGCGGCGGCAGCTACGGCTGGCAGGCGGCGGAAAAAGGCTATATTGGCATTTGCTGGACCAACTCCATCGCCGTAATGCCGCCGTGGGGCGCAAAAGAGTGTc                  <  5:149726/149‑1 (MQ=255)
                                                                                                                                        gTATTGGTCTGGTGGCACTACGTAATGCCAACCACTGGATGCGCGGCGGCAGCTACGGCTGGCAGGCGGCGGAAAAAGGCTATATTGGCATTTGc                                                                   <  1:169406/95‑1 (MQ=255)
                                                                                                                                        gTATTGGTCTGGTGGCACTACGTAATGCCAACCACTGGATGCGCGGCGGCAGCTACGGCTGGCAGGCGGCGGAAAAAGGCTATATTGGCATTTGc                                                                   >  2:169406/1‑95 (MQ=255)
                                                                                                                                        gTATTGGTCTGGTGGCACTACGTAATGCCAACCACTGGATGCGCGGCGGCAGCTACGGCTGGCAGGCGGCGGAAAAAGGCTATATTGGCATTTGCTGGACCAACTCCATCGCCGTAATGCCGCCGTGGGGCGCAAAAGAGTGTCGCATa             <  5:89934/149‑1 (MQ=255)
                                                                                                                                         tATTGGTCTGGTGGCACTACGTAATGCCAACCACTGGATGCGCGGCGGCAGCTACGGCTGGCAGGCGGCGGAAAACGGCTATATTGGCATTTGCTGGACCAACTCCATCGCCGTAATgc                                          <  8:149017/119‑1 (MQ=255)
                                                                                                                                         tATTGGTCTGGTGGCACTACGTAATGCCAACCACTGGATGCGCGGCGGCAGCTACGGCTGGCAGGCGGCGGAAAACGGCTATATTGGCATTTGCTGGACCAACTCCATCGCCGTAATgc                                          >  7:149017/1‑119 (MQ=255)
                                                                                                                                           ttGGGCTGGTGGCACTACGTAATGCCAACCACTGGATGCGCGGCGGCAGCTACGGCTGGCAGGCGGCGGAAAAAGGCTATATTGGCATTTGCTGGACCAACTCCATCGCCGTTATGCCGCCGTGGGGCGCAAAAGAGTGTCGCATAGGc          <  3:62478/149‑1 (MQ=255)
                                                                                                                                                   gtggCACTACGTAATGCCAACCACTGGATGCGCGGCGGCAGCTACGGCTGGCAGGCGGCGGAAAAAGGCTATATTGGCATTTGCTGGACCAACTCCATCGCCGTAATGCCGCCGTGGGGCGCAAAAGAGTGTCGCATAGGCACTAAccc  >  5:4709/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GATATCATTCCTGATGCCCAACCCAAACGTATAACCAGCCTCGGCGCAATTGAACAGTGGGACGCCCAGCGTTCGATCGGTAACCTGACAGCGAAAAAGATGATGGATCGCGCCATTGAACTGGCTGCCGATCACGGTATTGGTCTGGTGGCACTACGTAATGCCAACCACTGGATGCGCGGCGGCAGCTACGGCTGGCAGGCGGCGGAAAAAGGCTATATTGGCATTTGCTGGACCAACTCCATCGCCGTAATGCCGCCGTGGGGCGCAAAAGAGTGTCGCATAGGCACTAACCC  >  NC_000913/3742910‑3743205

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