Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,597,083 T→G 8.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 8.5% 125.9 / 4.5 47intergenic (‑78/+18)nemA/PP_3174FMN‑linked N‑ethylmaleimide reductase/oxidoreductase
Reads supporting (aligned to +/- strand):  ref base T (22/21);  new base G (2/2);  total (24/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.66e-01

TGAGGGGGGCCATGAATACGCGGTTGGGCAGGGTGAGTGGGCCGATTTGCAGCGGTTGCAAGAGTTTCATTGCAGCTCCAGAGGTGGAAAACCAATGGTGCGCACGATACCGGGGTGATTGGGGAAAGGGGGGACCGGGAGCCCGGCTCAGGGCTCCCCAATGCGTCAGCTGCGTGCCAGCCGCCCACTCACATAATGCGCGACGTCATTGAAGCCCGACGTGGACGCATGCCCTGGCGTTACCAACGAATCGATGAATGCCTCATCCTCCGCGGTAATCTTCGCC  >  NC_002947/3596936‑3597221
                                                                                                                                                   |                                                                                                                                          
tGAGGGGGGCCATGAATCCGCGGTTGGGCTGGGTGAGTGGGCAGATTTGCAGCGGTTGCAAGAGTTTAATTGCAGCTCCATAGGTGGAAAACCAATGGTGCGCACGATACCGGGGTGATTGGGGAAAGGGGGGACCGGGAGCCCGGCTc                                                                                                                                           <  8:340746/149‑1 (MQ=255)
  cggggggCCATGAATCCGCGGTTGGGCAGGGTGAGTGGGCCGATTGGCAGCGGTTGCAAGAGTTTCATTGCAGCTCCAGAGGTGGAAAACCAATGGTGCGCACGATACCGGGGTGATTGGGGAACGGGGGGACCGGGAGCCCGGCTCAg                                                                                                                                         <  8:96741/148‑1 (MQ=255)
      gggCCATGAATACGCGGTTGGGCAGGGTGAGTGGGCCGATTTGCAGCGGTTGCAAGAGTTTCATTGCAGCTCCAGAGGTGGAAAACCAATGGTGCGCACGATACCGGGGTGATTGGGGAAAGGGGGGACCGGGAGCCCGGCTCCGGGCt                                                                                                                                     >  2:86549/1‑149 (MQ=255)
         ccATGATTAGGGGGTGGGGCTGGTGTAGTGTGCCGTTTTGCATGGGTTGAAAGATTTTCATTGCAGCTCCAGAGGTGGATATCCAATGTTGGGAACGATACCGGGGTGATTGGGGAATGGGGGGACCGGGAGCCCTGAGCCGGGCTccc                                                                                                                                  <  6:222292/149‑1 (MQ=255)
          cttcaaTCCGCGGTTGGGCTGGGTGAGTGGGCCGATTTGCAGCGGTTGCAAGAGTTTCATTGCAGCTCCAGAGGTGGAAAACCAATGGTGCGCACGATACCGGGGTGATTGGGGAAAGGGGGGACCGGGAGCCCGGCTCAGGGCTcccc                                                                                                                                 <  6:50151/145‑1 (MQ=255)
          cATGAATACGCGGTTGGGCAGGGTGAGTGGGCCGATTTGCAGCGGTTGCAAGAGTTTCATTGCAGCTCCAGAGGTGGAAAACCAATGGTGCGCACGATACCGGGGTGATTGGGGAAAGGGGGGACCGGGAGCCCGGCTAAGGGCTccc                                                                                                                                  >  7:303241/1‑148 (MQ=255)
                tACGCGGTTGGGCAGGGTGAGTGGGCCGATTTGCAGCGGTTGCAAGAGTTTCATTGCAGCTCCAGAGGTGGAAAACCAATGGTGCGCACGATACCGGGGTGATTGGGGAAAGGGGGGACCGGGAGCCCGGCTCAGGGCTCCCCAATGc                                                                                                                            >  2:144702/1‑148 (MQ=255)
                 tcgcgGGCGGGCCGGGTGAGTGGGCCTATTTGCCGCGGTTGCATCAGTTTCATTGCAGCTCCATCGGTCTAAAACCAATGGTGCGCACGCTTCCGGGGTTATTGGGGAACGGGGGGACCGGGAGCCCGGCTCAGGGCTCCCCAATGCGt                                                                                                                          <  4:137963/148‑1 (MQ=255)
                   gcgGTTGGGCAGGGTGAGTGGGCCGATTTGCAGCGGTTGCAAGAGTTTCATTGCAGCTCCAGAGGTGGAAAACCAATGGTGCGCACGATACCGGGGTGATTGGGGAAAGGGGGGACCGGGAGCCCGGCTCAGGGCTCCCCAATGCGt                                                                                                                          >  3:29361/1‑147 (MQ=255)
                          ggCAGCGTGAGTGGGCCGATTTGCAGCGGTTGCAAGAGTTTCATTGCAGCTCCAGAGGTGGAAAACCAATGGTGCGCACGATACCGGGGTGATTGGGGAAAGGGGGGACCGGGAGCCCGGCTCAGGGCTCCCCAAAGCCTCCTCTGcat                                                                                                                 >  4:193859/1‑147 (MQ=255)
                                     ggggCCGTTTTGCAGCGGTTGCAAGAGTTTCATTGCAGCTCCAGAGGGGGAAAACCTATGGTGCGCACGATACCGGGGTGATTGGGGAAAGGGGGGACCGGGAGCCCGGCTCAGGGCTCCCCAATGCGTCAGCTGCGTGCCAgccgcc                                                                                                       <  4:212985/147‑1 (MQ=255)
                                             tttGCAGCGGTGCCAAGACTTTCATTGCAGCGCCAGAGGTGCAAAACCAATGGTCCGCACGATCCCGGGGTGATTGGGCAAAGGGCGCACCGGGAGCCCGGCTCAGGGCTCCCCAATGCGTCAGCTGCGTGCCAGCCGCCCACTcac                                                                                                <  7:31650/147‑1 (MQ=255)
                                                gccgCGGTTGCAAGAGTTCCATTGCAGCCCCAGAGGTGGAAACCCAATGGTGCGCACGATACCGGGGGGATTGGGGAAAGGGGGCACCGGGAGCCCGGCTCAGGGCTCCCCAATGCGTCAGCTGCGTGCCAGCCGCCCACTCACATaa                                                                                            <  5:370917/145‑1 (MQ=255)
                                                gCAGCGGTTGCAAGAGTTTCATTGCAGCTCCAGAGGTGGAAAACCAATGGTGCGCACGATACCGGGGTGATTGGGGAAAGGGGGGACCGGGAGCCCGGCTCAGGGCTCCCCAATGCGTCAGCTGCGTGCCAGCCGCCCACTCACATAAt                                                                                           >  3:86921/1‑149 (MQ=255)
                                                gCAGCGGTTGCAAGAGTTTCATTGCAGCTCCAGAGGTGGAAAACCAATGGTGCGCACGATACCGGGGTGATTGGGGAAAGGGGGGACCGGGAGCCCGGCTCAGGGCTCCCCAATGCGTCAGCTGCGTGCCAGACGCCCACTCACATAAt                                                                                           >  5:68677/1‑149 (MQ=255)
                                                         gCAAGAGTTTCATTGCAGCTCCAGAGGTGGAAAACCAATGGTGCGCACGATACCGGGGTGATTGGGGAAAGGGGGGACCGGGAGCCCGGCTCAGGGCTCCCCAATGCGTCAGCTGCGTGCCAGCCGCCCACTCACATAATGCGCGACGt                                                                                  >  1:223170/1‑149 (MQ=255)
                                                           aaGAGTTTCATTGCAGCTCCAGAGGTGGAAAACCAATGGTGCGCACGATACCGGGGTGATTGGGGAAAGGGGGGACCGGGAGCCCGGCTCAGGGCTCCCCAATGCGTCAGCTGCGTGCCAGCCGCCCACGCACATAATGCGCGACGTCa                                                                                >  8:268265/1‑149 (MQ=255)
                                                                    aTTGCAGCTCCAGAGGTGGAAAACCAATGGTGCGCACGATACCGGGGTGATTGGGGAAAGGGGGGACCGGGAGCCCGGCTCAGGGCTCCCCAATGCGTCAGCTGCGTGCCAGCCGCCCACTCACATAATGCGCGACGTCATTGAAGccc                                                                       >  8:56149/1‑149 (MQ=255)
                                                                     ttGCAGCTCCAGAGGTGGAAAACCAATGGTGCGCACGATACCGGGGTGATTGGGGAAAGGGGGGACCGGGAGCCCGGCTCAGGGCTCCCCAATGCGTCAGCTGCGTGCCAGCCGCCCACTCACATAATGCGCGACGTCATTGAAGccc                                                                       >  7:238286/1‑148 (MQ=255)
                                                                          gCTCCAGAGGTGGAAAACCAATGGTGCGCACGATACCGGGGTGATTGGGGAAAGGGGGGACCGGGAGCCCGGCTCAGGGCTCCCCAATGCGTCAGCTGCGTGCCAGCCGCCCACTCACATAATGCGCGACGTCATTGAAGCCCGACGTg                                                                 >  8:118529/1‑149 (MQ=255)
                                                                              cAGAGGTGGAAAACCAATGGTGCGCACGATACCGGGGTGATTGGGGAAAGGGGGGACCGGGAGCCCGGCTCAGGGCTCCCCAATGCGTCAGCTGCGTGCCAGCCGCCCACTCACATAATGCGCGACGTCATTGAAGCCCGACGTGGa                                                               <  6:100897/147‑1 (MQ=255)
                                                                               agagGTGGAAAACCAATGGTGCGCACGATACCGGGGTGATTGGGGAAAGGGGGGACCGGGAGCCCGGCTCAGGGCTCCCCAATGCGTCAGCTGCGTGCCAGCCGCCCACTCACATAATGCGCGACGTCATTGAAGCCCGACGTGGACg                                                             >  4:162779/1‑148 (MQ=255)
                                                                                    tGGAAAACCAATGGTGCGCACGATACCGGGGTGATTGGGGAAAGGGGGGACCGGGAGCCCGGCTCAGGGCTCCCCAATGCGTCAGCTGCGTGCCAGCCGCCCACTCACATAATGCGCGACGTCATTGAAGCCCGACGTGGACGCATGcc                                                       >  8:110597/1‑149 (MQ=255)
                                                                                    tGGAAAACCAATGGTGCGCACGATACCGGGGTGATTGGGGAAAGGGGGGACCGGGAGCCCGGCTCAGGGCTCCCCAATGCGTCAGCTGCGTGCCAGCCGCCCACTCACATAATGCGCGACGTCATTGAAGCCCGACGTGGACGCATGc                                                        <  8:303241/148‑1 (MQ=255)
                                                                                        aaaCCAATGGTGCGCACGATACCGGGGTGATTGGGGAAAGGGGGGACCGGGAGCCCGGCTCAGGGCTCCCCAATGCGTCAGCTGCGTGCCAGCCGCCCACTCACATAATGCGCGACGTCATTGAAGCCCGACGTGGACGCATGCCCTg                                                    <  1:144702/148‑1 (MQ=255)
                                                                                         aaCCAATGGTGCGCACGATACCGGGGTGATTGGGGAAAGGGGGGACCGGGAGCCCGGCTCAGGGCTCCCCAATGCGTCAGCTGCGTGCCAGCCGCCCACTCACATAATGCGCGACGTCATTGAAGCCCGACGTGGACGCATGCCCTGGc                                                  >  3:239937/1‑149 (MQ=255)
                                                                                            cAATGGTGCGCACGATACCGGGGTGATTGGGGAAAGGGGGGACCGGGAGCCCGGCTCAGGGCTCCCCAATGCGTCAGCTGCGTGCCAGCCGCCCACTCACATAATGCGCGACGTCATTGAAGCCCGACGTGGACGCATGCCCTGGCGtt                                               <  4:86921/149‑1 (MQ=255)
                                                                                                 gTGCGCACGATACCGGGGTGATTGGGGAAAGGGGGGACCGGGAGCCCGGCTCAGGGCTCCCCAATGCGTCAGCTGCGTGCCAGCCGCCCACTCACATAATGCGCGACGTCATTGAAGCCCGACGTGGACGCATGCCCTGGCGt                                                <  5:288496/143‑1 (MQ=255)
                                                                                                 gTGCGCACGATACCGGGGTGATTGGGGAAAGGGGGGACCGGGAGCCCGGCTCAGGGCTCCCCAATGCGTCAGCTGCGTGCCAGCCGCCCACTCACATAATGCGCGACGTCATTGAAGCCCGACGTGGACGCATGCCCTGGCGt                                                >  6:288496/1‑143 (MQ=255)
                                                                                                      cACGATACCGGGGTGATTGGGGAAAGGGGGGACCGGGAGCCCGGCTCAGGGCTCCCCAATGCGTCAGCTGCGTGCCAGCCGCCCACTCACATAATGCGCGACGTCATTGAAg                                                                          >  4:114359/1‑112 (MQ=255)
                                                                                                      cACGATACCGGGGTGATTGGGGAAAGGGGGGACCGGGAGCCCGGCTCAGGGCTCCCCAATGCGTCAGCTGCGTGCCAGCCGCCCACTCACATAATGCGCGACGTCATTGAAg                                                                          <  3:114359/112‑1 (MQ=255)
                                                                                                            aCCGGGGTTATTGGGGAATGGGGGGACCGGGAGCCCGGCTCAGGGCTCCCCAATGCGTCAGCTGCGTGCCAGCCGCCCACTCACATAATGCGCGACGTCATTGAAGCCCGACGTGGACGCATGCCCTGGCGTTACCAACGAATCGATg                                <  3:193859/148‑1 (MQ=255)
                                                                                                              cGGGGTGATTGGGGAAAGGGGGGACCGGGAGCCCGGCTCAGGGCTCCCCAATGCGTCAGCTGCGTGCCAGCCGCCCACTCACATAATGCGCGACGTCATTGAAGCCCGACGTGGACGCATGCCCTGGCGTTACCAACGAATCGATGAAt                             <  8:238286/149‑1 (MQ=255)
                                                                                                              cGGGGTGATTGGGGAAAGGGGGGACCGGGAGCCCGGCTCAGGGCTCCCCAATGCGTCAGCTGCGTGCCAGCAGCCCACTCACATAATGCGCGACGTCATTGAAGCCCGACGTGGACGCATGCCCTGGCGTTACCAACGAATCGATGAAt                             <  2:252600/149‑1 (MQ=255)
                                                                                                              cGGGGTCATTGGGGAAAGGGGGGACCGGGAGCCCGGCTCAGGGCTCCCCAATGCGTCAGCTGCGTGCCAGCCGCCCACTCACATAATGCGCGACGTCATTGAAGCCCGACGTGGACGCATGCCCTGGCGTTACCAACGAATCGATGAAt                             <  7:110312/149‑1 (MQ=255)
                                                                                                               ggggTGATTGGGGAAAGGGGGGACCGGGAGCCCGGCTCAGGGCTCCCCAATGCGTCAGCTGCGTGCCAGCCGCCCACTCACATAATGCGCGACGTCATTGAAGCCCGACGTGGACGCATGCCCTGGCGTTACCAACGAATCGATGAAt                             >  8:110312/1‑148 (MQ=255)
                                                                                                                 ggTGATTGGGGAAAGGGGGGACCGGGAGCCCGGCTCAGGGCTCCCCAATGCGTCAGCTGCGTGCCAGCCGCCCACTCACATAATGCGCGACGTCATTGAAGCCCGACGTGGACGCATGCCCTg                                                    >  3:284315/1‑123 (MQ=255)
                                                                                                                 ggTGATTGGGGAAAGGGGGGACCGGGAGCCCGGCTCAGGGCTCCCCAATGCGTCAGCTGCGTGCCAGCCGCCCACTCACATAATGCGCGACGTCATTGAAGCCCGACGTGGACGCATGCCCTg                                                    <  4:284315/123‑1 (MQ=255)
                                                                                                                     aTTGGGGAAAGGGGGGACCGGGAGCCCGGCTCAGGGCTCCCCAATGCGTCAGCTGCGTGCCAGCCGCCCACTCACATAATGCGCGACGTCATTGAAGCCCGACGTGGACGCATGCCCTGGCGTTACCAACGAATCGATGAATGCCTCa                       <  4:239937/148‑1 (MQ=255)
                                                                                                                     aTTGGGGAAACGGGGGACCGGGAGCCCGGCTCAGGGCTCCCCAATGCGTCAGCTGCGCGCCAGCCGCCCACTCACATAATGCGCGACGTCATTGAAGCCCGACGTGGACGCATGCCCTGGCGTTACCAACGAATCGATGAATGCCTCa                       <  4:29361/148‑1 (MQ=255)
                                                                                                                                 ggggACCGGGAGCCCGGCTCAGGGCTCCCCAATGAGTCAGCTGCGTGCCAGCCGCCCACTCACATAATGCGCGACGTCATTGAAGCCCGACGTGGACGCATGCCCTGGCGTTACCAACGAATCGATGAATGCCTCATCCTCCGCGGTaa          >  4:254729/1‑149 (MQ=255)
                                                                                                                                  gggACCGGGAGCCCGGCTCAGGGCTCCCCAATGCGTCAGCTGCGTGCCAGCCGCCCACTCACATAATGCGCGACGTCATTGAAGCCCGACGTGGACGCATGCCCTGGCGTTACCAACGAATCGATGAATGCCTCATCCTCCGCGGTAAt         <  7:110597/149‑1 (MQ=255)
                                                                                                                                  gcattggggAGCCCTGAGCCGGGCTCCCCAATGCGTCAGCTGCGTGCCAGCCGCCCACTCACATAATGCGCGACGTCATTGAAGCCCGACGTGGACGCATGCCCTGGCGTTACCAACGAATCGATGAATGCCTCATCCTCCGCGGTaa          >  4:228768/7‑148 (MQ=255)
                                                                                                                                       ggggAGCCCTGAGCCGGGCTCCCCAATGCGTCAGCTGCGTGCCAGCCGCCCACTCACATAATGCGCGACGTCATTGAAg                                                                          >  3:258325/2‑79 (MQ=38)
                                                                                                                                       ggggAGCCCTGAGCCGGGCTCCCCAATGCGTCAGCTGCGTGCCAGCCGCCCACTCACATAATGCGCGACGTCATTGAAg                                                                          <  4:258325/78‑1 (MQ=38)
                                                                                                                                         ggAGCCCGGCTCAGGGCTCCCCAATGCGTCAGCTGCGTGCCAGCCGCCCACTCACATAATGCGCGACGTCATTGAAGCCCGACGTGGACGCATGCCCTGGCGTTACCAACGAATCGATGAATGCCTCATCCTCCGCGGTAATCTTcgcc  >  7:273976/1‑149 (MQ=255)
                                                                                                                                         ggAGCCCGGCTCAGGGCTCCCCAATGCGTCAGCTGCGTGCCAGCCGCCCACTCACATAATGCGCGACGTCATTGAAGCCCGACGTGGACGCATGCCCTGGCGTTACCAACGAATCGATGAATGCCTCATCCTCCGCGGTAATCTTcgcc  >  1:322467/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                          
TGAGGGGGGCCATGAATACGCGGTTGGGCAGGGTGAGTGGGCCGATTTGCAGCGGTTGCAAGAGTTTCATTGCAGCTCCAGAGGTGGAAAACCAATGGTGCGCACGATACCGGGGTGATTGGGGAAAGGGGGGACCGGGAGCCCGGCTCAGGGCTCCCCAATGCGTCAGCTGCGTGCCAGCCGCCCACTCACATAATGCGCGACGTCATTGAAGCCCGACGTGGACGCATGCCCTGGCGTTACCAACGAATCGATGAATGCCTCATCCTCCGCGGTAATCTTCGCC  >  NC_002947/3596936‑3597221

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