Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,960,613 T→C 12.7% intergenic (‑229/‑100) mnxG ← / → PP_3491 manganese‑oxidizing multicopper oxidase/enoly‑coenzyme A hydratase/isomerase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,960,6130TC12.7% 112.9 / 7.6 47intergenic (‑229/‑100)mnxG/PP_3491manganese‑oxidizing multicopper oxidase/enoly‑coenzyme A hydratase/isomerase family protein
Reads supporting (aligned to +/- strand):  ref base T (20/21);  new base C (3/3);  total (23/24)
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.05e-01

ACTTGTACGATCAACCGCAACGGATTGGGGAGGATCCCCCAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGAAAATCAATAACAAAGCAGGGAGAACGCAATGACTGCACACGCTCATACCCCGGCAGCAACAGACCATG  >  NC_002947/3960466‑3960752
                                                                                                                                                   |                                                                                                                                           
aCTTGTACGATCAACCGCAACGGATTGGGGAGGATCCCCCAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTg                                                                                                                                            >  5:184268/1‑149 (MQ=255)
    gTACGATCAACCGCAACGGATTGGGGAGGATCCCCCAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCg                                                                                                                                         <  3:154957/148‑1 (MQ=255)
    gTACGATCAACCGCAACGGATTGGGGAGGATCCCCCAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTCCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGa                                                                                                                                        >  7:64611/1‑149 (MQ=255)
           cAACCGCAACGGATTGGGGAGGATCCCCCAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTTCTCCACCAGCCCGCGCTGGCCTCAAACAAGGTAAAAACGCCCGGTACCCATTCACAATTGTCGATTGTGc                                                                                                                                  >  1:247998/1‑148 (MQ=255)
               cGCAACGGATTGGGGAGGATCCCGCAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGc                                                                                                                             <  5:247729/149‑1 (MQ=255)
                gCAACGGATTGGGGAGGATCCCCCAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGc                                                                                                                             >  3:262929/1‑148 (MQ=255)
                       aTTGGGGAGGATCCCCCAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCg                                                                                                                      >  2:279428/1‑148 (MQ=255)
                          gcggaggaTCCCCCAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGGGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCAt                                                                                                                   <  2:243437/146‑1 (MQ=255)
                             gaggaTCCCCCAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGc                                                                                                                 >  7:26015/1‑147 (MQ=255)
                                gaTCCCCCAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGACGATTGTGCAAGGCAAATTCGCATGCCAGGa                                                                                                            >  7:212901/1‑149 (MQ=255)
                                         gcagGTTCCCCCGATTGTGGGGACGTAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGAt                                                                                                   <  2:269468/149‑1 (MQ=255)
                                            ccgtCCCCCGATTCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGAtcct                                                                                                >  8:248972/4‑149 (MQ=255)
                                                  cccGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCg                                                                                          <  1:279428/149‑1 (MQ=255)
                                                      aaTCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACt                                                                                      <  4:262929/149‑1 (MQ=255)
                                                           ggggACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGa                                                                                 <  5:104609/149‑1 (MQ=255)
                                                           ggggACGGAATCCAGTGGGGGCAAAGGCTGCTCCTCAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGa                                                                                 <  8:64611/149‑1 (MQ=255)
                                                             ggACGGAATCCCGTGGGGGCAAAGGCTACTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACt                                                                               >  1:182790/1‑149 (MQ=255)
                                                                 ggAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCAc                                                                           >  6:111615/1‑149 (MQ=255)
                                                                      cccGTGGGGGCAAAGGCTGCTCGACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCGGCGATCCTCCAGCGAACTCCAGACTTAACTTgg                                                                       >  7:205648/1‑148 (MQ=255)
                                                                       ccGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGa                                                                      <  7:248972/148‑1 (MQ=255)
                                                                              ggCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGAAAATCAAt                                                              <  2:182790/149‑1 (MQ=255)
                                                                                 aaaGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGCCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGAAAATCAATAAc                                                           <  8:212901/149‑1 (MQ=255)
                                                                                      ctgctCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGAAAATCAATAACAAAg                                                       >  3:319082/1‑148 (MQ=255)
                                                                                         ctCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGAAAATCAATAACAAAGCAgg                                                   <  8:26015/149‑1 (MQ=255)
                                                                                          tCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGAAAATCAATAACAAAGCAgg                                                   >  3:181228/1‑148 (MQ=255)
                                                                                             acaAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGAAAATCAATAACAAAGCAGGgag                                                >  3:29027/1‑148 (MQ=255)
                                                                                              caAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGAAAATCAATAACAAAGCAGGGAGaa                                              <  8:51716/149‑1 (MQ=255)
                                                                                                           gCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGAAAATCAATAACAAAGCAGGGAGAACGCAATGACTGca                                 >  7:229541/1‑149 (MQ=255)
                                                                                                            ccTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGGTTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGAAAATCAATAACAAAGCAGGGAGAACGc                                           <  7:219767/138‑1 (MQ=255)
                                                                                                            ccTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGAAAATCAATAACAAAGCAGGGAGAACGc                                           >  8:219767/1‑138 (MQ=255)
                                                                                                            ccTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGAAAATCAATAACAAAGCAGGGAGAACGCAATGACTGca                                 >  3:285658/1‑148 (MQ=255)
                                                                                                             cTCACACAAGGTAAATACGCACGGTAACCTTTCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCATCCAGAGAACTCCAGACTTCACTTGGAAAATCAATAACAAAGCAGGGAGAACGCAATGACTGcaca                               <  4:181228/149‑1 (MQ=255)
                                                                                                                      ggTATAATCGGACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGAAAATCAATAACAAAGCAGGGAGAACGCAATGACTGCACACGCTCATAc                      <  6:184268/149‑1 (MQ=255)
                                                                                                                        tAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGAAAATCAATAACAAAGCAGGGAGAACGCAATGACTGCACACGCTCATAccc                    <  2:247998/149‑1 (MQ=255)
                                                                                                                         aaaaaCGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGAAAATCAATAACAAAGCAGGGAGAACGCAATGACTGCACACGCTCATAcccc                   >  2:207804/1‑149 (MQ=255)
                                                                                                                             acggacggTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGAAAATCAATAACAAAGCAGGGAGAACGCAATGACTGCACACGCTCATACCCCGgca               <  4:167993/149‑1 (MQ=255)
                                                                                                                                  cGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCa                                                                                             >  3:227615/1‑66 (MQ=255)
                                                                                                                                  cGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCa                                                                                             <  4:227615/66‑1 (MQ=255)
                                                                                                                                        ccTTGCACAATCGACAATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGAAAATCAATAACAAAGCAGGGAGAACGCAATGACTGCACACGCTCATACCCCGGCAGCAACAGACCa    >  6:176210/1‑149 (MQ=255)
                                                                                                                                         cTTGCACAATCGACAATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGAAAATCAATAACAAAGCAGGGAGAACGCAATGACTGCACACGCTCATACCCCgg                 <  8:184184/135‑1 (MQ=255)
                                                                                                                                         cTTGCACAATCGACAATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGAAAATCAATAACAAAGCAGGGAGAACGCAATGACTGCACACGCTCATACCCCGGCAGCAACAGACCAt   <  5:176210/149‑1 (MQ=255)
                                                                                                                                         cTTGCACAATCGACAATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGAAAATCAATAACAAAGCAGGGAGAACGCAATGACTGCACACGCTCATACCCAgg                 >  7:184184/1‑135 (MQ=255)
                                                                                                                                          ttGCACAATCGACAATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGAAAATCAATAACAAAGc                                                      >  6:312208/1‑97 (MQ=255)
                                                                                                                                          ttGCACAATCGACAATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGAAAATCAATAACAAAGc                                                      <  5:312208/97‑1 (MQ=255)
                                                                                                                                            gCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGAAAATCAATAACAAAGCAGGGAGAACGCAATGACTGCACACGCTCATACCCCGGCAGCAACAGACCATg  <  8:205648/147‑1 (MQ=255)
                                                                                                                                              acaATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGAAAATCAAt                                                              <  6:3271/85‑1 (MQ=255)
                                                                                                                                              acaATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGAAAATCAAt                                                              >  5:3271/1‑85 (MQ=255)
                                                                                                                                                   |                                                                                                                                           
ACTTGTACGATCAACCGCAACGGATTGGGGAGGATCCCCCAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGAAAATCAATAACAAAGCAGGGAGAACGCAATGACTGCACACGCTCATACCCCGGCAGCAACAGACCATG  >  NC_002947/3960466‑3960752

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