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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0090483,824,8920TC100.0% 304.1 / NA 88A82A (GCT→GCCpyrEorotate phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (41/47);  total (41/47)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

CCGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGCCACCACAACCGCTGTGGCACTGGCGGAGCATCACGACCTGGACCTGCCGTACTGCTTTAACCGCAAAGAAGCAAAAGAC  >  AP009048/3824825‑3824958
                                                                   |                                                                  
ctggCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCcgtg                                                                 >  1:2479960/3‑71 (MQ=255)
ctggCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCcgtg                                                                 >  1:1594698/3‑71 (MQ=255)
   gCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTgg                                                                <  1:4354227/69‑1 (MQ=255)
   gCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGcg                                                              >  1:3135177/1‑70 (MQ=255)
   gCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGcg                                                              >  1:4158777/1‑70 (MQ=255)
     aTTGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCt                                                            >  1:23210/1‑71 (MQ=255)
     aTTGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCt                                                            <  1:4907681/71‑1 (MQ=255)
      ttGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTg                                                           >  1:5323551/1‑71 (MQ=255)
      ttGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTg                                                           >  1:1551228/1‑71 (MQ=255)
        gAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGc                                                         >  1:2680430/1‑71 (MQ=255)
         aGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGc                                                         <  1:3849895/70‑1 (MQ=255)
         aGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCg                                                        <  1:4745166/71‑1 (MQ=255)
           ttCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCg                                                        <  1:1974802/69‑1 (MQ=255)
            tCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAg                                                     >  1:1478316/1‑71 (MQ=255)
            tCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAg                                                     >  1:393734/1‑71 (MQ=255)
            tCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAg                                                     >  1:3995840/1‑71 (MQ=255)
             cGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAg                                                     >  1:2586625/1‑70 (MQ=255)
             cGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAg                                                     >  1:1752790/1‑70 (MQ=255)
               aTCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCa                                                   <  1:3516315/70‑1 (MQ=255)
                  tgctgTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCAt                                                  <  1:1837579/68‑1 (MQ=255)
                  tgctgTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCAt                                                  <  1:4464207/68‑1 (MQ=255)
                  tgctgTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCAt                                                  <  1:2563559/68‑1 (MQ=255)
                    ctgTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGcg                                                              <  1:1737040/54‑2 (MQ=255)
                    ctgTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCg                                                        <  1:3399705/60‑1 (MQ=255)
                    ctgTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCa                                                   <  1:2547187/65‑1 (MQ=255)
                    ctgTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCat                                               <  1:249203/69‑2 (MQ=255)
                    ctgTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCa                                                <  1:1972707/68‑1 (MQ=255)
                    ctgTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATga                                             >  1:3953469/1‑71 (MQ=255)
                    ctgTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATga                                             >  1:2561727/1‑71 (MQ=255)
                       tttGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATgacct                                          >  1:1559841/1‑71 (MQ=255)
                         tGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGa                                                      >  1:606778/1‑57 (MQ=255)
                          ggCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATgaccttg                                        <  1:4302801/70‑3 (MQ=255)
                          ggCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATgaccttg                                        <  1:4866765/70‑3 (MQ=255)
                           gCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCat                                               <  1:44452/62‑2 (MQ=255)
                           gCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATgaccttgac                                      <  1:1032726/71‑1 (MQ=255)
                            cccTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATgaccttg                                        <  1:176838/68‑3 (MQ=255)
                             ccTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCa                                                   <  1:4265021/56‑1 (MQ=255)
                              cTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATgaccttgacctg                                   >  1:1407351/1‑71 (MQ=255)
                              cTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATgaccttgacc                                     <  1:3064459/69‑1 (MQ=255)
                               tGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGc                                  >  1:236181/1‑71 (MQ=255)
                               tGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGc                                  >  1:1571272/1‑71 (MQ=255)
                                gCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTgg                                                                >  1:1292321/1‑40 (MQ=37)
                                gCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATgaccttgacctg                                   >  1:5343904/1‑69 (MQ=255)
                                gCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGcc                                 >  1:2630670/1‑71 (MQ=255)
                                 cTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATgaccttgac                                      >  1:5163541/1‑65 (MQ=255)
                                 cTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGc                                  <  1:4735984/69‑1 (MQ=255)
                                 cTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCg                                >  1:2393895/1‑71 (MQ=255)
                                 cTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCg                                <  1:2728394/71‑1 (MQ=255)
                                    aCAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCg                                                        <  1:3341729/44‑1 (MQ=25)
                                    aCAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGt                               <  1:2199800/69‑1 (MQ=255)
                                    aCAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTa                              <  1:3888244/70‑1 (MQ=255)
                                    aCAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTa                              <  1:1416437/70‑1 (MQ=255)
                                    aCAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTAc                             >  1:2931801/1‑71 (MQ=255)
                                     cAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACt                            <  1:3402405/71‑1 (MQ=255)
                                      aaaGGGATCCCGATTGCGACCACTACCGCCGTGGc                                                               <  1:267554/35‑1 (MQ=25)
                                      aaaGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACt                            >  1:665569/1‑70 (MQ=255)
                                      aaaGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTg                           <  1:1187243/71‑1 (MQ=255)
                                       aaGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGc                          >  1:249454/1‑71 (MQ=255)
                                        aGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACtt                           <  1:4805076/69‑2 (MQ=255)
                                        aGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCt                         >  1:5325971/1‑71 (MQ=255)
                                          ggATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTg                           >  1:1938844/1‑67 (MQ=255)
                                          ggATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCtt                        >  1:1682255/1‑70 (MQ=255)
                                           gATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCttt                       >  1:914967/1‑70 (MQ=255)
                                              cccGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCttt                       <  1:3829576/67‑1 (MQ=255)
                                              cccGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAAcc                   >  1:4560910/1‑71 (MQ=255)
                                              cccGATTGCGACCACTACCGCCGTGGCGCTAGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAAcc                   >  1:1966408/1‑71 (MQ=255)
                                              cccGATTGCGACCACTACCGCCGTGGCGCTAGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAAcc                   >  1:1972628/1‑71 (MQ=255)
                                               ccGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTa                      <  1:4371101/67‑1 (MQ=255)
                                                 gATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATgacctt                                         <  1:357017/46‑2 (MQ=19)
                                                 gATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGc                 >  1:5027929/1‑70 (MQ=255)
                                                 gATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGc                 >  1:1331533/1‑70 (MQ=255)
                                                 gATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGc                 >  1:2878316/1‑70 (MQ=255)
                                                 gATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGc                 >  1:4813155/1‑70 (MQ=255)
                                                  aTTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTAc                             <  1:1848147/57‑1 (MQ=255)
                                                       gaccacTACCGCCGTGGCGCTGGCGGAGCATCATgacctt                                         <  1:2273824/39‑2 (MQ=16)
                                                       gaccacTACCGCCGTGGCGCTGGCGGAGCATCATgac                                            <  1:3864486/36‑1 (MQ=16)
                                                         ccacTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCg                                <  1:2239891/47‑1 (MQ=25)
                                                         ccacTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCg                                <  1:4028579/47‑1 (MQ=25)
                                                         ccacTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAaaga            >  1:1097051/1‑67 (MQ=255)
                                                         ccacTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGc         >  1:400672/1‑70 (MQ=255)
                                                          cactACCGCCGTGGCGCTGGCGGAGCATCATgacctt                                         <  1:846051/37‑2 (MQ=16)
                                                          cactACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAaaga            >  1:3048149/1‑66 (MQ=255)
                                                          cactACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCaa       <  1:3644079/71‑1 (MQ=255)
                                                          cactACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCaa       <  1:1886298/71‑1 (MQ=255)
                                                            ctaCCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCaaaa     <  1:3964653/69‑1 (MQ=255)
                                                            ctaCCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCaa       <  1:2837012/67‑1 (MQ=255)
                                                              aCCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTa                              >  1:1973744/1‑44 (MQ=25)
                                                              aCCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTaa                     <  1:1569717/53‑1 (MQ=37)
                                                               ccgccgTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAAc                    >  1:4929748/1‑53 (MQ=37)
                                                               ccgccgTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGa   >  1:4187187/1‑70 (MQ=255)
                                                               ccgccgTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGAc  >  1:2466350/1‑71 (MQ=255)
                                                               ccgccgTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGAc  >  1:1156487/1‑71 (MQ=255)
                                                                cGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCt                         <  1:1488623/47‑1 (MQ=25)
                                                                cGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCCTACTGCTTTAACCGCAAAGAGGCaaaa     <  1:2610851/67‑1 (MQ=255)
                                                                   |                                                                  
CCGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGCCACCACAACCGCTGTGGCACTGGCGGAGCATCACGACCTGGACCTGCCGTACTGCTTTAACCGCAAAGAAGCAAAAGAC  >  AP009048/3824825‑3824958

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: