Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 649,648 A→G G281G (GGT→GGC citE ← citrate lyase, citryl‑ACP lyase (beta) subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913649,6480AG96.0% 83.7 / ‑6.1 25G281G (GGT→GGCcitEcitrate lyase, citryl‑ACP lyase (beta) subunit
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (14/10);  minor base C (0/1);  total (14/11)
Fisher's exact test for biased strand distribution p-value = 4.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCCAGGCCGCTACCCGTTCTTGTCGTTGAGATTGTTCAATTTTCTGCGTCATTTTGATTGCCTTATTCTTCGCGGATGCCGGAAAGTTCTGCACGGGAGAGCACCAGACGGGCGCGATCGATAACCGGACCGTCCACCATCTTGCCGTTCAGGGAAACCACGCCGAGGCCTTCGCGAGCGGCGGCTTCAGCGGCTTCTACGACGCGGCGGGCGTGATCCACTTCTTTCTGGGTCGGTGCGTAGAGGTTGTGCAGCAGATCAAT  >  NC_000913/649520‑649782
                                                                                                                                |                                                                                                                                      
tCCAGGCCGCTACCCGTTCTTGTCGTTGAGATTGTTCAATTTTCTGCGTCATTTTGATTGCCTTATTCTTCGCGGATGCCGGAAAGTTCTGCACGGGAGAGCACCAGACGGGCGCGATCGATAACCGGGCCGTccacca                                                                                                                              <  2:661644/139‑1 (MQ=255)
        gCTACCCGTTCTTGTCGTTGAGATTGTTCAATTTTCTGCGTCATTTTGATTGCCTTATTCTTCGCGGATGCCGGAAAGTTCTGCACGGGAGAGCACCAGACGGGCGCGATCGATAACCGGGCCGTCCACCATCTTGCCg                                                                                                                      >  2:810774/1‑139 (MQ=255)
         cTACCCGTTCTTGTCGTTGAGATTGTTCAATTTTCTGCGTCATTTTGATTGCCTTATTCTTCGCGGATGCCGGAAAGTTCTGCACGGGAGAGCACCAGACGGGCGCGATCGATAACCGGGCCGTCCACCATCTTGCCGt                                                                                                                     >  1:44509/1‑139 (MQ=255)
              cGTTCTTGTCGTTGAGATTGTTCAATTTTCTGCGTCATTTTGATTGCCTTATTCTTCGCGGATGCCGGAAAGTTCTGCACGGGAGAGCACCAGACGGGCGCGATCGATAACCGGGCCGTCCACCATCTTGCCGTTCAgg                                                                                                                <  1:810774/139‑1 (MQ=255)
                                            tGCGTCATTTTGATTGCCTTATTCTTCGCGGATGCCGGAAAGTTCTGCACGGGAGAGCACCAGACGGGCGCGATCGATAACCGGGCCGTCCACCATCTTGCCGTTCAGGGAAACCACGCCGAGGCCTTCGCGAgcggcg                                                                                  >  2:610942/1‑139 (MQ=255)
                                                                    ttcGCGGATTCCGGAAACTTCTGCACGTGAGAGCACCAGACGGGCCCGATCAATACCCGGCCCGTCCACCCTCCTGCCGTTCAGGGCAACCACGCCGAGGCCTTCGCGAGCTGCGGCTTCAGCGGCTTCTACGACgcgg                                                          <  2:615141/139‑1 (MQ=255)
                                                                      cgcgGATGCCGGAAAGTTCTGCACGGGAGAGCACCAGACGGGCGCGATCGATAACCGGGCCGTCCACCATCTTGCCGTTCAGGGAAACCACGCCGAGGCCTTCGCGAGCGGCGGCTTCAGCGGCTTCTACGACgcggcg                                                        >  2:364365/1‑139 (MQ=255)
                                                                           aTGCCGGAAAGTTCTGCACGGGAGAGCACCAGACGGGCGCGATCGATAACCGGGCCGTCCACCATCTTGCCGTTCAGGGAAACCACGCCGAGGCCTTCGCGAGCGGCGGCTTCAGCGGCTTCTACGACGCGGCGGGCGt                                                   >  2:419313/1‑139 (MQ=255)
                                                                                   aaGTTCTGCACGGGAGAGCACCAGACGGGCGCGATCGATAACCGGGCCGTCCACCATCTTGCCGTTCAGGGAAACCACGCCGAGGCCTTCGCGAGCGGCGGCTTCAGCGGCTTCTACGACGCGGCGGGCGTGATCCAct                                           <  1:364365/139‑1 (MQ=255)
                                                                                     gTTCTGCACGGGAGAGCACCAGACGGGCGCGATCGATAACCGGGCCGTCCACCATCTTGCCGTTCAGGGAAACCACGCCGAGGCCTTCGCGAGCGGCGGCTTCAGCGGCTTCTACGACGCGGCGGGCGTGATCCAct                                           >  1:370317/1‑137 (MQ=255)
                                                                                     gTTCTGCACGGGAGAGCACCAGACGGGCGCGATCGATAACCGGGCCGTCCACCATCTTGCCGTTCAGGGAAACCACGCCGAGGCCTTCGCGAGCGGCGGCTTCAGCGGCTTCTACGACGCGGCGGGCGTGATCCAct                                           <  2:370317/137‑1 (MQ=255)
                                                                                                     cACCAGACGGGCGCGATCGATAACCGGGCCGTCCACCATCTTGCCGTTCAGGGAAACCACGCCGAGGCCTTCGCGAgcg                                                                                     <  2:670909/79‑1 (MQ=255)
                                                                                                     cACCAGACGGGCGCGATCGATAACCGGGCCGTCCACCATCTTGCCGTTCAGGGAAACCACGCCGAGGCCTTCGCGAgcg                                                                                     >  1:670909/1‑79 (MQ=255)
                                                                                                          gACGGGCGCGATCGATAACCGGGCCGTCCACCATCTTGCCGTTCAGGGAAACCACGCCGAGGCCTTCGCGAGCGGCGGCTTCa                                                                            >  2:361827/1‑83 (MQ=255)
                                                                                                          gACGGGCGCGATCGATAACCGGGCCGTCCACCATCTTGCCGTTCAGGGAAACCACGCCGAGGCCTTCGCGAGCGGCGGCTTCa                                                                            <  1:361827/83‑1 (MQ=255)
                                                                                                          gACGGGCGCGATCGATAACCGGGCCGTCCACCATCTTGCCGTTCAGGGAAACCACGCCGAGGCCTTCGCGAGCGGCGGCTTCAGCGGCTTCTACGACGCGGCGGGCGTGATCCACTTCTTTCTGGGTc                               >  1:604863/1‑128 (MQ=255)
                                                                                                          gACGGGCGCGATCGATAACCGGGCCGTCCACCATCTTGCCGTTCAGGGAAACCACGCCGAGGCCTTCGCGAGCGGCGGCTTCAGCGGCTTCTACGACGCGGCGGGCGTGATCCACTTCTTTCTGGGTc                               <  2:604863/128‑1 (MQ=255)
                                                                                                          gACGGGCGCGATCGATAACCGGGCCGTCCACCATCTTGCCGTTCAGGGAAACCACGCCGAGGCCTTCGCGAGCGGCGGCTTCAGCGGCTTCTACGACGCGGCGGGCGTGATCCACTTCTTTCTGGGTCGGTGCGTagag                    >  1:628004/1‑139 (MQ=255)
                                                                                                          gACGGGCGCGATCGATAACCGGGCCGTCCACCATCTTGCCGTTCAGGGAAACCACGCCGAGGCCTTCGCGAGCGGCGGCTTCAGCGGCTTCTACGACGCGGCGGGCGTGATCCACTTCTTTCTGGGTCGGTGCGTagag                    >  1:150284/1‑139 (MQ=255)
                                                                                                             gggCGCGATCGATAACCGGGCCGTCCACCATCTTGCCGTTCAGGGAAACCACGCCg                                                                                                    >  2:249342/1‑56 (MQ=255)
                                                                                                             gggCGCGATCGATAACCGGGCCGTCCACCATCTTGCCGTTCAGGGAAACCACGCCg                                                                                                    <  1:249342/56‑1 (MQ=255)
                                                                                                              ggCGCGATCGATAACCGGGCCGTCCACCATCTTGCCGTTCAGGGAAACCACGCCGAGGCCTTCGCGAGCGGCGGCTTCAGCGGCTTCTACGACGCGGCGGGCGTGATCCACTTCTTTCTGGGTCGGTGCGTAGAGGTtg                >  1:790640/1‑139 (MQ=255)
                                                                                                                   gatcgatAACCGGGCCGTCCACCATCTTGCCGTTCAGGGAAACCACGCCGAGGCCTTCGCGAGCGGCGGCTTCAGCGGCTTCTACGACGCGGCGGGCGTGATCCACTTCTTTCTGGGTCGGTGCGTAGAGGTTGTgcag           <  2:790640/139‑1 (MQ=255)
                                                                                                                     tcgatAACCGGGCCGTCCACCATCTTGCCGTTCAGGGAAACCACGCCGAGGCCTTCGCGAGCGGCGGCTTCAGCGGCTTCTACGACGCGGCGGGCGTGATCCACTTCTTTCTGGGTCGGTGCGTAGAGGTTGTgcagca         >  1:558752/1‑139 (MQ=255)
                                                                                                                            ccGGGCCGTCCACCATCTTGCCGTTCAGGGAAACCACGCCGAGGCCTTCGCGAGCGGCGGCTTCAGCGGCTTCTACGACGCGGCGGGCGTGATCCACTTCTTTCTGGGTCGGTGCGTAGAGGTTGTGCAGCAGATCAAt  <  2:628004/139‑1 (MQ=255)
                                                                                                                                |                                                                                                                                      
TCCAGGCCGCTACCCGTTCTTGTCGTTGAGATTGTTCAATTTTCTGCGTCATTTTGATTGCCTTATTCTTCGCGGATGCCGGAAAGTTCTGCACGGGAGAGCACCAGACGGGCGCGATCGATAACCGGACCGTCCACCATCTTGCCGTTCAGGGAAACCACGCCGAGGCCTTCGCGAGCGGCGGCTTCAGCGGCTTCTACGACGCGGCGGGCGTGATCCACTTCTTTCTGGGTCGGTGCGTAGAGGTTGTGCAGCAGATCAAT  >  NC_000913/649520‑649782

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: