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

GAATGTGTTCTACCAGCGGCGCATGGCAATCACTCAGGGAGCTGGACGCCAGGGTCAGGTTTTTAAAGCCCATCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGAAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATCGCGCGCTTTTTCAGCCTGGTAGGTTTGCTTTGGCGAGTTCTGGAAAGCGGCAAGATC  >  NC_000913/649238‑649502
                                                                                                                                      |                                                                                                                                  
gAATGTGTTCTACCAGCGGCGCATGGCAATCACTCAGGGAGCTGGACGCCAGGGTCAGGTTTTTAAAGCCCATCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATgggagga                                                                                                                                <  2:814014/139‑1 (MQ=255)
  aTGTGTTCTACCAGCGGCGCATGGCAATCAATCAGGGAGCTGGACGCCAGGGTCAGGGTTTTAAAGCCCATCTTGGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGaaa                                                                                                                              >  2:714716/1‑139 (MQ=255)
            ccAGCGGCGCATGGCAATCACTCAGGGAGCTGGACGCCAGGGTCAGGTTTTTAAAGCCCATCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGcc                                                                                                                    >  2:695645/1‑139 (MQ=255)
             cAGCGGCGCATGGCAATCACTCAGGGAGCTGGACGCCAGGGTCAGGTTTTTAAAGCCCATCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCg                                                                                                                   <  1:105903/139‑1 (MQ=255)
             cAGCGGCGCATGGCAATCACTCAGGGAGCTGGACGCCAGGGTCAGGTTTTTAAAGCCCATCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCg                                                                                                                   <  1:695645/139‑1 (MQ=255)
             cAGCGGCGCATGGCAATCACTCAGGGAGCTGGACGCCAGGGTCAGGTTTTTAAAGCCCATCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCg                                                                                                                   <  1:139933/139‑1 (MQ=255)
               gcggcgCATGGCAATCACTCAGGGAGCTGGACGCCAGGGTCAGGTTTTTAAAGCCCATCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCg                                                                                                                   >  2:230447/1‑137 (MQ=255)
               gcggcgCATGGCAATCACTCAGGGAGCTGGACGCCAGGGTCAGGTTTTTAAAGCCCATCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCg                                                                                                                   <  1:230447/137‑1 (MQ=255)
               gcggcgCATGGCAATCACTCAGGGAGCTGGACGCCAGGGTCAGGTTTTTAAAGCCCATCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTc                                                                                                                 >  1:87215/1‑139 (MQ=255)
                    gcATGGCAATCACTCAGGGAGCTGGACGCCAGGGTCAGGTTTTTAAAGCCCATCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTa                                                                                                            <  2:124445/139‑1 (MQ=255)
                                       agttGGACGCCAGGGTCAGGTTTTTAAAGCCCATCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGc                                                                                         <  1:84902/136‑1 (MQ=255)
                                               gCCAGGGTCAGGTTTTTAAAGCCCATCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGaaa                                                                                                                              >  2:188362/1‑94 (MQ=255)
                                               gCCAGGGTCAGGTTTTTAAAGCCCATCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGaaa                                                                                                                              <  1:188362/94‑1 (MQ=255)
                                                               tAAAGCCCATCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCACTGAAAGCGTGATgggagg                                                                                                                                 <  1:524891/75‑1 (MQ=255)
                                                               tAAAGCCCATCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATgggagg                                                                                                                                 >  2:524891/1‑75 (MQ=255)
                                                               tAAAGCCCATCTTCGCGATGACGTCCATCACCATATTGACAGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCg                                                                                                                   <  1:640410/89‑1 (MQ=255)
                                                               tAAAGCCCATCTTCGCGATGACGTCCATCACCATATTGACAGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCg                                                                                                                   >  2:640410/1‑89 (MQ=255)
                                                                                  gACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGca                                                                       >  1:410543/1‑114 (MQ=255)
                                                                                  gACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGca                                                                       <  2:410543/114‑1 (MQ=255)
                                                                                                tcttGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCACACAGTTTGCGATCGCGCGCTTTTTCAg                                             >  2:503195/3‑126 (MQ=255)
                                                                                                tatTGACGGTCAGGTTACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATCGCGCGCTTTTTCAg                                             <  1:503195/126‑1 (MQ=255)
                                                                                                 atTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATCGCGCGCTTTTTCAGCCTGGTAGGTTTGc                               >  2:58545/1‑139 (MQ=255)
                                                                                                  ttGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCg                                                                                                                           <  2:727298/46‑1 (MQ=255)
                                                                                                  ttGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCg                                                                                                                           >  1:727298/1‑46 (MQ=255)
                                                                                                        gtcaggtcaCCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCg                                                                                                                   >  2:375221/1‑48 (MQ=255)
                                                                                                        gtcaggtcaCCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCg                                                                                                                   <  1:375221/48‑1 (MQ=255)
                                                                                                         tcaggtcaCCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGca                                                                       <  1:132488/91‑1 (MQ=255)
                                                                                                         tcaggtcaCCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGca                                                                       >  2:132488/1‑91 (MQ=255)
                                                                                                                         gAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATCGCGCGCTTTTTCAGCCTGGTAGGTTTGCTTTGGCGAGTTCTGGAAAGCGGCa       <  1:58545/139‑1 (MQ=255)
                                                                                                                              cGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATCGCGCGCTTTTTCAGCCTGGTAGGTTTGCTTTGGCGAGTTCTGGAAAGCGGCAAGATc  <  2:87215/139‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                  
GAATGTGTTCTACCAGCGGCGCATGGCAATCACTCAGGGAGCTGGACGCCAGGGTCAGGTTTTTAAAGCCCATCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGAAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATCGCGCGCTTTTTCAGCCTGGTAGGTTTGCTTTGGCGAGTTCTGGAAAGCGGCAAGATC  >  NC_000913/649238‑649502

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