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,7810TC98.0% 174.9 / ‑4.6 50D45D (GAT→GACpyrEorotate phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (33/16);  minor base A (0/1);  total (33/17)
Fisher's exact test for biased strand distribution p-value = 3.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

ACGCTGAAATCCGGGCGCAAAAGCCCCTATTTCTTCAACGCCGGGCTGTTTAATACCGGGCGCGATCTGGCACTGTTAGGCCGTTTTTACGCTGAAGCGTTGGTGGATTCCGGCATTGAGTTCGATCTGCTG  >  AP009048/3824716‑3824847
                                                                 |                                                                  
aCGCTGAAATCCGGGCGCAAAAGCCCCTATTTCTTCAACGCCGGGCTGTTTAATACCGGGCGCGACCTgg                                                                >  1:2770270/1‑70 (MQ=255)
   cTGAAATCCGGGCGCAAAAGCCCCTATTTCTTCAACGCCGGGCTGTTTAATACCGGGCGCGACCTGGCACt                                                            >  1:266599/1‑71 (MQ=255)
      aaaTCCGGGCGCAAAAGCCCCTATTTCTTCAACGCCGGGCTGTTTAATACCGGGAGCGACCTGGCACTGtt                                                         <  1:805537/71‑1 (MQ=255)
        aTCCGGGCGCAAAAGCCCCTATTTCTTCAACGCCGGGCTGTTTAATACCGGGCGCGACCTGGCACTGt                                                          <  1:411533/68‑1 (MQ=255)
        aTCCGGGCGCAAAAGCCCCTATTTCTTCAACGCCGGGCTGTTTAATACCGGGCGCGACCTGGCACTGTTAg                                                       >  1:377742/1‑71 (MQ=255)
         tCCGGGCGCAAAAGCCCCTATTTCTTCAACGCCGGGCTGTTTAATACCGGGCGCGACCTGGCACTGTTAgg                                                      >  1:4470076/1‑71 (MQ=255)
                  aaaaGCCCCTATTTCTTCAACGCCGGGCTGTTTAATACCGGGCGCGACCTGGCACTGTTAGGCCGt                                                  >  1:726458/1‑66 (MQ=255)
                    aaGCCCCTATTTCTTCAACGCCGGGCTGTTTAATACCGGGCGCGACCTGGCACTGTTAGGCCGtttt                                               <  1:3802045/67‑1 (MQ=255)
                    aaGCCCCTATTTCTTCAACGCCGGGCTGTTTAATACCGGGCGCGACCTGGCACTGTTAGGCCGTTTTTAc                                            >  1:2716165/1‑70 (MQ=255)
                    aaGCCCCTATTTCTTCAACGCCGGGCTGTTTAATACCGGGCGCGACCTGGCACTGTTAGGCCGTTTTTACg                                           >  1:1729139/1‑71 (MQ=255)
                     aGCCCCTATTTCTTCAACGCCGGGCTGTTTAATACCGGGCGCGACCTGGCACTg                                                           >  1:3010954/1‑54 (MQ=255)
                     aGCCCCTATTTCTTCAACGCCGGGCTGTTTAATACCGGGCGCGACCTGGCACTGTTAGGCCGTTTTTACGc                                          >  1:1436412/1‑71 (MQ=255)
                      gCCCCTATTTCTTCAACGCCGGGCTGTTTAATACCGGGCGCGACCTGGCACTGTTAGGCCGTTTTTACGCt                                         <  1:3266873/71‑1 (MQ=255)
                       ccccTATTTCTTCAACGCCGGGCTGTTTAATACCGGGCGCGACCTGGCACTGTTAGGCCGTTTTTACGCTg                                        <  1:2588312/71‑1 (MQ=255)
                       ccccTATTTCTTCAACGCCGGGCTGTTTAATACCGGGCGCGACCTGGCACTGTTAGGCCGTTTTTACGCTg                                        <  1:1666537/71‑1 (MQ=255)
                            aTTTCTTCAACGCCGGGCTGTTTAATACCGGGCGCGACCTGGCACTGTTAGGCCGtt                                                 <  1:3091536/57‑1 (MQ=255)
                            aTTTCTTCAACGCCGGGCTGTTTAATACCGGGCGCGACCTGGCACTGTTAGGCCGTTTTTACGCTGAAg                                     >  1:2481489/1‑69 (MQ=255)
                             tttCTTCAACGCCGGGCTGTTTAATACCGGGCGCGACCTGGCACTGTTAGGcc                                                    <  1:692942/53‑1 (MQ=255)
                             tttCTTCAACGCCGGGCTGTTTAATACCGGGCGCGACCTGGCACTGTTAGGCCGTTTTTACGCTGAAGcat                                  >  1:2384411/1‑69 (MQ=255)
                             tttCTTCAACGCCGGGCTGTTTAATACCGGGCGCGACCTGGCACTGTTAGGCCGTTTTTACGCTGAAGcat                                  >  1:2145598/1‑69 (MQ=255)
                             tttCTTCAACGCCGGGCTGTTTAATACCGGGCGCGACCTGGCACTGTTAGGCCGTTTTTACGCTGAAGcat                                  >  1:3692513/1‑69 (MQ=255)
                             tttCTTAAACGCCGGGCTGTTTAATACCGGGCGCGACCTGGCACTGTTAGGCCGTTTTTACGCTGAAGcat                                  >  1:661721/1‑69 (MQ=255)
                              ttcttcAACGCCGGGCTGTTTAATACCGGGCGCGACCTGGCACTGTTAGGCCGTTTTTACGCTGAAGCAtt                                 >  1:4243934/1‑71 (MQ=255)
                               tcttcAACGCCTGGCTGTTTAATACCGGGCGCGACCTGGCACTGTTAGGCCGt                                                  >  1:2876917/1‑53 (MQ=39)
                               tcttcAACGCCGGGCTGTTTAATACCGGGCGCGACCTGGCACTGTTAgg                                                      >  1:3900927/1‑49 (MQ=255)
                               tcttcAACGCCGGGCTGTTTAATACCGGGCGCGACCTGGCACTGTTAgg                                                      >  1:969908/1‑49 (MQ=255)
                                 ttcAACGCCGGGCTGTTTAATACCGGGCGCGACCTGGCACTGTTAGGCCGt                                                  <  1:4653301/51‑1 (MQ=255)
                                   cAACGCCGGGCTGTTTAATACCGGGCGCGACCTGGCACTGTTAGGCCGTTTTTACGCTGAAGca                                   >  1:4458785/1‑63 (MQ=255)
                                       gCCGGGCTGTTTAATACCGGGCGCGACCTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTc                        <  1:1859859/71‑1 (MQ=255)
                                       gCCGGGCTGTTTAATACCGGGCGCGACCTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTc                        >  1:178636/1‑71 (MQ=255)
                                        ccGGGCTGTTTAATACCGGGCGCGACCTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTct                       >  1:4940603/1‑70 (MQ=255)
                                            gCTGTTTAATACCGGGCGCGACCTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGc                    <  1:3154869/70‑1 (MQ=255)
                                                tttAATACCGGGCGCGACCTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCAtt                 <  1:3744749/69‑1 (MQ=255)
                                                  tAATACCGGGCGCGACCTGGCACTGTTAGGCCGtttt                                               >  1:4330220/1‑37 (MQ=255)
                                                  tAATACCGGGCGCGACCTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAg              >  1:1688792/1‑70 (MQ=255)
                                                  tAATACCGGGCGCGACCTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAg              <  1:3847495/70‑1 (MQ=255)
                                                  tAATACCGGGCGCGACCTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAg              >  1:5176121/1‑70 (MQ=255)
                                                  tAATACCGGGCGCGACCTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGt             >  1:311995/1‑71 (MQ=255)
                                                  tAATACCGGGCGCGACCTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGt             >  1:3877081/1‑71 (MQ=255)
                                                  tAATACCGGGCGCGACCTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGt             >  1:882993/1‑71 (MQ=255)
                                                  tAATACCGGGCGCGACCTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGt             >  1:4968628/1‑71 (MQ=255)
                                                   aaTACCGGGCGCGACCTGGCACTGTTAGGCCGGTTTTACGCTGAAGCATTGGTGGATTc                        >  1:2212865/1‑59 (MQ=255)
                                                     tACCGGGCGCGACCTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCg          <  1:3629610/71‑1 (MQ=255)
                                                     tACCGGGCGCGAACTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCg          <  1:3919408/71‑1 (MQ=255)
                                                        cGGGCGCGACCTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTct                       <  1:4345755/55‑2 (MQ=255)
                                                         gggCGCGACCTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAg              >  1:1735975/1‑63 (MQ=255)
                                                           gcgcGACCTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGATctgc    >  1:383580/1‑71 (MQ=255)
                                                           gcgcGACCTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGATctg     >  1:1015085/1‑70 (MQ=255)
                                                            cgcgACCTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGATctgct   >  1:2785846/1‑71 (MQ=255)
                                                             gcgACCTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTCGATTCTGGCATTGAGTTCGATctgctg  <  1:4800625/71‑1 (MQ=255)
                                                                 |                                                                  
ACGCTGAAATCCGGGCGCAAAAGCCCCTATTTCTTCAACGCCGGGCTGTTTAATACCGGGCGCGATCTGGCACTGTTAGGCCGTTTTTACGCTGAAGCGTTGGTGGATTCCGGCATTGAGTTCGATCTGCTG  >  AP009048/3824716‑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: