Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 3,824,781 T→C 100% D45D (GAT→GAC pyrE → orotate phosphoribosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0090483,824,7810TC100.0% 183.1 / NA 52D45D (GAT→GACpyrEorotate phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (33/19);  total (33/19)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TTACGCTGAAATCCGGGCGCAAAAGCCCCTATTTCTTCAACGCCGGGCTGTTTAATACCGGGCGCGATCTGGCACTGTTAGGCCGTTTTTACGCTGAAGCGTTGGTGGATTCCGGCATTGAGTTCGATCTGCTG  >  AP009048/3824714‑3824847
                                                                   |                                                                  
ttACGCTGAAATCCGGGCGCAAAAGCCCCTATTTCTTCAACGCCGGGCTGTTTAATACCGGGCGCGACCTg                                                                 >  1:4468875/1‑71 (MQ=255)
  aCGCTGAAATCCGGGCGCAAAAGCCCCTATTTCTTCAACGCCGGGCTGTTTAATACCGGGCGCGACCTgg                                                                >  1:816967/1‑70 (MQ=255)
   cGCTGAAATCCGGGCGCAAAAGCCCCTATTTCTTCAACGCCGGGCTGTTTAATACCGGGCGCGACCTGGc                                                               <  1:3553845/70‑1 (MQ=255)
   cGCTGAAATCCGGGCGCAAAAGCCCCTATTTCTTCAACGCCGGGCTGTTTAATACCGGGCGCGACCTGGCa                                                              <  1:2528581/71‑1 (MQ=255)
        aaaTCCGGGCGCAAAAGCCCCTATTTCTTCAACGCCGGGCTGTTTAATACCGGGCGCGACCTGGCACTg                                                           <  1:2753503/69‑1 (MQ=255)
        aaaTCCGGGCGCAAAAGCCCCTATTTCTTCAACGCCGGGCTGTTTAATACCGGGCGCGACCTGGCACTGtt                                                         <  1:232781/71‑1 (MQ=255)
           tCCGGGCGCAAAAGCCCCTATTTCTTCAACGCCGGGCTGTTTAATACCGGGCGCGACCt                                                                  >  1:463653/1‑59 (MQ=255)
                     aaaGCCCCTATTTCTTCAACGCCGGGCTGTTTAATACCGGGCGCGACCTGGCACTGTTAGGCCGTTTTTAc                                            <  1:1256218/71‑1 (MQ=255)
                     aaaGCCCCTATTTCTTCAACGCCGGGCTGTTTAATACCGGGCGCGACCTGGCACTGTTAGGCCGTTTTTAc                                            >  1:2491670/1‑71 (MQ=255)
                      aaGCCCCTATTTCTTCAACGCCGGGCTGTTTAATACCGGGCGCGACCTgg                                                                >  1:1151334/1‑50 (MQ=255)
                          cccTATTTCTTCAACGCCGGGCTGTTTAATACCGGGCGCGACCTGGCACTGTTAGGCCGtttt                                               <  1:2834804/63‑1 (MQ=255)
                            cTATTTCTTCAACGCCGGGCTGTTTAATACCGGGCGCGACCTGGCACTGTTAGGCCGTTTTTa                                             >  1:2985376/1‑63 (MQ=255)
                            cTATTTCTTCAACGCCGGGCTGTTTAATACCGGGCGCGACCTGGCACTGTTAGGCCGTTTTTACGCTTAAg                                     >  1:3955612/1‑71 (MQ=255)
                            cTATTTCTTCAACGCCGGGCTGTTTAATACCGGGCGCGACCTGGCACTGTTAGGCCGTTTTTACGCTGAAg                                     >  1:180854/1‑71 (MQ=255)
                              aTTTCTTCAACGCCGGGCTGTTTAATACCGGGCGCGACCTGGCACTGTTAGGCCGTTTTTACGCTGAAGca                                   >  1:3872667/1‑70 (MQ=255)
                              aTTTCTTCAACGCCGGGCTGTTTAATACCGGGCGCGACCTGGCACTGTTAGGCCGTTTTTACGCTGAAGca                                   >  1:2749721/1‑70 (MQ=255)
                              aTTTCTTCAACGCCGGGCTGTTTAATACCGGGCGCGACCTGGCACTGTTAGGCCGTTTTTACGCTGAAGca                                   >  1:2986863/1‑70 (MQ=255)
                                  cttcAACGCCGGGCTGTTTAATACCGGGCGCGACCTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATtgg                               >  1:3666996/1‑71 (MQ=255)
                                   ttcCACGCCGGGCTGTTTAATACCGGGCGCGACCTGGCACTGTTAGGCCGTTTTTACGCTGAAGcat                                  >  1:262544/1‑65 (MQ=255)
                                     cAACGCCGGGCTGTTTAATACCGGGCGCGACCTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATtgg                               <  1:2297095/68‑1 (MQ=255)
                                      aaCGCCGGGCTGTTTAATACCGGGCGCGACCTGGCACTGt                                                          >  1:1222344/1‑40 (MQ=255)
                                       aCGCCGGGCTGTTTAATACCGGGCGCGACCTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGAt                          >  1:4027081/1‑71 (MQ=255)
                                         gCCGGGCTGTTTAATACCGGGCGCGACCTGGCACTGtt                                                         >  1:4608994/1‑38 (MQ=255)
                                         gCCGGGCTGTTTAATACCGGGCGCGACCTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTc                        >  1:705966/1‑71 (MQ=255)
                                         gCCGGGCTGTTTAATACCGGGCGCGACCTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTc                        <  1:3542914/71‑1 (MQ=255)
                                              gCTGTTTAATACCGGGCGCGACCTGGCACTGTTATGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCa                   <  1:1043450/71‑1 (MQ=255)
                                              gCTGTTTAATACCGGGCGCGACCTGGCACTGTTAGTCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCa                   <  1:2895270/71‑1 (MQ=255)
                                               cTGTTTAATACCGGGCGCGACCTGGCACTGTTAGGCCGTTTTTACGCTGAAg                                     >  1:4426887/1‑52 (MQ=255)
                                                tGTTTAATACCGGGCGCGACCTGGCACTGTTAGGCCGTTTTTACGCTGAAGcat                                  <  1:2383614/54‑3 (MQ=255)
                                                tGTTTAATACCGGGCGCGACCTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCAtt                 >  1:3305842/1‑71 (MQ=255)
                                                  tttAATACCGGGCGCGACCTGGCACTGTTAGGCCGTTTTTACGCTGAAGca                                   <  1:3158301/51‑2 (MQ=255)
                                                   ttAATACCGGGCGCGACCTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAg              >  1:94957/1‑71 (MQ=255)
                                                    tAATACCGGGCGCGACCTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGt             >  1:3285295/1‑71 (MQ=255)
                                                    tAATACCGGGCGCGACCTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGt             >  1:4374833/1‑71 (MQ=255)
                                                     aaTACCGGGCGCGACCTGGCACTGTTAGGCCGTTTTTa                                             <  1:2217807/38‑1 (MQ=255)
                                                      aTACCGGGCGCGACCTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGtt            <  1:4281109/70‑1 (MQ=255)
                                                      aTACCGGGCGCGACCTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTc           >  1:206068/1‑71 (MQ=255)
                                                       tACCGGGCGCGACCTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGa               >  1:4438409/1‑66 (MQ=255)
                                                         ccGGGCGCGACCTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATtgg                               <  1:1567726/48‑1 (MQ=39)
                                                         ccGGGCGCGACCTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGa         >  1:2788190/1‑70 (MQ=255)
                                                         ccGGGCGCGACCTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGa         >  1:2782653/1‑70 (MQ=255)
                                                         ccGGGCGCGACCTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGAt        >  1:2158296/1‑71 (MQ=255)
                                                         ccGGGCGCGACCTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGAt        >  1:4070336/1‑71 (MQ=255)
                                                         ccGGGCGCGACCTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGAt        >  1:4139143/1‑71 (MQ=255)
                                                         ccGGGCGCGACCTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGAt        <  1:117839/71‑1 (MQ=255)
                                                         ccGGGCGCGACCTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGAt        <  1:3404113/71‑1 (MQ=255)
                                                            ggCGCGACCTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGATctg     >  1:2845918/1‑71 (MQ=255)
                                                            ggCGCGACCTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGATc       >  1:4364819/1‑69 (MQ=255)
                                                             gcgcGACCTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGAt        <  1:457459/67‑1 (MQ=255)
                                                              cgcgACCTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGa                           <  1:1256603/47‑1 (MQ=39)
                                                              cgcgACCTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGATctgc    >  1:1499247/1‑70 (MQ=255)
                                                               gcgACCTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGATctgctg  >  1:3535733/1‑71 (MQ=255)
                                                                   |                                                                  
TTACGCTGAAATCCGGGCGCAAAAGCCCCTATTTCTTCAACGCCGGGCTGTTTAATACCGGGCGCGATCTGGCACTGTTAGGCCGTTTTTACGCTGAAGCGTTGGTGGATTCCGGCATTGAGTTCGATCTGCTG  >  AP009048/3824714‑3824847

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