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% 121.5 / NA 34F67S (TTC→TCC) citFcitrate lyase, citrate‑ACP transferase (alpha) subunit
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (16/18);  total (16/18)

GCGAATGTGTTCTACCAGCGGCGCATGGCAATCACTCAGGGAGCTGGACGCCAGGGTCAGGTTTTTAAAGCCCATCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGAAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATCGCGCGCTTTTTCAGCCTGGTAGGTTTGCTTTGGCGAGTTCTGGAAAGCGGCAAGATC  >  NC_000913/649236‑649502
                                                                                                                                        |                                                                                                                                  
gCGAATGTGTTCTACCAGCGGCGCATGGCAATCACTCAGGGAGCTGGACGCCAGGGTCAGGTTTTTAAAGCACATCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATgggag                                                                                                                                  >  1:457457/1‑139 (MQ=255)
            tACCAGCGGCGCATGGCAATCACTCAGGGAGCTGGACGCCAGGGTCAGGTTTTTAAAGCCCATCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATg                                                                                                                      <  2:457457/139‑1 (MQ=255)
               cAGCGGGGCATGGCAATCACTGAGGGAGTTTGACGCCATGGTCAGGTTTTTAACGCCCATCTTCGCGATGATGTCCAACACCTTATTGGCGGTGAGGTCACCGTCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCg                                                                                                                   <  1:693103/139‑1 (MQ=255)
               cAGCGGCGCATGGCAATCACTCAGGGAGCTGGACGCCAGGGTCAGGTTTTTAAAGCCCATCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCg                                                                                                                   <  2:753646/139‑1 (MQ=255)
                                 aCTCAGGGAGCTGGACGCCAGGGTCAGGTTTTTAAAGCCCATCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGcga                                                                                                 <  1:750882/139‑1 (MQ=255)
                                 aCTCAGGGAGCTGGACGCCAGGGTCAGGTTTTTAAAGCCCATCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGcga                                                                                                 <  1:708995/139‑1 (MQ=255)
                                  cTCAGGGAGCTGGACGCCAGGGTCAGGTTTTTAAAGCCCATCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGcgac                                                                                                >  1:454042/1‑139 (MQ=255)
                                      gggAGCTGGACGCCAGGGTCAGGTTTTTAAAGCCCATCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAAt                                                                                            >  2:772047/1‑139 (MQ=255)
                                      gggAGCTGGACGCCAGGGTCAGGTTTTTAAAGCCCATCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAAt                                                                                            <  1:631148/139‑1 (MQ=255)
                                             ggACGCCAGGGTCAGGTTTTTAAAGCCCATCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATgggagga                                                                                                                                >  2:511140/1‑96 (MQ=255)
                                             ggACGCCAGGGTCAGGTTTTTAAAGCCCATCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATgggagga                                                                                                                                <  1:511140/96‑1 (MQ=255)
                                             ggACGCCAGGGTCAGGTTTTTAAAGCCCATCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGcttct                                                                                     >  1:431365/1‑139 (MQ=255)
                                                     gggTCAGGTTTTTAAAGCCCATCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGtt                                                                             >  2:595151/1‑139 (MQ=255)
                                                                 tAAAGCCCATCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGCTGGGCGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGttt                                                                 >  2:59492/1‑139 (MQ=255)
                                                                       ccATCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATc                                                           >  2:41492/1‑139 (MQ=255)
                                                                        cATCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATcg                                                          <  2:431365/139‑1 (MQ=255)
                                                                            ttCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATcgcgcg                                                      <  2:337555/139‑1 (MQ=255)
                                                                                  aTGACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATCGCGCGCttttt                                                <  1:59492/139‑1 (MQ=255)
                                                                                                     tGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATCGCGCGCTTTTTCAGCCTGGTAGGTTTGCtt                             >  1:105156/1‑139 (MQ=255)
                                                                                                       aCGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATCGCGCGCTTTTTCAGCCTGGTAGGTTTGCTTTg                           >  2:61003/1‑139 (MQ=255)
                                                                                                          gtcaggtcaCCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATCGCGCGCTTTTTCAGCCTGGTAGGTTTGCTTTGGCg                        >  2:757562/1‑139 (MQ=255)
                                                                                                          gtcaggtcaCCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATCGCGCGCTTTTTCAGCCTGGTAGGTTTGCTTTGGCg                        >  1:339518/1‑139 (MQ=255)
                                                                                                            caggtcaCCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATCGCGCGCTTTTTCAGCCTGGTAGGTTTGCTTTGGCGAg                      >  2:777768/1‑139 (MQ=255)
                                                                                                              ggtcaCCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGc                                                                                         >  1:745522/1‑70 (MQ=255)
                                                                                                              ggtcaCCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGc                                                                                         <  2:745522/70‑1 (MQ=255)
                                                                                                                    cGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATCGCGCGCTTTTTCAg                                             >  1:223842/1‑108 (MQ=255)
                                                                                                                    cGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATCGCGCGCTTTTTCAg                                             <  2:223842/108‑1 (MQ=255)
                                                                                                                    cGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATCGCGCGCTTTTTCAGCCTGGTAGGTTTGCTTTGGCGAGTTCTGGaa              <  1:595151/139‑1 (MQ=255)
                                                                                                                     gCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATCGCGCGCTTTTTCAGCCTGGTAGGTTTGCTTTGGCGAGTTCTGGaaa             >  1:635298/1‑139 (MQ=255)
                                                                                                                         cGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAAAGCGACGAATCGCTTCTTCCAgg                                                                               >  2:313283/1‑69 (MQ=255)
                                                                                                                         cGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAAAGCGACGAATCGCTTCTTCCAgg                                                                               <  1:313283/69‑1 (MQ=255)
                                                                                                                           gAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATCGCGCGCTTTTTCAGCCTGGTAGGTTTGCTTTGGCGAGTTCTGGAAAGCGGCa       <  2:454042/139‑1 (MQ=255)
                                                                                                                           gAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATCGCGCGCTTTTTCAGCCTGGTAGGTTTGCTTTGGCGAGTTCTGGAAAGCGGCa       <  2:105156/139‑1 (MQ=255)
                                                                                                                             aaGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATCGCGCGCTTTTTCAGCCTGGTAGGTTTGCTTTGGCGAGTTCTGGAAAGCGGCAAg     <  1:777768/139‑1 (MQ=255)
                                                                                                                                cGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATCGCGCGCTTTTTCAGCCTGGTAGGTTTGCTTTGGCGAGTTCTGGAAAGCGGCAAGATc  <  1:757562/139‑1 (MQ=255)
                                                                                                                                 gTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATCGCGCGCTTTTTCAGCCt                                          >  2:531374/1‑98 (MQ=255)
                                                                                                                                 gTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATCGCGCGCTTTTTCAGCCt                                          <  1:531374/98‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                  
GCGAATGTGTTCTACCAGCGGCGCATGGCAATCACTCAGGGAGCTGGACGCCAGGGTCAGGTTTTTAAAGCCCATCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGAAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATCGCGCGCTTTTTCAGCCTGGTAGGTTTGCTTTGGCGAGTTCTGGAAAGCGGCAAGATC  >  NC_000913/649236‑649502

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