Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 790,155 A→T 11.9% intergenic (+43/+17) gdhA → / ← PP_0676 NADP‑dependent glutamate dehydrogenase/transcription elongation factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947790,1550AT11.9% 107.5 / 6.7 42intergenic (+43/+17)gdhA/PP_0676NADP‑dependent glutamate dehydrogenase/transcription elongation factor
Reads supporting (aligned to +/- strand):  ref base A (17/20);  new base T (2/3);  total (19/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 = 9.75e-01

ACTACGGTGAAGAGGCTGATGGCCGGATCAACTACGTCAAGGGGGCGAACATCGCGGGCTTCGTGAAAGTGGCCGATGCGATGCTGGCTCAGGGCGTGGTCTGATCTGAGACCGTAGGGGGCCGCTTTGCGTCGCGATGGGCTGCAAAGCAGCCCCAGGGCCCTCAACCCTCGCCCAGAATCTCAACCACCTCTATCACCTGTTCCCCCACCGGCCGCCGCCAGACCACCTCATCCCCCGGCCCAGCCCCGAGCAAAGCCCGCCCCAATGGCGAACCCCAGTTGATCAACCC  >  NC_002947/790009‑790300
                                                                                                                                                  |                                                                                                                                                 
aCTACGGTGAAGAGGCTGATGGCCGGATCAACTACGTCAAGGGGGCGAACATCGCGGGCTTCGTGAAAGTGGCCGATGCGATGCTGGCTCAGGGCGTGGTCTGATCTGAGACCGTAGGGGGCCGCTTTGCGTCGCGATGGGCTGCAAag                                                                                                                                                 <  1:306260/149‑1 (MQ=255)
    cGGTGAAGAGGCTGATGGCCGGATCAACTACGTCAAGGGGGCGAACATCGCGGGCTTCGTGAAAGTGGCCGATGCGATGCTGGCTCAGGGCGTGGTCTGATCTGAGACCGTAGGGGGCCGCTTTGCGTCGCGATGGGCTGCAAagcagc                                                                                                                                             >  5:182325/1‑149 (MQ=255)
        gAAGAGGCTGATGGCCGGATCAACTACGTCAAGGGGGCGAACATCGCGGGCTTCGTGAAAGTGGCCGATGCGATGCTGGCTCAGGGCGTGGTCTGATCTGAGACCGTAGGGGGCCGCTTTGCGTCGCGATGGGCTGCAAAGCAGccc                                                                                                                                           >  4:136672/1‑147 (MQ=255)
          agagGCTGATGGCCGGATCAACTACGTCAAGGGGGCGAACATCGCGGGCTTCGTGAAAGTGGCCGATGCGATGCTGGCTCAGGGCGTGGTCTGATCTGAGACCGTAGGGGGCCGCTTTGCGTCGCGATGGGCTGCTTTGCAGcccatc                                                                                                                                        <  3:237100/148‑4 (MQ=255)
             ggCTGATGGCCGGATCAAGTACGTCAAGGGGGCGAACATCGCGGGCTTCGTGAAAGTGGCCGATGCGATGCTGGCTCAGGGCGTGGTCTGATCTGAGACCGTAGGGGGCCGCTTTGCGTCGCGATGGGCTGCAAAGCAGCCCCAGGGcc                                                                                                                                    <  6:182325/149‑1 (MQ=255)
             ggCTGATGGCCGGATCAACTACGTCAAGGGGGCGAACATCGCGGGCTTCGTGAAAGTGGCCGATGCGATGCTGGCTCAGGGCGTGGTCTGATCTGAGACCGTAGGGGGCCGCTTTGCGTCGCGATGGGCTGCAAAGCAGCCCCAGGGcc                                                                                                                                    <  3:136672/149‑1 (MQ=255)
                       cGGATCAACTACGTCAAGGGGGCGAACATCGCGGGCTTCGTGAAAGTGGCCGATGCGATGCTGGCTCAGGGCGTGGTCTGATCTGAGACCGTAGGGGGCCGCTTTGCGTCGCGATGGGCTGCTTTGCAGcccatc                                                                                                                                        >  7:172347/1‑132 (MQ=255)
                       cGGATCAACTACGTCAAGGGGGCGAACATCGCGGGCTTCGTGAAAGTGGCCGATGCGATGCTGGCTCAGGGCGTGGTCTGATCTGAGACCGTAGGGGGCCGCTTTGCGTCGCGATGGGCTGCTTTGCAGcccatc                                                                                                                                        <  8:172347/135‑4 (MQ=255)
                        ggATCAACTACGTCAAGGGGGCGAACATCGCGGGCTTCGTGAAAGTGGCCGATGCGATGCTGGCTCAGGGCGTGGTCTGATCTGAGACCGTAGGGGGCCGCTTTGCGTCGCGATGGGCTGCAAAGCAGCCCCAGGGCCCTCAACCCTCg                                                                                                                         <  2:117709/149‑1 (MQ=255)
                             aaCTACGTCAAGGGGGCGAACATCGCGGGCTTCGTGAAAGTGGCCGATGCGATGCTGGCTCAGGGCGTGGTCTGATCTGAGACCGTAGGGGGCCGCTTTGCGTCGCGATGGGCTGCAAAGCAGCCCCAggg                                                                                                                                      >  8:262905/1‑131 (MQ=255)
                             aaCTACGTCAAGGGGGCGAACATCGCGGGCTTCGTGAAAGTGGCCGATGCGATGCTGGCTCAGGGCGTGGTCTGATCTGAGACCGTAGGGGGCCGCTTTGCGTCGCGATGGGCTGCAAAGCAGCCCCAggg                                                                                                                                      <  7:262905/131‑1 (MQ=255)
                              aCTACGTCAAGGGGGCGAACATCGCGGGCTTCGTGAAAGTGGCCGATGCGATGCTGGCTCAGGGCGTGGTCTGATCTGAGACCGTAGGGGGCCGCTTTGCGTCGCGATGGGCTGCTTTGCAGcccat                                                                                                                                         <  1:144529/127‑3 (MQ=255)
                              aCTACGTCAAGGGGGCGAACATCGCGGGCTTCGTGAAAGTGGCCGATGCGATGCTGGCTCAGGGCGTGGTCTGATCTGAGACCGTAGGGGGCCGCTTTGCGTCGCGATGGGCTGCTTTGCAGcccat                                                                                                                                         >  2:144529/1‑125 (MQ=255)
                                  cGTCAAGGGGGCGAACATCGCGGGCTTCGTGAAAGTGGCCGATGCGATGCTGGCTCAGGGCGTGGTCTGATCTGAGACCGTAGGGGGCCGCTTTGCGTCGCGATGGGCTGCAAAGCAGCCCCAGGGCCCTCAACCCTCGCCCAGAAtct                                                                                                               >  7:210108/1‑149 (MQ=255)
                                  cGTCAAGGGGGCGAACATCGCGGGCTTCGTGAAAGTGGCCGATGCGATGCTGGCTCAGGGCGTGGTCTGATCTGAGACCGTAGGGGGCCGCTTTGCGTCGCGATGGGCTGCAAAGCAGCCCCAGGGCCCTCAACCCTCGCCCAGAAtc                                                                                                                <  7:208370/148‑1 (MQ=255)
                                     cAAGGGGGCGAACATCGCGGGCTTCGTGAAAGTGGCCGATGCGATGCTGGCTCAGGGCGTGGTCTGATCTGAGACCGTAGGGGGCCGCTTTGCGTCGCGATGGGCTGCAAAGCAGCCCCAGGGCCCTCAACCCTCGCCCAGCATCTCaa                                                                                                            >  3:13334/1‑149 (MQ=255)
                                       aGGGGGCGAACATCGCGGGCTTCGTGAAAGTGGCCGATGCGATGGTGGCTCAGGGCGTGGTCTGATCTGAGACCGTAGGGGGCCGCTTTGCGTCGCGATGGGCTGCAAAGCAGCCCCAGGGCCCTCAACCCTCGCCCAGAATCTCAacc                                                                                                          <  8:210108/149‑1 (MQ=255)
                                            gCGAACATCGCGGGCTTCGTGAAAGTGGCCGATGCGATGCTGGCTCAGGGCGTGGTCTGATCTGAGACCGTAGGGGGCCGCTTTGCGTCGCGATGGGCTGCAAAGCAGCCCCAGGGCCCTCAACCCTCGCCCAGAATCTCAACCACct                                                                                                      >  2:283655/1‑148 (MQ=255)
                                                        ggCTTCGTGAAAGTGGCCGATGCTATGCTTGCTCAGGGCGTGGTCTGATCTGAGACCGTAGGGGGCCGCTTTGCGTCGCGATGGGCTGCAAAGCAGCGCCAGGGCCCTCAACCCTCGCCCAGAATCTCAACCACCTCTATCACCTGTTc                                                                                         <  1:57432/149‑1 (MQ=255)
                                                         gCTTCGTGAAAGTGGCCGATGCGATGCTGGCTCAGGGCGTGGTCTGATCTGAGACCGTAGGGGGCCGCTTTGCGTCGCGATGGGCTGCAAAGCAGCCCCAGGGCCCTCAACCCTCGCCCAGAATCTCAACCACCTCTATCACCTGTTcc                                                                                        >  5:220542/1‑149 (MQ=255)
                                                                 aaaGTGGCCGATGCGATGCTGGCTCAGGGCGTGGTCTGATCTGAGACCGTAGGGGGCCGCTTTGCGTCGCGATGGGCTGCAAAGCAGCCCCAGGGCCCTCAACCCTCGCCCAGAATCTCAACCACctct                                                                                                    <  4:216252/129‑1 (MQ=255)
                                                                 aaaGTGGCCGATGCGATGCTGGCTCAGGGCGTGGTCTGATCTGAGACCGTAGGGGGCCGCTTTGCGTCGCGATGGGCTGCAAAGCAGCCCCAGGGCCCTCAACCCTCGCCCAGAATCTCAACCACctct                                                                                                    >  3:216252/1‑129 (MQ=255)
                                                                 aaaGTGGCCGATGCGATGCTGGCTCAGGGCGTGGTCTGATCTGAGACCGTAGGGGGCCGCTTTGCGTCGCGATGGGCTGCAAAGCAGCCCCAGGGCCCTCAACCCTCGCCCAGAATCTCAACCACCTCTATCACCTGTTCCCCCACCgg                                                                                >  1:72359/1‑149 (MQ=255)
                                                                 aaaGTGGCCGATGCGATGCTGGCTCAGGGCGTGGTCTGATCTGAGACCGTAGGGGGCCGCTTTGCGTCGCGATGGGCTGCAAAGCAGCCCCAGGGCCCTCAACCCTCGCCCAGAATCTCAACCACCTCTATCACCTGTTCCCCCACCgg                                                                                <  1:283655/149‑1 (MQ=255)
                                                                    gTGGCCGATGCGATGCTGGCTCAGGGCGTGGTCTGATCTGAGACCGTAGGGGGCCGCTTTGCGTCGCGATGGGCTGCAAAGCAGCCCCAggg                                                                                                                                      <  5:220929/92‑1 (MQ=255)
                                                                    gTGGCCGATGCGATGCTGGCTCAGGGCGTGGTCTGATCTGAGACCGTAGGGGGCCGCTTTGCGTCGCGATGGGCTGCAAAGCAGCCCCAggg                                                                                                                                      >  6:220929/1‑92 (MQ=255)
                                                                                 tgcTGGCTCAGGGCGTGGTCTGATCTGAGACCGTAGCGGGCCGCTTTGCGTCGCGATGGGCAGCAAAGCACCCCCAGGGCCAACAACCCTCGCCCAGAATCTCAACCACCTATAAAGCCAGTTCACCCACCCGCCGCCGCAAGaccac                                                                 >  4:30336/1‑148 (MQ=255)
                                                                                          aGGGCGTGGTCTGATCTGAGACCGTAGGGGGCCGCTTTGCGTCGCGATGGGCTGCAAAGCAGCCCCAGGGCCCTCAACCCTCGCCCAg                                                                                                                    >  2:257116/1‑88 (MQ=255)
                                                                                          aGGGCGTGGTCTGATCTGAGACCGTAGGGGGCCGCTTTGCGTCGCGATGGGCTGCAAAGCAGCCCCAGGGCCCTCAACCCTCGCCCAg                                                                                                                    <  1:257116/88‑1 (MQ=255)
                                                                                             gCGTGGTCTGATCTGAGACCGTAGGGGGCCGCTTTGCGTCGCGATGGGCTGCAAAGCAGCCCCAGGGCCCTCAACCCTCGCCCAGAATCTCAACCACCTCTATCACCTGTTCCCCCACCGGCCGCCGCCAGACCACCTCATCCCCCGGc                                                    <  2:72359/149‑1 (MQ=255)
                                                                                                     tgatctgaGACCGTAGGGGGCCGCTTTGCGTCGCGATGGGCTGCAAAGCAGCCCCAGGGCCCTCAACCCTCGCCCAGAATCTCAACCACCTCTATCACCTGTTCCCCCACCGGCCGCCGCCAGaccac                                                                 >  7:3971/1‑128 (MQ=255)
                                                                                                     tgatctgaGACCGTAGGGGGCCGCTTTGCGTCGCGATGGGCTGCAAAGCAGCCCCAGGGCCCTCAACCCTCGCCCAGAATCTCAACCACCTCTATCACCTGTTCCCCCACCGGCCGCCGCCAGaccac                                                                 <  8:3971/128‑1 (MQ=255)
                                                                                                     tgatctgaGACCGTAGGGGGCCGCTTTGCGTCGCGATGGGCTGCAAAGCAGCCCCAGGGCCCTCAACCCTCGCCCAGAATCTCAACCACCTCTATCACCTGTTCCCCCACCGGCCGCCGCCAGACCACCTCATCCCCCGGCCCAGcccc                                            >  4:114728/1‑149 (MQ=255)
                                                                                                          tgaGACCGTAGGGGGCCGCTTTGCGTCGCGATGGGCTGCAAAGCAGCCCCAGGGCCCTCAACCCTCGCCCAGAATCTCAACCACCTCTATCACCTGTTCCCCCACCCGCCCCCGCCAGACCACCACCTCCCCCCGCCCAGcccccacc                                        >  6:142470/1‑144 (MQ=255)
                                                                                                                cGTAGGGGGCCGCTTTGCGTCGCGATGGGCTGCAAAGCAGCCCCAGGGCCCTCAACCCTCGCACAGAATCTCAACCACCTCTATCACCTGTTCCCCCACCGGCCGCCGCCAGACCACCTCATCCCCCGGcc                                                   >  4:104493/1‑131 (MQ=255)
                                                                                                                cGTAGGGGGCCGCTTTGCGTCGCGATGGGCTGCAAAGCAGCCCCAGGGCCCTCAACCCTCGCACAGAATCTCAACCACCTCTATCACCTGTTCCCCCACCGGCCGCCGCCAGACCACCTCATCCCCCGGcc                                                   <  3:104493/131‑1 (MQ=255)
                                                                                                                  tAGGGGGCCGCTTTGCGTCGCGATGGGCTGCAAAGCAGCCCCAGGGCCCTCAACCCTCGCCCAGAATCTCAACCACCTCTATCACCTGTTCCCCCACCGGCCGCCGCCAGACCACCTCATCCCCCGGCCCAGCCCCGAGCAAAgcccgc                               <  6:220542/149‑1 (MQ=255)
                                                                                                                                 cGTCGCGATGGGCTGCAAAGCAGCCCCAGGGCCCTCAACCCTCGCCCAGAATCTCAACCACCTCTATCACCTGTTCCCCCACCGGCCGCCGCCAGACCACCTCATCCCCCGGCCCAGCCCCGAGCAAAGCCCGCCCCAATGGCGAAccc                <  1:143975/149‑1 (MQ=255)
                                                                                                                                  gTCGCGATGGGCTGCAAAGCAGCCCCAGGGCCCTCAACCCTCGCCCAGAATCTCAACCACCTCTATCACCTGTTCCCCCCCCGGCCGCCGCCAGACCCCCTCATCCCCCGGCCCACCCCCGAGCAAACCCCGCCCCAATGGCGCAcccc               >  1:287169/1‑149 (MQ=255)
                                                                                                                                             cTGCAAAGCAGCCCCAGGGCCCTCAACCCTCGCCCAGAATCTCAACCACCTCTATCACCTGTTCCCCCACCGGCCGCCGCCAGACCACCTCATCCCCCGGCCCAGCCCCGAGCAAAGCCCGCCCCAATGGCGAACCCCAGTTGATCAAc    <  4:13334/149‑1 (MQ=255)
                                                                                                                                                cAAAGCAGCCCCAGGGCCCTCAACCCTCGCCCAGAATCTCAACCACCTCTATCACCTGTTCCCCCACCGGCCGCCGCCAGACCACCTCATCCCCCGGCCCAGCCCCGAGCAAAGCCCGCCCCAATGGCGAACCCCAGTTGATCAAccc  <  3:30336/148‑1 (MQ=255)
                                                                                                                                                cAAAGCAGCCCCAGGGCCCTCAACCCTCGCCCAGAATCTCAACCACCTCTATCACCTGTTCCCCCACCGGCCGCCGCCAGACCACCTCATCCCCCGGCCCAGCCCCGAGCAAAGCCCGCCCCAATGGCGAACCCCAGTTGATCAAcc   <  6:62294/147‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
ACTACGGTGAAGAGGCTGATGGCCGGATCAACTACGTCAAGGGGGCGAACATCGCGGGCTTCGTGAAAGTGGCCGATGCGATGCTGGCTCAGGGCGTGGTCTGATCTGAGACCGTAGGGGGCCGCTTTGCGTCGCGATGGGCTGCAAAGCAGCCCCAGGGCCCTCAACCCTCGCCCAGAATCTCAACCACCTCTATCACCTGTTCCCCCACCGGCCGCCGCCAGACCACCTCATCCCCCGGCCCAGCCCCGAGCAAAGCCCGCCCCAATGGCGAACCCCAGTTGATCAACCC  >  NC_002947/790009‑790300

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: