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,6480AG98.0% 169.2 / ‑5.9 50G281G (GGT→GGCcitEcitrate lyase, citryl‑ACP lyase (beta) subunit
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (25/24);  total (26/24)
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

GCTACCCGTTCTTGTCGTTGAGATTGTTCAATTTTCTGCGTCATTTTGATTGCCTTATTCTTCGCGGATGCCGGAAAGTTCTGCACGGGAGAGCACCAGACGGGCGCGATCGATAACCGGACCGTCCACCATCTTGCCGTTCAGGGAAACCACGCCGAGGCCTTCGCGAGCGGCGGCTTCAGCGGCTTCTACGACGCGGCGGGCGTGATCCACTTCTTTCTGGGTCGGTGCGTAGAGGTTGTGCAGCAGATCAATCTGA  >  NC_000913/649528‑649786
                                                                                                                        |                                                                                                                                          
gCTACCCGTTCTTGTCGTTGAGATTGTTCAATTTTCTGCGTCATTTTGATTGCCTTATTCTTCGCGGATGCCGGAAAGTTCTGCACGGGAGAGCACCAGACGGGCGCGATCGATAACCGGGCCGTCCACCATCTTGCCg                                                                                                                          >  1:461089/1‑139 (MQ=255)
 cTACCCGTTCTTGTCGTTGAGATTGTTCAATTTTCTGCGTCATTTTGATTGCCTTATTCTTCGCGGATGCCGGAAAGTTCTGCACGGGAGAGCACCAGACGGGCGCGATCGATAACCGGGCCGTCCACCATCTTGCCGt                                                                                                                         >  2:795326/1‑139 (MQ=255)
       gTTCTTGTCGTTGAGATTGTTCAATTTTCTGCGTCATTTTGATTGCCTTATTCTTCGCGGATGCCGGAAAGTTCTGCACGGGAGAGCACCAGACGGGCGCGATCGATAACCGGGCCGTCCACCATCTTGCCGTTCAggg                                                                                                                   >  2:519907/1‑139 (MQ=255)
                          ttCAATTTTCTGCGTCATTTTGATTGCCTTATTCTTCGCGGATGCCGGAAAGTTCTGCACGGGAGAGCACCAGACGGGCGCGATCGATAACCGGGCCGTCCACCATCTTGCCGTTCAGGGAAACCACGCCGAGGcc                                                                                                   <  1:210494/136‑1 (MQ=255)
                          ttCAATTTTCTGCGTCATTTTGATTGCCTTATTCTTCGCGGATGCCGGAAAGTTCTGCACGGGAGAGCACCAGACGGGCGCGATCGATAACCGGGCCGTCCACCATCTTGCCGTTCAGGGAAACCACGCCGAGGcc                                                                                                   >  2:210494/1‑136 (MQ=255)
                           tCAATTTTCTGCGTCATTTTGATTGCCTTATTCTTCGCGGATGCCGGAAAGTTCTGCACGGGAGAGCACCAGACGGGCGCGATCGATAACCGGGCCGTCCACCATCTTGCCGTTCAGGGAAACCACGCCGAGGCCTTcg                                                                                               >  2:1639395/1‑139 (MQ=255)
                               ttttCTGCGTCATTTTGATTGCCTTATTCTTCGCGGATGCCGGAAAGTTCTGCACGGGAGAGCACCAGACGGGCGCGATCGATAACCGGGCCGTCCACCATCTTGCCGTTCAGGGAAACCACGCCGAGGCCTTCGCGAg                                                                                           >  2:1524163/1‑139 (MQ=255)
                                    tGCGTCATTTTGATTGCCTTATTCTTCGCGGATGCCGGAAAGTTCTGCACGGGAGAGCACCAGACGGGCGCGATCGATAACCGGGCCGTCCACCATCTTGCCGTTCAGGGAAACCACGCCGAGGCCTTCGCGAgcggcg                                                                                      >  1:273911/1‑139 (MQ=255)
                                         cATTTTGATTGCCTTATTCTTCGCGGATGCCGGAAAGTTCTGCACGGGAGAGCACCAGACGGGCGCGATCGATAACCGGGCCGTCAACCATCTTGCCGTTCAGGGAAACCACGCCGAGGCCTTCGCGAGCGGCGGCTTc                                                                                 >  2:1794170/1‑139 (MQ=255)
                                            tttGATTGCCTTATTCTTCGCGGATGCCGGAAAGTTCTGCACGGGAGAGCACCAGACGGGCGCGATCGATAACCGGGCCGTCCACCATCTTGCCGTTCAGGGAAACCACGCCGAGGCCTTCGCGAGCGGCGGCTTCAGc                                                                              <  2:1021664/139‑1 (MQ=255)
                                            tttGATTGCCTTATTCTTCGCGGATGCCGGAAAGTTCTGCACGGGAGAGCACCAGACGGGCGCGATCGATAACCGGGCCGTCCACCATCTTGCCGTTCAGGGAAACCACGCCGAGGCCTTCGCGAGCGGCGGCTTCAGc                                                                              <  1:795326/139‑1 (MQ=255)
                                                     cTTATTCTTCGCGGATGCCGGAAAGTTCTGCACGGGAGAGCACCAGACGGGCGCGATCGATAACCGGGCCGTCCACCATCTTGCCGTTCAGGGAAACCACGCCGAGGCCTTCGCGAGCGGCGGCTTCAGCGGCTTCTac                                                                     >  1:335343/1‑139 (MQ=255)
                                                            ttcGCGGATGCCGGAAAGTTCTGCACGGGAGAGCACCAGACGGGCGCGATCGATAACCGGGCCGTCCACCATCTTGCCGTTCAGGGAAACCACGCCGAGGCCTTCGCGAGCGGCGGCTTCAGCGGCTTCTACGACgcgg                                                              <  2:273911/139‑1 (MQ=255)
                                                              cgcgGATGCCGGAATGTTCTGCACGGGAGAGCACCCGAAGGGCGCGATCGATAACCGGACCGTCCACCATCTTGGCGGTCAGGAGAACCACGCCCAGGCCAACGCAAGCGGCGGCGTCAGCGGCTTCTACGGCgcgacg                                                            >  1:1243587/1‑139 (MQ=255)
                                                              cgcgGATGCCGGAAAGTTCTGCACGGGAGAGCACCAGACGGGCGCGATCGATAACCGGGCCGTcca                                                                                                                                     >  2:1754552/1‑66 (MQ=255)
                                                              cgcgGATGCCGGAAAGTTCTGCACGGGAGAGCACCAGACGGGCGCGATCGATAACCGGGCCGTcca                                                                                                                                     <  1:1754552/66‑1 (MQ=255)
                                                              cgcgGATGCCGGAAAGTTCTGCACGGGAGAGCACCAGACGGGCGCGATCGATAACCGGGCCGTCCACCATCTTGCCGTTCAGGGAAACCACGcc                                                                                                         >  2:574162/1‑94 (MQ=255)
                                                              cgcgGATGCCGGAAAGTTCTGCACGGGAGAGCACCAGACGGGCGCGATCGATAACCGGGCCGTCCACCATCTTGCCGTTCAGGGAAACCACGcc                                                                                                         <  1:574162/94‑1 (MQ=255)
                                                              cgcgGATGCCGGAAAGTTCTGCACGGGAGAGCACCAGACGGGCGCGATCGATAACCGGGCCGTCCACCATCTTGCCGTTCAGGGAAACCACGCCGAg                                                                                                      <  1:423531/97‑1 (MQ=255)
                                                              cgcgGATGCCGGAAAGTTCTGCACGGGAGAGCACCAGACGGGCGCGATCGATAACCGGGCCGTCCACCATCTTGCCGTTCAGGGAAACCACGCCGAg                                                                                                      >  2:423531/1‑97 (MQ=255)
                                                              cgcgGATGCCGGAAAGTTCTGCACGGGAGAGCACCAGACGGGCGCGATCGATAACCGGGCCGTCCACCATCTTGCCGTTCAGGGAAACCACGCCGAGGCCTTCGCGAGCGGCGGCTTCAGCgg                                                                            >  1:276674/1‑123 (MQ=255)
                                                              cgcgGATGCCGGAAAGTTCTGCACGGGAGAGCACCAGACGGGCGCGATCGATAACCGGGCCGTCCACCATCTTGCCGTTCAGGGAAACCACGCCGAGGCCTTCGCGAGCGGCGGCTTCAGCgg                                                                            <  2:276674/123‑1 (MQ=255)
                                                              cgcgGATGCCGGAAAGTTCTGCACGGGAGAGCACCAGACGGGCGCGATCGATAACCGGGCCGTCCACCATCTTGCCGTTCAGGGAAACCACGCCGAGGCCTTCGCGAGCGGCGGCTTCAGCGGCTTCTACGACgcggcg                                                            >  1:64545/1‑139 (MQ=255)
                                                              cgcgGATGCCGGAAAGTTCTGCACGGGAGAGCACCAGACGGGCGCGATCGATAACCGGGCCGTCCACCATCTTGCCGTTCAGGGAAACCACGCCGAGGCCTTCGCGAGCGGCGGCTTCAGCGGCTTCTACGACgcggcg                                                            >  1:1137343/1‑139 (MQ=255)
                                                              cgcgGATGCCGGAAAGTTCTGCACGGGAGAGCACCAGACGGGCGCGATCGATAACCGGGCCGTCCACCATCTTGCCGTTCAGGGAAACCACGCCGAGGCCTTCGCGAGCGGCGGCTTCAGCGGCTTCTACGACgcggcg                                                            >  2:1381791/1‑139 (MQ=255)
                                                               gcgGATGCCGGAAAGTTCTGCACGGGAGAGCACCAGACGGGCGCGATCGATAACCGGGCCGTCCACCATCTTGCCGTTCAGGGAAACCACGCCGAGGCCTTCGCGAGCGGCGGCTTCAGCGGCTTCTACGACGcggcgg                                                           >  1:83893/1‑139 (MQ=255)
                                                                 ggATGCCGGAAAGTTCTGCACGGGAGAGCACCAGACGGGCGCGATCGATAACCGGGCCGTCCACCATCTTGCCGTTCAGGGAAACCACGCCGAGGCCTTCGCGAGCGGCGGCTTCAGCGGCTTCTACGACGCggcgggc                                                         <  1:1369285/139‑1 (MQ=255)
                                                                    tGCCGGAAAGTTCTGCACGGGAGAGCACCAGACGGGCGCGATCGATAACCGGGCCGTCCACCATCTTGCCGTTCAGGGAAACCACGCCGAGGCCTTCGCGAGCGGCGGCTTCAGCGGCTTCTACGACGCGGCGGGCGTg                                                      <  2:83893/139‑1 (MQ=255)
                                                                      ccGGAAAGTTCTGCACGGGAGAGCACCAGACGGGCGCGATCGATAACCGGGCCGTCCACCATCTTGCCGTTCAGGGAAACCACGCCGAGGCCTTCGCGAGCGGCGGCTTCAGCGGCTTCTACGACGCGGCGGGCGTGAt                                                    <  1:1473300/139‑1 (MQ=255)
                                                                           aaGTTCTGCACGGGAGAGCACCAGACGGGCGCGATCGATAACCGGGCCGTCCACCATCTTGCCGTTCAGGGAAACCACGCCGAGGCCTTCGCGAGCGGCGGCTTCAGCGGCTTCTACGACGCGGCGGGCGTGATCCAct                                               <  1:1648438/139‑1 (MQ=255)
                                                                                cTGCACGGGAGAGCACCAGACGGGCGCGATCGATAACCGGGCCGTc                                                                                                                                       <  1:156508/46‑1 (MQ=255)
                                                                                cTGCACGGGAGAGCACCAGACGGGCGCGATCGATAACCGGGCCGTc                                                                                                                                       >  2:156508/1‑46 (MQ=255)
                                                                                 tGCACGGGAGAGCACCAGACGGGCGCGATCGATAACCGGGCCGTCCACCATCTTGCCGTTCAGGGAAACCACGCCGAGGCCTTCGCGAGCGGCGGCTTCAGCGGCTTCTACGACGCGGCGGGCGTGATCCACttctttc                                         <  1:112134/139‑1 (MQ=255)
                                                                                   cACGGGAGAGCACCAGACGGGCGCGATCGATAACCGGGCCGTCCACCATCTTGCCGTTCAGGGAAACCACGCCGAGGCCTTCGCGAGCGGCGGCTTCAGCGGCTTCTACGACGCGGCGGGCGTGATCCACTTCTTTCTg                                       <  1:1639395/139‑1 (MQ=255)
                                                                                      gggAGAGCACCAGACGGGCGCGATCGATAACCGGGCCGTCCACCATCTTGCCGTTCAGGGAAACCACGCCGAGGCCTTCGCGAGCGGCGGCTTCAGCGGCTTCTACGACGCGGCGGGCGTGATCCACTTCTTTCTGGGt                                    <  2:461089/139‑1 (MQ=255)
                                                                                        gagaGCACCAGACGGGCGCGATCGATAACCGGGCCGTCCACCATCTTGCCGTTCAGGGAAACCACGCCGAGGCCTTCGCGAGCGGCGGCTTCAGCGGCTTCTacgacg                                                                 >  1:188535/1‑108 (MQ=255)
                                                                                        gagaGCACCAGACGGGCGCGATCGATAACCGGGCCGTCCACCATCTTGCCGTTCAGGGAAACCACGCCGAGGCCTTCGCGAGCGGCGGCTTCAGCGGCTTCTacgacg                                                                 <  2:188535/108‑1 (MQ=255)
                                                                                           agCACCAGACGGGCGCGATCGATAACCGGGCCGTCCACCATCTTGCCGTTCAGGGAAACCACGCCGAGGCCTTCGCGAGCGGCGGCTTCAGCGGCTTCTACGACGCggcgggcg                                                        <  2:1821581/114‑1 (MQ=255)
                                                                                           agCACCAGACGGGCGCGATCGATAACCGGGCCGTCCACCATCTTGCCGTTCAGGGAAACCACGCCGAGGCCTTCGCGAGCGGCGGCTTCAGCGGCTTCTACGACGCggcgggcg                                                        >  1:1821581/1‑114 (MQ=255)
                                                                                             cACCAGACGGGCGCGATCGATAACCGGGCCGTCCACCATCTTGCCGTTCAGGGAAACCACGCCGAGGCCTTCGCGAGCGGCGGCTTCAGCGGCTTCTACGACGCGGCGGGCGTGATCCACTTCTTTCTGGGTCGGTGCg                             >  2:1495622/1‑139 (MQ=255)
                                                                                               ccAGACGGGCGCGATCGATAACCGGGCCGTCCACCATCTTGCCGTTCAGGGAAACCACGCCGAGGCCTTCGCGAGCGGCGGCTTCAGCGGCTTCTACGACGCGGCGGGCGTGATCCACTTCTTTCTGGGTCGGTGCGTa                           >  2:695653/1‑139 (MQ=255)
                                                                                                  gACGGGCGCGATCGATAACCGGGCCGTCCACCATCTTGCCGt                                                                                                                         <  1:817406/42‑1 (MQ=255)
                                                                                                  gACGGGCGCGATCGATAACCGGGCCGTCCACCATCTTGCCGt                                                                                                                         >  2:817406/1‑42 (MQ=255)
                                                                                                  gACGGGCGCGATCGATAACCGGGCCGTCCACCATCTTGCCGTTCAg                                                                                                                     >  1:1277266/1‑46 (MQ=255)
                                                                                                  gACGGGCGCGATCGATAACCGGGCCGTCCACCATCTTGCCGTTCAg                                                                                                                     <  2:1277266/46‑1 (MQ=255)
                                                                                                  gACGGGCGCGATCGATAACCGGGCCGTCCACCATCTTGCCGTTCAGGGAAACCACGCCGAGGCCTTCGCGAGCGGCGGCTTCAGCGGCTTCTACGACGCGGCGGGCGTGATc                                                   <  1:752805/112‑1 (MQ=255)
                                                                                                  gACGGGCGCGATCGATAACCGGGCCGTCCACCATCTTGCCGTTCAGGGAAACCACGCCGAGGCCTTCGCGAGCGGCGGCTTCAGCGGCTTCTACGACGCGGCGGGCGTGATc                                                   >  2:752805/1‑112 (MQ=255)
                                                                                                          cgatcgatAACCGGGCCGTCCACCATCTTGCCGTTCAGGGAAACCACGCCGAGGCCTTCGCGAGCGGCGGCTTCAGCGGCTTCTACGACTCGGCGGGCGTGATCCACTTCTTTCTGGGTCGGTGCGTAGAGGTTGTgca                <  1:1794170/139‑1 (MQ=255)
                                                                                                           gatcgatAACCGGGCCGTCCACCATCTTGCCGTTCAGGGAAACCACGCCGAGGCCTTCGCGAGCGGCGGCTTCAGCGGCTTCTACGACGCGGCGGGCGTGATCCACTTCTTTCTGGGTCGGTGCGTAGAGGTTGTgcag               <  1:1381791/139‑1 (MQ=255)
                                                                                                               gatAACCGGGCCGTCCACCATCTTGCCGTTCAGGGAAACCACGCCGAGGCCTTCGCGAGCGGCGGCTTCAGCGGCTTCTACGACGCGGCGGGCGTGATCCACTTCTTTCTGGGTCGGTGCGTAGAGGTTGTGCAGCAGa           <  1:695653/139‑1 (MQ=255)
                                                                                                                        aCCGTCCACCATCTTGCCGTTCAGGGAAACCACGCCGAGGCCTTCGCGAGCGGCGGCTTCAGCGGCTTCTACGACGCGGCGGGCGTGATCCACTTCTTTCTGGGTCGGTGCGTAGAGGTTGTGCAGCAGATCAATCTGa  >  2:50903/1‑139 (MQ=255)
                                                                                                                        |                                                                                                                                          
GCTACCCGTTCTTGTCGTTGAGATTGTTCAATTTTCTGCGTCATTTTGATTGCCTTATTCTTCGCGGATGCCGGAAAGTTCTGCACGGGAGAGCACCAGACGGGCGCGATCGATAACCGGACCGTCCACCATCTTGCCGTTCAGGGAAACCACGCCGAGGCCTTCGCGAGCGGCGGCTTCAGCGGCTTCTACGACGCGGCGGGCGTGATCCACTTCTTTCTGGGTCGGTGCGTAGAGGTTGTGCAGCAGATCAATCTGA  >  NC_000913/649528‑649786

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