Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 3,824,814 G→A 100% A56A (GCG→GCA pyrE → orotate phosphoribosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0090483,824,8140GA100.0% 228.8 / ‑6.7 78A56A (GCG→GCApyrEorotate phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (44/33);  minor base T (0/1);  total (44/34)
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.

CTTCAACGCCGGGCTGTTTAATACCGGGCGCGATCTGGCACTGTTAGGCCGTTTTTACGCTGAAGCGTTGGTGGATTCCGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGC  >  AP009048/3824748‑3824879
                                                                  |                                                                 
cttcAACGCCGGGCTGTTTAATACCGGGCGCGACCTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATtgg                                                               >  1:3666996/1‑71 (MQ=255)
   cAACGCCGGGCTGTTTAATACCGGGCGCGACCTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATtgg                                                               <  1:2297095/68‑1 (MQ=255)
     aCGCCGGGCTGTTTAATACCGGGCGCGACCTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGAt                                                          >  1:4027081/1‑71 (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)
              tGTTTAATACCGGGCGCGACCTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCAtt                                                 >  1:3305842/1‑71 (MQ=255)
                 ttAATACCGGGCGCGACCTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAg                                              >  1:94957/1‑71 (MQ=255)
                  tAATACCGGGCGCGACCTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGt                                             >  1:4374833/1‑71 (MQ=255)
                  tAATACCGGGCGCGACCTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGt                                             >  1:3285295/1‑71 (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:117839/71‑1 (MQ=255)
                       ccGGGCGCGACCTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGAt                                        >  1:4070336/1‑71 (MQ=255)
                       ccGGGCGCGACCTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGAt                                        >  1:2158296/1‑71 (MQ=255)
                       ccGGGCGCGACCTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGAt                                        >  1:4139143/1‑71 (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)
                                 ccTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGATCTGCTGTTTg                              >  1:1984867/2‑71 (MQ=255)
                                 ccTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGATCTGCTGTTTg                              >  1:3240178/2‑71 (MQ=255)
                                 ccTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGATCTGCTGTTTg                              >  1:3858881/2‑71 (MQ=255)
                                   tGGCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGATCTGCTGtt                                >  1:3754040/1‑67 (MQ=255)
                                    ggCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGATCTGCTGTTTgg                             >  1:1094897/1‑69 (MQ=255)
                                    ggCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCAgg                                                     <  1:324476/45‑1 (MQ=38)
                                     gCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGATCTGCTGt                                 >  1:2279917/1‑64 (MQ=255)
                                     gCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGATCTGCTGTTTGGccc                          >  1:1005310/1‑71 (MQ=255)
                                     gCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGATCTGCTGTTTGGc                            <  1:822305/69‑1 (MQ=255)
                                      cACTGTTAGGCCGTTTTTACGCTGAAGCTTTGGTGGATTCTGGCATTGAGTTCGATCTGCTGTTTGGCCCt                         <  1:2045673/71‑1 (MQ=255)
                                       aCTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTc                                                        <  1:3119344/39‑1 (MQ=255)
                                       aCTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGATCTGCTGttt                               >  1:2542665/1‑64 (MQ=255)
                                       aCTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGATCTGCTGTTTGGCCCTg                        >  1:2469317/1‑71 (MQ=255)
                                       aCTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGCCATTGAGTTCGATCTGCTGTTTGGCCCTg                        <  1:2698616/71‑1 (MQ=255)
                                        cTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCAtt                                                 <  1:4439147/45‑1 (MQ=38)
                                         tGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGATCTGCTGTTTGGcc                           >  1:3743949/1‑66 (MQ=255)
                                         tGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCt                      >  1:4525788/1‑71 (MQ=255)
                                         tGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCt                      >  1:1626164/1‑71 (MQ=255)
                                         tGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCt                      >  1:3912772/1‑71 (MQ=255)
                                          gTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGATCTGCTGt                                 <  1:3533710/59‑1 (MQ=255)
                                           ttAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGAt                                        >  1:1678535/1‑51 (MQ=255)
                                           ttAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGATCTGCTGTTTGGCCCTGc                       >  1:3627518/1‑68 (MQ=255)
                                           atAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTgg                                                     >  1:4006884/2‑38 (MQ=37)
                                            tAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCAtt                                                 >  1:408364/1‑41 (MQ=38)
                                            tAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTa                    >  1:3367309/1‑70 (MQ=255)
                                            tAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTa                    >  1:2458549/1‑70 (MQ=255)
                                             aGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCg                                          <  1:4322680/47‑1 (MQ=39)
                                             aGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTAc                   >  1:1565679/1‑70 (MQ=255)
                                              ggCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACa                  <  1:4333436/70‑1 (MQ=255)
                                               gCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGt                                             <  1:1034522/42‑1 (MQ=38)
                                               gCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACaaa                <  1:2424188/71‑1 (MQ=255)
                                                 cGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAgg              >  1:14063/1‑71 (MQ=255)
                                                  gTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGAt                                        >  1:1697544/1‑44 (MQ=38)
                                                  gTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAggg             <  1:1536608/71‑1 (MQ=255)
                                                  gTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAggg             <  1:343236/71‑1 (MQ=255)
                                                    ttttACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGATCTGCTGTTTGGCCCTGc                       <  1:941929/59‑1 (MQ=255)
                                                     tttACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAg               >  1:2174700/1‑66 (MQ=255)
                                                     tttACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGAt           <  1:4330933/70‑1 (MQ=255)
                                                       tACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATcc         <  1:2487395/70‑1 (MQ=255)
                                                        aCGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGATCTGCTGTTTg                              >  1:941734/1‑48 (MQ=39)
                                                        aCGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATccc        <  1:3647606/70‑1 (MQ=255)
                                                        aCGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATccc        >  1:3612353/1‑70 (MQ=255)
                                                         cGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGTATCCCGa      <  1:2673844/71‑1 (MQ=255)
                                                         cGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGa      <  1:535655/71‑1 (MQ=255)
                                                          gCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGa      <  1:1258216/70‑1 (MQ=255)
                                                          gCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGAt     <  1:2179273/71‑1 (MQ=255)
                                                            tGAAGCATTGGTGGATTCTGGCATTGAGTTCGATCTGCTGTTTgg                             <  1:1583131/45‑1 (MQ=38)
                                                              aaGCATTGGTGGATTCTGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGc  >  1:1197125/1‑70 (MQ=255)
                                                              aaGCATTGGTGGATTCTGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGc  <  1:4045436/70‑1 (MQ=255)
                                                              aaGCATTGGTGGATTCTGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGc  <  1:1900282/70‑1 (MQ=255)
                                                                  |                                                                 
CTTCAACGCCGGGCTGTTTAATACCGGGCGCGATCTGGCACTGTTAGGCCGTTTTTACGCTGAAGCGTTGGTGGATTCCGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGC  >  AP009048/3824748‑3824879

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: