Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 649,372 A→G F67S (TTC→TCC)  citF ← citrate lyase, citrate‑ACP transferase (alpha) subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913649,3720AG98.0% 176.2 / ‑4.6 51F67S (TTC→TCC) citFcitrate lyase, citrate‑ACP transferase (alpha) subunit
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (22/28);  total (23/28)
Fisher's exact test for biased strand distribution p-value = 4.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.37e-01

TTCTACCAGCGGCGCATGGCAATCACTCAGGGAGCTGGACGCCAGGGTCAGGTTTTTAAAGCCCATCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGAAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATCGCGCGCTTTTTCAGCCTGGTAGGTTTGCTTTGGCGAGTTCTGGAAAGCGGCA  >  NC_000913/649245‑649497
                                                                                                                               |                                                                                                                             
ttCTACCAGCGGCGCATGGCAATCACTCAGGGAGCTGGACGCCAGGGTCAGGTTTTTAAAGCCCATCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTc                                                                                                                    <  2:867419/139‑1 (MQ=255)
 tCTACCAGCGGCGCATGGCAATCACTCAGGGAGCTGGACGCCAGGGTCAGGTTTTTAAAGCCCATCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCa                                                                                                                   >  1:1124952/1‑139 (MQ=255)
 tCTACCAGCGGCGCATGGCAATCACTCAGGGAGCTGGACGCCAGGGTCAGGTTTTTAAAGCCCATCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCa                                                                                                                   <  2:491481/139‑1 (MQ=255)
   tACCAGCGGCGCATGGCAATCACTCAGGGAGCTGGACGCCAGGGTCAGGTTTTTAAAGCCCATCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATg                                                                                                                 <  2:329975/139‑1 (MQ=255)
   tACCAGCGGCGCATGGCAATCACTCAGGGAGCTGGACGCCAGGGTCAGGTTTTTAAAGCCCATCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATg                                                                                                                 <  1:449287/139‑1 (MQ=255)
      cAGCGGCGCATGGCAATCACTCAGGGAGCTGGACGCCAGGGTCAGGTTTTTAAAGCCCATCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCg                                                                                                              >  1:828105/1‑139 (MQ=255)
      cAGCGGCGCATGGCAATCACTCAGGGAGCTGGACGCCAGGGTCAGGTTTTTAAAGCCCATCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCg                                                                                                              <  2:598349/139‑1 (MQ=255)
      cAGCGGCGCATGGCAATCACTCAGGGAGCTGGACGCCAGGGTCAGGTTTTTAAAGCCCATCTTCGCGATGACGTCCATCACCAGATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGTAACCGTCATGCCg                                                                                                              <  1:823060/139‑1 (MQ=255)
        gcggcgCATGGCAATCACTCAGGGAGCTGGACGCCAGGGTCAGGTTTTTAAAGCCCATCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTc                                                                                                            >  2:50220/1‑139 (MQ=255)
              cATGGCAATCACTCAGGGAGCTGGACGCCAGGGTCAGGTTTTTAAAGCCCATCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTaa                                                                                                      >  1:1074888/1‑139 (MQ=255)
                 cgCAATCACTCAGGGAGCTGGACGCCAGGGTCAGGTTTTTAAAGCCCATCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAAcc                                                                                                   >  1:583488/2‑139 (MQ=255)
                      tCACTCAGGGAGCTGGACGCCAGGGTCAGGTTTTTAAAGCCCATCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGc                                                                                              >  1:1069002/1‑139 (MQ=255)
                      tCACTCAGGGAGCTGGACGCCAGGGTCAGGTTTTTAAAGCCCATCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGc                                                                                              <  2:1069002/139‑1 (MQ=255)
                         cTCAGGGAGCTGGACGCCAGGGTCAGGTTTTTAAAGCCCATCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCg                                                                                                              <  2:614904/120‑1 (MQ=255)
                         cTCAGGGAGCTGGACGCCAGGGTCAGGTTTTTAAAGCCCATCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCg                                                                                                              >  1:614904/1‑120 (MQ=255)
                             gggAGCTGGACGCCAGGGTCAGGTTTTTAAAGCCCATCTTCGCGATGACGTCCATCACCATATTGACGGTCCGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAAt                                                                                       >  1:30790/1‑139 (MQ=255)
                                        gCCAGGGTCAGGTTTTTAAAGCCCATCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCa                                                                            >  2:402337/1‑139 (MQ=255)
                                          cAGGGTCAGGTTTTTAAAGCCCATCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCg                                                                                                                      <  2:647592/95‑1 (MQ=255)
                                          cAGGGTCAGGTTTTTAAAGCCCATCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCg                                                                                                                      >  1:647592/1‑95 (MQ=255)
                                                 aGGTTTTTAAAGCCCATCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGgcgc                                                                   <  2:493534/139‑1 (MQ=255)
                                                  ggTTTTTAAAGCCCATCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGca                                                                  <  1:50220/139‑1 (MQ=255)
                                                             cccATCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGAt                                                       >  2:965500/1‑139 (MQ=255)
                                                              ccATCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCg                                                                                                              >  1:1180306/1‑83 (MQ=255)
                                                              ccATCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCg                                                                                                              <  2:1180306/83‑1 (MQ=255)
                                                                aTCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATcgc                                                    >  1:138673/1‑139 (MQ=255)
                                                                   ttCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATcgcgcg                                                 <  1:258378/139‑1 (MQ=255)
                                                                        gATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATCGCGCGCtttt                                            <  2:1124952/139‑1 (MQ=255)
                                                                        gATGACGATCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCCGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATCGCGCGCtttt                                            <  1:740062/139‑1 (MQ=255)
                                                                              gTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATCGCGCGCTTTTTCAGcc                                      <  2:583488/139‑1 (MQ=255)
                                                                               tCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATCGCGCGCTTTTTCAGCCt                                     <  1:402337/139‑1 (MQ=255)
                                                                                     caatatggtgaggtcaggtcaCCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATcg                                                     >  2:629157/12‑117 (MQ=255)
                                                                                     caatatggtgaggtcaggtcaCCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATcg                                                     <  1:629157/106‑1 (MQ=255)
                                                                                      ccATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATCGCGCGCTTTTTCAGCCTGGTAGGt                              >  2:411027/1‑139 (MQ=255)
                                                                                                ggtcaggtcaCCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCg                                                                                                              >  1:1135073/1‑49 (MQ=255)
                                                                                                ggtcaggtcaCCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCg                                                                                                              <  2:1135073/49‑1 (MQ=255)
                                                                                                ggtcaggtcaCCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATCGCGCGCTTTTTCAGCCTGGTAGGTTTGCTTTGGc                    <  2:1074888/139‑1 (MQ=255)
                                                                                                ggtcaggtcaCCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATCGCGCGCTTTTTCAGCCTGGTAGGTTTGCTTTGGc                    <  2:285590/139‑1 (MQ=255)
                                                                                                 gtcaggtcaCCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATCGCGCGCTTTTTCAGCCTGGTAGGTTTGCTTTGGCg                   >  1:78249/1‑139 (MQ=255)
                                                                                                   caggtcaCCGCCACGGAAAGCGTGATGGAAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATCGCGCGCTTTTTCAGCCTGGTAGGTTTGCTTTGGCGAg                 >  1:250859/1‑139 (MQ=255)
                                                                                                    aggtcaCCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCg                                                                                                              >  1:331529/1‑45 (MQ=255)
                                                                                                    aggtcaCCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCg                                                                                                              <  2:331529/45‑1 (MQ=255)
                                                                                                     ggtcaCCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATCGCGCGCTTTTTCAGCCTGGTAGGTTTGCTTTGGCGAGtt               >  1:973161/1‑139 (MQ=255)
                                                                                                         aCCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATCGCGCGCTTTTTCAGCCTGGTAGGTTTGCTTTGGCGAGTTCTgg           <  2:78249/139‑1 (MQ=255)
                                                                                                           cGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATCGCGCGCTTTTTCAGCCTAGTAGGTCTGCTTTGGCGAGTTCTTGaa         <  2:828105/139‑1 (MQ=255)
                                                                                                                cGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATcgcgcg                                                 <  1:780736/94‑1 (MQ=255)
                                                                                                                cGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATcgcgcg                                                 >  2:780736/1‑94 (MQ=255)
                                                                                                                cGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATcgcgcg                                                 <  2:779151/94‑1 (MQ=255)
                                                                                                                cGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATcgcgcg                                                 >  1:779151/1‑94 (MQ=255)
                                                                                                                cGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATCGCGCGCTTTTTCAGCCTGGTAGGTTTGCTTTGGCGAGTTCTGGAAAGCgg    >  2:708340/1‑139 (MQ=255)
                                                                                                                  gAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATCGCGCGCTTTTTCAGCCTGGTAGGTTTGCTTTGGCGAGTTCTGGAAAGCGGCa  <  1:411027/139‑1 (MQ=255)
                                                                                                                  gAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATCGCGCGCTTTTTCAGCCTGGTAGGTTTGCTTTGGCGAGTTCTGGAAAGCGGCa  <  1:965500/139‑1 (MQ=255)
                                                                                                                               |                                                                                                                             
TTCTACCAGCGGCGCATGGCAATCACTCAGGGAGCTGGACGCCAGGGTCAGGTTTTTAAAGCCCATCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGAAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATCGCGCGCTTTTTCAGCCTGGTAGGTTTGCTTTGGCGAGTTCTGGAAAGCGGCA  >  NC_000913/649245‑649497

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