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,3720AG100.0% 223.4 / NA 62F67S (TTC→TCC) citFcitrate lyase, citrate‑ACP transferase (alpha) subunit
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (25/37);  total (25/37)

TACCAGCGGCGCATGGCAATCACTCAGGGAGCTGGACGCCAGGGTCAGGTTTTTAAAGCCCATCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGAAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATCGCGCGCTTTTTCAGCCTGGTAGGTTTGCTTTGGCGAGTTCTGGAAAGCGGCA  >  NC_000913/649248‑649497
                                                                                                                            |                                                                                                                             
tACCAGCGGCGCATGGCAATCACTCAGGGAGCTGGACGCCAGGGTCAGGTTTTTAAAGCCCATCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATg                                                                                                                 <  2:250991/139‑1 (MQ=255)
tACCAGCGGCGCATGGCAATCACTCAGGGAGCTGGACGCCAGGGTCAGGTTTTTAAAGCCCATCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATg                                                                                                                 <  2:175394/139‑1 (MQ=255)
tACCAGCGGCGCATGGCAATCACTCAGGGAGCTGGACGCCAGGGTCAGGTTTTTAAAGCCCATCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATg                                                                                                                 <  2:139291/139‑1 (MQ=255)
tACCAGCGGCGCATGGCAATCACTCAGGGAGCTGGACGCCAGGGTCAGGTTTTTAAAGCCCATCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATg                                                                                                                 <  2:456655/139‑1 (MQ=255)
  ccAGCGGCGCATGGCAATCACTCAGGGAGCTGGACGCCAGGGTCAGGTTTTTAAAGCCCATCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGcc                                                                                                               <  1:676438/139‑1 (MQ=255)
   cAGCGGCGCATGGCAATCACTCAGGGAGCTGGACGCCAGGGTCAGGTTTTTAAAGCCCATCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCg                                                                                                              <  1:876156/139‑1 (MQ=255)
       ggcgCATGGCAATCACTCAGGGAGCTGGACGCCAGGGTCAGGTTTTTAAAGCCCATCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCt                                                                                                          <  2:1293053/139‑1 (MQ=255)
        gcgcATGGCAATCACTCAGGGAGCTGGACGCCAGGGTCAGGTTTTTAAAGCCCATCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTg                                                                                                         <  1:411886/139‑1 (MQ=255)
                   tCACTCAGGGAGCTGGACGCCAGGGTCAGGTTTTTAAAGCCCATCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGc                                                                                              >  1:1240574/1‑139 (MQ=255)
                   tCACTCAGGGAGCTGGACGCCAGGGTCAGGTTTTTAAAGCCCATCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGc                                                                                              >  1:377909/1‑139 (MQ=255)
                         aGGGAGCTGGACGCCAGGGTCAGGTTTTTAAAGCCCATCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGaa                                                                                        <  1:555562/139‑1 (MQ=255)
                             aGCTGGACGCCAGGGTCAGGTTTTTAAAGCCCATCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGc                                                                                    <  2:481303/139‑1 (MQ=255)
                              gCTGGACGCCAGGGTCAGGTTTTTAAAGCCCATCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCg                                                                                                              <  2:683096/112‑1 (MQ=255)
                              gCTGGACGCCAGGGTCAGGTTTTTAAAGCCCATCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCg                                                                                                              >  1:683096/1‑112 (MQ=255)
                                       cAGGGTCAGGTTTTTAAAGCCCATCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAgg                                                                          <  1:683468/139‑1 (MQ=255)
                                          ggTCAGGTTTTTAAAGCCCATCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTg                                                                       <  1:642070/139‑1 (MQ=255)
                                              aGGTTTTTAAAGCCCATCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGgcgc                                                                   <  2:399525/139‑1 (MQ=255)
                                               ggTTTTTAAAGCCCATCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTaaa                                                                                                     <  2:411393/104‑1 (MQ=255)
                                               ggTTTTTAAAGCCCATCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTaaa                                                                                                     >  1:411393/1‑104 (MQ=255)
                                               ggTTTTTAAAGCCCATCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGca                                                                  <  2:520929/139‑1 (MQ=255)
                                               ggTTTTTAAAGCCCATCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGca                                                                  <  1:1158464/139‑1 (MQ=255)
                                               ggTTTTTAAAGCCCATCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGca                                                                  <  2:129402/139‑1 (MQ=255)
                                               ggTTTTTAAAGCCCATCTTCGCGATGACGTCCAGCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGca                                                                  <  1:134341/139‑1 (MQ=255)
                                                          cccATCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGAt                                                       >  1:25972/1‑139 (MQ=255)
                                                             aTCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATcgc                                                    >  2:1000567/1‑139 (MQ=255)
                                                               cTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATcgcgc                                                  <  1:1000567/139‑1 (MQ=255)
                                                                    cgATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGca                                                                  >  2:612036/1‑118 (MQ=255)
                                                                    cgATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGca                                                                  <  1:612036/118‑1 (MQ=255)
                                                                    cgATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATCGCGCGCttt                                             >  2:306032/1‑139 (MQ=255)
                                                                       tGACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGTTTCTTCCAgg                                                                          >  1:789227/1‑107 (MQ=255)
                                                                       tGACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGTTTCTTCCAgg                                                                          <  2:789227/107‑1 (MQ=255)
                                                                            tCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGc                                                                                    <  2:1188513/92‑1 (MQ=255)
                                                                            tCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGc                                                                                    >  1:1188513/1‑92 (MQ=255)
                                                                              cATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGca                                                                  >  1:187292/1‑108 (MQ=255)
                                                                              cATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGca                                                                  <  2:187292/108‑1 (MQ=255)
                                                                                   ccATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATCGCGCGCTTTTTCAGCCTGGTAGGt                              >  1:392979/1‑139 (MQ=255)
                                                                                      tatTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGca                                                                  <  2:659602/100‑1 (MQ=255)
                                                                                      tatTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGca                                                                  >  1:659602/1‑100 (MQ=255)
                                                                                        ttGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGgcg                                                                    >  2:532429/1‑96 (MQ=255)
                                                                                        ttGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGgcg                                                                    <  1:532429/96‑1 (MQ=255)
                                                                                         tGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATCGCGCGCTTTTTCAGCCTGGTAGGTTTGCtt                        >  2:589334/1‑139 (MQ=255)
                                                                                           aCGGTCATGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATCGCGCGCTTTTTCAGCCTGGTAGGTTTGCTTTg                      >  2:3672/1‑139 (MQ=255)
                                                                                           aCGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAGCCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATCGCGCGCTTTTTCAGCCTGGTAGGTTTGCTTTg                      >  1:1093534/1‑139 (MQ=255)
                                                                                              gtcaggtcaCCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATCGCGCGCTTTTTCAGCCTGGTAGGTTTGCTTTGGc                    >  1:438313/1‑138 (MQ=255)
                                                                                              gtcaggtcaCCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATCGCGCGCTTTTTCAGCCTGGTAGGTTTGCTTTGGc                    <  2:438313/138‑1 (MQ=255)
                                                                                               tcaggtcaCCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATCGCGCGCTTTTTCAGCCTGGTAGGTTTGCTTTGGc                    >  2:109799/1‑137 (MQ=255)
                                                                                               tcaggtcaCCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATCGCGCGCTTTTTCAGCCTGGTAGGTTTGCTTTGGc                    <  1:109799/137‑1 (MQ=255)
                                                                                                caggtcaCCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGc                                                                                    <  2:889344/72‑1 (MQ=255)
                                                                                                caggtcaCCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGc                                                                                    >  1:889344/1‑72 (MQ=255)
                                                                                                caggtcaCCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGca                                                                  <  2:1104293/90‑1 (MQ=255)
                                                                                                caggtcaCCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGca                                                                  >  1:1104293/1‑90 (MQ=255)
                                                                                                 aggtcaCCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCg                                                                                                              <  2:1142569/45‑1 (MQ=255)
                                                                                                 aggtcaCCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCg                                                                                                              >  1:1142569/1‑45 (MQ=255)
                                                                                                        cGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGca                                                                  <  2:506358/82‑1 (MQ=255)
                                                                                                        cGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGca                                                                  >  1:506358/1‑82 (MQ=255)
                                                                                                             cGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTc                                                                                                            <  2:582483/35‑1 (MQ=255)
                                                                                                             cGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTc                                                                                                            >  1:582483/1‑35 (MQ=255)
                                                                                                             cGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATCGCGCGCTTTTTCAGCCTGGTAGGTTTGCTTTGGCGAGTTCTGGAAAGCgg    >  2:1228824/1‑139 (MQ=255)
                                                                                                               gAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAGCCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATCGCGCGCTTTTTCAGCCTGGTAGGTTTGCTTTGGCGAGTTCTGGAAAGCGGCa  <  2:1093534/139‑1 (MQ=255)
                                                                                                               gAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATCGCGCGCTTTTTCAGCCTGGTAGGTTTGCTTTGGCGAGTTCTGGAAAGCGGCa  <  2:25972/139‑1 (MQ=255)
                                                                                                               gAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATCGCGCGCTTTTTCAGCCTGGTAGGTTTGCTTTGGCGAGTTCTGGAAAGCGGCa  <  2:502046/139‑1 (MQ=255)
                                                                                                               gAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATCGCGCGCTTTTTCAGCCTGGTAGGTTTGCTTTGGCGAGTTCTGGAAAGCGGCa  <  1:3672/139‑1 (MQ=255)
                                                                                                                            |                                                                                                                             
TACCAGCGGCGCATGGCAATCACTCAGGGAGCTGGACGCCAGGGTCAGGTTTTTAAAGCCCATCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGAAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATCGCGCGCTTTTTCAGCCTGGTAGGTTTGCTTTGGCGAGTTCTGGAAAGCGGCA  >  NC_000913/649248‑649497

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: