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% 100.2 / NA 33intergenic (+53/+700)gltA/yeiWcitrate synthase/putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base G (14/18);  total (14/19)
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

AAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTTATTTTTCTGGTGCCGGCCTCTTCGCGGC‑CTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCA  >  NC_002947/4740689‑4740951
                                                                                                                               |                                                                                                                                        
aaGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCTTTGCCCTGAAAGACCGCAAGTAAGCTGAAGTAGCCGAACGAAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTcgcg                                                                                                                     >  6:42946/1‑148 (MQ=255)
       gCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGC‑CTTGc                                                                                                             <  6:125622/148‑1 (MQ=255)
              gcgcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGC‑CTTGCCCGCGaa                                                                                                      <  1:156101/148‑1 (MQ=255)
              gcgcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGC‑CTTGCCCGCGaa                                                                                                      <  2:45144/148‑1 (MQ=255)
               cgcCAGCGGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGC‑CTTGCCCGCGAag                                                                                                     <  8:74013/148‑1 (MQ=255)
                tccAGCTGTACACCGGCGAGCAGAAGCGCGAAATCGTTTCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGC‑CTTGCCCGCGAaga                                                                                                    <  8:128942/147‑1 (MQ=255)
                gcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGG‑CAAGGCCGCGAaga                                                                                                    <  5:82770/148‑1 (MQ=38)
                gcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGC‑CTTGCCCGCGAaga                                                                                                    <  2:46186/148‑1 (MQ=255)
                gcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGC‑CTTGCCCGCGAaga                                                                                                    <  7:14403/148‑1 (MQ=255)
                gcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCGCTTCGCGGC‑CTTGCCCGCGAaga                                                                                                    <  3:15621/148‑1 (MQ=255)
                   aGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGC‑CTTGCCCGCGAAGAGGc                                                                                                 <  2:62519/148‑1 (MQ=255)
                      tGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGCCCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCCCATAGCGGC‑CCTGCCCCCGACGCGGCCgg                                                                                              >  5:114338/1‑148 (MQ=255)
                         acacCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGC‑CTTGCCCGCGAAGAGGCCGGc                                                                                             >  7:41031/1‑146 (MQ=255)
                                           gcgACATGGTGGCCCTGAACGACCCAAATTACGCTGCAGTAGCCGAACGCAATACGGGTGCCCTCGGGCAGCTTTTTTTTTTTGTCCTGG‑GCCGGCCTCTTCGCGGC‑CTTGCCCGCGAATAGGCCGGCACAGACACCACAAAACCTCa                                                                         <  8:30034/148‑1 (MQ=255)
                                                      gCCCTGAAAGACGGCATGTAAGGTGCAGTAGCCTAACGCAAGAAGGGTGCCGTGGGGCAG‑CTTTTTTTATTTGTCTTGTGCCGGCCTCTGCGCGGC‑TTTGCCCGCTAATAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTcc                                                               <  8:21874/147‑1 (MQ=255)
                                                          tGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTAGCGGC‑CTTGCCCCCGAAGAGGCCGGCACAGACCCCACAAAACCTCAGGTTTTATccgccc                                                           >  8:84809/1‑147 (MQ=255)
                                                                  cGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGC‑CTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCa                                                  >  6:661/1‑148 (MQ=255)
                                                                  cGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGC‑CTTGCCAGCCAAGAGGCCGGCCACGACACCACAAAACCTCAGGTTTTCTACGCACCcgcacgca                                                  >  6:7769/1‑148 (MQ=255)
                                                                        tAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGC‑CTTGCCCGCGCAGCGGCCGGCCCAGCCACCCCAAAACCTCACTTTTTATCCCCCCGCGCACCCAGCCCtt                                            >  1:152334/1‑148 (MQ=255)
                                                                           gcgccAGTCGCCCACCGCACGACGGCGGCCCGCGGGCAGCC‑TTTTTTATTTGTCGTGGGCCGGCCGCTTCGCGGC‑CTGGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCa                                         <  8:119189/144‑1 (MQ=255)
                                                                           gCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCGGGCCTCTTCGCGGC‑CTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCa                                         <  7:120876/148‑1 (MQ=255)
                                                                             tccAGTAGCCGAACGAAAGACGGCTGCCCGCGGGCAGCC‑TTTTTAATTTGCCCGGGGCCGGCCTCTGCGCGGC‑CTGGCCCGCTAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAgg                                       <  5:12601/146‑1 (MQ=255)
                                                                               cAGTAGCCGCACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGC‑CGTGCCCGCGACGCGGCCGGCCCAGACCCCACAAACCCGCAGGTTTTCTCCGCCCGCTCCCGCAGCCCTTTGAgggg                                     >  7:9562/1‑147 (MQ=255)
                                                                               cAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGC‑CATGCCCGCGAAGAGGCCGGCCCAGACACCACAAAACCTCAGTTTTTATACGCCCGCTCACGCAGCCCGTTCAGGGt                                     >  8:159540/1‑148 (MQ=255)
                                                                               cAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCCCGGC‑CTTGCCCGCGAATAGGCCCGCACAGACACCACAAAACCGCAGTTTTTCTACGCCCGCTCACGCAGCCCTTTCAGGGt                                     >  2:116238/1‑148 (MQ=255)
                                                                                 gTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTGCGCGGC‑CTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTAt                                   <  7:24613/148‑1 (MQ=255)
                                                                                           taaaaaagGCTGCCCGAGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTAGCGGC‑CTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTACGCCCCCTCACCCAGCCCGTTCAGGGTATTGAACGgcg                          >  2:161141/6‑147 (MQ=255)
                                                                                            cAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGC‑CTTGCCCGCTAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCa                        <  2:5804/148‑1 (MQ=255)
                                                                                              aaaagGCTGCCCGAGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGC‑CTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCAt                       >  8:152185/3‑147 (MQ=255)
                                                                                                 agGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGC‑CTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATccac                   <  7:53518/148‑1 (MQ=255)
                                                                                                          tCGGGCCGCC‑TTTTTTATTTGTCCTGGGCCGGCCTCTTCGCGGC‑GTGGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACtt          <  5:107265/148‑1 (MQ=255)
                                                                                                                gctgttttttttttttTGCGGGGCCCGCCTCGGCGGGGCGCTGGCCCGCTTCGCGCCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACttgttg      <  6:126973/144‑1 (MQ=255)
                                                                                                                   ccttttttATTTGTCCTGTGCCGGCCTCTTCGCGGC‑CTTGCCCGCGAAGAGGCCGGCACAGCCACCACCAAACCTCAGTTTTTCTCCGCCCGCTCCCGCAGCCCGTTCAGGGGATTGAACGGCGCATCCACCACGAACTTGTTGGcca  >  7:88125/3‑148 (MQ=255)
                                                                                                                   ccttttttATTTGTCCTGTGCCGGCATCTTCGCGGC‑CTTGCCCGCCAAGAGGCCGGCACAGCCCCCCCCACACCCCAGGTTTTGTCCGGCCGCGCACGCAGCCCTGTCAGGGTATTGAACGGCGCATCCACCCAGAACTGGTTGGcca  >  6:179419/3‑148 (MQ=255)
                                                                                                                               |                                                                                                                                        
AAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTTATTTTTCTGGTGCCGGCCTCTTCGCGGC‑CTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCA  >  NC_002947/4740689‑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: ‑