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.94.6% 154.5 / 5.8 38intergenic (+42/+711)gltA/yeiWcitrate synthase/putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base T (0/2);  new base . (20/15);  total (20/18)
Fisher's exact test for biased strand distribution p-value = 2.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.60e-01

ACTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTTATTTTTCTGGTGCCGGCCTCTTCGCGGCC‑TTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTT  >  NC_002947/4740660‑4740943
                                                                                                                                                 |                                                                                                                                           
aCTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCtt                                                                                                                                            <  4:246846/147‑1 (MQ=255)
  tGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCtttt                                                                                                                                          <  5:334021/147‑1 (MQ=255)
     aaagaaaTGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCtttttta                                                                                                                                       >  3:372537/1‑146 (MQ=255)
     aaagaaaTGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCtttttta                                                                                                                                       <  2:304197/147‑2 (MQ=255)
     aaagaaaTGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCtttttta                                                                                                                                       <  4:379555/147‑2 (MQ=255)
        gaaaTGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCttttttattt                                                                                                                                    >  7:92763/1‑143 (MQ=255)
        gaaaTGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCttttttattt                                                                                                                                    >  5:80109/1‑143 (MQ=255)
        gaaaTGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGAAGTAGACGACAGCAAGAAGGCTGCCCTCGGGCAGCCtattttattt                                                                                                                                    >  8:199680/1‑143 (MQ=255)
          aaTGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCttttttatttg                                                                                                                                   >  3:318930/1‑141 (MQ=255)
                    aGCCCGTACAAGATCGGTCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCgg                                                                                                                       <  8:317178/147‑1 (MQ=255)
                    aGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGcc                                                                                                                         <  6:80109/145‑1 (MQ=255)
                    aGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCgg                                                                                                                       <  3:178464/147‑1 (MQ=255)
                       ccGTACAAGATCGGCCGTCCGCGCCAGCTATACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGAGCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCct                                                                                                                    >  8:172372/1‑147 (MQ=255)
                              aGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGGAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTcgcg                                                                                                             <  7:199680/147‑1 (MQ=255)
                                  cGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCC‑t                                                                                                        <  5:1611/147‑1 (MQ=255)
                                    gCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCC‑TTg                                                                                                      >  6:110365/1‑147 (MQ=255)
                                           gcgcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAACTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCC‑TTGCCcgcg                                                                                                <  4:295183/146‑1 (MQ=255)
                                            cgcCAGCTGTACACCGGAGATCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGATCGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCC‑TTGCCCGCGaa                                                                                              <  8:257101/147‑1 (MQ=255)
                                                   tGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCC‑TTGCCCGCGAAGAGGCCg                                                                                       >  2:88314/1‑147 (MQ=255)
                                                   tGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCATCGCGGCC‑TTTCCCGCGAAGCGGCCg                                                                                       >  7:153729/1‑147 (MQ=255)
                                                    gTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCACTTCGCCGCC‑GTGCCCGCGAAGACGCCgg                                                                                      >  5:221586/1‑147 (MQ=255)
                                                      acacCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGAGCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCC‑TTGCCCGCGAAGAGGCCGGca                                                                                    <  7:172372/147‑1 (MQ=255)
                                                           ggCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCC‑TTGCCCGCGAAGAGGCCGGCACAGac                                                                               >  6:323583/1‑147 (MQ=255)
                                                              gagcgtaaGCGCGACATGGTTGCCCTCAACGCCCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCC‑TTGCCCGCTAAGAGGCCGGCACAGAcacc                                                                            <  2:240325/141‑1 (MQ=255)
                                                                gCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCACTTGGCGGCC‑TGTGCCGCGAAGAGGCCGGGACAGccacca                                                                           >  4:10235/1‑146 (MQ=255)
                                                                 cAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCttttttat                                                                                                                                      <  2:214653/88‑3 (MQ=255)
                                                                 cAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCttttttat                                                                                                                                      >  1:214653/1‑86 (MQ=255)
                                                                           aCATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCA‑TTGCCCCCGAAGAGGCCGGCACAGACACCACAAAACCTAAgg                                                               >  7:285379/1‑146 (MQ=255)
                                                                                          aacaCCCCAAGTAAGCTCCAGTAGCCGAACCCAACAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGCCCTGTGCCGGCCTCTTCGCGCCC‑TTGCCCGCTAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTccgcccg                                                  <  8:92763/142‑1 (MQ=255)
                                                                                          aaGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCC‑TTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTccgcccg                                                  <  5:168232/145‑1 (MQ=255)
                                                                                          aaGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCACATAGCCGCCATGGCCCGCGAA‑CGGCCGCCCCACACACCACAAAAACAGAATTTTTATACGCCCCCt                                                >  8:312385/1‑147 (MQ=255)
                                                                                                gCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCC‑TTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCa                                          >  4:138033/1‑147 (MQ=255)
                                                                                                gCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCC‑TTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTATACGCCCCCACACGCa                                          >  4:155874/1‑147 (MQ=255)
                                                                                                  aaGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCC‑TTGCCCGCGAACAGGCCGGCACAGACACAACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAg                                         >  1:45023/1‑146 (MQ=255)
                                                                                                    ttAAGCTGCATTAGCCGACCCCAACACGGCTGCCCTCGGGCAGCCTTTTTTT‑TTTGTCCTGTGCCGGCCTCTTCGCGGCC‑TTGCCCGCTATGAGGCCGGc                                                                                     <  2:286973/99‑1 (MQ=25)
                                                                                                    gTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTAGCGGCC‑TTGCCCGCGAAGAGGCCGGc                                                                                     >  1:286973/1‑100 (MQ=255)
                                                                                                        gcgccAGTTGCCGAACGCAATATGGGGGCCCTGGGGCAGCC‑TTTTTTATTTGTCCTGTGCGGGCCTTTTCGCGGCT‑TTGCCCGCTATGTGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTc                                  <  6:221586/143‑1 (MQ=255)
                                                                                                         ctcaaGTACCCTAACGAAAGAAGGCTGCCCTCGGGCAG‑CATTTTTTATTTGCCCTGGGCCGCCCTCTTCGCGCCC‑TTGCCCGCTACGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCa                                 <  3:21291/143‑1 (MQ=255)
                                                                                                         cTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCC‑TTGCCCCCGAAGAGGCCGGCACAGACCCCACCAAACCTCAGTTTTTCTCCGCCCGCGCACGCAGCCCTTTCa                                 >  5:314369/1‑147 (MQ=255)
                                                                                                          tgaaGTAGCCGAAGCCAAGACGGCTGCCCTCGGCCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCC‑TTGCCCGAGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAg                                <  3:322874/144‑1 (MQ=255)
                                                                                                          gtcAGTACGCGAACGGAAGAAGCCTGCCCTCGTGCTGCCTTTTTTT‑TTTGTCCTGTGCCGGCCGCTTCGCGGCC‑TTGCCCGCGATTAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAg                                <  7:70127/145‑1 (MQ=255)
                                                                                                            cAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCC‑TTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAgg                               >  3:265768/1‑146 (MQ=255)
                                                                                                            cAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGAAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCC‑TTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAggg                              <  3:236011/147‑1 (MQ=255)
                                                                                                                           agaagGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCC‑TTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCAt               >  3:281478/1‑147 (MQ=255)
                                                                                                                              agGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCC‑TTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATcca            >  1:256314/1‑147 (MQ=255)
                                                                                                                                gCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCC‑TTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATccacc          >  2:62234/1‑147 (MQ=255)
                                                                                                                                     ccTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCC‑TTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGaa     >  6:354193/1‑147 (MQ=255)
                                                                                                                                        agggCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCC‑ATGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACtt  >  2:41382/2‑147 (MQ=255)
                                                                                                                                          ggCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCC‑TTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACtt  <  5:110365/145‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                           
ACTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTTATTTTTCTGGTGCCGGCCTCTTCGCGGCC‑TTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTT  >  NC_002947/4740660‑4740943

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