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.90.5% 144.5 / 11.8 42intergenic (+42/+711)gltA/yeiWcitrate synthase/putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base T (0/4);  new base . (20/18);  total (20/22)
Fisher's exact test for biased strand distribution p-value = 1.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

AAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTTATTTTTCTGGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTG  >  NC_002947/4740669‑4740947
                                                                                                                                        |                                                                                                                                               
aaaTGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTTGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCttttttatttgtc                                                                                                                                     <  1:307819/149‑8 (MQ=255)
   tGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCtt                                                                                                                                                >  8:63718/1‑135 (MQ=255)
   tGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCtt                                                                                                                                                <  7:63718/135‑1 (MQ=255)
          cAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGc                                                                                                                          >  5:263177/1‑149 (MQ=255)
          cAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCACAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTCTTTATTTGTCCTGTGCCGGc                                                                                                                          <  4:111140/149‑1 (MQ=255)
           aGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGc                                                                                                                          >  4:411417/1‑148 (MQ=255)
               cGTACAAGATCGGCCGTCCGCACCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCttttt                                                                                                                                             <  7:226696/126‑1 (MQ=255)
               cGTACAAGATCGGCCGTCCGCACCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCttttt                                                                                                                                             >  8:226696/1‑126 (MQ=255)
                gTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCttttt                                                                                                                                             >  7:275347/1‑125 (MQ=255)
                gTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCttttt                                                                                                                                             <  8:275347/125‑1 (MQ=255)
                      gATCGGCCGTCCTCGCCAGCTGTACTCCGGCGAGCAGAAGCGCGACATCGTTGCGCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGcc                                                                                                             <  3:38660/149‑1 (MQ=255)
                      gATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCttttttat                                                                                                                                          <  7:124958/122‑3 (MQ=255)
                      gATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCttttttat                                                                                                                                          >  8:124958/1‑120 (MQ=255)
                       aTCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGACAGACCGCAAGTAAGCTGCAGTAGCCGACCGCAAGAAGGCTGCCCTCGGGCAGCCGTTTTTATTTGTCCTGGGCCGGCCTCTTCGCGGCCt                                                                                                            >  3:342758/1‑149 (MQ=255)
                       aTCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCt                                                                                                            >  6:430180/1‑149 (MQ=255)
                        tCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCtt                                                                                                           <  5:124189/149‑1 (MQ=255)
                          ggCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCACTTCGCGGCCtgtc                                                                                                         >  3:130348/1‑146 (MQ=255)
                                   cgcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAaga                                                                                                <  2:279996/149‑1 (MQ=255)
                                    gcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAagag                                                                                               >  4:246569/1‑149 (MQ=255)
                                      cAGCTGTACACCGGCGAGCAGTAGAGCGACATCGTTGCCCTGAAAGACCGAAAGTAAGCTGCAGTAGCCGAACGCAATAAGTCTGCCCTCGGGCAGCCTTTTTTTTTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCGCGCGAAGAGGc                                                                                             <  4:342758/149‑1 (MQ=255)
                                               acCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGa                                                                                    <  7:125568/149‑1 (MQ=255)
                                                          gacgcgcGACCTCGTGGCCCTGAAACACCTCACTTACGCGGCATCAGCCTACCCCAACAGGGCTCCCCTCGGCATGCTTTTTTTATTTGCCCTGGGCCGGCCTCTTCGCGGCCTTGCCCGCTATGTGGCCGGCACAGACACCACAAAAc                                                                         <  4:130348/146‑1 (MQ=255)
                                                                 gACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCACATACGCGCCATTCCCCGCAAAGAGCCCGGACACGCCCCACAAAAccaaaatt                                                                  >  2:361180/1‑143 (MQ=255)
                                                                      cGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGGCCGCGAaga                                                                                                >  7:396546/1‑114 (MQ=37)
                                                                      cGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGGCCGCGAaga                                                                                                <  8:396546/114‑1 (MQ=37)
                                                                      cGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAagag                                                                                               <  4:412773/115‑1 (MQ=255)
                                                                      cGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAagag                                                                                               >  3:412773/1‑115 (MQ=255)
                                                                      cGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTAGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCATCGCGGCCTTGCCCGCGAAGACGCCGGCACAGACACCACAAAACCGCAGGTTTtc                                                              >  5:402927/1‑148 (MQ=255)
                                                                       gTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTACCGGCCTTTCCCGCGCAGAGGCCGGCACCGACCCCCCAAAACCGCAGGTTTtatc                                                            >  6:364524/1‑149 (MQ=255)
                                                                        ttGCCCTCAAAGACCGCAAGTAAGCTGCAGTCGCCCAACCCAACAAGGCTGCCCTCGGCCAGCCTTTTTTATTTGTCCTGGGCCGGCCTCTTCGCGGCCTTGCCCGCTAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTcc                                                           <  4:374778/149‑1 (MQ=255)
                                                                             cTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTACTccgcc                                                        >  7:373840/1‑147 (MQ=255)
                                                                             cTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTACTccgcc                                                        <  8:373840/147‑1 (MQ=255)
                                                                             cTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCACTTAGCGGCCTTGCCCGCGCAGAGGCCGGCACAGCCACCACAAAACCACACTTTTTATacgcccg                                                      >  1:180492/1‑149 (MQ=255)
                                                                              tGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGGCCGCGAagag                                                                                               <  1:332831/107‑1 (MQ=25)
                                                                                       gCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACCCCACAAAACCTCAGTTTTTCTACGCCCGCGAACGCAGc                                            >  2:297777/1‑149 (MQ=255)
                                                                                        cAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAaga                                                                                                >  1:383707/1‑96 (MQ=38)
                                                                                        cAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAaga                                                                                                <  2:383707/96‑1 (MQ=38)
                                                                                         gaGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAagag                                                                                               >  7:14847/2‑96 (MQ=38)
                                                                                         gaGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAagag                                                                                               <  8:14847/95‑1 (MQ=38)
                                                                                              aGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTc                                      <  6:337899/148‑1 (MQ=255)
                                                                                                cTGCCGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGa                                                                                    <  5:433174/100‑1 (MQ=37)
                                                                                                cTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGa                                                                                    >  6:433174/1‑100 (MQ=38)
                                                                                                 tGCAGTAGCCGAACGCAAGAAGGCTGCCCTCTGGAAGCCTTTTTTATTTGTCCTGTTCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCATTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAggg                                  <  5:430180/149‑1 (MQ=255)
                                                                                                        gCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGa                            <  3:246569/148‑1 (MQ=255)
                                                                                                                 acgaagGCTCCCCTCGGGCAGCCTTTTTTTTTTGTCCTGTGCCGCCCTCTTCGCGGCCTTGCCCACGAAGAGGCCGGCACAGACACCACAAAACCTCAg                                                                    <  2:183369/97‑1 (MQ=37)
                                                                                                                 aagaagGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCGTGGCCCCGAAGAGGCCGGGACAGACACCACAAAACCTCAg                                                                    >  1:183369/1‑99 (MQ=25)
                                                                                                                        cTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACg           <  3:411417/149‑1 (MQ=255)
                                                                                                                              tCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGcc                                           >  6:332313/1‑111 (MQ=38)
                                                                                                                              tCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGcc                                           <  5:332313/111‑1 (MQ=38)
                                                                                                                              tCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGAAAAAACAAAACCTCAGTTTTTATACGCCCCCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCCTCCACAACGAACttg     >  1:344823/1‑149 (MQ=255)
                                                                                                                                 ggCAGCCTTTTTTTTTTGTCCGGTCCCGGCCTCTGCGCGGCCTTGCCGGCGAAGTGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACttgttg  <  6:402927/149‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                               
AAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTTATTTTTCTGGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTG  >  NC_002947/4740669‑4740947

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: