Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 3,825,072 C→T 100% A142A (GCC→GCT pyrE → orotate phosphoribosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0090483,825,0720CT100.0% 234.3 / NA 78A142A (GCC→GCTpyrEorotate phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (33/45);  total (33/45)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GTAATGCTGGTAGATGATGTGATCACCGCCGGAACGGCGATTCGCGAGTCGATGGAGATTATTCAGGCCAATGGCGCGACGCTTGCTGGCGTGTTGATTTCGCTCGATCGTCAGGAACGCGGGCGCGGCGAGATTT  >  AP009048/3825004‑3825139
                                                                    |                                                                   
gTAATGCTGGTAGATGATGTGATCACCGCCGGAACGGCGATTCGCGAATCGATGGAGATTATTCAGGCTaa                                                                   >  1:451702/1‑71 (MQ=255)
gTAATGCTGGTAGATGATGTGATCACCGCCGGAACGGCGATTCGCGAATCGATGGAGATTATTCAGGCTaa                                                                   <  1:291911/71‑1 (MQ=255)
 tAATGCTGGTAGATGATGTGATCACCGCCGGAACGGCGATTCGCGAATCGATGGAGATTATTCAGGCTaa                                                                   <  1:2912019/70‑1 (MQ=255)
  aaTGCTGGTAGATGATGTGATCACCGCCGGAACGGCGATTCGCGAATCGATGGAGATTATTCAGGCTAATg                                                                 >  1:406701/1‑71 (MQ=255)
    tGCTGGTAGATGATGTGATCACCGCCGGAACGGCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGgc                                                               <  1:1275995/71‑1 (MQ=255)
         gTAGATGATGTGATCACCGCCGGAACGGCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGAc                                                          >  1:2804411/1‑71 (MQ=255)
         gTAGATGATGTGATCACCGCCGGAACGGCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGAc                                                          <  1:3705377/71‑1 (MQ=255)
          tAGATGATGTGATCACCGCCGGAACGGCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGAc                                                          >  1:1504258/1‑70 (MQ=255)
              tgatgTGATCACCGCCGGAACGGCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGcgt                                                      >  1:1320742/1‑68 (MQ=255)
              tgatgTGATCACCGCCGGAACGGCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTg                                                     <  1:549396/71‑1 (MQ=255)
                 tgtgATCACCGCCGGAACGGCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGcgcg                                                            <  1:2105981/61‑1 (MQ=255)
                 tgtgATCACCGCCGGAACGGCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTg                                                  <  1:2417389/71‑1 (MQ=255)
                 tgtgATCACCGCCGGAACGGCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTg                                                  >  1:3300849/1‑71 (MQ=255)
                  gtgATCACCGCCGGAACGGCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTg                                                     <  1:885950/67‑1 (MQ=255)
                   tgATCACCGCCGGAACGGCGATTCGCGAATCGATGGAGATTATTCAGGCTaa                                                                   <  1:3738126/52‑1 (MQ=39)
                   tgATCACCGCCGGAACGGCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGATGGc                                                >  1:4149700/1‑71 (MQ=255)
                      tCACCGCCGGAACGGCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCgta                                             >  1:1649644/1‑70 (MQ=255)
                      tCACCGCCGGAACGGCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGATGGCgta                                             <  1:1124544/71‑2 (MQ=255)
                       cACCGCCGGAACGGCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCgtat                                            >  1:1586172/1‑69 (MQ=255)
                        aCCGCCGGAACGGCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTAtt                                           >  1:3431158/1‑71 (MQ=255)
                          cgccgGAACGGCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGa                                         >  1:212467/1‑71 (MQ=255)
                            ccgGAACGGCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGCCGCTTGCTgg                                                 >  1:2632190/1‑61 (MQ=255)
                              ggAACGGCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTc                                     <  1:2888825/71‑1 (MQ=255)
                                aaCGGCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGc                                   <  1:1641416/71‑1 (MQ=255)
                                aaCGGCGATTCGCGAATCGATGGAGATTATTCAGGATAATGGCGCGa                                                           >  1:2647512/1‑47 (MQ=37)
                                 aCGGCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCt                                                   >  1:1367202/1‑54 (MQ=39)
                                 aCGGCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCt                                  >  1:3202809/1‑71 (MQ=255)
                                  cGGCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCtt                                 <  1:4522355/71‑2 (MQ=255)
                                  cGGCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCtt                                 >  1:1341393/1‑70 (MQ=255)
                                  cGGCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCtt                                 <  1:4295501/71‑2 (MQ=255)
                                   ggCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGc                                   <  1:162912/68‑1 (MQ=255)
                                   ggCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCtt                                 >  1:1754405/1‑69 (MQ=255)
                                    gCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCt                                  >  1:646567/1‑68 (MQ=255)
                                    gCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGa                               >  1:3304369/1‑71 (MQ=255)
                                    gCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGa                               >  1:2788483/1‑71 (MQ=255)
                                     cGATTCGCGAATCGATGGAGATTATTCAGGCTAATgg                                                                <  1:1714816/37‑1 (MQ=37)
                                     cGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCt                                  >  1:3030462/1‑67 (MQ=255)
                                     cGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGa                               >  1:4403973/1‑70 (MQ=255)
                                      gATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATc                             >  1:2647793/1‑71 (MQ=255)
                                        ttCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTg                                                     <  1:1324806/45‑1 (MQ=38)
                                        ttCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGt                           <  1:3742535/71‑1 (MQ=255)
                                        ttCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGt                           <  1:2687081/71‑1 (MQ=255)
                                        ttCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGt                           <  1:517736/71‑1 (MQ=255)
                                         tCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGt                           >  1:4349746/1‑70 (MQ=255)
                                          cgcgAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTc                                     >  1:866565/1‑59 (MQ=255)
                                          cgcgAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCg                            <  1:1977238/68‑1 (MQ=255)
                                          cgcgAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTc                          <  1:2946924/70‑1 (MQ=255)
                                             gaatCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAgg                       >  1:4062546/4‑70 (MQ=255)
                                             gaatCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGa                      <  1:4452387/68‑1 (MQ=255)
                                               atCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAAc                    <  1:1278153/70‑1 (MQ=255)
                                                tCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGa                                         >  1:2281312/1‑49 (MQ=38)
                                                tCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAAcg                   <  1:4420731/71‑1 (MQ=255)
                                                 cGATTGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAAc                    <  1:715837/69‑1 (MQ=255)
                                                  gATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAAcgg                  <  1:4599494/70‑2 (MQ=255)
                                                  gATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAAcgg                  <  1:629530/70‑2 (MQ=255)
                                                   aTGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTAtt                                           <  1:3686640/44‑1 (MQ=38)
                                                     ggAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGgc              >  1:3189050/1‑71 (MQ=255)
                                                      gagaTTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTTATCGTCAGGAACGGGGgcg             <  1:32197/71‑1 (MQ=255)
                                                      gagaTTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGgc              <  1:927143/70‑1 (MQ=255)
                                                       agaTTATTCAGGCTAATGGCGCGACGCTTGCTGGCg                                               <  1:4340365/36‑1 (MQ=255)
                                                       agaTTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGgcgc            >  1:1673456/1‑71 (MQ=255)
                                                       agaTTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGgcgc            <  1:3545027/71‑1 (MQ=255)
                                                        gaTTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCttg                                >  1:1876353/1‑48 (MQ=255)
                                                        gaTTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACgggg                >  1:1208203/1‑66 (MQ=255)
                                                         attattCAGGCTAATGGCGCGACGCTTGCTGGCgta                                             >  1:3690063/1‑35 (MQ=38)
                                                          ttattCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGt                           <  1:3819662/53‑1 (MQ=255)
                                                          ttattCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGt                           <  1:2147224/53‑1 (MQ=255)
                                                          ttattCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCgcggc         <  1:2763703/71‑1 (MQ=255)
                                                           tattCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGa                                         >  1:1484997/1‑38 (MQ=37)
                                                           tattCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCg                            <  1:425581/51‑1 (MQ=255)
                                                             ttCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCgag      <  1:1555226/71‑1 (MQ=255)
                                                               cAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCa                         <  1:4210280/50‑1 (MQ=255)
                                                                aGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCg                                    <  1:1783954/38‑1 (MQ=37)
                                                                aGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGAtt   <  1:1617157/71‑1 (MQ=255)
                                                                 ggCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGAttt  <  1:4522319/71‑1 (MQ=255)
                                                                 ggCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGAttt  <  1:3627377/71‑1 (MQ=255)
                                                                 ggCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTCAGGAACGGGGGCGCGGCGAGAttt  <  1:3087504/71‑1 (MQ=255)
                                                                 ggCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGGCAGGATCGGGGGCGCGGCGAGAt    <  1:620104/69‑1 (MQ=255)
                                                                    |                                                                   
GTAATGCTGGTAGATGATGTGATCACCGCCGGAACGGCGATTCGCGAGTCGATGGAGATTATTCAGGCCAATGGCGCGACGCTTGCTGGCGTGTTGATTTCGCTCGATCGTCAGGAACGCGGGCGCGGCGAGATTT  >  AP009048/3825004‑3825139

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: