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.92.3% 87.9 / 8.5 28intergenic (+42/+711)gltA/yeiWcitrate synthase/putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base T (0/2);  new base . (7/17);  total (7/21)
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

GAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTTATTTTTCTGGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGG  >  NC_002947/4740664‑4740948
                                                                                                                                             |                                                                                                                                               
gaaagaaaTGCTCTCCAGCCCATACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCtttttta                                                                                                                                           <  1:136096/148‑2 (MQ=255)
   agaaaTGCTCTCCAGCCCGTACAAGATCGGCCGTCAGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGATGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCttttttatttg                                                                                                                                       >  1:23076/1‑144 (MQ=255)
    gaaaTGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCttttttatttgt                                                                                                                                      >  3:108615/1‑143 (MQ=255)
     aaaTGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGCGCAGCCttttttatttgtc                                                                                                                                     >  1:153928/1‑142 (MQ=255)
             tcCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCg                                                                                                                            >  2:175712/1‑149 (MQ=255)
               cAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGcc                                                                                                                             <  2:23076/146‑1 (MQ=255)
               cAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCg                                                                                                                            <  3:91144/147‑1 (MQ=255)
               cAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCg                                                                                                                            <  4:1425/147‑1 (MQ=255)
                  cccGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCctc                                                                                                                       >  5:136019/1‑149 (MQ=255)
                      tACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCctaatcg                                                                                                                   >  6:148144/1‑144 (MQ=255)
                        cAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTcgcg                                                                                                                 <  3:104876/149‑1 (MQ=255)
                        cAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTcgc                                                                                                                  <  6:31721/148‑1 (MQ=255)
                         aaGATCGGCCGGCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCgg                                                                                                                <  5:107484/149‑1 (MQ=255)
                          aGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGc                                                                                                               <  7:112636/149‑1 (MQ=255)
                                      cgcgCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGaa                                                                                                   <  2:153870/149‑1 (MQ=255)
                                      cgcgCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGaa                                                                                                   <  1:175712/149‑1 (MQ=255)
                                            aGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGcc                                                                                             >  4:44599/1‑149 (MQ=255)
                                                                  gcgcTACATCGTTGCCCTGAAAGACCGAAAGTAAGCTGCAGTAGCCGAACGAAAGATGGCTGCCCTCGGGCAG‑CGTTTTTTATTTGTCGTGTGCCGGCCTCTTCGGGGCCTTGCCCGCTAAGAGGCCGGCACAGACACCACAAAACCTc                                                                       <  2:109088/149‑1 (MQ=255)
                                                                     cgACATCGTTGCCCTGAGAGACCGCAAGTAAGCTGCAGTAGCCGAACCCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGt                                                                    <  8:75690/149‑1 (MQ=255)
                                                                         ttggtggCCCTGAAACCCCGCATGTAAGCTGCAGTAGCCGAACGCAAGACGGCTGCCCTCGGGCTGCGTTTTTTT‑TTTGTCCTGTGCCGGCGTCTTCGCGTCCTGGCCGGCTAAGAGGCCGGCACAGACACCACAAAACCTCAGttttt                                                                <  2:107061/143‑1 (MQ=255)
                                                                                   tGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTAGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTccgcccg                                                       >  2:1821/1‑148 (MQ=255)
                                                                                            gCAAGTAAGCTCCAGACGCCGAACCCAAGACGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGCCCGGGGCCGGCCGCTGCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGc                                             <  3:151990/149‑1 (MQ=255)
                                                                                                    gctgctGAAGCCGAACGCAAGAAGGCTGCCAGGGGGCAGCC‑TTTTTTATTTGACCTGTGCCGCCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAg                                     <  5:148144/149‑1 (MQ=255)
                                                                                                     cTGCAGTCGCCGACCGCAAGACGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCGGTGCCGGCCTCTTCGCGGCGTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAgg                                    <  6:136019/149‑1 (MQ=255)
                                                                                                     cTGCAGTAGCCCAACGCAAGACGGCTGCCCTCGGGCTGCC‑TTTTTTATTTGTCCTGTTCCGGCCTCTTCGCGGCCTTGCCCGCTAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAgg                                    <  6:135334/149‑1 (MQ=255)
                                                                                                         aGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCGGCGAAGCCGCCGGCACAGGCCCCACACCACCTCAGGTTTTCGCCGCCCGCTCCCGCAGCCCTTTCCGGGGAt                                >  7:32849/1‑149 (MQ=255)
                                                                                                          gTAGCCGACCGCACGAAGGCTGCCCTGGGGCAGCC‑TTTTTTATTGGTCCGGGGCCGGCCTCTTCGCGGCCTCGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTAtt                               <  6:150291/149‑1 (MQ=255)
                                                                                                                 cccGCAAGAATGCGGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTGCGCGGCCGTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGgc                        <  7:113651/147‑1 (MQ=255)
                                                                                                                   cGCAAGAAGGCTGCCCTCGGGCAG‑CGTTTTTTATTTGGCCTGGGCCGGCCTCTGCGCGGCCTGGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGgcgc                      <  5:169876/149‑1 (MQ=255)
                                                                                                                                     gggCTGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTGGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACttg      <  4:108615/147‑1 (MQ=255)
                                                                                                                                     gggCAG‑CTTTTTTTTTTTGGCCTGTGCCGGCCTTTTCGCGGCCTTGCCCGCTATGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACttgtt    <  4:161010/149‑1 (MQ=255)
                                                                                                                                       gCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTgg  >  8:1617/1‑149 (MQ=255)
                                                                                                                                             |                                                                                                                                               
GAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTTATTTTTCTGGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGG  >  NC_002947/4740664‑4740948

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