Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 126,947 G→T 100% R418L (CGT→CTT)  aceF → pyruvate dehydrogenase, dihydrolipoyltransacetylase component E2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913126,9470GT100.0% 100.3 / NA 33R418L (CGT→CTT) aceFpyruvate dehydrogenase, dihydrolipoyltransacetylase component E2
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (22/11);  total (22/11)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

AAACGTGCAGAAGCAGCTCCGGCAGCGACTGGCGGTGGTATCCCTGGCATGCTGCCGTGGCCGAAGGTGGACTTCAGCAAGTTTGGTGAAATCGAAGAAGTGGAACTGGGCCGCATCCAGAAAATCTCTGGTGCGAACCTGAGCCGTAACTGGGTAATGATCCCGCATGTTACTCACTTCGACAAAACCGATATCACCGAGTTGGAAGCGTTCCGTAAACAGCAGAACGAAGAAGCGGCGAAACGTAAGCTGGATGTGAAGATCACCCCGGTTGTCTTC  >  NC_000913/126802‑127080
                                                                                                                                                 |                                                                                                                                     
aaaCGTGCAGAAGCAGCCCCGGCAGCGACTGGCGGTGGTATCCCTGGCATGCTGCCGTGGCCGAAGGTGGACTTCAGCAAGTTTGGTGAAATCGAAGAAGTGGAACTGGGCCGCATCCAGAAAATCTCTGGTGCGAACCTGAGCCTTaa                                                                                                                                    <  2:234139/149‑1 (MQ=255)
 aaCGTGCAGAAGCAGCTCCGGCAGCGACTGGCGGTGGTATCCCTGGCATGCTGCCGTGGCCGAAGGTGGACTTCAGCAAGTTTGGTGAAATCGAAGAAGTGGAACTGGGCCGGATCCAGAAAATCTCTGGTGCGAACCTGAGCCTTAAc                                                                                                                                   >  1:732044/1‑149 (MQ=255)
                               gCGGTGGTATCCCTGGCATGCTGCCGTGGCCGAAGGTGGACTTCAGCAAGTTTGGTGAAATCGAAGAAGTGGAACTGGGCCGCATCCAGAAAATCTCTGGTGCGAACCTGAGGCTTAACTGGGTAGTGGTCCCGCATGTTAGTccgttg                                                                                                     >  1:230264/1‑144 (MQ=255)
                                    ggtATCCCTGGCATGCTGCCGTGGCCGAAGGTGGACTTCAGCAAGTTTGGTGAAATCGAAGAAGTGGAACTGGGCCGCATCCAGAAAATCTCTGGTGCGAACCTGAGCCTTAACTGGGTAATGATCCCGCATGTTACTCACTTCGACaa                                                                                                >  1:106828/1‑149 (MQ=255)
                                    ggtATCCCTGGCATGCTGCCGTGGCCGAAGGTGGACTTCAGCAAGTTTGGTGAAATCGAAGAAGTGGAACTGGGCCGCATCCAGAAAATCTCTGGTGCGAACCTGAGCCTTAACTGGGTAATGATCCCGCATGTTACTCACTTCGACaa                                                                                                >  1:276726/1‑149 (MQ=255)
                                       aTCCCTGGCATGCTGCCGTGGCCGAAGGTGGACTTCAGCAAGTTTGGTGGAATCGAAGAAGTGGAACTGGGCCGCATCCAGAAAATCTCTGGTGCGAACCTGAGCCTTAACTGGGTAATGATCCCGCATGTTACTCACTTCGACAAAAc                                                                                             >  1:158121/1‑149 (MQ=255)
                                       aTCCCTGGCATGCTGCCGTGGCCGAAGGTGGACTTCAGCAAGTTTGGTGAAATCGAAGAAGTGGAACTGGGCCGCATCCAGAAAATCTCTGGTGCGAACCTGAGCCTTAACTGGGTAATGATCCCGCATGTTACTCACTTCGACAAAAc                                                                                             <  2:16409/149‑1 (MQ=255)
                                              gCATGCTGCCGTGGCCGAAGGTGGACTTCAGCAAGTTTGGTGAAATGGAAGAAGTGGAACTGGGCCGCATCCAGAAAATCTCTGGTGCGAACCTGAGCCTTAACTGGGTAATGATCCCGCATGTTCCTCACTTAGACAAACCAGCTATc                                                                                      >  2:154797/1‑149 (MQ=255)
                                              gCATGCTGCCGTGGCCGAAGGTGGACTTCAGCAAGTTTGGTGAAATCGAAGAAGTGGAACTGGGCCGCATCCAGAAAATCTCTGGTGCGAACCTGAGCCTTAACTGGGTAATGATCCCGCATGTTACTCACTTCGACAAAACCGATATc                                                                                      >  2:777819/1‑149 (MQ=255)
                                                 tgctgcCGTGGCCGAAGGTGGACTTCAGCAAGTTTGGTGAAATCGAAGAAGTGGAACTGGGCCGCATCCAGAAAATCTCTGGTGCGAACCTGAGCCTTAACTGGGTAATGATCCCGCATGTTACTGACTTCGACAAAACCGATATCAcc                                                                                   >  2:259173/1‑149 (MQ=255)
                                                 tgctgcCGTGGCCGAAGGTGGACTTCAGCAAGTTTGGTGAAATCGAAGAAGTGGAACTGGGCCGCATCCAGAAAATCTCTGGTGCGAACCTGAGCCTTAACTGGGTAATGATACCGCATGTTACTCACTTCGACAAAACCGATATCAcc                                                                                   >  2:259158/1‑149 (MQ=255)
                                                     gcCGTGGCCGAAGGTGGACTTCAGCAAGTTTGGTGACATCGAAGAAGTGGAACTGGGCCGCAAACAGAAAAACTCTGGTGCGAACATGAGTCTTAACTGGGTAATGATGCCGCATGTTACTCACTTCGACAAAACCGATATGACCGAGt                                                                               >  2:70952/1‑149 (MQ=255)
                                                        gTGGCCGAAGGTGGACTTCAGCAAGTTTGGTGAAATCGAAGAAGTGGAACTGGGCCGCATCCAGAAAATCTCTGGTGCGAACCTGAGCCTTAACTGGGTAATGATCCCGCATGTTACTCACTTAGACAAAACCGATATCACCGAGTTgg                                                                            >  1:71715/1‑149 (MQ=255)
                                                             cGAAGGTGGACTTCAGCAAGTTTGGTGAAATCGAAGAAGTGGAACTGGGCCGCATCCAGAAAATCTCTGGTGCGAACCTGAGCCTTAACTGGGTAATGATCCCGCATGTTACTCACTTCGACAAAACCGATCTCACCGAGTTGGagccg                                                                       >  1:178055/1‑145 (MQ=255)
                                                             cGAAGGTGGACTTCAGCAAGTTTGGTGAAATCGAAGAAGTGGAACTGGGCCGCATCCAGAAAATCTCTGGTGCGAACCTGAGCCTTAACTGGGTAATGATCCCGCATGTTACTCACTTCGACAAAACCGATATCACCGAGTTGGAAGCg                                                                       >  2:525429/1‑149 (MQ=255)
                                                                 ggTGGACTTCAGCAAGTTTGGTGAAATCGAAGAAGTGGAACTGGGCCGCATCCAGAAAATCTCTGGTGCGAACCTGAGCCTTAACTGGGTAATGATCCCGCATGTTACTCACTTCGACAAAACCGATATCACCGAGTTGGAAGCGTTcc                                                                   <  1:790927/149‑1 (MQ=255)
                                                                        ttCAGCAAGTTTGGTGAAATCGAAGAAGTGGAACTGGGCCGCATCCAGAAAATCTCTGGTGCGAACCTGAGCCTTAACTGGGTAATGATCCCGCATGTTACTCACTTCGACAAAACCGATATCACCGCGTTGGAAGCGTTCCGTAAAca                                                            >  2:419498/1‑149 (MQ=255)
                                                                         gcAGCAAGTTTGGTAAAATCGAAGAAGTGGAAGTGGGCCGCATCCAGAAAACCTCTGGTGCGAACCTGAGCCTTAACTGGGTAATGATCCCGCATGTTACTCACTTCGACAAAACCGATATCACCGAGTTGGAAGCGTTCCGTAAAcag                                                           <  1:212676/148‑1 (MQ=255)
                                                                              aaGTTTGGTGAAATCGAAGAAGTGGAACTGGGCCGCATCCAGAAAATCTCTGGTGCGAACCTGAGCCTTAACTGGGTAATGATCCCGCATGTTACTCACTTCGACAAAACCGATATCACCGAGTTGGAAGCGTTCCGTAAACAGCAGaa                                                      <  1:737622/149‑1 (MQ=255)
                                                                                    ggggAAATCGAAGAAGTGGCCCTGGCCCGCCTCCACACAATCTCGAGTGCCACCATGATGCTTAGCTGGGTTAGGAGCCGGCAGGGTACTCACTTCGACAAATCCGATATCACCGAGTTGGAAGCGTTCCGTAAACAGCAGAACgaaga                                                <  2:666975/146‑1 (MQ=255)
                                                                                    ggTGAAATCGAAGAAGTGGAACTGGGCCGCATCCAGAAAATCTCTGGTGCGAACCTGAGCCTTAACTGGGTAATGATCCCGCATGTTACTCACTTCGACAAAACCGATATCACCGCGTTGGAAGCGTTCCGTAAACAGCAGAACgaaga                                                >  2:301895/1‑149 (MQ=255)
                                                                                         aaTCGAAGAAGTGGAACTGGGCCGCATCCAGAAAATCTCTGGTGCGAACCTGAGCCTTAACTGGGTAATGATCCCGCATGTTACTCACTTCGACAAAACCGATATCACCGAGTTGGAAGCGTTCCGTAAACAGCAGAACGAAGAAgcgg                                           >  2:501466/1‑149 (MQ=255)
                                                                                          aTCGAAGAAGTGGAACTGGGCCGCATCCAGAAAATCTCGGGTGCGAACCTGAGCCTTAACTGGGTAATGATCCCGCATGTTACTCACTTCGACAAAACCGATATCACCGAGTTGGAAGCGTTCCGTAAACAGCAGAACGAAGAAgcggc                                          <  1:140278/149‑1 (MQ=255)
                                                                                            cGAAGAAGTGGAACTGGGCCGCATCCAGAAAATCTCTGGTGCGAACCTGAGCCTTAACTGGGTAATGATCCCGCATGTTACTCACTTCGACAAAACCGATATCACCGAGTTGGAAGCGTTCCGTAGACAGCAGAACGAAGAAGCAGCGa                                        >  1:318411/1‑149 (MQ=255)
                                                                                              aagaagTGGAACTGGGCCGCATCCAGAAAATCTCTGGTGCGAACCTGAGCCTTAACTGGGTAATGATCCCGCATGTTACTCACTTCGACAAAACCGATATCACCGAGTTGGAAGCGTTCCGTAAACAGCAGAACGAAGAAGCGGCGaaa                                      >  1:641246/1‑149 (MQ=255)
                                                                                                     ggAACTGGGCCGCATCCAGAAAATCTCTGGTGCGAACCTGAGCCTTAACTGGGTAATGATCCCGCATGTTACTCACTTCGACAAAACCGATATCACCGAGTTGGAAGCGTTCCGTAAACAGCAGAACGAAGAAGCGGCGAAACGTAAGc                               <  2:125934/149‑1 (MQ=255)
                                                                                                     ggAACGGGGCCGCATCCAGAAAATCTCTGGTGCGAACCTGAGCCTTAACTGGGTAATGATCCCGCATGTTACTCACTTCGACAAAACCGATATCACCGAGTTGGAAGCGTTCCGTAAACAGCAGAACGAAGAAGCGGCGAAACGTAAGc                               <  2:64180/149‑1 (MQ=255)
                                                                                                          tGGGCCGCCTCCAGAAAATCTCTGGTGCGAACCTGAGCCTTAACTGGGTAATGATCCCGCATGTTACTCACTTCGACAAAACCGATATCACCGAGTTGGAAGCGTTCCGTAAACAGCAGAACGAAGAAGCGGCGAAACGTAAGCTGGAt                          <  1:650893/149‑1 (MQ=255)
                                                                                                           gggCCGCATCCAGAAAATCTCTGGTGCGAACCTGAGACTTAACTGGGTAATGATCCCGCATGTTAGTCACTTAGACAAAACGGATATCACCGAGTTGGAAGCGTTCCGTGAACAGCAGAACGAAGAAGCGGGGAAACGTAAGCTGGAtg                         >  2:415818/1‑149 (MQ=255)
                                                                                                                 cATCCAGAAAATCTCTGGTGCGAACCTGGGCCTTAACTGGGTAATGATCCCGCATGTTACTCACTTCGACAAAACCGATATCACCGAGTTGGAAGCGTTCCGTAAACAGCAGAACGAAGAAGCGGCGAAACGTAAGCTGGATGTGAAGa                   <  2:244922/149‑1 (MQ=255)
                                                                                                                 cATCCAGAAAATCTCTGGTGCGAACCTGAGCCTTAACTGGGTAATGATCCCGCATGTTACTCACTTCGACAAAACCGATATCACCGAGTTGGAAGCGTTCCGTAAACAGCAGAACGAAGAAGCGGCGAAACGTAAGCTGGATGTGAAgg                   >  1:118160/1‑148 (MQ=255)
                                                                                                                      aGAAAATCTCTGGTGCGAACCTGAGCCTTAACTGGGTAATGATCCCGCATGTTACTCACTTCGACAAAACCGATATCACCGAGTTGGAAGCGTTCCGTAAACCGCAGAACGAAGAAGCGGCGAAACGTCAGCTGCACGTGAAGCGCAcc              >  2:31544/1‑149 (MQ=255)
                                                                                                                                  gTGCGAACCTGAGCCTTAACTGGGTAATGATCCCGCATGTTACTCACTTCGACAAAACCGATATCACCGAGTTGGAAGCGTTCCGTAAACAGCAGAACGAAGAAGCGGCGAAACGTAAGCTGGATGTGAAGATCACCCCGGTTGTCtgc  >  2:417834/1‑147 (MQ=255)
                                                                                                                                                 |                                                                                                                                     
AAACGTGCAGAAGCAGCTCCGGCAGCGACTGGCGGTGGTATCCCTGGCATGCTGCCGTGGCCGAAGGTGGACTTCAGCAAGTTTGGTGAAATCGAAGAAGTGGAACTGGGCCGCATCCAGAAAATCTCTGGTGCGAACCTGAGCCGTAACTGGGTAATGATCCCGCATGTTACTCACTTCGACAAAACCGATATCACCGAGTTGGAAGCGTTCCGTAAACAGCAGAACGAAGAAGCGGCGAAACGTAAGCTGGATGTGAAGATCACCCCGGTTGTCTTC  >  NC_000913/126802‑127080

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