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.83.9% 91.8 / 15.5 31intergenic (+42/+711)gltA/yeiWcitrate synthase/putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base T (0/5);  new base . (18/8);  total (18/13)
Fisher's exact test for biased strand distribution p-value = 7.57e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTTATTTTTCTGGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCA  >  NC_002947/4740661‑4740936
                                                                                                                                                |                                                                                                                                    
cTGGAAAGTAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGGCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCtttt                                                                                                                                   <  2:162698/148‑1 (MQ=255)
cTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCtttt                                                                                                                                   <  1:144137/148‑1 (MQ=255)
   gaaagaaaTGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCtttttta                                                                                                                                >  6:70044/1‑147 (MQ=255)
    aaagaaaTGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCtttttta                                                                                                                                <  8:56230/147‑2 (MQ=255)
    aaagaaaTGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCGGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCtttttt                                                                                                                                 <  3:154361/146‑1 (MQ=255)
    aaagaaaTGCTCTCCAGCCCGTACAAGAGCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCGAGGGCAGCcttcttgc                                                                                                                                <  4:167801/147‑3 (MQ=255)
     aagaaaTGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTCGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCttttttatt                                                                                                                              >  4:23738/1‑145 (MQ=255)
        aaaTGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAACTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCttttttatttgt                                                                                                                           >  3:100616/1‑142 (MQ=255)
                   aGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGc                                                                                                               <  5:182793/148‑1 (MQ=255)
                         tACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTc                                                                                                         <  8:133603/148‑1 (MQ=255)
                            aaGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTcgcg                                                                                                      >  8:1697/1‑148 (MQ=255)
                                   gCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGc                                                                                              >  7:155455/1‑148 (MQ=255)
                                           cgcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAag                                                                                      <  8:155455/148‑1 (MQ=255)
                                            gcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAaga                                                                                     <  3:103879/148‑1 (MQ=255)
                                               aGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCCCCGCCTTGCCCGCGAAGAGGc                                                                                  >  7:168772/1‑148 (MQ=255)
                                                    tACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGcc                                                                                 >  2:29429/1‑144 (MQ=255)
                                                     acacCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCg                                                                                <  1:29429/144‑1 (MQ=255)
                                                                 aGAAGCGCTACTTGGTTGCCCTTAATCACCGCATGTATGGTGCAGTCGCCGAACCCATTTCGGCGGCCGTGGGGCAGCTTTTTTTTTTTGTCCTTTGCGGGGCTCTTGGCGGCGTTGCCCGTTAAGAGGCCGGCACAGACACCACaaa                                                                <  4:136371/148‑1 (MQ=255)
                                                                   aaGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGCACAGGCCGGCACCGCCACCACaaaa                                                               >  6:121174/1‑147 (MQ=255)
                                                                       gctaCATGGTTGCCGTGAATGACCGAATGTAAGCTGCTGTACCCGAACGCAACAAGGGGGCCCTCGTTCTGCTTTTTTTATTTGTCCTGTGGCGGCCTTTTCGGGGCCTTGCCCGCTATTAGGCCGGCACAGACACCACAAAACCTCa                                                          <  8:168772/145‑1 (MQ=255)
                                                                       gcgACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTAGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCCCAAAACCTCa                                                          >  6:79559/1‑148 (MQ=255)
                                                                         gACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCGCTTAGCGGCCTTTGCCGCGCACAGGGCGGGCCACGCCCCACAACACCGCAgg                                                        >  5:116978/1‑147 (MQ=255)
                                                                          aCATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCACGAAGAGGCCGGCACAGACACCACAAAACCTCAGttt                                                      <  4:113536/148‑1 (MQ=255)
                                                                          aCATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCACGAAGAGGCCGGCACAGACACCACAAAACCTCAGttt                                                      <  4:111058/148‑1 (MQ=255)
                                                                                     cTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGCAGAGGCCGGCACAGACCCCACAAAACCTCAGTTTTTCTccgccc                                            >  3:70602/1‑148 (MQ=255)
                                                                                     cTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTAGCGGCCGTGGCCGCGCAGCGGCCGGCACAGCCACCACAAAACCTCAGTTTTTGTccgccc                                            >  5:105279/1‑148 (MQ=255)
                                                                                     cTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCGTCGCCGCCGTGCCCGCGCAGAGGGCGGCCCCGACACCACAAAACCCCCCGTTTTGTccgccc                                            >  5:15537/1‑148 (MQ=255)
                                                                                          aGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCAAATAGCGGCCGTTGCCCCGCACAAGCCGGCCCCCACACCACAAAAACACAATTTTTTTtcgcccccgca                                       >  2:100331/1‑143 (MQ=255)
                                                                                           gACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGCAGAGGCCGGGCCAGACACCCCAAACCCTCAGTTTTTGTGCGCCCGCGCAc                                      >  5:109326/1‑148 (MQ=255)
                                                                                              cccAAGTACGCGCCTGTCGCCGAACGAACGAAGGCGGCCCTCGGGCAGCTTTTTTTATTTGCCCTGTGCCGCCCTGTCGGCGGCCTGGCCCGCTAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCa                                   <  2:133302/146‑1 (MQ=255)
                                                                                                      aGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCCCAGACCCCACACAACCTCAGTTTTTCTCCGCCCGCGCCCGCAGCCCtg                             >  7:45233/1‑145 (MQ=255)
                                                                                                          gCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTAGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCGCAGTTTTTATCCGCCCGCTCACGCACCCCGTTCAgg                        >  4:2997/1‑147 (MQ=255)
                                                                                                                 ccGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGACGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGaa                <  7:1697/148‑1 (MQ=255)
                                                                                                                   gAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACg              >  4:146367/1‑148 (MQ=255)
                                                                                                                       gCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGgcgc          >  8:56583/1‑148 (MQ=255)
                                                                                                                       gCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCACTTAGCGGCCTTGCCCGCGAAGAGGCCGGCCCAGACACCACAAAACCTCAGTTTTTATACCCCCCCTCACGCACCCCATTCAGGGTATTTAACGgggc          >  4:162838/1‑148 (MQ=255)
                                                                                                                          atatgggtGCCCTCGGGCAGCTTTTTTTATTTGTCGTGTGCCGGCCTCTTCGCGGCCTTGCCCGCTAAGTGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATc       <  2:48926/141‑1 (MQ=255)
                                                                                                                           aaagGCTGCCCGAGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCGTGCCCGCGCAGCGGGCGGCCCCGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCCCGCAGCCCTTTGAgg                        >  7:73438/2‑130 (MQ=255)
                                                                                                                               gCTGCCCTCGGGCAGCATTTTTTTTTTGTCCTGTGCCGGCCGCTTCGCGGCCTTGCCCGCTAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATccacca  <  8:39462/148‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                    
CTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTTATTTTTCTGGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCA  >  NC_002947/4740661‑4740936

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: