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% 177.7 / NA 50D45D (GAT→GACpyrEorotate phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (28/22);  total (28/22)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GCTGAAATCCGGGCGCAAAAGCCCCTATTTCTTCAACGCCGGGCTGTTTAATACCGGGCGCGATCTGGCACTGTTAGGCCGTTTTTACGCTGAAGCGTTGGTGGATTCCGGCATTGAGTTCGATCTGCTG  >  AP009048/3824718‑3824847
                                                               |                                                                  
gCTGAAATCCGGGCGCAAAAGCCCCTATTTCTTCAACGCCGGGCTGTTTAATACCGGGCGCGACCTGGCAc                                                             >  1:2853046/1‑71 (MQ=255)
 cTGAAATCCGGGCGCAAAAGCCCCTATTTCTTCAACGCCGGGCTGTTTAATACCGGGCGCGACCTGGCACt                                                            <  1:3532907/71‑1 (MQ=255)
    aaaTCCGGGCGCAAAAGCCCCTATTTCTTCAACGCCGGGCTGTTTAATACCGGGCGCGACCTGGCACTg                                                           <  1:2766199/69‑1 (MQ=255)
    aaaTCCGGGCGCAAAAGCCCCTATTTCTTCAACGCCGGGCTGTTTAATACCGGGCGCGACCTGGCACTGt                                                          <  1:248542/70‑1 (MQ=255)
             cgcAAAAGCCCCTATTTCTTCAACGCCGGGCTGTTTAATACCGGGCGCGACCTGGCACTGTTAGGCCGttt                                                <  1:1626794/71‑1 (MQ=255)
                aaaaGCCCCTATTTCTTCAACGCCGGGCTGTTTAATACCGGGCGCGACCTGGc                                                               >  1:2508111/1‑53 (MQ=255)
                 aaaGCCCCTATTTCTTCAACGCCGGGCTGTTTAATACCGGGCGCGACCTGGCACTGtt                                                         <  1:273040/58‑1 (MQ=255)
                  aaGCCCCTATTTCTTCAACGCAGGGCTGTTTAATACCGGGCGCGACCTGGCACTGTTAGGCCGTTTTTACg                                           >  1:44868/1‑71 (MQ=255)
                      cccTATTTCTTCAACGCCGGGCTGTTTAATACCGGGCGCGACCTGGCACTGTTa                                                        <  1:3501968/54‑1 (MQ=255)
                        cTATTTCTTCAACGCCGGGCTGTTTAATACCGGGCGCGACCTGGCACTGTTAGGCCGTTTTTACGCTGAAg                                     >  1:2978393/1‑71 (MQ=255)
                         tATTTCTTCAACGCCGGGCTGTTTAATACCGGGCGCGACCTGGCACTGTTAGGCCGTTTTTACGCTGAAg                                     >  1:193811/1‑70 (MQ=255)
                          aTTTCTTCAACGCCGGGCTGTTTAATACCGGGCGCGACCTGGCa                                                              >  1:3778467/1‑44 (MQ=255)
                          aTTTCTTCAACGCCGGGCTGTTTAATACCGGGCGCGACCTGGCACTGTTAGGCCGTTTTTACGCTGAAGca                                   >  1:229556/1‑70 (MQ=255)
                           tttCTTCAACGCCGGGCTGTTTAATACCGGGCGCGACCTGGCACTGTTAGGCCGTTTTTACGCTGAAGcat                                  >  1:1927876/1‑69 (MQ=255)
                            ttcttcAACGCCGGGCTGTTTAATACCGGGCGCGACCTGGCACTGTTAGGCCGTTTTTa                                             <  1:393783/59‑1 (MQ=255)
                            ttcttcAACGCCGGGCTGTTTAATACCGGGCGCGACCTGGCACTGTTAGGCCGTTTTTACGCt                                         >  1:3157166/1‑63 (MQ=255)
                             tcttcAACGCCGGGCTGTTTAATACCGGGCGCGACCTGGCACTGTTa                                                        >  1:2377689/1‑47 (MQ=255)
                             tcttcAACGCCGGGCTGTTTAATACCGGGCGCGACCTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATtg                                >  1:140208/1‑71 (MQ=255)
                              cttcAACGCCGGGCTGTTTAATACCGGGCGCGACCTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATtgg                               >  1:17139/1‑71 (MQ=255)
                              cttcAACGCCGGGCTGTTTAATACCGGGCGCGACCTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATtg                                <  1:325807/70‑1 (MQ=255)
                                 cAACGCCGGGCTGTTTAATACCGGGCGCGACCTGGCACTGTTAgg                                                      >  1:4415041/1‑45 (MQ=255)
                                  aaCGCCGGGCTGTTTAATACCGGGCGCGACCTGGCACTGTTa                                                        >  1:4857274/1‑42 (MQ=255)
                                    cGCCGGGCTGTTTAATACCGGGCGAGACCTGGCACTGTTAGGCCGTTTTTa                                             >  1:129727/1‑51 (MQ=255)
                                       cGGGCTGTTTAATACCGGGCGCGACCTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTctg                      <  1:4334646/71‑3 (MQ=255)
                                       cGGGCTGTTTAATACCGGGCGCGACCTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTctg                      >  1:3756746/1‑69 (MQ=255)
                                        gggCTGTTTAATACCGGGCGCGACCTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGGGGATTctg                      <  1:4182921/70‑3 (MQ=255)
                                          gCTGTTTAATACCGGGCGCGACCTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCa                   <  1:4169760/71‑1 (MQ=255)
                                          gCTGTTTAATACCGGGCGCGACCTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCa                   <  1:2590442/71‑1 (MQ=255)
                                            tGTTTAATACCGGGCGCGACCTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCAtt                 >  1:1805446/1‑71 (MQ=255)
                                             gTTTAATACCGGGCGCGACCTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCAtt                 >  1:2233035/1‑70 (MQ=255)
                                              tttAATACCGGGCGCGACCTTGCACTGTTAGGCCGttttt                                              <  1:3520864/40‑1 (MQ=37)
                                              tttAATACCGGGCGCGACCTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGa                           <  1:2857493/59‑1 (MQ=255)
                                              tttAATACCGGGCGCGACCTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCAtt                 >  1:192887/1‑69 (MQ=255)
                                               ttAATACCGGGCGCGACCTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATtg                                <  1:1903015/53‑1 (MQ=255)
                                               ttAATACCGGGCGCGACCTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAg              >  1:2538604/1‑71 (MQ=255)
                                                tAATACCGGGCGCGACCTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAg              <  1:4300840/70‑1 (MQ=255)
                                                tAATACCGGGCGCGACCTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGt             >  1:712476/1‑71 (MQ=255)
                                                tAATACCGGGCGCGACCTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGt             >  1:4848177/1‑71 (MQ=255)
                                                 aaTACCGGGCGCGACCTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGtt            >  1:4675492/1‑71 (MQ=255)
                                                  aTACCGGGCGCGACCTGGCACTGTTAGGCCGTTTTTACGCTGAAGCAtt                                 <  1:2773460/49‑1 (MQ=39)
                                                   tACCGGGCGCGACCTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCAtt                 <  1:2431110/64‑1 (MQ=255)
                                                   tACCGGGCGCGACCTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTc           <  1:2098439/70‑1 (MQ=255)
                                                     ccGGGCGCGACCTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCAtt                 >  1:4772664/1‑62 (MQ=255)
                                                      cGGGCGCGACCTGGCACTGTTAGGCCGTTTTTACGCTGAAg                                     <  1:2334007/41‑1 (MQ=255)
                                                      cGGGCGCGACCTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTct                       <  1:1926450/55‑2 (MQ=255)
                                                       gggCGCGACCTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGtt            >  1:2985307/1‑65 (MQ=255)
                                                       gggCGCGACCTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGATct      >  1:2087091/1‑71 (MQ=255)
                                                        ggCGCGACCTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGATct      >  1:3848493/1‑70 (MQ=255)
                                                           gcgACCTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGATctgctg  >  1:4080204/1‑71 (MQ=255)
                                                            cgACCTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCAtt                 <  1:3678462/55‑1 (MQ=255)
                                                               |                                                                  
GCTGAAATCCGGGCGCAAAAGCCCCTATTTCTTCAACGCCGGGCTGTTTAATACCGGGCGCGATCTGGCACTGTTAGGCCGTTTTTACGCTGAAGCGTTGGTGGATTCCGGCATTGAGTTCGATCTGCTG  >  AP009048/3824718‑3824847

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