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% 297.3 / NA 99D186D (GAC→GATpyrEorotate phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (47/52);  total (47/52)
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
                                                                    |                                                                   
aTTTCGGCGATTCAGGAAGTTGAGCGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAGAGATCt                                                                   >  1:425622/1‑71 (MQ=255)
 tttCGGCGATTCAGGAAGTTGAGCGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTg                                                                  <  1:767453/71‑1 (MQ=255)
   tCGGCGATTCAGGAAGTTGAGCGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGa                                                                 <  1:108477/70‑1 (MQ=255)
   tCGGCGATTCAGGAAGTTGAGCGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGAt                                                                >  1:140280/1‑71 (MQ=255)
    cGGCGATTCAGGAAGTTGAGCGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGAtt                                                               >  1:5812721/1‑71 (MQ=255)
    cGGCGATTCAGGAAGTTGAGCGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGAtt                                                               >  1:3693224/1‑71 (MQ=255)
       cGATTCAGGAAGTTGAGCGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCt                                                                   <  1:1837525/64‑1 (MQ=255)
        gATTCAGGAAGTTGAGCGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGcc                                                            >  1:6827631/1‑69 (MQ=255)
           tCAGGAAGTTGAGCGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTatc                                                        >  1:6846132/1‑69 (MQ=255)
              ggAAGTTGAGCGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCt                                                                   <  1:2276646/57‑1 (MQ=255)
              ggAAGTTGAGCGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGc                                                             <  1:7660005/63‑1 (MQ=255)
              ggAAGTTGAGCGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTgg                                                     <  1:1550208/71‑1 (MQ=255)
              ggAAGTTGAGCGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTg                                                      >  1:7330236/1‑70 (MQ=255)
                aaGTTGAGCGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTgg                                                     >  1:6722452/1‑69 (MQ=255)
                 aGTTGAGCGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAag                                                  <  1:2094922/71‑1 (MQ=255)
                  gTTGAGCGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTgg                                                     <  1:5848457/67‑1 (MQ=255)
                  gTTGAGCGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGaa                                                   >  1:802478/1‑69 (MQ=255)
                  gTTGAGCGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAaga                                                 <  1:251676/71‑1 (MQ=255)
                  gTTGAGCGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAaga                                                 <  1:3269838/71‑1 (MQ=255)
                  gTTGAGCGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAaga                                                 <  1:3594807/71‑1 (MQ=255)
                  gTTGAGCGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAaga                                                 <  1:145593/71‑1 (MQ=255)
                   ttGAGCGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAagag                                                <  1:832222/71‑1 (MQ=255)
                    tGAGCGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGa                                                                 >  1:5283011/1‑53 (MQ=255)
                    tGAGCGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGa                                                                 >  1:796249/1‑53 (MQ=255)
                    tGAGCGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAagag                                                >  1:7545473/1‑70 (MQ=255)
                     gAGCGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAaga                                                 >  1:404399/1‑68 (MQ=255)
                      aGCGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAg                                             >  1:2813475/1‑71 (MQ=255)
                       gCGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAgaga                                               >  1:3147174/1‑68 (MQ=255)
                       gCGTGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGc                                            <  1:786453/71‑1 (MQ=255)
                          tGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGcc                                           <  1:7353509/69‑1 (MQ=255)
                          tGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGcc                                           <  1:5879748/69‑1 (MQ=255)
                          tGATTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCgg                                         <  1:1606917/71‑1 (MQ=255)
                            aTTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCgg                                         >  1:7330531/1‑69 (MQ=255)
                            aTTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCgg                                         >  1:7036496/1‑69 (MQ=255)
                            aTTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCgg                                         >  1:5849133/1‑69 (MQ=255)
                            aTTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGaa                                       <  1:5910878/71‑1 (MQ=255)
                            aTTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGaa                                       <  1:5783270/71‑1 (MQ=255)
                            aTTACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGa                                        <  1:5559905/70‑1 (MQ=255)
                             ttACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGAt                                                                >  1:1561286/1‑45 (MQ=255)
                             ttACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGAt                                                                >  1:5444733/1‑45 (MQ=255)
                             ttACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGAt                                                                >  1:5773352/1‑45 (MQ=255)
                             ttACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGAt                                                                >  1:4384149/1‑45 (MQ=255)
                              tACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGAt                                                                >  1:6275930/1‑44 (MQ=255)
                              tACAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAAt                                     >  1:3717963/1‑71 (MQ=255)
                               aCAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAagag                                                >  1:1338930/1‑59 (MQ=38)
                                cAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATgg                                   >  1:3776052/1‑71 (MQ=255)
                                cAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATg                                    <  1:6926011/70‑1 (MQ=255)
                                cAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATg                                    <  1:4790229/70‑1 (MQ=255)
                                cAACTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATg                                    <  1:987778/70‑1 (MQ=255)
                                 aaCTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGc                                            <  1:554199/61‑1 (MQ=255)
                                 aaCTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGc                                  <  1:1412897/71‑1 (MQ=255)
                                   cTGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCg                                 >  1:3510420/1‑70 (MQ=255)
                                    tGCAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCgg                                >  1:5284581/1‑70 (MQ=255)
                                      cAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCg                                 <  1:3258414/67‑1 (MQ=255)
                                      cAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCg                                 <  1:7676805/67‑1 (MQ=255)
                                      cAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGaa                              >  1:3668861/1‑70 (MQ=255)
                                      cAAAGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGa                               <  1:2694892/69‑1 (MQ=255)
                                      cAAAGTGATCTCTATCATCACCCTGAAAGATATGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCgg                                >  1:6504479/1‑68 (MQ=255)
                                       aaaGTGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACa                            <  1:2316765/71‑1 (MQ=255)
                                          gTGATCTCTATCATCACCCTGAAAGATCTGATGGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATCt                         >  1:6465145/1‑71 (MQ=255)
                                           tGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGc                                            >  1:4897106/1‑51 (MQ=38)
                                           tGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGc                                            >  1:5288726/1‑51 (MQ=38)
                                           tGATCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGc                                            >  1:2667087/1‑51 (MQ=38)
                                             aTCTCTATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAAt                                     <  1:1772639/56‑1 (MQ=38)
                                                tctATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGaa                                              <  1:5017695/44‑1 (MQ=25)
                                                tctATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATCTggcggc                   >  1:3900297/1‑71 (MQ=255)
                                                 ctATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATCTGgcggcg                  >  1:2687971/1‑71 (MQ=255)
                                                  tATCATCACCCTGAAAGATCTGATTGCCTATCTGGa                                                    <  1:6962896/36‑1 (MQ=21)
                                                  tATCATCACCCTGAAAGATCTGATTGCCTATCTGGa                                                    <  1:1290945/36‑1 (MQ=21)
                                                  tATCATCACCCTGAAAGATCTGATTGCCTATCTGGa                                                    <  1:5909923/36‑1 (MQ=21)
                                                  tATCATCACCCTGAAAGATCTGATTGCCTATCTGGa                                                    <  1:2956585/36‑1 (MQ=21)
                                                  tATCATCACCCTGAAAGATCTGATTGCCTATCTGGa                                                    <  1:3878447/36‑1 (MQ=21)
                                                  tATCATCACCCTGAAAGATCTGATTGCCTATCTGGa                                                    <  1:5738160/36‑1 (MQ=21)
                                                  tATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATCTggcggc                   <  1:2485798/69‑1 (MQ=39)
                                                  tATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATCTggcggc                   <  1:2654715/69‑1 (MQ=39)
                                                  tATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATCTggcggc                   <  1:7011387/69‑1 (MQ=39)
                                                  tATCATCACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATCTGGcggcgg                 >  1:6820695/1‑71 (MQ=255)
                                                   atcatcACCCTGAAAGATCTGATTGCCTATCTGGa                                                    <  1:6046047/35‑1 (MQ=21)
                                                   atcatcACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAg                                             >  1:5645982/1‑42 (MQ=25)
                                                   atcatcACCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAg                                             >  1:2507891/1‑42 (MQ=25)
                                                    tcatcaCCCTGAAAGATCTGATTGCCTATCTGTAAGAGAAGCCGGAAATGGCGGAACATCTGGCGGCGGtt               <  1:2948899/71‑1 (MQ=255)
                                                     catcaCCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGc                                  >  1:3481267/1‑51 (MQ=38)
                                                      atcaCCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATCTGGCGGCGGtt               >  1:1660342/1‑69 (MQ=39)
                                                      atcaCCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATCTGGCGGCGGTTaa             <  1:4819452/71‑1 (MQ=255)
                                                       tcaCCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATCTGGCGGCGGTTAAg            <  1:6820222/71‑1 (MQ=255)
                                                       tcaCCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATCTGGCGGCGGTTAAg            <  1:5303332/71‑1 (MQ=255)
                                                        caCCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAAt                                     <  1:4128047/45‑1 (MQ=37)
                                                        caCCCTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATCTggcgg                    <  1:3603088/62‑1 (MQ=38)
                                                           ccTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAAc                             <  1:2717707/50‑1 (MQ=38)
                                                            cTGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATCTGGCGGCGGTTAAgg           <  1:3597620/67‑1 (MQ=39)
                                                             tGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACAt                           >  1:3660055/1‑50 (MQ=38)
                                                             tGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATCTGGCGGCGGTTAAGGCCTAt      >  1:5741753/1‑71 (MQ=255)
                                                             tGAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATCTGGCGGCGGTTAAGGCCTAt      >  1:7101017/1‑71 (MQ=255)
                                                             agAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATCTGGCGCCGGTTaa             >  1:5323970/2‑64 (MQ=255)
                                                              gAAAGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATCTGGCGGCGGtt               <  1:371193/61‑1 (MQ=38)
                                                               aaaGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAAc                             >  1:417316/1‑46 (MQ=37)
                                                               aaaGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATCTGGCGGCGGTTAAGGCCTATcg    >  1:7479453/1‑71 (MQ=255)
                                                                 aGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATCTGGCGGCGGTTAAGGCCTATcgcg  <  1:1863989/71‑1 (MQ=255)
                                                                 aGATCTGATTGCCTATCTGGAAGAGAAGCCGGAAATGGCGGAACATCTGGCGGCGGTTAAGGCCTATcgcg  <  1:6358385/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: