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,8190TC97.0% 96.1 / NA 33intergenic (+56/+697)gltA/yeiWcitrate synthase/putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (19/13);  minor base G (0/1);  total (19/14)
Fisher's exact test for biased strand distribution p-value = 4.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*NC_0029474,740,8200GT100.0% 109.7 / NA 33intergenic (+57/+696)gltA/yeiWcitrate synthase/putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (19/14);  total (19/14)

AGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTTATTTTTCTGGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCA  >  NC_002947/4740680‑4740936
                                                                                                                                           ||                                                                                                                     
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)
                                                                                                                            ccttttttATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCGGCCCTTTCAg                         >  6:163889/3‑110 (MQ=38)
                                                                                                                             cttttttATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCGGCCCTTTCAgg                        <  5:163889/109‑1 (MQ=39)
                                                                                                                                           ||                                                                                                                     
AGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTTATTTTTCTGGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCA  >  NC_002947/4740680‑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: