Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,012,599 T→C 7.5% intergenic (+206/+15) yqeA → / ← yqeB putative amino acid kinase/XdhC‑CoxI family protein with NAD(P)‑binding Rossman fold

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,012,5990TC 7.5% 172.2 / 3.6 53intergenic (+206/+15)yqeA/yqeBputative amino acid kinase/XdhC‑CoxI family protein with NAD(P)‑binding Rossman fold
Reads supporting (aligned to +/- strand):  ref base T (25/24);  new base C (2/2);  total (27/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

ATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATGCATCAACATCATTAACGCCTCAAGTACGCCGCCGCCAATCGCCCGGGCTTTATCAGAAACGCTGGTGAAA  >  NC_000913/3012508‑3012728
                                                                                           |                                                                                                                                 
aTAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATc                                                                                                                                  <  4:2214/93‑1 (MQ=255)
aTAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATc                                                                                                                                  >  3:2214/1‑93 (MQ=255)
      tGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCActt                                                                                                 >  1:198625/1‑120 (MQ=255)
      tGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCActt                                                                                                 <  2:198625/120‑1 (MQ=255)
          aTGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCActtct                                                                                               >  4:266877/1‑118 (MQ=255)
          aTGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCActtct                                                                                               <  3:266877/118‑1 (MQ=255)
           tGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATGCATCAAca                                                                >  3:986194/1‑148 (MQ=255)
              aCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATGCATCAAcatca                                                             >  2:871039/1‑148 (MQ=255)
                 tGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATGCATCAACATCATTa                                                          <  1:871039/148‑1 (MQ=255)
                       ccTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCAt                                                                                                                     >  4:282293/1‑83 (MQ=255)
                       ccTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCAt                                                                                                                     <  3:282293/83‑1 (MQ=255)
                       ccTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCAt                                                                                                                     >  4:171527/1‑83 (MQ=255)
                       ccTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCAt                                                                                                                     <  3:171527/83‑1 (MQ=255)
                         tGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCAt                                                                                                                                   >  4:1039998/1‑67 (MQ=255)
                         tGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCAt                                                                                                                                   <  3:1039998/67‑1 (MQ=255)
                           aTAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCg                                                                                        >  2:67985/1‑108 (MQ=255)
                           aTAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCg                                                                                        <  1:67985/108‑1 (MQ=255)
                           aTAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACg                                                                                >  3:1078550/1‑116 (MQ=255)
                           aTAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACg                                                                                <  4:1078550/116‑1 (MQ=255)
                             aaGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTg                                                                                           >  2:1127931/1‑103 (MQ=255)
                             aaGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTg                                                                                           <  1:1127931/103‑1 (MQ=255)
                                     cgcATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTa                                                                                                                  >  3:1114656/1‑72 (MQ=255)
                                     cgcATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTa                                                                                                                  <  4:1114656/72‑1 (MQ=255)
                                        aTCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATGCAt                                                                     >  3:1122983/1‑114 (MQ=255)
                                        aTCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATGCAt                                                                     <  4:1122983/114‑1 (MQ=255)
                                         tCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATGCATCAACATCATTAACGCCTCAAGTAcgccg                                         <  2:633797/141‑1 (MQ=255)
                                         tCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATGCATCAACATCATTAACGCCTCAAGTAcgccg                                         >  1:633797/1‑141 (MQ=255)
                                             gCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGccc                                                                             >  2:873087/1‑101 (MQ=255)
                                             gCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGccc                                                                             <  1:873087/101‑1 (MQ=255)
                                              cAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTg                                                                           <  4:418249/102‑1 (MQ=255)
                                              cAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTg                                                                           >  3:418249/1‑102 (MQ=255)
                                              cAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATGCATCAACATCATTAACGCCTCAAGTACGCCGCCGCCAATCGcc                             >  3:1121917/1‑148 (MQ=255)
                                              cAATTTGGCGTTTGCCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATGCATCAACATCATTAACGCCTCAAGTACGCCGCCGCCAATCGcc                             >  2:227486/1‑148 (MQ=255)
                                                         ttGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATGCATCAACATCATTAACGCCTCAAGTACGCCGCCGCCAATCGCCCGGGCTTTATc                  <  4:1121917/148‑1 (MQ=255)
                                                               tcaGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATGCATCAACATCATTAACGCCTCAAGTAcg                                            <  1:123281/116‑1 (MQ=255)
                                                               tcaGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATGCATCAACATCATTAACGCCTCAAGTAcg                                            >  2:123281/1‑116 (MQ=255)
                                                                  gTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATGCATCAACATCATTAACGCCTCAAGTACGCCGCCGCCAATCGCCCGGGCTTTATCAGAAACGCt         <  1:227486/148‑1 (MQ=255)
                                                                      cAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATGCATCAACATCATTAACGCCTCAAGTACGCCGCCGCCAATCGCCCGGGCTTTATCAGAAACGCt         <  4:898652/144‑1 (MQ=255)
                                                                      cAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATGCATCAACATCATTAACGCCTCAAGTACGCCGCCGCCAATCGCCCGGGCTTTATCAGAAACGCt         >  3:898652/1‑144 (MQ=255)
                                                                      cAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATGCATCAACATCATTAACGCCTCAAGTACGCCGCCGCCAATCGCCCGGGCTTTATCAGAAACGCt         >  3:356022/1‑144 (MQ=255)
                                                                      cAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATGCATCAACATCATTAACGCCTCAAGTACGCCGCCGCCAATCGCCCGGGCTTTATCAGAAACGCt         <  4:356022/144‑1 (MQ=255)
                                                                        ggCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATGCATCAACATCATTAACGCCTCAAGTACGCCGCCGCCAATCGCCCGGGCTTTATCAGAAACGCTGGTGaa   >  1:502820/1‑148 (MQ=255)
                                                                         gCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATGCATCAACATCATTAACGCCTCAAGTACGCCGCCGCCAATCGCCCGGGCTTTATCAGAAACGCTGGTGaaa  <  2:502820/148‑1 (MQ=255)
                                                                            gggtgggtAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATGCATCAACATCATTAACGCCTCAAGTACGCCGCCGCCAATCGCCCGGGCTTTATCAGAAACg           <  3:1053941/136‑1 (MQ=255)
                                                                            gggtgggtAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATGCATCAACATCATTAACGCCTCAAGTACGCCGCCGCCAATCGCCCGGGCTTTATCAGAAACg           >  4:1053941/1‑136 (MQ=255)
                                                                                    ttacccACCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATGCATCAACATCATTAACGCCTCAAGTACGCCGCCGCCAATCGCCCggg                         <  2:290594/110‑1 (MQ=255)
                                                                                    ttacccACCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATGCATCAACATCATTAACGCCTCAAGTACGCCGCCGCCAATCGCCCggg                         >  1:290594/5‑114 (MQ=255)
                                                                                    ttacccACCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATGCATCAACATCATTAACGCCTCAAGTACGCCGCCGCCAATCGCCCggg                         >  1:310125/5‑114 (MQ=255)
                                                                                    ttacccACCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATGCATCAACATCATTAACGCCTCAAGTACGCCGCCGCCAATCGCCCggg                         <  2:310125/110‑1 (MQ=255)
                                                                                    aaTACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATGCATCAACATCATTAACGCCTCAAGTa                                              >  4:682893/1‑93 (MQ=255)
                                                                                    aaTACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATGCATCAACATCATTAACGCCTCAAGTa                                              <  3:682893/93‑1 (MQ=255)
                                                                                      tACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATGCATCAACATCATTAACGCCTCAAGTACGCCGCCGCCAATCGCCCGGGCTTTATCAGAAACGCTGGTGaa   <  3:603215/134‑1 (MQ=255)
                                                                                      tACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATGCATCAACATCATTAACGCCTCAAGTACGCCGCCGCCAATCGCCCGGGCTTTATCAGAAACGCTGGTGaa   >  4:603215/1‑134 (MQ=255)
                                                                                        ccATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTg                                                                           >  3:1032041/1‑60 (MQ=255)
                                                                                        ccATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTg                                                                           <  4:1032041/60‑1 (MQ=255)
                                                                                           |                                                                                                                                 
ATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATGCATCAACATCATTAACGCCTCAAGTACGCCGCCGCCAATCGCCCGGGCTTTATCAGAAACGCTGGTGAAA  >  NC_000913/3012508‑3012728

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: