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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0090483,825,2040CT100.0% 223.5 / NA 74D186D (GAC→GATpyrEorotate phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (35/39);  total (35/39)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

ATTTCGGCGATTCAGGAAGTTGAGCGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGACCTGATTGCTTACCTGGAAGAGAAGCCGGAAATGGCGGAACATCTGGCGGCGGTTAAGGCCTATCGCG  >  AP009048/3825136‑3825271
                                                                    |                                                                   
aTTTCGGCGATTCAGGAAGTTGAGCGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCt                                                                   <  1:2598789/71‑1 (MQ=255)
aTTTCGGCGATTCAGGAAGTTGAGCGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCt                                                                   <  1:3585972/71‑1 (MQ=255)
 tttCGGCGATTCAGGAAGTTGAGCGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTg                                                                  <  1:2166777/71‑1 (MQ=255)
 tttCGGCGATTCAGGAAGTTGAGCGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTg                                                                  <  1:4019323/71‑1 (MQ=255)
  ttCGGCGATTCAGGAAGTTGAGCGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTg                                                                  >  1:298980/1‑70 (MQ=255)
   tCGGCGATTCAGGAAGTTGAGCGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGAt                                                                >  1:1097539/1‑71 (MQ=255)
     ggCGATTCAGGAAGTTGAGCGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTg                                                              <  1:4272197/71‑1 (MQ=255)
       cGATTCAGGAAGTTGAGCGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGcc                                                            >  1:4007931/1‑70 (MQ=255)
         aTTCAGGAAGTTGAGCGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGa                                                                 <  1:4953974/64‑1 (MQ=255)
         aTTCAGGAAGTTGAGCGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGAtt                                                               <  1:207257/66‑1 (MQ=255)
            cAGGAAGTTGAGCGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTatc                                                        >  1:3468556/1‑68 (MQ=255)
              ggAAGTTGAGCGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTgg                                                     <  1:1694183/71‑1 (MQ=255)
              ggAAGTTGAGCGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTg                                                      >  1:4322960/1‑70 (MQ=255)
                 aGTTGAGCGTGTTTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAag                                                  <  1:1135267/71‑1 (MQ=255)
                 aGTTGAGCGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCt                                                                   >  1:3544279/1‑54 (MQ=255)
                 aGTTGAGCGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCt                                                                   >  1:3922505/1‑54 (MQ=255)
                   ttGAGCGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAagag                                                <  1:2282409/71‑1 (MQ=255)
                    tGAGCGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAagag                                                >  1:2661754/1‑70 (MQ=255)
                    tGAGCGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAag                                                  >  1:2354713/1‑68 (MQ=255)
                     gAGCGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAagag                                                >  1:4377842/1‑69 (MQ=255)
                     gAGCGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAagag                                                >  1:2241810/1‑69 (MQ=255)
                        cGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGcc                                           <  1:1646192/71‑1 (MQ=255)
                        cGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGcc                                           <  1:67620/71‑1 (MQ=255)
                         gTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAgaga                                               >  1:4258595/1‑66 (MQ=255)
                             ttACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGaaa                                      <  1:2615440/71‑1 (MQ=255)
                             ttACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGaa                                       <  1:3539447/70‑1 (MQ=255)
                              tACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGa                                                                 >  1:2911706/1‑43 (MQ=255)
                              tACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAgaga                                               >  1:2594038/1‑61 (MQ=255)
                              tACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAAt                                     >  1:4200201/1‑71 (MQ=255)
                              tACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAAt                                     >  1:3902058/1‑71 (MQ=255)
                               aCAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGa                                        >  1:2874556/1‑67 (MQ=255)
                               aCAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAAt                                     <  1:1318402/70‑1 (MQ=255)
                                cAACTGCTAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATg                                    <  1:3585312/70‑1 (MQ=255)
                                cAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATgg                                   <  1:4744873/71‑1 (MQ=255)
                                cAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATgg                                   <  1:1022560/71‑1 (MQ=255)
                                cAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATgg                                   <  1:1742235/71‑1 (MQ=255)
                                 aaCTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCgg                                         <  1:1021976/64‑1 (MQ=255)
                                  aCTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCg                                 <  1:1868467/71‑1 (MQ=255)
                                   cTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCg                                 >  1:4104293/1‑70 (MQ=255)
                                    tGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTgg                                                     <  1:4575471/49‑1 (MQ=38)
                                    tGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTgg                                                     <  1:44249/49‑1 (MQ=38)
                                     gCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTat                                                         >  1:4789264/1‑43 (MQ=38)
                                      cataGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGa                               <  1:756402/66‑1 (MQ=255)
                                      cAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGa                               <  1:1369398/69‑1 (MQ=255)
                                       aaaGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGa                                        <  1:3925388/59‑1 (MQ=38)
                                       aaaGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAAc                             <  1:4320331/70‑1 (MQ=255)
                                         aGTGATCTCTATCATCACCCTGAATGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATc                          >  1:5670/1‑71 (MQ=255)
                                         aGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCTGGAAATGGCGGAACAt                           <  1:2347012/70‑1 (MQ=255)
                                           tGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACAt                           >  1:1097074/1‑68 (MQ=255)
                                           tGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATCTg                        >  1:1179415/1‑71 (MQ=255)
                                              tctcTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCCGAACAt                           <  1:12102/65‑1 (MQ=255)
                                                tctATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATCTggcgg                    >  1:3988191/1‑70 (MQ=255)
                                                 ctATCATCACCATGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACAt                           >  1:1806508/1‑62 (MQ=255)
                                                  tATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATCt                         <  1:1672203/63‑1 (MQ=255)
                                                   atcatcACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATCTGGCGGCTGt                >  1:3085364/1‑71 (MQ=255)
                                                    tcatcaCCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATCTGGCGGCGGtt               <  1:2607877/71‑1 (MQ=255)
                                                      atcaCCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATCTGGCGGCGGTTa              <  1:4129509/70‑1 (MQ=255)
                                                      atcaCCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATCTGGCGGCGGTTa              <  1:5146238/70‑1 (MQ=255)
                                                       tcaCCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATCTggcggc                   >  1:3511601/1‑64 (MQ=38)
                                                         aCCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGaa                                       >  1:2558676/1‑42 (MQ=25)
                                                           ccTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATg                                    <  1:2986861/43‑1 (MQ=25)
                                                            cTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCg                                 >  1:4448692/1‑45 (MQ=37)
                                                            cTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATCTGGCGGCGGTTAAGGCCTa       <  1:3513343/71‑1 (MQ=255)
                                                              gAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATCTGGCGGCGGTTAAGGCCTATc     >  1:3041606/1‑71 (MQ=255)
                                                               aaaGATCTGATTGCCTATCTGGAAGAGAAGCCGTAAATGGCGGAAc                             >  1:448469/1‑46 (MQ=25)
                                                               aaaGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATCTGGCGGCGGTTaa             >  1:96517/1‑62 (MQ=38)
                                                               aaaGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATCTGGCGGCGGTTAAGGCCTATcg    >  1:1456536/1‑71 (MQ=255)
                                                               aaaGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATCTGGCGGCGGTTAAGGCCTATcg    >  1:1414763/1‑71 (MQ=255)
                                                               aaaGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATCTGGCGGCGGTTAAGGCCTATc     <  1:4844396/70‑1 (MQ=255)
                                                                aaGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATCTGGCGGCGGTTAAg            >  1:3734219/1‑62 (MQ=38)
                                                                aaGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATCTGGCGGCGGTTAAGGCCTATcgc   <  1:62046/71‑1 (MQ=255)
                                                                aaGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATCTGGCGGCGGTTAAGGCCTATcgc   <  1:2182828/71‑1 (MQ=255)
                                                                 aGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATCTGGCGGCGGTTAAGGCCt        >  1:4268967/1‑65 (MQ=39)
                                                                 aGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATCTGGCGGCGGTTAAGGCCTATcgcg  <  1:3339955/71‑1 (MQ=255)
                                                                    |                                                                   
ATTTCGGCGATTCAGGAAGTTGAGCGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGACCTGATTGCTTACCTGGAAGAGAAGCCGGAAATGGCGGAACATCTGGCGGCGGTTAAGGCCTATCGCG  >  AP009048/3825136‑3825271

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: