Predicted mutation
evidence seq id position mutation annotation gene description
RA AM260479 245,288 T→G L247V (TTG→GTG)  betA2 → choline dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AM260479245,2880TG100.0% 108.0 / NA 29L247V (TTG→GTG) betA2choline dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (18/11);  total (18/11)

ACCGGGCCTGCAATGGCGGCCGCGCGCACCTGCGGGTGCTGACCGATACGCAGGCGCTGCGCATAGTCTTCGAAGGCCGCCATGCGGCCGGGGTGCAGTTGGTGCGCGAGGGCCACCGGCAACTGCTGCGCGCGCGCCGCGACGTGATCGTGTGCGCGGGCACGTTCGGCTCGCCGCAGCTGCTGATGGTGTCGGGCGTCGGCCCCGCCGCGCA  >  AM260479/245190‑245403
                                                                                                  |                                                                                                                   
aCCGGGCCTGCAATGGCGGCCGCGCGCACCTGCGGGTGCTGACCGATACGCAGGCGCTGCGCATAGTCTTCGAAGGCCGCCATGCGGCCGGGGTGCAGGTGGTGCGCGAGGGCCACCGGCAACTGCTgcgcgcg                                                                                  >  2:437716/1‑134 (MQ=255)
        tGCAATGGCGGCCGCGCGCACCTGCGGGTGCTGACCGATACGCAGGCGCTGCGCATAGTCTTCGAAGGCCGCCATGCGGCCGGGGTGCAGGTGGTGCGCGAGGGCCACCGGCAACTGCTGCGCGCGCGCCGCGa                                                                          >  1:1045647/1‑134 (MQ=255)
                           aCCTGCGGGTGCTGACCGATACGCAGGCGCTGCGCATAGTCTTCGAAGGCCGCCATGCGGCCGGGGTGCAGGTGGTGCGCGAGGGCCACCGGCAACTGCTGCGCGCGCGCCGCGACGTGATCGTGTGCGCGGGc                                                       <  2:1045647/134‑1 (MQ=255)
                                    tGCTGACCGATACGCAGGCGCTGCGCATAGTCTTCGAAGGCCGCCATGCGGCCGGGGTGCAGGTGGTGCGCGAGGGCCACCGGCAACTGCTgcgcgcgc                                                                                 >  2:262591/1‑99 (MQ=255)
                                    tGCTGACCGATACGCAGGCGCTGCGCATAGTCTTCGAAGGCCGCCATGCGGCCGGGGTGCAGGTGGTGCGCGAGGGCCACCGGCAACTGCTgcgcgcgc                                                                                 <  1:262591/99‑1 (MQ=255)
                                         aCCGATACGCAGGCGCTGCGCATAGTCTTCGAAGGCCGCCATGCGGCCGGGGTGCAGGTGGTGCGCGAGGGCCACCGGCAACTGCTGCGCGCGCGCCGCGACGTGATCGTGTGCGCGGGCACGTTCGGCTcgc                                          <  1:437716/133‑1 (MQ=255)
                                                 gCAGGCGCTGCGCATAGTCTTCGAAGGCCGCCATGCGGCCGGGGTGCAGGTGGTGCGCGAgg                                                                                                         >  2:1030275/1‑62 (MQ=255)
                                                 gCAGGCGCTGCGCATAGTCTTCGAAGGCCGCCATGCGGCCGGGGTGCAGGTGGTGCGCGAgg                                                                                                         <  1:1030275/62‑1 (MQ=255)
                                                 gCAGGCGCTGCGCATAGTCTTCGAAGGCCGCCATGCGGCCGGGGTGCAGGTGGTGCGCGAGGGCCACCGGCAACTGCTGCGCGCGCGCCGCGACGTGATCGTGTGCGCGGGCACGTTCGGCTCGCCGCAgctgc                                 >  2:585572/1‑134 (MQ=255)
                                                  cAGGCGCTGCGCATAGTCTTCGAAGGCCGCCATGCGGCCGGGGTGCAGGTGGTGCGCGAGGGCCACCGGCAACTGCTGCGCGCGCGCCGCGACGTGATCGTGTGCGCGGGCACGTTCGGCTCGCCGCAgctgct                                >  1:186345/1‑134 (MQ=255)
                                                  cAGGCGCTGCGCATAGTCTTCGAAGGCCGCCATGCGGCCGGGGTGCAGGTGGTGCGCGAGGGCCACCGGCAACTGCTGCGCGCGCGCCGCGACGTGATCGTGTGCGCGGGCACGTTCGGCTCGCCGCAgctgct                                >  1:271497/1‑134 (MQ=255)
                                                       gcTGCGCATAGTCTTCGAAGGCAGCCATGCGGCCGGGGTGCAGGTGGTGCGCGAGGGCCACCGGCAACTGCTGCGCGCGCGCCGCGACGTGATCGTGTGCGCGGGCACGGTCGGCTCGCCGCAGCTGCTGATgg                           >  1:129349/1‑134 (MQ=255)
                                                            gcATAGTCTTCGAAGGCCGCCATGCGGCCGGGGTGCAGGTGGTGCGCGAGGGCCACCGGCAACTGCTgcgcgcgcgc                                                                               <  1:749687/77‑1 (MQ=255)
                                                            gcATAGTCTTCGAAGGCCGCCATGCGGCCGGGGTGCAGGTGGTGCGCGAGGGCCACCGGCAACTGCTgcgcgcgcgc                                                                               >  2:749687/1‑77 (MQ=255)
                                                            gcATAGTCTTCGAAGGCCGCCATGCGGCCGGGGTGCAGGTGGTGCGCGAGGGCCACCGGCAACTGCTGCGCGCGCGCCGCGACGTGATCgtg                                                                >  2:185064/1‑92 (MQ=255)
                                                            gcATAGTCTTCGAAGGCCGCCATGCGGCCGGGGTGCAGGTGGTGCGCGAGGGCCACCGGCAACTGCTGCGCGCGCGCCGCGACGTGATCgtg                                                                <  1:185064/92‑1 (MQ=255)
                                                               tAGTCTTCGAAGGCCGCCATGCGGCCGGGGTGCAGGTGGTGCGCGAGGGCCACCGGCAACTGCTGCGCGCGCGCCGCGACGTGATCGTGTGCGCGGGCACGTTCGGCTCGCCGCAGCTGCTGATGGTGTCGGGc                   >  1:33191/1‑134 (MQ=255)
                                                                 gTCTTCGAAGGCCGCCATGCGGCCGGGGTGCAGGTGGTGCGCGAGGGCCACCGGCAACTGCTGCGCGCGCGCCGCGACGTGATCGTGTGCGCGGGCACGTTCGGCTCGCCGCAGCTGCTGATGGTGTCGGGCGt                 >  2:558761/1‑134 (MQ=255)
                                                                         aGGCCGCCATGCGGCCGGGGTGCAGGTGGTGCGCGAGGGc                                                                                                       >  1:112403/1‑40 (MQ=255)
                                                                         aGGCCGCCATGCGGCCGGGGTGCAGGTGGTGCGCGAGGGc                                                                                                       <  2:112403/40‑1 (MQ=255)
                                                                         aGGCCGCCATGCGGCCGGGGTGCAGGTGGTGCGCGAGGGCCACCGGCAACTGCTgcgcgcgc                                                                                 <  1:403446/62‑1 (MQ=255)
                                                                         aGGCCGCCATGCGGCCGGGGTGCAGGTGGTGCGCGAGGGCCACCGGCAACTGCTgcgcgcgc                                                                                 >  2:403446/1‑62 (MQ=255)
                                                                         aGGCCGCCATGCGGCCGGGGTGCAGGTGGTGCGCGAGGGCCACCGGCAACTGCTGCGCGCGCGCCGCGa                                                                          <  1:2293/69‑1 (MQ=255)
                                                                         aGGCCGCCATGCGGCCGGGGTGCAGGTGGTGCGCGAGGGCCACCGGCAACTGCTGCGCGCGCGCCGCGa                                                                          >  2:2293/1‑69 (MQ=255)
                                                                            ccgccATGCGGCCGGGGTGCAGGTGGTGCGCGAGGGCCACCGGCAACTGCTGCGCGCGCGCCGCGACGTGATCGTGTGCGCGGGCACGTTCGGCTCGCCGCAGCTGCTGATGGTGTCGGGCGTCGGCCccgccg      <  1:929347/134‑1 (MQ=255)
                                                                                cATGCGGCCGGGGTGCAGGTGGTGCGCGAGGGCCACCGGCAACTGCTGCGCGCGCGCCGCGACGTGATCGTGTGCGCGGGCACGTTCGGCTCGCCGCAGCTGCTGATGGTGTCGGGCGTCGGCCCCGCCGCGCa  >  2:20772/1‑134 (MQ=255)
                                                                                cATGCGGCCGGGGTGCAGGTGGGGCGCGAGGGCCACCGGCAACTGCTGCGCGCGCGCCGCGACGTGATCGTGTGCGCGGGCACGTTCGGCTCGCCGCAGCTGCTGATGGTGTCGGGCGTCGGCCCCGCCGCGCa  >  2:922245/1‑134 (MQ=255)
                                                                                  tGCGGCCGGGGTGCAGGTGGTGCGCGAGGGCCACCGGCAACTGCTGCGCGCGCGCCGCGACGTGATCGtgtg                                                              >  2:414972/1‑72 (MQ=255)
                                                                                  tGCGGCCGGGGTGCAGGTGGTGCGCGAGGGCCACCGGCAACTGCTGCGCGCGCGCCGCGACGTGATCGtgtg                                                              <  1:414972/72‑1 (MQ=255)
                                                                                                  |                                                                                                                   
ACCGGGCCTGCAATGGCGGCCGCGCGCACCTGCGGGTGCTGACCGATACGCAGGCGCTGCGCATAGTCTTCGAAGGCCGCCATGCGGCCGGGGTGCAGTTGGTGCGCGAGGGCCACCGGCAACTGCTGCGCGCGCGCCGCGACGTGATCGTGTGCGCGGGCACGTTCGGCTCGCCGCAGCTGCTGATGGTGTCGGGCGTCGGCCCCGCCGCGCA  >  AM260479/245190‑245403

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