Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 649,964 +G coding (527/909 nt) 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,9601.G97.8% 162.0 / ‑5.0 47P177P (CCG→CCCcitEcitrate lyase, citryl‑ACP lyase (beta) subunit
Reads supporting (aligned to +/- strand):  ref base . (0/1);  major base G (23/22);  minor base T (1/0);  total (24/23)
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

GTTAGCGTCGGAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGAAGCGTGAGCGATTTCCACTGCGCG  >  NC_000913/649864‑650055
                                                                                                 |                                                                                               
gttAGCGTCGGAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCggg                                                                                              <  1:191402/100‑1 (MQ=255)
    gCGTCGGAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCgggggag                                                                                          >  1:302837/1‑97 (MQ=255)
    gCGTCGGAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCgggggag                                                                                          <  1:1485987/100‑4 (MQ=255)
             tAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGTGTTCTGTa                                                                                  <  2:694264/100‑1 (MQ=255)
             tAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTa                                                                                  <  1:1041825/100‑1 (MQ=255)
             tAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTa                                                                                  <  1:1581193/100‑1 (MQ=255)
             tAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGGTCAGTTCCTTCCGGGGGAGCGTTCTGTa                                                                                  <  1:630022/100‑1 (MQ=255)
               gACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACg                                                                                >  1:641235/1‑100 (MQ=255)
               gACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACg                                                                                >  1:855890/1‑100 (MQ=255)
               gACGGTATCGAACGCCTGAATACCCGAAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACg                                                                                >  1:961098/1‑100 (MQ=255)
                aCGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGc                                                                               >  1:834862/1‑100 (MQ=255)
                 cGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCa                                                                              >  1:193049/1‑100 (MQ=255)
                   ttATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAgg                                                                            <  1:1154052/99‑1 (MQ=255)
                    tATCGAACGCCTGAATACCCGCAGCGCGCGCGGCCTGCAAAATGGAACAGCGTTCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAgga                                                                           >  2:1154047/1‑99 (MQ=255)
                     aTCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGtt                                                                          >  2:363180/1‑100 (MQ=255)
                       cGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTgc                                                                        >  2:297722/1‑100 (MQ=255)
                       cGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTgc                                                                        >  1:1533432/1‑100 (MQ=255)
                       cGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTgc                                                                        >  2:518102/1‑100 (MQ=255)
                          aCGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGca                                                                     <  2:746824/100‑1 (MQ=255)
                          aCGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGca                                                                     <  1:45607/100‑1 (MQ=255)
                           cGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGcac                                                                    <  1:1533618/100‑1 (MQ=255)
                           cGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGcac                                                                    <  1:386610/100‑1 (MQ=255)
                              cTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATa                                                                 >  1:1185052/1‑100 (MQ=255)
                              cTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATa                                                                 <  1:1232173/100‑1 (MQ=255)
                                  aTACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGtct                                                             <  2:672185/100‑1 (MQ=255)
                                     cccGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTCCGAACAGCACTTCAGTTCCTTCCTGCGCCGCGTTCTGTACGCAGGTTGCACACATAGtcttct                                                          >  1:1245321/1‑100 (MQ=255)
                                       cGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTgn                                                        <  2:993852/100‑2 (MQ=255)
                                                 gcgGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCg                                              <  1:1465669/100‑1 (MQ=255)
                                                  ctgCCTGGAAAATTGAGCAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGc                                               <  1:683602/95‑1 (MQ=255)
                                                  atctcTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGa                                             <  1:892130/96‑1 (MQ=255)
                                                      ntGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGaagat                                         >  2:892125/2‑96 (MQ=255)
                                                      cTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATAcc                                         <  1:429608/100‑1 (MQ=255)
                                                      cTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATAcc                                         <  2:84988/100‑1 (MQ=255)
                                                       tGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCg                                        >  1:809026/1‑100 (MQ=255)
                                                       tGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTGCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCg                                        <  2:593163/100‑1 (MQ=255)
                                                        gCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGa                                       >  1:500019/1‑100 (MQ=255)
                                                         cAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGAt                                      >  1:577398/1‑100 (MQ=255)
                                                          aaaaTGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATc                                     <  2:1305309/100‑1 (MQ=255)
                                                                   cAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCg                            >  1:1223420/1‑100 (MQ=255)
                                                                    aGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCgg                           <  2:809021/100‑1 (MQ=255)
                                                                     gCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGa                          >  1:278828/1‑100 (MQ=255)
                                                                       gTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGAAg                        >  1:1322665/1‑100 (MQ=255)
                                                                       gTGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGAAg                        >  1:807916/1‑100 (MQ=255)
                                                                        tGCGAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGAAGc                       <  1:423873/100‑1 (MQ=255)
                                                                           gAACAGCAGTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCTAGGGCGATACCGATCAAACGCTCGGAAGCGTg                    >  1:294135/1‑100 (MQ=255)
                                                                               agcagTTCAGTTCCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGAAGCGTGAGCg                <  2:502118/100‑1 (MQ=255)
                                                                                   gatcagttcCTTCAGGGGGAGCGTTCTGTACACAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGAAGCGTGAGCGAttt            >  1:1065270/3‑100 (MQ=255)
                                                                                     tcagttcCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACACAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGAAGCGTGAGCGATTTcc          <  2:855885/100‑1 (MQ=255)
                                                                                      cagttcCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGAAGCGTGAGCGATTTCCa         >  1:619894/1‑100 (MQ=255)
                                                                                      cagttcCTTCCGGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGAAGCGTGAGCGATTTCCa         >  2:1087103/1‑100 (MQ=255)
                                                                                         ttccttccgGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGAAGCGTGAGCGATTTCCACTg      <  1:251548/100‑1 (MQ=255)
                                                                                            cttccgGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGAAGCGTGAGCGATTTCCACTgcgc   <  1:93430/100‑1 (MQ=255)
                                                                                            cttccgGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGAAGCGTGAGCGATTTCCACTgcgc   >  2:855058/1‑100 (MQ=255)
                                                                                            cttccGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGAAGCGTGAGCGATTTCCACTGcgcg  >  1:1403262/1‑100 (MQ=255)
                                                                                                 |                                                                                               
GTTAGCGTCGGAATAGACGGTATCGAACGCCTGAATACCCGCAGAGCGCGCGGCCTGCAAAATGGAACAGCGTGCGAACAGCAGTTCAGTTCCTTCCGGGGAGCGTTCTGTACGCAGGTTGCGCACATAGTCTTCTGCACCGAGGGCGATACCGATCAAACGCTCGGAAGCGTGAGCGATTTCCACTGCGCG  >  NC_000913/649864‑650055

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: