Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,841,978 T→A 12.4% intergenic (‑114/+138) yqhD ← / ← PP_2493 NAD(P)‑dependent alcohol dehydrogenase/RadC family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,841,9780TA12.4% 92.0 / 5.6 34intergenic (‑114/+138)yqhD/PP_2493NAD(P)‑dependent alcohol dehydrogenase/RadC family protein
Reads supporting (aligned to +/- strand):  ref base T (15/14);  new base A (2/2);  total (18/16)
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.39e-01

TTCGCGTGGGGTTGTGGAAGTCAAAATTGAGCATAGATTTCTCTACTCCTTAAAAGTAGACCAGTCGTATATTTTTTGGTGAAATTTTTATTGCAGAAGCTCGTCGGCGCGCATTGCGCACCAAATAGCTGAAGATTGGTTAGAAGCTTGCTTCTAAAGGCAACCAAGTTGGTGTGGCATTCGTCTGCATTCCTGGTCCCGAGCCAAAGGCGTTGCCGCCAAGCAGCGTAAGCTCCGCGCGATGTCGAGGGTTGTTGCTAGGCACTGATGATAGGGCATCGGGAGTCATAGGTTTA  >  NC_002947/2841831‑2842126
                                                                                                                                                   |                                                                                                                                                    
ttCGCGTGGGGTTGTGGAAGTCAAAATTGAGCATAGATTTCTCTACTCCTTAAAAGTAGACCAGTCGTATATTTTTTGGTGAAATTTTTATTGCAGAAGCTCGTCGGCGCGCATTGCGCACCAAATAGCTGAAGATTGGTTAGAAgcttg                                                                                                                                                    <  1:1131742/150‑1 (MQ=255)
      tGGGGTTGTGGAAGTCAAAATTGAGCATAGATTTCTCTACTCCTTAAAAGTAGACCAGTCGTATATTTTTTGGTGAAATTTTTATTGCAGAAGCTCGTCGGCGCGCATTGCGCACCAAATAGCTGAAGATTGGTTAGAAGCTTGCTTCTa                                                                                                                                              >  1:609079/1‑150 (MQ=255)
       ggggTTGTGGAAGTCAAAATTGAGCATAGATTTCTCTACTCCTTAAAAGTAGACCAGTCGTATATTTTTTGGTGAAATTTTTATTGCAGAAGCTCGTCGGCGCGCATTGCGCACCAAATAGCTGAAGATTGGTTAGAAGCTTGCTTCTaa                                                                                                                                             >  1:583896/1‑150 (MQ=255)
         ggTTGTGGAAGTCAAAATTGAGCATAGATTTCTCTACTCCTTAAAAGTAGACCAGTCGTATATTTTTTGGTGAAATTTTTATTGCAGAAGCTCGTCGGCGCGCATTGCGCACCAAATAGCTGAAGATTGGTTAGAAGCTTGCTTCTAAAg                                                                                                                                           <  1:878306/150‑1 (MQ=255)
           ttGTGGAAGTCAAAATTGAGCATAGATTTCTCTACTCCTTAAAAGTAGACCAGTCGTATATTTTTTGGTGAAATTTTTATTGCAGAAGCTCGTCGGCGCGCATTGCGCACCAAATAGCTGAAGATTGGTTAGAAGCAAGCTTCTaaccaa                                                                                                                                         <  2:384318/150‑5 (MQ=255)
              tgGAAGTCAAAATTGAGCATAGATTTCTCTACTCCTTAAAAGTAGACCAGTCGTATATTTTTTGGTGAAATTTTTATTGCAGAAGCTCGTCGGCGCGCATTGCGCACCAAATAGCTGAAGATTGGTTAGAAGCTTGCTTCTAAAGGCAAc                                                                                                                                      >  2:602683/1‑150 (MQ=255)
                                            aCTCCTTAAAAGTAGACCAGTCGTATATTTTTTGGTGAAATTTTTATTGCAGAAGCTCGTCGGCGCGCATTGCGCACCAAATAGCTGAAGATTGGTTAGAAGCTTGCTTCTAAAGGCAACCAAGTTGGTGTGGCATTCGTCTGCATTCCt                                                                                                        >  2:1019567/1‑150 (MQ=255)
                                                        tAGACCAGTCGTATATTTTTTGGTGAAATTTTTATTGCAGAAGCTCGTCGGCGCGCATTGCGCACCAAATAGCTGAAGATTGGTTAGAAGCTTGCTTCTaaa                                                                                                                                            <  1:1120199/102‑1 (MQ=255)
                                                        tAGACCAGTCGTATATTTTTTGGTGAAATTTTTATTGCAGAAGCTCGTCGGCGCGCATTGCGCACCAAATAGCTGAAGATTGGTTAGAAGCTTGCTTCTaaa                                                                                                                                            >  2:1120199/1‑102 (MQ=255)
                                                           aCCAGTCGTATATTTTTTGGTGAAATTTTTATTGCAGAAGCTCGTCGGCGCGCATTGCGCACCAAATAGCTGAAGATTGGTTAGAAgctt                                                                                                                                                     <  2:986214/90‑1 (MQ=255)
                                                           aCCAGTCGTATATTTTTTGGTGAAATTTTTATTGCAGAAGCTCGTCGGCGCGCATTGCGCACCAAATAGCTGAAGATTGGTTAGAAgctt                                                                                                                                                     >  1:986214/1‑90 (MQ=255)
                                                           aCCAGTCGTATATTTTTTGGTGAAATTTTTATTGCAGAAGCTCGTCGACGCGCATTGCGCACCAAATAGCTGAAGATTGGTTAGAAGCTTGCTTCTAAAGGCAACCAAGTTGGTGTGGCATTCGTCTGCATTCCTGGTCCCGATCCAAAg                                                                                         >  2:175681/1‑150 (MQ=255)
                                                                tCGTATATTTTTTGGTGAAATTTTTATTGCAGAAGCTCGTCGGCGCGCATTGCGCACCAAATAGCTGAAGATTGGTTAGAAGCTTGCTTCTAAAGGCAACCAAGTTGGTGTGGCATTCGTCTGCATTCCTGGTCCCGAGCCAAAGGCGtt                                                                                    <  2:714969/150‑1 (MQ=255)
                                                                             ggTGAAATTTTTATTGCAGAAGCTCGTCGGCGCGCATTGCGCACCAAATAGCTGAAGATTGGTTAGAAGCTTGCTTCTAAAGGCAACCAAGTTGGTGTGGCATTCGTCTGCATTCCTGGTCCCGAGCCAAAGGCGTTGCCGCCAagcagc                                                                       <  2:702887/150‑1 (MQ=255)
                                                                             ggTGAAATTTTTATTGCAGAAGCTCGTCGGCGCGCATTGCGCACCAAATAGCTGAAGATTGGTTAGAAGCTTGCTTCTAAAGGCAACCAAGTTGGTGTGGCATTCGTCTGCATTCCTGGTCCCGAGCCAAAGGCGTTGCCGCCAagcagc                                                                       <  2:583896/150‑1 (MQ=255)
                                                                                gAAATTTTTATTGCAGAAGCTCGTCGGCGCGCATTGCGCACCAAATAGCTGAAGATTGGTTAGAAGCTTGCTTCTAAAGGCAACCAAGTTGGTGTGGCATTCGTCTGCATTCCTGGTCCCGAGCCAAAGGCGTTGCCGCCAAGCAGCGTa                                                                    >  1:845438/1‑150 (MQ=255)
                                                                                    tttttATTGCAGAAGCTCGTCGGCGCGCATTGCGCACCAAATAGCTGAAGATTGGTTAGAAGCTTGCTTCTAAAGGCAACCAAGTTGGTGTGGCATTCGTCTGCATTCCTGGTCCCGAGCCAAAGGCGTTGCCGCCAAGCAGCGTAAGCt                                                                >  1:1008852/1‑150 (MQ=255)
                                                                                       tgaTTGCTGAAGCTCGTCGGCGCGCATTGGGCACCAAATAGCTGAAGATTGGTTAGAAGCTTGCTTCTAAAGGGAACCAAGTTGGTGTGGCATTCGTCTGCACTGCTTCTCGCTAGCCAAAAGCGTTGCCGCCAAGCAGCGTAAGCTccc                                                             <  2:18408/148‑2 (MQ=255)
                                                                                             cAGAAGCTCGTCGGCGCGCATTGCGCACCAAATAGCTGAAGATTGGTTAGAAGCTTGCTTCTAAAGGCAACCAAGTTGGTGTGGCATTCGTCTGCATTCCTGGTCCCGAGCCAAAGGCGTTGCCGCCAAGCAGCGTAAGCTCCGCGCGAt                                                       >  2:825960/1‑150 (MQ=255)
                                                                                                aaGCTCGTCGGCGCGCATTGCGCACCAAATAGCTGAAGATTGGTTAGAAGCTTGCTTCTAAAGGCAACCAAGTTGGTGTGGCATTCGTCTGCATTCCTGGTCCCGAGCCAAAGGCGTTGCCGCCAAGCAGCGTAAGCTCCGCGCGATGTc                                                    <  1:825960/150‑1 (MQ=255)
                                                                                                       tcgGCGCGCATTGCGCACCAAATAGCTGAAGATTGGTTAGAAGCTTGCTTCTAAAGGCAACCAAGTTGGTGTGGCATTCGTCTGCATTCCTGGTCCCGAGCCAAAGGCGTTGCCGCCAagcagc                                                                       <  2:214040/124‑1 (MQ=255)
                                                                                                       tcgGCGCGCATTGCGCACCAAATAGCTGAAGATTGGTTAGAAGCTTGCTTCTAAAGGCAACCAAGTTGGTGTGGCATTCGTCTGCATTCCTGGTCCCGAGCCAAAGGCGTTGCCGCCAagcagc                                                                       >  1:214040/1‑124 (MQ=255)
                                                                                                          gcgcgcATTGCGCACCAAATAGCTGAAGATTGGTTAGAAGCTTGCTTCTAAAGGCAACCAAGTTGGTGTGGCATTCGTCTGCATTCCTGGTCCCGAGCCAAAGGCGTTGCCGCCAAGCAGCGTAAGCTCcgc                                                            <  2:331529/132‑1 (MQ=255)
                                                                                                          gcgcgcATTGCGCACCAAATAGCTGAAGATTGGTTAGAAGCTTGCTTCTAAAGGCAACCAAGTTGGTGTGGCATTCGTCTGCATTCCTGGTCCCGAGCCAAAGGCGTTGCCGCCAAGCAGCGTAAGCTCcgc                                                            >  1:331529/1‑132 (MQ=255)
                                                                                                             cgcATTGCGCACCAAATAGCTGAAGATTGGTTAGAAGCTTGCTTCTAAAGGCAACCAAGTTGGTGTGGCATTCGTCTGCATTCCTGGTCCCGAGCCAAAGGCGTTGCCGCCAAGCAGCGTAAGCTCCGCGCGATGTCGAGGGTTGTTGCt                                       <  1:602683/150‑1 (MQ=255)
                                                                                                               cATTGCGCACCAAATAGCTGAAGATTGGTTAGAAGCTTGCTTCTAAAGGCAACCAAGTTGGTGTGGCATTCGTCTGCATTCCTGGTCCCGAGCCAAAGGCGTTGCCGCCAAGCAGCGTAAGCTCCGCGCGATGTCGAGGGTTGTTGCTAg                                     >  2:183802/1‑150 (MQ=255)
                                                                                                                aTTGCGCACCAAATAGCTGAAGATTGGTTAGAAGCTTGCTTCTACAGGCTGGTGAGTTGGTGTGGCATTCGTCTGCATTCCTGGTCCCGAGCCAAAGGCGTGGCCGCCAAGCAGCGTAAGCTCCGCGCGATGTCGAGGGTTGTTGCTAgg                                    >  1:936759/1‑150 (MQ=255)
                                                                                                                aTTGCGCACCAAATAGCTGAAGATTGGTTAGAAGCTTGCTTCTAAAGGCAACCAAGTTGGTGTGGCATTCGTCTGCATTCCTGGTCCCGAGCCAAAGGCGTTGCCGCCAAGCAGCGTAAGCTCCGCGCGATGTCGAGGgttg                                            >  2:468961/1‑142 (MQ=255)
                                                                                                                aTTGCGCACCAAATAGCTGAAGATTGGTTAGAAGCTTGCTTCTAAAGGCAACCAAGTTGGTGTGGCATTCGTCTGCATTCCTGGTCCCGAGCCAAAGGCGTTGCCGCCAAGCAGCGTAAGCTCCGCGCGATGTCGAGGgttg                                            <  1:468961/142‑1 (MQ=255)
                                                                                                                     gcACCAAATAGCTGAAGATTGGTTAGATCGCTGCTTCTAAAGGCAACCAAGTTGGTGTGGCATTCGTCTGCATTCCTGGTCCCGAGCCAAAGGCGTTGCCGCCAAGCAGCGTAAGCTCCGCGCGATGTCGAGGGTTGTTGCTAGGCACtg                               >  1:937938/1‑150 (MQ=255)
                                                                                                                      cACCAAATAGCTGAAGATTGGTTAGAAGCTTGCTTCTAAAGGCAACCAAGTTGGTGTGGCATTCGTCTGCATTCCTGGTCCCGAGCCAAAGGCGTTGCCGCCAAGCAGCGTAAGCTCCGCGCGATGTCGAGGGTTGTTGCTAGGCACtga                              <  2:609079/150‑1 (MQ=255)
                                                                                                                         cAAATAGCTGAAGATTGGTTAGAAGCTTGCTTCTAAAGGCAACCAAGTTGGTGTGGCATTCGTCTGCATTCCTGGTCCCGAGCCAAAGGCGTTGCCGCCAAGCAGCGTAAGCTCCGCGCGATGTCGAGGGTTGTTGCTAGGCACtgatga                           <  2:1114243/150‑1 (MQ=255)
                                                                                                                            aTAGCTGAAGATTGGTTAGAAGCTTGCTTCTAAAGGCAACCAAGTTGGTGTGGCATTCGTCTGCATTCCTGGTCCCGAGCCAAAGGCGTTGCCGCCAAGCAGCGTAAGCTCCGCGCGATGTCGAGGGTTGTTGCTAGGCACTGATGATAg                        >  2:87692/1‑150 (MQ=255)
                                                                                                                                     ttgcctttAGAAGCAAGCTTCTAAAGGCAACCAAGTTGGTGTGGCATTCGTCTGCATTCCTGGTCCCGAGCCAAAGGCGTTGCCGCCAAGCAGCGTAAGCTCCGCGCGATGTCGAGGGTTGTTGCTAGGCACTGATGATAGGGCATCggg               >  2:606683/7‑150 (MQ=255)
                                                                                                                                      aTTGGTTAGAAGCTTGCTTCTAAAGGCAACCAAGTTGGTGTGGCATTCGTCTGCATTCCTGGTCCCGAGCCAAAGGCGTTGCCGCCAAGCAGCGTAAGCTCCGCGCGATGTCGAGGGTTGTTGCTAGGCACTGATGATAGGGCATCGGGa              <  2:845438/150‑1 (MQ=255)
                                                                                                                                         ctttAGAAGCAAGCTTCTAAAGGCAACCAAGTTGGTGTGGCATTCGTCTGCATTCCTGGTCCCGAGCCAAAGGCGTTGCCGCCAAGCAGCGTAAGCTCCGCGCGATGTCGAg                                                 <  1:302968/110‑1 (MQ=255)
                                                                                                                                         ctttAGAAGCAAGCTTCTAAAGGCAACCAAGTTGGTGTGGCATTCGTCTGCATTCCTGGTCCCGAGCCAAAGGCGTTGCCGCCAAGCAGCGTAAGCTCCGCGCGATGTCGAg                                                 >  2:302968/3‑112 (MQ=255)
                                                                                                                                                  cttgcttCTAAAGGCAACCAAGTTGGTGTGGCATTCGTCTGCATTCCTGGTCCCGAGCCAAAGGCGTTGCCGCCAAGCAGCGTAAGCTCCGCGCGATGTCGAGGGTTGTTGCTAGGCACTGATGATAGGGCATCGGGAGTCATAGGTTTa  >  1:864434/1‑150 (MQ=255)
                                                                                                                                                   ttgcttCTAAAGGCAACCAAGTTGGTGTGGCATTCGTCTGCATTCCTGGTCCCGa                                                                                                >  2:821228/1‑55 (MQ=255)
                                                                                                                                                   ttgcttCTAAAGGCAACCAAGTTGGTGTGGCATTCGTCTGCATTCCTGGTCCCGa                                                                                                <  1:821228/55‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
TTCGCGTGGGGTTGTGGAAGTCAAAATTGAGCATAGATTTCTCTACTCCTTAAAAGTAGACCAGTCGTATATTTTTTGGTGAAATTTTTATTGCAGAAGCTCGTCGGCGCGCATTGCGCACCAAATAGCTGAAGATTGGTTAGAAGCTTGCTTCTAAAGGCAACCAAGTTGGTGTGGCATTCGTCTGCATTCCTGGTCCCGAGCCAAAGGCGTTGCCGCCAAGCAGCGTAAGCTCCGCGCGATGTCGAGGGTTGTTGCTAGGCACTGATGATAGGGCATCGGGAGTCATAGGTTTA  >  NC_002947/2841831‑2842126

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