Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_002947 4,740,811 (T)7→6 intergenic (+48/+705) 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,8050T.90.5% 58.1 / 9.4 24intergenic (+42/+711)gltA/yeiWcitrate synthase/putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base T (0/2);  major base . (13/6);  minor base G (0/2);  total (13/11)
Fisher's exact test for biased strand distribution p-value = 1.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTTATTTTTCTGGTGCCGGCCT‑CTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTT  >  NC_002947/4740661‑4740946
                                                                                                                                                |                                                                                                                                              
cTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAAAAGGCTGCCCGAGGGCAGCCttct                                                                                                                                             <  3:92886/148‑3 (MQ=255)
   gaaagaaaTGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCtttttta                                                                                                                                          >  5:93903/1‑147 (MQ=255)
    aaagaaaTGCTCTCCAGCCCGTACAAGATCGGCCGTCCTCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCttttttat                                                                                                                                         <  8:84502/148‑3 (MQ=255)
    aaagaaaTGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCtttttt                                                                                                                                           <  3:57347/146‑1 (MQ=255)
     aagaaaTGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCttttttatt                                                                                                                                        >  6:169512/1‑145 (MQ=255)
       gaaaTGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCttttttatttg                                                                                                                                      >  1:122042/1‑143 (MQ=255)
         aaTGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCttttttatttgtc                                                                                                                                    <  2:51756/148‑8 (MQ=255)
         aaTGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCttttttatttgtc                                                                                                                                    <  5:185465/148‑8 (MQ=255)
                       cGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCct‑ct                                                                                                                    >  6:93630/1‑148 (MQ=255)
                            aaGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCT‑CTTcgcg                                                                                                               >  4:154004/1‑148 (MQ=255)
                              gATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCT‑CTTCGCGGc                                                                                                             <  5:42406/148‑1 (MQ=255)
                                tCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCT‑CTTCGCGGCCt                                                                                                           >  5:12390/1‑148 (MQ=255)
                                    ccGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCT‑CTTCGCGGCCTTGcc                                                                                                       <  6:149831/148‑1 (MQ=255)
                                            gcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCT‑CTTCGCGGGCAAGGCCGCGAaga                                                                                               <  8:88914/148‑1 (MQ=38)
                                            gcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCT‑CTTCGCGGCCTTGCCCGCGAaga                                                                                               <  6:12390/148‑1 (MQ=255)
                                             ccAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCT‑CTTCGCGGCCTTGCCCGCGAagag                                                                                              <  3:6893/148‑1 (MQ=255)
                                                   gTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCT‑CTTCGCGGCCTTGCCCGCGAAGAGGCCGGc                                                                                        >  6:173223/1‑148 (MQ=255)
                                                         cGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCT‑CATAGCCGCCATTCCCCCGAATAGGCCCGCACAGac                                                                                  >  6:123639/1‑148 (MQ=255)
                                                                cAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCA‑CATAGCCGCCGTGCCCGCGAAGACGGCGGCACACAccccacc                                                                            >  2:186798/1‑146 (MQ=255)
                                                                     gcgcGACATCGTTGCCCTGAAAGACCGCAAGTATGCTGCTGTAGCCGACCGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGGGCCGGCCT‑CTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCt                                                                      <  6:93903/148‑1 (MQ=255)
                                                                              cGTTGCCCTTAAAGACCGCAAGTAAGCTGCAGTAGCCCAACGCAAGAAGGCTGCCCTCGGGCAG‑CATTTTTTATTTGTCGTGTGCCGGCCT‑CTTCGCGGCCTTGCCCGCTAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTtc                                                             <  6:191646/148‑1 (MQ=255)
                                                                                  gCCCTGAAAGACCGCAAGTACGCTGCAGTAGCCTAACGCAATACGGCTGCCCGGGGGCAG‑CGTTTTTTATTGGTCCTGTGCCGGCCT‑CTTCGCGGCCTTGCCCGCTAATAGGCCGGCACAGACACCACAAAACCTCAGTTTTtct                                                            <  2:129129/145‑1 (MQ=255)
                                                                                      tGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCT‑CTTAGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACCAAACCTCAGTTTTTCTccgccc                                                      >  4:31828/1‑147 (MQ=255)
                                                                                            aCCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGGCTAATTGCCGGCCTTGCCC‑CGAAGTGGGCGGCACCAACAACACAAAAACCAAGGTTTTCTaagcccccgacac                                               >  4:2614/1‑141 (MQ=255)
                                                                                                 aaGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑ATTTTTATTTGTCCTGTGCCGGACT‑CTTCGCGGACTTTCACGCGAAGAGGCCGGCACAGCCACCACAAAACCTCAGTTTTTATCCGCCCCCTCACGCAGcc                                          >  4:146588/1‑148 (MQ=255)
                                                                                                        cTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCT‑CTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAg                                   >  2:4307/1‑148 (MQ=255)
                                                                                                         ggCAGTACCCCAACGCAAGACGGCGGCCCTCGGGCAG‑CGTTTTTTATTTGCCCGGTGCCGGCCG‑CTGCGCGGCCTTGCCCGCTAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAgg                                  <  5:169512/147‑1 (MQ=255)
                                                                                                                   gcaCGCATGACGGGTGCCGTGGGGCTGCTTTTTTTTTTTGGCCTGG‑GCCGGCGG‑CTGCGCGGCCTGGCCCGGGATGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACg                        <  2:122042/146‑1 (MQ=255)
                                                                                                                       gCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCT‑ATTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGCCACCCCAAAACCTCAGTTTTTCTCCGCCCGCAAACGCAGCCCATTCAGGGTATTGAACGgcgc                    >  6:56293/1‑148 (MQ=38)
                                                                                                                          agaagGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCT‑CTTCGCGGCCTTGCCCGCGCAGAGGCCGGCCCAGACACCACAAAACCTCAGTGTTTCTCCGCCCGCTCACGCCGCCCGTTCAGGGTATTGAACGGCGCATc                 >  8:170902/1‑148 (MQ=255)
                                                                                                                                 tGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGACGGCCT‑ATTAGCGGCCTTGCCCGCGAATAGGACGGCACACACAAAAAAAAAACAAAGTTTTTATCACCCACCCCACGCAACCCATTCAGGGTATTGAACGGCGCATCcaccac           >  8:23773/1‑147 (MQ=255)
                                                                                                                                          ggaagcTTTTTTTATTGGGCCTGGGCCGGCCT‑CTTCGGGGCCTGGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACttgtt  <  3:8519/143‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                              
CTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTTATTTTTCTGGTGCCGGCCT‑CTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTT  >  NC_002947/4740661‑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: ‑