Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 647,863 A→G 8.4% M311T (ATG→ACG)  citF ← citrate lyase, citrate‑ACP transferase (alpha) subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP009048647,8630AG 8.4% 194.0 / 8.5 71M311T (ATG→ACG) citFcitrate lyase, citrate‑ACP transferase (alpha) subunit
Reads supporting (aligned to +/- strand):  ref base A (33/32);  new base G (2/4);  total (35/36)
Fisher's exact test for biased strand distribution p-value = 6.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCAGCATGGCTGTCAAAGCTCTGCACATCCAGCAGTTTGCGGATCAGACCTTTTTCGTGCAGGTCAACCATCGTCGCGGTAATACCGCCAAGGGCGAAGTCGGCGCGAATATCGCGGCTACGCATTTTGTCTTCCAG  >  AP009048/647793‑647930
                                                                      |                                                                   
cgcAGCATGGCTGTCAAAGCTCTGCACATCCAGCAGTTTGCGGATCAGACCTTTTTCGTGCAGGTCAACCa                                                                     <  1:3946051/71‑1 (MQ=255)
cgcAGCATGGCTGTCAAAGCTCTGCACATCCAGCAGTTTGCGGATCAGACCTTTTTCGTGCAGGTCAACCa                                                                     <  1:1652496/71‑1 (MQ=255)
     caTGGCTGTCAAAGCTCTGCACATCCAGCAGTTTGCGGATCAGACCTTTTTCGTGCAGGTCAACCAtcgtc                                                                <  1:3649155/71‑1 (MQ=255)
       tGGCTGTCAAAGCTCTGCACATCCAGCAGTTTGCGGATCAGACCTTTTTCGTGCAGGTCAACCATCGTcgc                                                              <  1:891202/71‑1 (MQ=255)
       tGGCTGTCAAAGCTCTGCACATCCAGCAGTTTGCGGATCAGACCTTTTTCGTGCAGGTCAACCATCGTcgc                                                              >  1:2347317/1‑71 (MQ=255)
         gCTGTCAAAGCTCTGCACATCCAGCAGTTTGCGGATCAGACCTTTTTCGTGCAGGTCAACCAtc                                                                   <  1:744589/64‑1 (MQ=255)
             tCAAAGCTCTGCACATCCAGCAGTTTGCGGATCAGACCTTTTTCGTGCAGGTCAACCATCGTCGCGGTaa                                                         <  1:3403736/70‑1 (MQ=255)
               aaaGCTCTGCACATCCAGCAGTTTGCGGATCAGACCTTTTTCGTGCAGGTCAACCATCGTCGCGGTaa                                                         <  1:1985416/68‑1 (MQ=255)
               aaaGCTCTGCACATCCAGCAGTTTGCGGATCAGACCTTTTTCGTGCAGGTCAACCATCGTCGCGGTAATa                                                       <  1:3954781/70‑1 (MQ=255)
               aaaGCTCTGCACATCCAGCAGTTTGCGGATCAGACCTTTTTCGTGCAGGTCAACCATCGTCGCGGTAATAc                                                      >  1:1636460/1‑71 (MQ=255)
               aaaGCTCTGCACATCCAGCAGTTTGCGGATCAGACCTTTTTCGTGCAGGTCAACCATCGTCGCGGTAATAc                                                      <  1:533607/71‑1 (MQ=255)
                 aGCTCTGCACATCCAGCAGTTTGCGGATCAGACCTTTTTCGTGCAGGTCAACCATCGTcgc                                                              <  1:4298971/61‑1 (MQ=255)
                   ctctGCACATCCAGCAGTTTGCGGATCAGACCTTTTTCGTGCAGGTCAACCATCGTCGCGGTAATACCGcc                                                  <  1:2025050/71‑1 (MQ=255)
                        cacaTCCAGCAGTTTGCGGATCAGACCTTTTTCGTGCAGGTCAACCGTCGTCGCGGTAATACCGCCaa                                                <  1:2238542/68‑1 (MQ=255)
                         acaTCCAGCAGTTTGCGGATCAGACCTTTTTCGTGCAGGTCAACCATCGTCGCGGTAATACCg                                                    >  1:400496/1‑63 (MQ=255)
                         acaTCCAGCAGTTTGCGGATCAGACCTTTTTCGTGCAGGTCAACCATCGTCGCGGTAATACCGCCAAGGGc                                            <  1:633915/71‑1 (MQ=255)
                            tCCAGCAGTTTGCGGATCAGACCTTTTTCGTGCAGGTCAACCGTCGTCGCGGTAATACCGCCAAGGGCg                                           >  1:717512/1‑69 (MQ=255)
                            tCCAGCAGTTTGCGGATCAGACCTTTTTCGTGCAGGTCAACCATCGTCGCGGTAAt                                                        <  1:135720/56‑1 (MQ=255)
                             ccAGCAGTTTGCGGATCAGACCTTTTTCGTGCAGGTCAACCATCGTCGCGGTAATACCGCCAAGGGCGAAg                                        >  1:3401269/1‑71 (MQ=255)
                             ccAGCAGTTTGCGGATCAGACCTTTTTCGTGCAGGTCAACCATCGTCGCGGTAATACCGCCAAGGGCGAAg                                        >  1:797009/1‑71 (MQ=255)
                              cagcagTTTGCGGATCAGACCTTTTTCGTGCAGGTCAACCATCGTCGCGGTAATACCGCCAAGGGCGaa                                         >  1:73125/1‑69 (MQ=255)
                              cagcagTTTGCGGATCAGACCTTTTTCGTGCAGGTCAACCATCGTCGCGGTAATACCGCCAAGGGCGAAGt                                       <  1:1128217/71‑1 (MQ=255)
                               agcagTTTGCGGATCAGACCTTTTTCGTGCAGGTCAACCATCGTCGCGGTAATACCGcc                                                  <  1:1426101/59‑1 (MQ=255)
                                 cagTTTGCGGATCAGACCTTTTTCGTGCAGGTCAACCATCGTCGCGGTAAt                                                        <  1:2843634/51‑1 (MQ=255)
                                  agTTTGCGGATCAGACCTTTTTCGTGCAGGTCAACCATCGTCGCGGTAATACCGCCAAg                                               <  1:3897581/59‑1 (MQ=255)
                                    tttGCGGATCAGACCTTTTTCGTGCAGGTCAACCAtcgtc                                                                <  1:112205/40‑1 (MQ=255)
                                    tttGCGGATCAGACCTTTTTCGTGCAGGTCAACCATCGTCGCgg                                                            <  1:3926353/44‑1 (MQ=255)
                                    tttGCGGATCAGACCTTTTTCGTGCAGGTCAACCATCGTCGCGGTAATACCGCCAAGGGCGAAGTCGgcgc                                 <  1:4370977/71‑1 (MQ=255)
                                    tttGAGGATCAGACCTTTTTCGTGCAGGTCAACCATCGTCGCGGTAATACCGCCAAGGGCGAAGTCGgcgc                                 >  1:1879312/1‑71 (MQ=255)
                                      tGCGGATCAGACCTTTTTCGTGCAGGTCAACCAtc                                                                   >  1:104599/1‑35 (MQ=255)
                                      tGCGGATCAGACCTTTTTCGTGCAGGTCAACCATCGTCGCgg                                                            >  1:776938/1‑42 (MQ=255)
                                      tGCGGATCAGACCTTTTTCGTGCAGGTCAACCATCGTCGCGGTAATACCGCCaa                                                >  1:1051816/1‑54 (MQ=255)
                                         ggATCAGACCTTTTTCGTGCAGGTCAACCATCGTCGCGGTAATACCGCCAAGGGCGAAGTCGGCGCGAat                             >  1:3627927/1‑70 (MQ=255)
                                          gATCAGACCTTTTTCGTGCAGGTCAACCATCGTCGCGGTaa                                                         <  1:174543/41‑1 (MQ=255)
                                           aTCAGACCTTTTTCGTGCAGGTCAACCATCGTCGCGGTAATACCGCCAAGGGCGAAGTCg                                     >  1:2214548/1‑60 (MQ=255)
                                            tCAGACCTTTTTCGTGCAGGTCAACCATCGTCGCGGTAATACCGCCaa                                                >  1:4394330/1‑48 (MQ=255)
                                            tCAGACCTTTTTCGTGCAGGTCAACCATCGTCGCGGTAATACCGCCAAGGGCGaa                                         >  1:91139/1‑55 (MQ=255)
                                            tCAGACCTTTTTCGTGCAGGTCAACCATCGTCGCGGTAATACCGCCAAGGGCGAAGTCGGCGCGAATATcg                         <  1:40304/71‑1 (MQ=255)
                                             cAGACCTTTTTCGTGCAGGTCAACCATCGTCGCGGTAATACCGCCAAGGGCg                                           <  1:2551821/52‑1 (MQ=255)
                                             cAGACCTTTTTCGTGCAGGTCAACCATCGTCGCGGTAATACCGCCAAGGGCGAAGTCGgcgc                                 >  1:3708974/1‑62 (MQ=255)
                                             cAGACCTTTTTCGTGCAGGTCAACCATCGTCGCGGTAATACCGCCAAGGGCGAAGTCGGCGCGAATATcgc                        >  1:3287860/1‑71 (MQ=255)
                                              aGACCTTTTTCGTGCAGGTCAACCGTCGTCGCGGTAATACCGCCAAGGGCGAAGTCGGCGCGAATATcgc                        <  1:2923624/70‑1 (MQ=255)
                                              aGACCTTTTTCGTGCAGGTCAACCATCGTCGCGGt                                                           <  1:3621924/35‑1 (MQ=255)
                                              aGACCTTTTTCGTGCAGGTCAACCATCGTCGCGGTAATACCGCCAAGGGCGAAGTCGGCGCGAATATcgcg                       >  1:2218218/1‑71 (MQ=255)
                                              aGACCTTTTTCGTGCAGGTCAACCATCGTCGCGGTAATACCGCCAAGGGCGAAGTCGGCGCGAATATcgcg                       <  1:2888333/71‑1 (MQ=255)
                                              aGACCTTTTTCGTGCAGGTCAACCATCGTCGCGGTAATACCGCCAAGGGCGAAGTCGGCGCGAATATcgc                        <  1:276234/70‑1 (MQ=255)
                                              aGACCTTTTTCGTGCAGGTCAACCATCGTCGCGGTAATACCGCCAAGGGCGAAGTCGGCGCGAATATcgc                        <  1:4493736/70‑1 (MQ=255)
                                                aCCTTTTTCGTGCAGGTCAACCATCGTCGCGGTAATACCGCCa                                                 <  1:1431666/43‑1 (MQ=255)
                                                aCCTTTTTCGTGCAGGTCAACCATCGTCGCGGTAATACCGCCa                                                 <  1:1281382/43‑1 (MQ=255)
                                                aCCTTTTTCGTGCAGGTCAACCATCGTCGCGGTAATACCGCCAAGGGCGAAGTCGGCGCGAATATCGCGGc                     >  1:4160137/1‑71 (MQ=255)
                                                aCCTTTTTCGTGCAGGTCAACCATCGTCGCGGTAATACCGCCAAGGGCGAAGTCGGCGCGAATATCGCGGc                     >  1:4165694/1‑71 (MQ=255)
                                                 ccTTTTTCGTGCAGGTCAACCGTCGTCGCGGTAATACCGCCAAGGGCGAAGTCGGCGCGAATATCGCGGCt                    <  1:3146755/71‑1 (MQ=255)
                                                 ccTTTTTCGTGCAGGTCAACCATCGTCGCGGTAATACCGCCAAGCGCGAAGTCGGCGCGAATATcgccgc                     >  1:553210/1‑70 (MQ=255)
                                                  cTTTTTCGTGCAGGTCAACCATCGTCGCGGTAATACCGCCAAGGGCGAAGTCGGCGCGAATATCGCGGCTa                   >  1:1172715/1‑71 (MQ=255)
                                                  cTTTTTCGTGCAGGTCAACCATCGTCGCGGTAATACCGCCAAGGGCGAAGTCGGCGCGAATATCGCGGCTa                   >  1:4197639/1‑71 (MQ=255)
                                                  cTTTTTCGTGCAGGTCAACCATCGTCGCGGTAATACCGCCAAGGGCGAAGTCGGCGCGAATATCGCGGCTa                   >  1:4558070/1‑71 (MQ=255)
                                                    ttttCGTGCAGGTCAACCATCGTCGCGGTAATACCGCCAAGGGCGAAGTCGGCGCGAATATcg                         >  1:240188/1‑63 (MQ=255)
                                                    ttttCGTGCAGGTCAACAATCGTCGCGGTAATACCGCCAAGGGCGAAGTCGgcgc                                 >  1:425744/1‑55 (MQ=255)
                                                     tttCGTGCAGGTCAACCATCGTCGCGGTAATACCGCCAAgg                                              <  1:3897486/41‑1 (MQ=255)
                                                     tttCGTGCAGGTCAACCATCGTCGCGGTAATACCGCCAAGGGCGAAGTCGGCGCGAATATCGCGGCTACg                 >  1:2506116/1‑70 (MQ=255)
                                                     tttCGTGCAGGTCAACCATCGTCGCGGTAATACCGCCAAGGGCGAAGTCGGCGCGAATATCGCGGCTACg                 >  1:634324/1‑70 (MQ=255)
                                                     tttCGTGCAGGTCAACCATCGTCGCGGTAATACCGCCAAGGGCGAAGTCGGCGCGAATATCGCGGCTACg                 >  1:1163493/1‑70 (MQ=255)
                                                      ttCGTGCAGGTCAACCATCGTCGCGGTAATACCGCCAAGGGCGAAGTCGGCGCGAATATCGCGGCTACGc                <  1:3719397/70‑1 (MQ=255)
                                                        cGTGCAGGTCAACCATCGTCGCGGTAATACCGCCAAGGGCGAAGTCGGCGCGAATATCGCGGc                     >  1:682620/1‑63 (MQ=255)
                                                        cGTGCAGGTCAACCATCGTCGCGGTAATACCGCCAAGGGCGAAGTCGGCGCGAATATCGCGGCTACGCAt              <  1:1142035/70‑1 (MQ=255)
                                                             aGGTCAACCGTCGTCGCGGTAATACCGCCAAGGGCGAAGTCGGCGCGAATATCGCGGCTACGCATTTTGTc        >  1:4591857/1‑71 (MQ=255)
                                                               gTCAACCGTCGTCGCGGTAATACCGCCAAGGGCGAAGTCGGCGCGAATATCGCGGCTACGCATTTTGTCt       <  1:1873568/70‑1 (MQ=255)
                                                               gTCAACCATCGTCGCGGTAATACCGCCAAGGGCGAAGTCGGCGCGAATATCGCGGCTACGCa               >  1:4187581/1‑62 (MQ=255)
                                                                tCAACCATCGTCGCGGTAATACCGCCAAGGGCGAAGTCGGCGCGAATATCGCGGCTAc                  <  1:1181524/58‑1 (MQ=255)
                                                                  aaCCATCGTCGCGGTAATACCGCCAAGGGCGAAGTCGGCGCGAATATCGCgg                      <  1:2904085/52‑1 (MQ=255)
                                                                  aaCCATCGTCGCGGTAATACCGCCAAGGGCGAAGTCGGCGCGAATATCGCGGCTACGCATTTTGTCTTc     >  1:81215/1‑69 (MQ=255)
                                                                    ccATCGTCGCTGTAATACCGCCAAGGGCGAAGTCGGCGCGAATATCGCGGCTACGCATTTTGTCTTCCAg  >  1:1768004/1‑70 (MQ=255)
                                                                    ccATCGTCGCGGTAATACCGCCAAGGGCGAAGTCGGCGCGAATATCGCGGCTACGCATTTTGTCTTCCAg  >  1:1469857/1‑70 (MQ=255)
                                                                      |                                                                   
CGCAGCATGGCTGTCAAAGCTCTGCACATCCAGCAGTTTGCGGATCAGACCTTTTTCGTGCAGGTCAACCATCGTCGCGGTAATACCGCCAAGGGCGAAGTCGGCGCGAATATCGCGGCTACGCATTTTGTCTTCCAG  >  AP009048/647793‑647930

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