Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_002947 4,740,819 2 bp→CT intergenic (+56/+696) gltA → / ← yeiW citrate synthase/putative oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,740,8190TC95.5% 56.8 / ‑4.1 22intergenic (+56/+697)gltA/yeiWcitrate synthase/putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (12/9);  total (12/10)
Fisher's exact test for biased strand distribution p-value = 4.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.28e-01
*NC_0029474,740,8200GT95.5% 64.2 / ‑2.9 22intergenic (+57/+696)gltA/yeiWcitrate synthase/putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base T (12/9);  minor base . (0/1);  total (12/10)
Fisher's exact test for biased strand distribution p-value = 4.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.49e-01

GCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTTATTTTTCTGGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTT  >  NC_002947/4740681‑4740946
                                                                                                                                          ||                                                                                                                              
gCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCTGTGCCGGCct                                                                                                                      <  6:173339/148‑1 (MQ=255)
 cccGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCct                                                                                                                      >  3:161003/1‑148 (MQ=255)
    gTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCGCtt                                                                                                                   >  6:95838/1‑148 (MQ=255)
     tACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCAGCCTCTTc                                                                                                                  >  6:93374/1‑148 (MQ=255)
              ggCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTg                                                                                                         >  3:159212/1‑148 (MQ=255)
                      gcgcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGaa                                                                                                 <  8:120730/148‑1 (MQ=255)
                         ccAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAagag                                                                                              <  2:12148/148‑1 (MQ=255)
                         ccAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAagag                                                                                              <  2:9130/148‑1 (MQ=255)
                           aGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGGCCGCGAAGAGGc                                                                                            <  4:61037/148‑1 (MQ=38)
                                tACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCGTCGCGGCCTTGCCCGCGCATcggccggc                                                                                        >  5:18140/1‑147 (MQ=255)
                                   acCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGCGGCCGGcaccg                                                                                    >  7:117620/1‑146 (MQ=255)
                                                   gcgACATCGTGGCCCTGAAAGACCGCAAGTACGCTGCAGTAGCCGAACGCAAGACGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGCGCCAGCCTCTTCGCGGCCTTGCCCGCTAAGAGGCCGGCACAGACACCACAAAACCTCa                                                                    <  2:43809/148‑1 (MQ=255)
                                                     gACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCAGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCATGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGt                                                                  >  8:158910/1‑148 (MQ=39)
                                                                 cTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCATCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGCCACCACCAAACCTCAGTTTTTCTacgccc                                                      >  6:2330/1‑148 (MQ=255)
                                                                          cGCAAGTACGCTTCAGTAGCCTAACGCAACAAGGCTGCCCTCGGGATGCGTTTTTTATTTGGCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCa                                             <  1:3432/148‑1 (MQ=255)
                                                                          cGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTAGCGGGCTTGGCCGCGAAGAGGCCGGCACAGAAACCACAAAAACTAAGTTTTTCTCAGACCGCTGACGCa                                             >  2:133907/1‑148 (MQ=37)
                                                                                     tGCAGTAGCCCAACGCAACACGGCGGCCCTCGGGCTGCTTTTTTTATTTGTCCTGGGCCGGCCGCTTCGCGGCCTGGCCCGCTAAGCGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAgg                                  <  3:164309/148‑1 (MQ=255)
                                                                                        aGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGCAGAGGCCGGCCCCGACCCCACAAAACCGCAGTTTTTGTCCGCCCGCTCCCGCCGCCCGTTCAggggt                               >  5:142807/1‑146 (MQ=255)
                                                                                         gTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACCGACCCCCCAAAACCGCAGTTTTTCTCCGCCCGCTCACGCAGACCTTTCAGGGGAt                              >  7:94699/1‑148 (MQ=255)
                                                                                             ccGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGCAGCGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGaa                          >  3:81055/1‑148 (MQ=255)
                                                                                                 ccGCAACACGGCGGCCCTCGGCCAGCGTTTTTTATTTGCCCTGGCCCCGCCTCTGCGCGGCCTGGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGgc                      <  5:93374/147‑1 (MQ=255)
                                                                                                                     ggCAGCTTTTTTTTTTTGTCCTTTGCCGGCCTGTTGGCGGCGTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACttgtt  <  8:23975/148‑1 (MQ=255)
                                                                                                                                          ||                                                                                                                              
GCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTTATTTTTCTGGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTT  >  NC_002947/4740681‑4740946

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: