Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 3,824,919 G→T 100% L91L (CTG→CTT pyrE → orotate phosphoribosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0090483,824,9190GT98.6% 207.5 / ‑4.1 70L91L (CTG→CTTpyrEorotate phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base T (33/36);  minor base A (0/1);  total (33/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

CCCTGCTTACAAAGGGATCCCGATTGCCACCACAACCGCTGTGGCACTGGCGGAGCATCACGACCTGGACCTGCCGTACTGCTTTAACCGCAAAGAAGCAAAAGACCACGGTGAAGGCGGCAATCTGGTTGGT  >  AP009048/3824853‑3824985
                                                                  |                                                                  
cccTGCTTACAAAGTGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATgaccttgacc                                                                >  1:1593928/1‑71 (MQ=255)
cccTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATgaccttgacc                                                                >  1:3472502/1‑71 (MQ=255)
cccTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATgaccttga                                                                  >  1:3717598/1‑69 (MQ=255)
 ccTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCAtgacctgacctg                                                              <  1:2464907/71‑1 (MQ=255)
  cTGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATgaccttgacctg                                                              >  1:1797215/1‑71 (MQ=255)
   tGCTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGc                                                             >  1:3930200/1‑71 (MQ=255)
     cTTACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCg                                                           >  1:3303739/1‑71 (MQ=255)
      ttACAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCg                                                           >  1:2299084/1‑70 (MQ=255)
        aCAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGt                                                          <  1:766191/69‑1 (MQ=255)
        aCAAAGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTa                                                         <  1:1022044/70‑1 (MQ=255)
          aaaGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACt                                                       >  1:730349/1‑70 (MQ=255)
          aaaGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTg                                                      >  1:5003074/1‑71 (MQ=255)
            aGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATgaccttgac                                                                 >  1:4270914/1‑58 (MQ=255)
            aGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATgaccttgac                                                                 >  1:1026169/1‑58 (MQ=255)
            aGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCt                                                    >  1:3842773/1‑71 (MQ=255)
            aGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCt                                                    >  1:3138110/1‑71 (MQ=255)
            aGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCt                                                    >  1:2808094/1‑71 (MQ=255)
            aGGGATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCt                                                    >  1:2981782/1‑71 (MQ=255)
              ggATCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCttt                                                  <  1:2266032/71‑1 (MQ=255)
                 tCCCGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAAc                                               <  1:4258885/71‑1 (MQ=255)
                  cccGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTTCCGTACTGCTTTAAcc                                              >  1:5101386/1‑71 (MQ=255)
                   ccGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATgaccttgacctg                                                              >  1:1268162/1‑54 (MQ=25)
                   ccGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCg                                             <  1:2970445/71‑1 (MQ=255)
                    cGATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAAcc                                              >  1:616703/1‑69 (MQ=255)
                     gATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGc                                            >  1:1465747/1‑70 (MQ=255)
                     gATTGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCa                                           <  1:3379486/71‑1 (MQ=255)
                     gATTGCGACCACTACCGCCGTGCCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGc                                            >  1:4450107/1‑70 (MQ=255)
                        tGCGACCACTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCa                                           >  1:3155027/1‑68 (MQ=255)
                         gcgaccacTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCaaa                                         >  1:2370290/4‑69 (MQ=255)
                           gaccacTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAaagagg                                     <  1:3842048/70‑3 (MQ=255)
                            accacTACCGCCGTGGCGCTGGCGGAGCATCCTGACCTTGACCTGCCGTACTGCTTTAACCGCAaagagg                                     >  1:4199562/1‑68 (MQ=255)
                            accacTACCGCCGTGGCGCTGGCGGAGCATCATgaccttga                                                                  >  1:2749798/1‑41 (MQ=16)
                            accacTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGc                                    <  1:2627195/71‑1 (MQ=255)
                            accacTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGc                                    <  1:221514/71‑1 (MQ=255)
                             ccacTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGc                                    >  1:1566546/1‑70 (MQ=255)
                             ccacTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGc                                    >  1:1800212/1‑70 (MQ=255)
                             ccacTACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGc                                    >  1:2139405/1‑70 (MQ=255)
                              cactACCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCaa                                  <  1:2069658/71‑1 (MQ=255)
                                ctaCCGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCttt                                                  >  1:252124/3‑53 (MQ=37)
                                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)
                                    cGCCGTGGCGCTGGCGGAGCATCATGACCTTGACCTGCCGTa                                                         <  1:1004759/42‑1 (MQ=25)
                                    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)
                                             gcTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCttt                                                  <  1:3198096/39‑1 (MQ=38)
                                              cTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCaa                                  <  1:3149869/55‑1 (MQ=255)
                                              cTGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAgg                  >  1:5121200/1‑71 (MQ=255)
                                               tGGCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACgg                        >  1:1463398/1‑64 (MQ=255)
                                                 gCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTa                                                 <  1:2216593/37‑1 (MQ=38)
                                                 gCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAggggg               <  1:2219535/71‑1 (MQ=255)
                                                 gCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAggcgg               <  1:900494/71‑1 (MQ=255)
                                                 gCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAggcgg               <  1:3402987/71‑1 (MQ=255)
                                                 gCGGAGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAggcgg               <  1:1043984/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:446042/71‑1 (MQ=255)
                                                     aGCATCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAGGCGGCAAt           >  1:4546732/1‑71 (MQ=255)
                                                       catcatGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAGGCGGCa             >  1:4457915/1‑67 (MQ=255)
                                                       catcatGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTGAAGGCGGCAATCt         <  1:1691485/71‑1 (MQ=255)
                                                         tCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTg                      >  1:4838218/1‑56 (MQ=255)
                                                         tCATGACCTTGACCTGCCGTACTGCTTTAACCGCAAAGAGGCAAAAGACCACGGTg                      <  1:1041366/56‑1 (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)
                                                                  |                                                                  
CCCTGCTTACAAAGGGATCCCGATTGCCACCACAACCGCTGTGGCACTGGCGGAGCATCACGACCTGGACCTGCCGTACTGCTTTAACCGCAAAGAAGCAAAAGACCACGGTGAAGGCGGCAATCTGGTTGGT  >  AP009048/3824853‑3824985

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: