Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 3,824,949 A→G 100% E101E (GAA→GAG pyrE → orotate phosphoribosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0090483,824,9490AG100.0% 313.5 / NA 90E101E (GAA→GAGpyrEorotate phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (47/43);  total (47/43)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

ACCACAACCGCTGTGGCACTGGCGGAGCATCACGACCTGGACCTGCCGTACTGCTTTAACCGCAAAGAAGCAAAAGACCACGGTGAAGGCGGCAATCTGGTTGGTAGCGCGTTACAAGGACGCGTAATGCTGGTA  >  AP009048/3824881‑3825015
                                                                    |                                                                  
accacTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGc                                                                  <  1:221514/71‑1 (MQ=255)
accacTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGc                                                                  <  1:2627195/71‑1 (MQ=255)
 ccacTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGc                                                                  >  1:1566546/1‑70 (MQ=255)
 ccacTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGc                                                                  >  1:2139405/1‑70 (MQ=255)
 ccacTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGc                                                                  >  1:1800212/1‑70 (MQ=255)
  cactACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCaa                                                                <  1:2069658/71‑1 (MQ=255)
    ctaCCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCaaaa                                                              <  1:2502336/69‑1 (MQ=255)
    ctaCCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCaaaa                                                              <  1:4517601/69‑1 (MQ=255)
    ctaCCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCaaaa                                                              <  1:501875/69‑1 (MQ=255)
    ctaCCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCaa                                                                >  1:1627547/3‑69 (MQ=255)
     taCCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCaa                                                                <  1:4923786/67‑1 (MQ=255)
     taCCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCaa                                                                <  1:2607802/67‑1 (MQ=255)
       ccgccgTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAATAGGCAAAAGa                                                            >  1:1178309/1‑70 (MQ=255)
       ccgccgTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAg                                                             <  1:1888605/69‑1 (MQ=255)
       ccgccgTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGAc                                                           >  1:140365/1‑71 (MQ=255)
       ccgccgTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGAc                                                           >  1:1450498/1‑71 (MQ=255)
        cGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGAcc                                                          <  1:916455/71‑1 (MQ=255)
        cGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGAc                                                           >  1:3291487/1‑70 (MQ=255)
               gcgcTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTg                                                    <  1:4913456/67‑1 (MQ=255)
               gcgcTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGa                                                   <  1:4902635/68‑1 (MQ=255)
                cgcTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGaa                                                  <  1:4943401/69‑1 (MQ=255)
                  cTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCaa                                                                <  1:3149869/55‑1 (MQ=255)
                  cTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAgg                                                >  1:5121200/1‑71 (MQ=255)
                   tGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACgg                                                      >  1:1463398/1‑64 (MQ=255)
                     gCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAggggg                                             <  1:2219535/71‑1 (MQ=255)
                     gCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAggcgg                                             <  1:3402987/71‑1 (MQ=255)
                     gCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAggcgg                                             <  1:1043984/71‑1 (MQ=255)
                     gCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAggcgg                                             <  1:900494/71‑1 (MQ=255)
                     gCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAggcgg                                             >  1:4742852/1‑71 (MQ=255)
                     gCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAggcg                                              <  1:3770226/70‑1 (MQ=255)
                     gCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAggc                                               <  1:1582570/69‑1 (MQ=255)
                     gCGGAGCATCATGACCTAGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAggc                                               <  1:2548262/69‑1 (MQ=255)
                       ggAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCaa                                                                <  1:2627005/50‑1 (MQ=255)
                        gAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAGGCGGCaa                                          >  1:2876614/1‑71 (MQ=255)
                         aGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGa                                                            >  1:4535494/1‑52 (MQ=255)
                         aGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAGGCGGCAAt                                         >  1:4546732/1‑71 (MQ=255)
                         aGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAGGCGGCAAt                                         <  1:446042/71‑1 (MQ=255)
                           catcatGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAGGCGGCa                                           >  1:4457915/1‑67 (MQ=255)
                           catcatGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAGGCGGCAATCt                                       <  1:1691485/71‑1 (MQ=255)
                             tCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTg                                                    <  1:1041366/56‑1 (MQ=255)
                             tCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTg                                                    >  1:4838218/1‑56 (MQ=255)
                             tCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGa                                                   <  1:180695/57‑1 (MQ=255)
                             tCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAGGCGGCAATCtgg                                     <  1:227643/71‑1 (MQ=255)
                              catgaccttgacctgCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAGGCGGCAATCtggt                                    >  1:4395169/4‑71 (MQ=255)
                               atgacctttacctgCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAGGCGGCAATCtggt                                    >  1:2120980/3‑70 (MQ=255)
                               atgaccttgacctgCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAGGCGGCAATCtcgtt                                   >  1:4886440/3‑71 (MQ=255)
                                tgaccttgacctgCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAGGCGGCAATCtggttg                                  >  1:2256188/2‑71 (MQ=255)
                                 gaccttgacctgCCGTACTGCTTTAACCGCAAAGAGGCaaaa                                                              >  1:253125/1‑42 (MQ=38)
                                 gaccttgacctgCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAGGCGGCAATCtg                                      <  1:2881793/66‑1 (MQ=255)
                                 gaccttgacctgCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAGGGGCAATCtggttg                                  <  1:5089260/69‑1 (MQ=255)
                                  accttgacctgCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAGGCGGCAATCtggttggt                                <  1:889905/71‑1 (MQ=255)
                                  accttgacctgCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAGGCGGCAATCtggttgg                                 >  1:397848/1‑70 (MQ=255)
                                  accttgacatgCCGTACTGCTTTAACCGCAAAGAGGCAAAAGAACACGGTGAAGGCGGCAATc                                        >  1:1086625/1‑63 (MQ=255)
                                   ccttgacctgCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAGGCGGCAATCtggttgg                                 >  1:3929839/1‑69 (MQ=255)
                                     ttgacctgCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAGGCGGCAATCTGGTTGGTAgc                             >  1:170851/3‑71 (MQ=255)
                                      tgacctgCCGTACTGCTTTAACCGCAAAGAGGCAAAAGa                                                            >  1:2074985/2‑39 (MQ=39)
                                      tgacctgCCGTACTGCTTTAACCGCAAAGAGGCAAAAGAc                                                           >  1:145748/2‑40 (MQ=39)
                                      tgacctgCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAGGCGGCAATCtggttg                                  >  1:4649574/2‑65 (MQ=255)
                                       gacctgCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAGGCGGCAATCTGGTTGGTAgcg                            >  1:2079696/1‑70 (MQ=255)
                                       gacctgCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAGGCGGCAATCTGGTTGGTAgc                             <  1:452906/69‑1 (MQ=255)
                                          ctgCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAGGCGGCAATCTGGTTGGTAGcgcg                          <  1:2150866/69‑1 (MQ=255)
                                           tgCCGTACTGCTTTAACCTCAAAGAGGCAAAAGACCACGGTGAAGGCGGCAATCTGGTTGGTAGCGCGTTa                       >  1:3501470/1‑71 (MQ=255)
                                            gCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAGGCGGCAAt                                         >  1:4123186/1‑52 (MQ=255)
                                            gCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAGGCGGCAATCTGGTTGGTAGCGCGTTa                       >  1:119720/1‑70 (MQ=255)
                                             ccGTACTGCTTTAACCGCGAAGAGGCAAAAGACCACGGTGAAGGCGGCAATCtggt                                    <  1:604092/56‑1 (MQ=39)
                                             ccGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAGGCGGCAATCtggtt                                   >  1:2108556/1‑57 (MQ=255)
                                             ccGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAGGCGGCAATCTGGTTGGTAGCGCGTTAc                      <  1:132851/70‑1 (MQ=255)
                                             ccGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAGGCGGCAATCTGGTTGGTAGCGCGTTAc                      <  1:1332762/70‑1 (MQ=255)
                                             ccGCACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAGGCGGCAATCTGGTTGGTAGCGCGTTAc                      <  1:475048/70‑1 (MQ=255)
                                              cGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAGGCGGCAATCTGGTTGGTAGCGCGTTAc                      >  1:668589/1‑69 (MQ=255)
                                              cGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAGGCGGCAATCTGGTTGGTAGCGCGTTACaa                    <  1:567512/71‑1 (MQ=255)
                                                  cTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAGGCGGCAATCTGGTTGGTAGCGCGTTACAAGGAc                <  1:4216026/71‑1 (MQ=255)
                                                     cTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAGGCGGCAATCTGGTTGGTAGCGCGTTACAAGGAcgcg             >  1:4184009/1‑71 (MQ=255)
                                                     cTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAGGCGGCAATCTGGTTGGTAGCGCGTTACAAGGAcgcg             >  1:1396090/1‑71 (MQ=255)
                                                      tttAACCGCAAAGAGGCAAAAGACCACGGTGAAGGCGGCAATCTGGTTGGTa                               >  1:4000468/1‑52 (MQ=255)
                                                       ttAACCGCAAAGAGGCAAAAGACCACGGTGAAGGCGGCAATCTGGTTGGTAGCGCGTTACAAGGACGCGt            >  1:2364000/1‑70 (MQ=255)
                                                           ccGCAAAGAGGCAAAAGACCACGGTGAAGGCGGCAATCTGGTTGGTAGCGCGTTACAAGGACGCGTAATGc       <  1:3445523/71‑1 (MQ=255)
                                                            cGCAAAGAGGCAAAAGACCACGGTGAAGGCGGCAATCTGGTTGGTAGCGCGTTACAAGGACGCGTaa          >  1:1679190/1‑67 (MQ=255)
                                                            cGCAAAGAGGCAAAAGACCACGGTGAAGGCGGCAATCTGGTTGGTAGCGCGTTACAAGGACGCGTAATGCt      <  1:2695623/71‑1 (MQ=255)
                                                             gCAAAGAGGCAAAAGACCACGGTGAAGGCGGCAATCTGGTTGGTAGCGCGTTACAAg                   <  1:429569/57‑1 (MQ=255)
                                                             gCAAAGAGGCAAAAGACCACGGTGAAGGCGGCAATCTGGTTGGTAGCGCGTTACAAGGACGCGTAATGCTg     >  1:1162595/1‑71 (MQ=255)
                                                             gCAAAGAGGCAAAAGACCACGGTGAAGGCGGCAATCTGGTTGGTAGCGCGTTACAAGGACGCGTAATGCTg     >  1:3545647/1‑71 (MQ=255)
                                                              cAAAGAGGCAAAAGACCACGGTGAAGGCGGCAATCTGGTTGGTAGCGCGTTACAAGGAcg               >  1:4287101/1‑60 (MQ=255)
                                                              cAAAGAGGCAAAAGACCACGGTGAAGGCGGCAATCTGGTTGGTAGCGCGTTACAAGGAcg               >  1:67410/1‑60 (MQ=255)
                                                              cAAAGAGGCAAAAGACCACGGTGAAGGCCGCAATCtggtt                                   >  1:2336008/1‑40 (MQ=38)
                                                              cAAAGAGGCAAAAGACACGGTGAAGGCGGCAATCTGGTTGGTAGCGCGTTACAAGGAcg               >  1:5012814/1‑59 (MQ=255)
                                                               aaaGAGGCAAAAGACCACGGTGAAGGCGGCAATCtggtt                                   >  1:1058590/1‑39 (MQ=255)
                                                               aaaGAGGCAAAAGACCACGGTGAAGGCGGCAATCTGGTTGGTAGCGCGTTACAAGGACGCGTAATGCt      <  1:4933647/68‑1 (MQ=255)
                                                               aaaGAGGCAAAAGACCACGGTGAAGGCGGCAATCTGGTTGGTAGCGCGTTACAAGGACGCGTAATGCTgg    >  1:2645290/1‑70 (MQ=255)
                                                               aaaGAGGCAAAAGACCACGGTGAAGGCGGCAATCTGGTTGGTAGCGCGTTACAAGGACGCGTAATGCTGGt   <  1:3028187/71‑1 (MQ=255)
                                                                aagaggCAAAAGACCACGGTGAAGGCGGCAATCtggttggt                                >  1:2501278/1‑41 (MQ=255)
                                                                aagaggCAAAAGACCACGGTGAAGGCGGCAATCTGGTTGGTAGCGCGTTACAAGGACGCGTAATGCTgg    >  1:3652836/1‑69 (MQ=255)
                                                                aagaggCAAAAGACCACGGTGAAGGCGGCAATCTGGTTGGTAGCGCGTTACAAGGACGCGTAATGCTgg    >  1:4429002/1‑69 (MQ=255)
                                                                aagaggCAAAAGACCACGGTGAAGGCGGCAATCTGGTTGGTAGCGCGTTACAAGGACGCGTAATGCTGGTa  <  1:3845159/71‑1 (MQ=255)
                                                                 agaggCAAAAGACCACGGTGAAGGCGGCAATCTGGTTGGTAGCGCGTTACAAGGACGCGTAATGCTgg    >  1:3552173/1‑68 (MQ=255)
                                                                    |                                                                  
ACCACAACCGCTGTGGCACTGGCGGAGCATCACGACCTGGACCTGCCGTACTGCTTTAACCGCAAAGAAGCAAAAGACCACGGTGAAGGCGGCAATCTGGTTGGTAGCGCGTTACAAGGACGCGTAATGCTGGTA  >  AP009048/3824881‑3825015

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: