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,8190TC96.8% 89.6 / ‑4.2 31intergenic (+56/+697)gltA/yeiWcitrate synthase/putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (15/15);  total (15/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.56e-01
*NC_0029474,740,8200GT96.7% 97.3 / NA 31intergenic (+57/+696)gltA/yeiWcitrate synthase/putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base T (15/14);  total (15/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTTATTTTTCTGGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAA  >  NC_002947/4740680‑4740940
                                                                                                                                           ||                                                                                                                        
aGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCg                                                                                                                    <  3:128431/146‑1 (MQ=255)
aGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGc                                                                                                                  <  3:159631/148‑1 (MQ=255)
              cGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCtt                                                                                                    <  8:105381/148‑1 (MQ=255)
                gCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCCGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGc                                                                                                  >  1:14389/1‑148 (MQ=255)
                gCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGc                                                                                                  <  1:120130/148‑1 (MQ=255)
                 ccGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCACGCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCACATAGCGGccattcc                                                                                                 >  7:138155/1‑143 (MQ=255)
                  cGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGACCGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGccc                                                                                                >  1:36264/1‑148 (MQ=255)
                       gcgcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGaa                                                                                           <  6:108377/148‑1 (MQ=255)
                        cgcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCT‑GTGCCGGCCTCTTCGCGGGCAAGGCCGCGAaga                                                                                         <  7:137037/148‑1 (MQ=38)
                         gtcAGCTGTACACAGGCGAGCAGAAGCGCGTCATCGTTGCCCTGAAAGACCTCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGGCCGCGAaga                                                                                         <  6:52181/146‑1 (MQ=37)
                         gcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGGCCGCGAaga                                                                                         <  3:148837/148‑1 (MQ=38)
                         gcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGGCCGCGAaga                                                                                         <  8:34843/148‑1 (MQ=38)
                          ccAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGGCCGCGAagag                                                                                        <  7:14310/148‑1 (MQ=38)
                          ccAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCAGGCCACTTCGCGGCCTTGCCCGCGAagag                                                                                        >  6:48626/1‑148 (MQ=255)
                           cAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGCCGAgg                                                                                       >  1:53934/1‑148 (MQ=255)
                              cTGTACACCGGCGAGCAGAAGCGCGACAGCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCACGAAGGCTGCCCTCGGGCAGCC‑GTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGTCCg                                                                                    <  6:84100/148‑1 (MQ=255)
                                                    gctaCATCGTTGCCCTGAAAGACCGAAAGTAAGCTGCAGTAGCCTAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCa                                                              <  8:72506/145‑1 (MQ=255)
                                                       aCATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCATGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGtt                                                           >  2:36145/1‑148 (MQ=39)
                                                        cATCGTGGCCCTGAAAGACCGCAAGTAAGCGGCAGAACCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTGCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGttt                                                          <  7:24062/148‑1 (MQ=255)
                                                              tGCCCTGAAAGACCGCAAGTACGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCGGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTcc                                                    <  3:160926/148‑1 (MQ=255)
                                                                 ccTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCATAGCGCCCTTTCCCGCGAATAGGCCGGCACAGCAACCACAAAACCTCAGgttttaaacgcc                                                 >  8:44276/1‑141 (MQ=255)
                                                                   tGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGCAGCGGCCCGCACCGACCCCACAAAACCTCAGTTTTTCTccgcccg                                               >  6:57717/1‑148 (MQ=255)
                                                                   tGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCGCATCGCCGCCTTGCCCGCGCAGAGGCCGGCACAGACCACACAAAACCTCATGTTTTCTccccccg                                               >  5:6637/1‑148 (MQ=255)
                                                                                     cTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTAGCGGCCTTGCCCGCGAAGAGGCCGGCCCAGCCCCCACACAACCTCCGTTTTTCTCCGCCCGCTCCCGCCGCCCTTTCAg                             >  4:125048/1‑148 (MQ=255)
                                                                                      tGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCGCTTAGCGGCCGTGGCCGCGCAGAGGGCGGCCCCGCCCCCCCAAAACCCCAGGTTTTGTCCGCCCGCGCCCGCCCCCCCTTGAgg                            >  1:134639/1‑148 (MQ=255)
                                                                                             gCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCGCTTAGCGGCCTTGGCCGCGAAGAGGCCGGCACCGCCACCCCAAAACCACAGTTTTTGTCCGCCCCCCCACGCAGCCCGTTACGGGTGttta                     >  4:88253/1‑146 (MQ=255)
                                                                                                    gCAAGACGGCTGCCCTCGGCCAGCC‑TTTTTTATTTGTCCTGTGCCGCCCTCTTCGCGCCCTGGCCCGCTAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGgcgc              <  2:125330/148‑1 (MQ=255)
                                                                                                      aaaaagGCTGCCCGAGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGACTTGCCCGCGCAGAGGCCGGCACAGACACCACAAACCCGCAGTTGTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGgcgcct            >  4:148499/4‑146 (MQ=255)
                                                                                                          agGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTACGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATcca         >  5:26234/1‑147 (MQ=255)
                                                                                                              tGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTAGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTATACGCCCGCTCACGCAGCCCTTTAAGGGTATTGAACGGCGCATCcaccac     >  4:67579/1‑147 (MQ=255)
                                                                                                                 ccTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCGCTGCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGaa  <  7:72049/147‑1 (MQ=255)
                                                                                                                                           ||                                                                                                                        
AGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTTATTTTTCTGGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAA  >  NC_002947/4740680‑4740940

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