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% 102.6 / NA 27L247V (TTG→GTG) betA2choline dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (18/9);  total (18/9)

GTGCCTACCTGCACGGTGGCAACAAGCGCGACCGGGCCTGCAATGGCGGCCGCGCGCACCTGCGGGTGCTGACCGATACGCAGGCGCTGCGCATAGTCTTCGAAGGCCGCCATGCGGCCGGGGTGCAGTTGGTGCGCGAGGGCCACCGGCAACTGCTGCGCGCGCGCCGCGACGTGATCGTGTGCGCGGGCACGTTCGGCTCGCCGCAGCTGCTGATGGTGTCGGGCGTCGGCCCCGCCGCGCACTTGCG  >  AM260479/245160‑245409
                                                                                                                                |                                                                                                                         
gTGCCTACCTGCACGGTGGCAACAAGCGCGACCGGGCCTGCAATGGCGGCCGCGCGCACCTGCGGGTGCTGACCGATACGCAGGCGCTGCGCATAGTCTTCGAAGGCCGCCATGCGGCCGGGGTGCAGgtggtg                                                                                                                      <  1:557111/134‑1 (MQ=255)
    cTACCTGCACGGTGGCAACAAGCGCGACCGGGCCTGCAATGGCGGCCGCGCGCACCTGCGGGTGCTGACCGATACGCAGGCGCTGCGCATAGTCTTCGAAGGCCGCCATGCGGCCGGGGTGCAGGTGGTGcgcg                                                                                                                  <  1:636611/134‑1 (MQ=255)
                tGGCAACAAGCGCGACCGGGCCTGCAATGGCGGCCGCGCGCACCTGCGGGTGCTGACCGATACGCAGGCGCTGCGCATAGTCTTCGAAGGCCGCCATGCGGCCGGGGTGCAGGTGGTGCGCGAGGGCCACCgg                                                                                                       >  2:652339/1‑133 (MQ=255)
                    aacaaGCGCGACCGGGCCTGCAATGGCGGCCGCGCGCACCTGCGGGTGCTGACCGATACGCAGGCGCTGCGCATAGTCTTCGAAGGCCGCCATGCGGCCGGGGTGCAGGTGGTGCGCGAGGGc                                                                                                             >  2:994743/1‑123 (MQ=255)
                    aacaaGCGCGACCGGGCCTGCAATGGCGGCCGCGCGCACCTGCGGGTGCTGACCGATACGCAGGCGCTGCGCATAGTCTTCGAAGGCCGCCATGCGGCCGGGGTGCAGGTGGTGCGCGAGGGc                                                                                                             <  1:994743/123‑1 (MQ=255)
                          cgcgACCGGGCCTGCAATGGCGGCCGCGCGCACCTGCGGGTGCTGACCGATACGCAGGCGCTGCGCATAGTCTTCGAAGGCCGCCATGCGGCCGGGGTGCAGGTGGTGCGCGAGGGCCACCGGCAACTGCTgcg                                                                                            <  2:76606/134‑1 (MQ=255)
                              aCCGGGCCTGCAATGGCGGCCGCGCGCACCTGCGGGTGCTGACCGATACGCAGGCGCTGCGCATAGTCTTCGAAGGCCGCCATGCGGCCGGGGTGCAGGTGGTGCGCGAGGGCCACCGGCAACTGCTgcgcgcg                                                                                        >  2:948714/1‑134 (MQ=255)
                                    ggccgcaTGGCGGCCGCGCGCACCTGCGGGTGCTGACCGATACGCAGGCGCTGCGCATAGTCTTCGAAGGCCGCCATGCGGCCGGGGTGCAGGTGGTGCGCGAGGGCCACCGGCAACTGCTGCGCGCGcgccgc                                                                                  >  1:4313/7‑134 (MQ=255)
                                                                       aCCGATACGCAGGCGCTGCGCATAGTCTTCGAAGGCCGCCATGCGGCCGGGGTGCAGGTGGTGCGCGAGGGCCACCGGCAACTGCTGCGCGCGCGCCGCGACGTGATCGTGTGCGCGGGCACGTTCGGCTcgcc                                               >  2:362876/1‑134 (MQ=255)
                                                                         cGATACGCAGGCGCTGCGCATAGTCTTCGAAGGCCGCCATGCGGCCGGGGTGCAGGTGGTGCGCGAGGGCCACCGGCAACTGCTGCGCGCGCGCCGCGACGTGATCGTGTGCGCGGGCACGTTCGGCTcgccg                                              >  1:373578/1‑133 (MQ=255)
                                                                            tACGCAGGCGCTGCGCATAGTCTTCGAAGGCCGCCATGCGGCCGGGGTGCAGGTGGTGCGCGAGGGCCACCGGCAACTGCTGCGCGCGCGCCGCGACGTGATCGTGTGCGCGGGCACGTTCGGCTCGCCGCAgc                                          >  2:889239/1‑134 (MQ=255)
                                                                               gCAGGCGCTGCGCATAGTCTTCGAAGGCCGCCATGCGGCCGGGGTGCAGGTGGTGCGCGAGGGCCACCGGCAACTGCTGCGCGCGCGCCGCGACGTGATCGTGTGCGCGGGCACGTTCGGCTCGCCGCAgctgc                                       >  2:39643/1‑134 (MQ=255)
                                                                                       tGCGCATAGTCTTCGAAGGCCGCCATGCGGCCGGGGTGCAGGTGGTGCGCGAGGGCCACCGGCAACTGCTGCGCGCGCGCCGCGACGTGATCGTGTGCGCGGGCACGTTCGGCTCGCCGCAGCTGCTGATGgtg                               >  2:914179/1‑134 (MQ=255)
                                                                                             tAGTCTTCGAAGGCCGCCATGCGGCCGGGGTGCAGGTGGTGCGCGAGGGCCACCGGCAACTGCTGCGCGCGCGCCGCGACGt                                                                             <  1:272565/82‑1 (MQ=255)
                                                                                             tAGTCTTCGAAGGCCGCCATGCGGCCGGGGTGCAGGTGGTGCGCGAGGGCCACCGGCAACTGCTGCGCGCGCGCCGCGACGt                                                                             >  2:272565/1‑82 (MQ=255)
                                                                                                       aGGCCGCCATGCGGCCGGGGTGCAGGTGGTGCGCGAGGGCCACCGGCAACTGCTGCGCGCGCGCCGCGACGTGATCGTGTGCGCGGGCACGTTCGGCTCGCCGCAGCTGCTGATGGTGTCGGGCGTCGGCCccg               <  2:4313/134‑1 (MQ=255)
                                                                                                       aGGCCGCCATGCGGCCGGGGTGCAGGTGGTGCGCGAGGGCCACCGGCAACTGCTGCGCGCGCGCCGCGACGTGATCGTGTGCGCGGGCACGTTCGGCTCGCCGCAGCTGCTGATGGTGTCGGGCGTCGGCCccg               >  1:299061/1‑134 (MQ=255)
                                                                                                        ggCCGCCATGCGGCCGGGGTGCAGGTGGTGCGCGAGGGCCACCGGCAACTGCTgcgcgcgcgc                                                                                     <  2:586758/63‑1 (MQ=255)
                                                                                                        ggCCGCCATGCGGCCGGGGTGCAGGTGGTGCGCGAGGGCCACCGGCAACTGCTgcgcgcgcgc                                                                                     >  1:586758/1‑63 (MQ=255)
                                                                                                        ggCCGCCATGCGGCCGGGGTGAAGGTGGTGCGCGAGGGCTACCGGAAACCGCTGCGCGCGCGCCGCGACGTGATCGTGTGCGCGGGCACGTTCGGCTCGCCGCAGCTGCTGATGGTGTCGGGCGTCGGCCccgc              <  1:706112/134‑1 (MQ=255)
                                                                                                              cATGCGGCCGGGGTGCAGGTGGTGCGCGAGGGCCACCGGCAACTGCTGCGCGCGCGCCGCGACGTGATCGTGTGCGCGGGCACGTTCGGCTCGCCGCAGCTGCTGATGGTGTCGGGCGTCGGCCCCGCcgcgta        >  1:308299/1‑132 (MQ=255)
                                                                                                              cATGCGGCCGGGGTGCAGGTGGTGCGCGAGGGCCACCGGCAACTGCTGCGCGCGCGCCGCGACGTGATCGTGTGCGCGGGCACGTTCGGCTCGCCGCAGCTGCTGATGGTGTCGGGCGTCGGCCCCGCCGCGCa        >  2:395150/1‑134 (MQ=255)
                                                                                                              cATGCGGCCGGGGTGCAGGTGGTGCGCGAGGGCCACCGGCAACTGCTGCGCGCGCGCCGCGACGTGATCGTGTGCGCGGGCACGTTCGGCTCGCCGCAGCTGCTGATGGTGTCGGGCGTCGGCCCCGCCGCGCa        >  1:776817/1‑134 (MQ=255)
                                                                                                               aTGCGGCCGGGGTGCAGGTGGTGCGCGAGGGCCACCGGCAACTGCTGCGCGCGCGCCGCGACGTGATCGTGTGCGCGGGCACGTTCGGCTCGCCGCAGCTGCTGATGGTGTCGGGCGTCGGCCCCGCCGCGCAc       >  2:350490/1‑134 (MQ=255)
                                                                                                                tGCGGCCGGGGTGCAGGTGGTGCGCGAGGGCCACCGGCAACTGCTGCGCGCGCGCCGCGAc                                                                               >  2:108643/1‑61 (MQ=255)
                                                                                                                tGCGGCCGGGGTGCAGGTGGTGCGCGAGGGCCACCGGCAACTGCTGCGCGCGCGCCGCGAc                                                                               <  1:108643/61‑1 (MQ=255)
                                                                                                                  cGGCCGGGGTGCAGGTGGTGCGCGAGGGCCACCGGCAACTGCTGCGCGCGCGCCGCGACGTGATCGTGTGCGCGGGCACGTTCGGCTCGCCGCAGCTGCTGATGGTGTCGGGCGTCGGCCCCGCCGCGCACTTg    <  1:914179/134‑1 (MQ=255)
                                                                                                                    gCCGGGGTGCAGGTGGTGCGCGAGGGCCACCGGCAACTGCTGCGCGCGCGCCGCGACGTGATCGTGTGCGCGGGCACGTTCGGCTCGCCGCAGCTGCTGATGGTGTCGGGCGTCGGCCCCGCCGCGCACTTgcg  >  1:615233/1‑134 (MQ=255)
                                                                                                                                |                                                                                                                         
GTGCCTACCTGCACGGTGGCAACAAGCGCGACCGGGCCTGCAATGGCGGCCGCGCGCACCTGCGGGTGCTGACCGATACGCAGGCGCTGCGCATAGTCTTCGAAGGCCGCCATGCGGCCGGGGTGCAGTTGGTGCGCGAGGGCCACCGGCAACTGCTGCGCGCGCGCCGCGACGTGATCGTGTGCGCGGGCACGTTCGGCTCGCCGCAGCTGCTGATGGTGTCGGGCGTCGGCCCCGCCGCGCACTTGCG  >  AM260479/245160‑245409

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