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.0% 68.5 / ‑4.5 26intergenic (+56/+697)gltA/yeiWcitrate synthase/putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (13/11);  total (13/13)
Fisher's exact test for biased strand distribution p-value = 4.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.20e-01
*NC_0029474,740,8200GT92.3% 72.2 / ‑5.5 26intergenic (+57/+696)gltA/yeiWcitrate synthase/putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base G (0/2);  new base T (13/11);  total (13/13)
Fisher's exact test for biased strand distribution p-value = 4.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTTATTTTTCTGGTGCCGGCCT‑CTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACC  >  NC_002947/4740684‑4740953
                                                                                                                                       ||                                                                                                                                      
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)
                                                                                                                        ccttttttATTTGTCCTGTGCCGGCCT‑CTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGcca    <  5:93630/146‑1 (MQ=255)
                                                                                                                          ttttttATTTGTCCTGTGCCGGCCT‑CTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGccacc  <  3:45694/148‑1 (MQ=255)
                                                                                                                                       ||                                                                                                                                      
CGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTTATTTTTCTGGTGCCGGCCT‑CTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACC  >  NC_002947/4740684‑4740953

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