Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_002947 4,740,816 T→G intergenic (+53/+700) 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,8160TG97.0% 98.6 / NA 33intergenic (+53/+700)gltA/yeiWcitrate synthase/putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base G (13/19);  total (13/20)
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

AGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTTATTTTTCTGGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCA  >  NC_002947/4740680‑4740951
                                                                                                                                        |                                                                                                                                       
aGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCg                                                                                                                               <  3:73756/146‑1 (MQ=255)
                       gcgcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCtt                                                                                                                        >  2:84274/1‑130 (MQ=255)
                       gcgcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGaa                                                                                                      <  4:77585/148‑1 (MQ=255)
                        cgcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCT‑GTGCCGGCCTCTTCGCGGGCAAGGCCGCGAaga                                                                                                    <  2:18662/148‑1 (MQ=38)
                        cgcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTc                                                                                                                       <  1:84274/130‑1 (MQ=255)
                         gcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGGCCGCGAaga                                                                                                    <  1:4595/148‑1 (MQ=38)
                         gcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAaga                                                                                                    <  8:50153/148‑1 (MQ=255)
                             gCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCATCGCGGCCTTGCCCGCGAAGCgggc                                                                                                >  5:158765/1‑146 (MQ=255)
                                               gAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCACTTAGCGGGCTTTGCCGCGAAGCGGGCGGCACAGACACAACaaaa                                                                              >  1:67700/1‑148 (MQ=255)
                                                    gcgACATGGTTGCCCTTAAACACCGAATTTACGCTGCAGTAGCCTAACCCAATACGGGTGCCCGGGGGCTG‑CTTTTTTTATTTGTCGTGTGCGGGCCTCTTCGCGGCCTTGCCCGCTATGAGGCCGGCACAGACACCACAAAACCTCa                                                                         <  6:158765/148‑1 (MQ=255)
                                                    gcgACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGGCGCTTTGCGGCCGTGGCGGCGGAGGGGGCGGCCGAGCCCCCCCaaaaaccac                                                                         >  1:163488/1‑143 (MQ=255)
                                                           cGTTGCCCTGAAAGACCGCTAGTAAGCTGCAGTCGCGGAACGCAATAAGGCTGCCCTCGGGCTGCCTTTTTTTTTTGGGCT‑GTGCGGGCCTCTTCGCGGCGTTGCCCGGTAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTtc                                                                  <  4:55377/148‑1 (MQ=255)
                                                                  cTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCAGGCCTCTTCGCGGCCTTGACCGCGAAGAGGCCGGCACAGACACCACAAAACCTCACTTATTCTacgccc                                                           >  6:172653/1‑148 (MQ=255)
                                                                     aaaGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCt                                                        >  4:40246/1‑148 (MQ=255)
                                                                        gACCGCAAGTAAGCTGCAGTCGCCGAACGAATGAAGGCTGCCCTCGGGCAGCCTTTTTTT‑TTTGTCGTGTGCCGGCCTCTTGGCGGCGTTGCCCGCTAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCAc                                                     <  6:178821/148‑1 (MQ=255)
                                                                        gACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGGGAAGAGGCCGGCACAGACACCACAAAACCTCAGGGTTTGTCCGCCCCCTccc                                                     >  5:352/1‑146 (MQ=255)
                                                                                    gctgctTTAGCCGAACGAATGATGGCTGCCCTCGGGCAG‑CTTTTTTTTTTTTTCGTGTGCCGGCCTCTTCGGGGCCTTGCCCGCTAATTGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCa                                         <  6:352/148‑1 (MQ=255)
                                                                                    gCTGGAGTAGCCCAACGAAAGACGGGGGCCGTGGGGCTG‑GTTTTTTTTTTTGTCGTGGGCCGGCCTCTGCGCGGCCTTGCCGGGTATGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCa                                         <  4:102344/148‑1 (MQ=255)
                                                                                    gCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACCCCACaaa                                                                               >  6:130986/1‑110 (MQ=38)
                                                                                     cggCAGTCGCCCACCCCACGACGGCGGCCCTCGGCCAGCC‑TTTTTTATTTGGCCTGGCCCGGCCCCTCCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACaaaa                                                                              <  5:130986/108‑1 (MQ=16)
                                                                                      tGCAGTACCCGAACGCAAGACGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTGGCCCGCTAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAgg                                       <  4:112458/148‑1 (MQ=255)
                                                                                        cAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAggg                                      >  1:14669/1‑147 (MQ=255)
                                                                                         aGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCTAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGt                                     <  2:14669/147‑1 (MQ=255)
                                                                                         aGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAggt                                      <  3:40246/146‑2 (MQ=255)
                                                                                           tAGCCGACCGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTCCCCGCTAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTAtt                                  <  2:54900/148‑1 (MQ=255)
                                                                                                    gCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGgcg                          >  2:116779/1‑147 (MQ=255)
                                                                                                      aagaagGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTAGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCAt                       >  4:45311/1‑148 (MQ=255)
                                                                                                         aagGCTGCCCGAGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCGCACGCCGCCCTTTCAGGGTATTGAACGGCGCATccc                    >  7:146082/1‑147 (MQ=255)
                                                                                                                 ccTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTATTCGAGGCCTTGCCCCCGAAGAGGCCGGCACAGACACCACAAAACCTAAGTTTTTATCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCAACCACCACGAAc            >  2:127506/1‑148 (MQ=255)
                                                                                                                       gCAGCC‑TTTTTAATTTGTCCGGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACttgttg      <  4:72530/148‑1 (MQ=255)
                                                                                                                       gCAGCC‑TTTTTAATTTGCCTGGGGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACttgttg      <  6:179600/148‑1 (MQ=255)
                                                                                                                        gcagcgttttttATTTGTCCTGTGCCGGCCTCTTGGCGGCCTTGCCCGGGAATAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACttgttg      <  3:82556/142‑1 (MQ=255)
                                                                                                                            cTTTTTTTTTGGGTCTGGGGCCGCCGGCTGCGCGGCGTGCCCCGGCACGCGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGcca  <  2:73828/148‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                       
AGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTTATTTTTCTGGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCA  >  NC_002947/4740680‑4740951

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑