Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,597,083 T→G 9.5% intergenic (‑78/+18) nemA ← / ← PP_3174 FMN‑linked N‑ethylmaleimide reductase/oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,597,0830TG 9.5% 116.4 / 3.0 40intergenic (‑78/+18)nemA/PP_3174FMN‑linked N‑ethylmaleimide reductase/oxidoreductase
Reads supporting (aligned to +/- strand):  ref base T (20/16);  new base G (2/2);  total (22/18)
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.51e-01

CCATGAATACGCGGTTGGGCAGGGTGAGTGGGCCGATTTGCAGCGGTTGCAAGAGTTTCATTGCAGCTCCAGAGGTGGAAAACCAATGGTGCGCACGATACCGGGGTGATTGGGGAAAGGGGGGACCGGGAGCCCGGCTCAGGGCTCCCCAATGCGTCAGCTGCGTGCCAGCCGCCCACTCACATAATGCGCGACGTCATTGAAGCCCGACGTGGACGCATGCCCTGGCGTTACCAACGAATCGATGAATGCCTCATCCTCCGCGGTAATCTTCGCCGCCAATG  >  NC_002947/3596945‑3597228
                                                                                                                                          |                                                                                                                                                 
ccATGAATCCGCGGTGGGGCTGGGGGAGTGGGCCGATTTGCTGCGGTTGCAAGTGTTTCATTGCAGCTCCAGAGGTGTAAAACCAATGGTGCGCACGATACCGGGGTGATTGGGGAAAGGGGGGACCGGGAGCCCTGAGCCGGGCTccc                                                                                                                                         <  2:176489/149‑1 (MQ=255)
       tACGCGGTTGGGCAGGGTGAGTGGGCCGATTTGCAGCGGTTGCAAGAGTTTCATTGCAGCTCCAGAGGTGGAAAACCAATGGTGCGCACGATACCGGGGTGATTGGGGAAAGGGGGGACCGGGAGCCCGGCTCAGGGCTCCCCAATGCg                                                                                                                                  >  5:305356/1‑149 (MQ=255)
          gcgGTTGGGCAGGGTGAGTGGGCCGATTTGCAGCGGTTGCAAGAGTTTCATTGCAGCTCCAGAGGTGGAAAACCAATGGTGCGCACGATACCGGGGTGATTGGGGAAAGGGGGGACCGGGAGCCCGGCTCAGGGCTCCCCAATGCGTcc                                                                                                                               >  7:318742/1‑148 (MQ=255)
                  accgggTTAGTGGGCGTATTTGCTGCGGTTGCAAGAGTTTCATTGCAGCTCCAGAGTTGGAAAACGAATGGTGCGCACGATACCGGGGTGATTGGGGAAAGGGGGGACCGGGAGCCCGGCTCAGGGCTCCCCAATGCGTCAGCTGCGt                                                                                                                        <  4:270287/145‑1 (MQ=255)
                          aGTGGGCCGATTTGCAGCGGTTGCAAGAGTTTCATTGCAGCTCCAGAGGTGGAAAACCAATGGTGCGCACGATACCGGGGTGATTGGGGAAAGGGGGGACCGGGAGCCCGGCTCAGGGCTCCCCAATGCGTCATCTGCGTGCCAgccgc                                                                                                               >  5:297642/1‑149 (MQ=255)
                            ttggCCGTTTTTCAGCGGTTGCAAGATTTTCATTGCAGCTCCAGAGGTGGAAAACCAATGGTGCGCAGGATACGGGGGTGATTGGGGAAAGGGGGGACCGGGAGCCCGGCTCAGGTCTCCCCAATGCGTCAGCTGCGTGCCAgccgcc                                                                                                              <  2:281299/146‑1 (MQ=255)
                            tGGGCCGATTTGCAGCGGTTGCAAGAGTTTCATTGCAGCTCCAGAGGTGGAAAACCAATGGTGCGCACGATACCGGGGTGATTGGGGAAAGGGGGGACCGGGAGCCCGGCTCAGGGCTCCCCAATGCGTCAGCTGCGTGCCAgccgcc                                                                                                              >  3:135422/1‑148 (MQ=255)
                                 cGATTTGCAGCGGTTGCAAGAGTTTCATTGGAGCTCCAGTGGTGTAAAACCAATGGTGCGCACTATACCGGGGTGATTGGGGAAAGGGGGGACCGGGAGCCCGGCTCAGGGCTCCCCAATGCGTCAGCTGCGTGCCAGCCGCCCACTca                                                                                                        <  4:72000/149‑1 (MQ=255)
                                           cGGTTGCAAGAGTTTCATTGCAGCTCCAGAGGTGGAAAACCAATGGTGCGCACGATACCGGGGTGATTGGGGAAAGGGGGGACCGGGAGCCCGGCTCAGGGCTCCCCAATGCGTCAGCTGCGTGCCAGCCGCCCACTCACATAATgcgc                                                                                              <  4:135422/149‑1 (MQ=255)
                                                    gagTTTCATTGCAGCTCCAGAGGTGGAAAACCAATGGTGCGCACGATACCGGGGTGATTGGGGAAAGGGGGGACCGGGAGCCCGGCTCAGGGCTCCCCAATGCGTCAGCTGCGTGCCAGCCGCCCACTCACATAATGCGCGACGTCAtt                                                                                     >  4:91601/1‑149 (MQ=255)
                                                      gTTTCATTGCAGCTCCAGAGGTGGAAAACCAATGGTGCGCACGATACCGGGGTGTTTGGGGAAAGGGGGGACCGGGAGCCCGGCTCAGGGCTCCCCAATGCGTCAGCTGCGTGCCAGCCGCCCACTCACATAATGCGCGACGTCATTGa                                                                                   <  3:91601/149‑1 (MQ=255)
                                                       tttCATTGCAGCTCCAGAGGTGGAAAACCAATGGTGCGCACGATACCGGGGTGATTGGGGAAAGGGGGGACCGGGAGCCCGGCTCAGGGCTCCCCAATGCGTCAGCTGCGTGCCAGCCGCCCACTCACAt                                                                                                     <  3:200060/130‑1 (MQ=255)
                                                       tttCATTGCAGCTCCAGAGGTGGAAAACCAATGGTGCGCACGATACCGGGGTGATTGGGGAAAGGGGGGACCGGGAGCCCGGCTCAGGGCTCCCCAATGCGTCAGCTGCGTGCCAGCCGCCCACTCACAt                                                                                                     >  4:200060/1‑130 (MQ=255)
                                                        ttCATTGCAGCTCCAGAGGTGGAAAACCAATGGTGCGCACGATACCGGGGTGATTGGGGAAAGGGGGGACCGGGAGCCCGGCTCAGGGCTCCCCAATGCGTCAGCTGCGTGCCAGCCGCCCACTCACATAATGCGCGACGTCATTGAAg                                                                                 >  5:60463/1‑149 (MQ=255)
                                                          cATTGCAGCTCCAGCGGTGGAAAACCAAGGGTGCGCACGATACCGGGGTGATTGGGGAACGGGGGGACCGGGAGCCCGGCTCAGGGCTCCCCAATGCGTCAGCTGCGTGCCAGCCGCCCACTCACATAATGCGCGACGTCATTGAAGc                                                                                <  3:284438/148‑1 (MQ=255)
                                                           aTTGCAGCTCCAGAGGTGGAAAACCAATGGTGCGCACGATACCGGGGTGATTGGGGAAAGGGGGGACCGGGAGCCCGGCTCAGGGCTCCCCAATGCGTCAGCTGCGTGCCAGCCGCACACTCACACAATGCGCGACGTCATTGAAGccc                                                                              >  5:329581/1‑149 (MQ=255)
                                                           aTTGCACCTCCAGAGGTGGAAACCCAATGGTGCGCACGATACCGGGGTTATTGGGGAAAGGGGGGACCGGGAGCCCGGCTCAGGGCTCCCCAATGCGTCAGCTGCGTGCCAGCCGCCCACTCACATAATGCGCGACGTCATTGAAGccc                                                                              <  6:305356/149‑1 (MQ=255)
                                                            ttGCAGCTCCAGAGGTGGAAAACCAATGGTGCGCACGATACCGGGGTGATTGGGGAAAGGGGGGACCGGGAGCCCGGCTCAGGGCTCCCCAATGCGTCAGCTGCGTGCCAGCCGCCCACTCACATAATGCGCGACGTCATTGAAGCCCg                                                                             >  6:183214/1‑149 (MQ=255)
                                                               cAGCTCCAGCTGTGGAAAACCAATGGTGCGCACGATACCGGGGTGTTTGGGGAAAGGGGGGACCGGGAGCCCGGCTCAGGGCTCCCCAATGCGTCAGCTGCGTGCCAGCCGCCCACTCACATAATGCGCGACGTCATAGAAGCCCGAc                                                                           <  8:116882/148‑1 (MQ=255)
                                                                    cccgcggTGGAAAACCAATGGTGCGCACGATACCGGGGTTATTGGGGAAAGGGGGGACCGGGAGCCCGGCTCAGGGCTCCCCAATGCGTCAGCTGCGTGCCAGCCGCCCACTCACATAATgc                                                                                                <  2:324327/117‑1 (MQ=255)
                                                                    ccAGAGGTGGAAAACCAATGGTGCGCACGATACCGGGGTGATTGGGGAAAGGGGGGACCGGGAGCCCGGCTCAGGGCTCCCCAATGCGTCAGCTGCGTGCCAGCCGCCCACTCACATAATgc                                                                                                >  1:324327/1‑122 (MQ=255)
                                                                      agagGTGGAAAACCAATGGTGCGCACGATACCGGGGTGATTGGGGAAAGGGGGGACCGGGAGCCCGGCTCAGGGCTCCCCAATGCGTCAGCTGCGTGCCAGCCGCCCACTCACATAATGCGCGACGTCATTGAAGCCCGACGTGGACg                                                                    >  1:167597/1‑148 (MQ=255)
                                                                                aaCCAATGGTGCGCACGATACCGGGGTGATTGGGGAACGGGGGGACCGGGAGCCCGGCTCAGGGCTCCCCAATGCGTCAGCTGCGTGCCAGCCGCCCACTCACATAATGCGCGACGTCATTGAAGCCCGACGTGGACGCATGCCCTGGc                                                         <  5:157678/149‑1 (MQ=255)
                                                                                    aaTGGTGCGCACGATACCGGGGTGATTGGGGAAAGGGGGGACCGGGAGCCCGGCTCAGGGCTCCCCAATGCGTCAGCTGCGTGCCAGCCGCCCACTCACATAATGCGCGACGTCATTGAAGCCCGACGTGGACGCATGCCCTGGCGTTa                                                     >  2:145009/1‑149 (MQ=255)
                                                                                               cGATACCGGGGTGATTTGGGAAAGGGGGTACCGGGAGCCCGGCGCAGGGCTCCCCAATGCGTCCGCTGCGTGCCAGCCGCCCACTCACACAATGCGCGACGTCATTGAAGCCCGACGTGGACGCATGCCCTGGCGTTACCAACGAATCg                                          <  6:329581/149‑1 (MQ=255)
                                                                                                  tACCGGGGTGATTGGGGAAAGGGGGGACCGGGAGCCCGGCTCAGGGCTCCCCAATGCGTCAGCTGCGTGCCAGCCGCCCACTCACATAATGCGCGACGTCATTGAAGCCCGACGTGGCCGCATGCCCTGGCGTTACCAACGAATCGAt                                        >  8:326700/1‑148 (MQ=255)
                                                                                                    ccGGGGTGATTGGGGAAAGGGGGGACCGGGAGCCCGGCTCAGGGCTCCCCAATGCGTCAGCTGCGTGCCAGCCGCCCACTCACATAATGCGCGACGTCATTGAAGCCCGACGTGGACGCATGCCCTGGCGTTACCAACGAATCGATGa                                      <  1:145009/148‑1 (MQ=255)
                                                                                                         gTGATTGGGGAAAGGGGGGACCGGGAGCCCGGCTCAGGGCTCCCCAATGCGTCAGCTGCGTGCCAGCCGCCCACTCACATAATGCGCGACGTCATTGAAGCCCGACGTGGACGCATGCCCTGGCGTTACCAACGAATCGATGAATGCCt                                <  5:183214/149‑1 (MQ=255)
                                                                                                                         gggACCGGGAGCCCGGCTCAGGGCTCCCCAATGCGTCAGCTGCGTGCCAGCCGCCCACTCACATAATGCGCTACGTCATTGAAGCCCGACGTGGACGCATGCCCTGGc                                                         <  7:224635/108‑1 (MQ=255)
                                                                                                                         gggACCGGGAGCCCGGCTCAGGGCTCCCCAATGCGTCAGCTGCGTGCCAGCCGCCCACTCACAATATGCGCGACGTCATTGAAGCCCGACGTGGACGCATGCCCTGGc                                                         >  8:224635/1‑108 (MQ=255)
                                                                                                                          ggACCGGGAGCCCGGCTCAGGGCTCCCCAATGCGTCAGCTGCGTGCCAGCCTCCCACTCACATAATGCGCGACGTCATTGAAGCCCGACGTGGACGCATGCCCTGGCGTTACCAAc                                                <  1:143808/116‑1 (MQ=255)
                                                                                                                          ggACCGGGAGCCCGGCTCAGGGCTCCCCAATGCGTCAGCTGCGTGCCAGCCGCCCACTCACATAATGCGCGACGTCATTGAAGCCCGACGTGGACGCATGCCCTGGCGTTACCAACGAATCGATGAATGCCTCATCCTCCGCGGTAATc               <  7:326700/149‑1 (MQ=255)
                                                                                                                          ggACCGGGAGCCCGGCTCAGGGCTCCCCAATGCGTCAGCTGCGTGCCAGCCCCCCACTCACATAATGCGCGACGTCATTGAAGCCCGACGTGGACGCATGCCCTGGCGTTACCAAc                                                >  2:143808/1‑116 (MQ=255)
                                                                                                                               gggAGCCCTGAGCCGGGCTCCCCAATGCGTCAGCTGCGTGCCAGCCGCCCACTCACATAATGCGCGACGTCATTGAAGCCCGACGTGGACGCATGCCCTGGCGTTACCAACGAATCGATGAATGCCTCATCCTCCGCGGTAATCTTcgc          >  3:124199/1‑149 (MQ=255)
                                                                                                                               gggAGCCCTGAGCCGGGCTCCCCAATGCGTCAGCTGCGTGCCAGCCGCCCACTCACATAATGCGCGACGTCATTGAAGCCCGACGTGGACGCATGCCCTGGCGTTACCAACGAATCGATGAATGCCTCATCCTCCGCGGTAATCTTcgc          >  8:275045/1‑149 (MQ=255)
                                                                                                                               gggAGCCCGGCTCAGGGCTCCCCAATGCGTCAGCTGCGTGCCAGCCGCCCACTCACATAATGCGCGACGTCATTGAAGCCCGACGTGGACGCATGCCCTGGCGTTACCAACGAATCGATGCATGCCTCATCCTCCGCGGTAATCTTcgc          >  6:81698/1‑149 (MQ=255)
                                                                                                                                ggAGCCCGGCTCAGGGCTCCCCAATGCGTCAGCTGCGTGCCAGCCGCCCACTCACATAATGCGCGACGTCATTGAAGCCCGCCGTGGACGCATGCCCTGGCGTTACCAACGAATCGATGAATGCCTCATCCTCCGCGGTAATCTTcgcc         >  3:230720/1‑149 (MQ=255)
                                                                                                                                    cccGGCTCAGGGCTCCCCAATGCGTCAGCTGCGTGCCAGCCGCCCACTCACATAATGCGCGACGTCATTGAAGCCCGACGTGGACGCATGCCCTGGCGTTACCAACGAATCGATGAATGCCTCATCCTCCGCGGTAATCTTCgccgcc      >  2:47016/1‑148 (MQ=255)
                                                                                                                                      cGGCTCAGGGCTCCCCAATGCGTCAGCTGCGTGCCAGCCGCCCACTCACATAATGCGCGACGTCATTGAAGCCCGACGTGGACGCATGCCCTGGCGTTACCAACGAATCGATGAATGCCTCATCCTCCGCGGTAATCTTCGCCGCCAAt   >  5:117470/1‑149 (MQ=255)
                                                                                                                                       ggCTCAGGGCTCCCCAATGCGTCAGCTGCGTGCCAGCCGCCCACTCACATAATGCGCGACGTCATTGAAGCCCGACGTGGACGCATGCCCTGGCGTTACCAACGAATCGATGAATGCCTCATCCTCCGCGGTAATCTTCGCCGCCAATg  >  6:116187/1‑149 (MQ=255)
                                                                                                                                          |                                                                                                                                                 
CCATGAATACGCGGTTGGGCAGGGTGAGTGGGCCGATTTGCAGCGGTTGCAAGAGTTTCATTGCAGCTCCAGAGGTGGAAAACCAATGGTGCGCACGATACCGGGGTGATTGGGGAAAGGGGGGACCGGGAGCCCGGCTCAGGGCTCCCCAATGCGTCAGCTGCGTGCCAGCCGCCCACTCACATAATGCGCGACGTCATTGAAGCCCGACGTGGACGCATGCCCTGGCGTTACCAACGAATCGATGAATGCCTCATCCTCCGCGGTAATCTTCGCCGCCAATG  >  NC_002947/3596945‑3597228

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