Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 3,824,826 C→T 100% S60S (TCC→TCT pyrE → orotate phosphoribosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0090483,824,8260CT98.9% 261.4 / ‑4.2 87S60S (TCC→TCTpyrEorotate phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (48/38);  minor base A (0/1);  total (48/39)
Fisher's exact test for biased strand distribution p-value = 4.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-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.

GCTGTTTAATACCGGGCGCGATCTGGCACTGTTAGGCCGTTTTTACGCTGAAGCGTTGGTGGATTCCGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGCCACCACAACCGCT  >  AP009048/3824760‑3824892
                                                                  |                                                                  
gCTGTTTAATACCGGGCGCGACCTGGCACTGTTATGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCa                                                                <  1:1043450/71‑1 (MQ=255)
gCTGTTTAATACCGGGCGCGACCTGGCACTGTTAGTCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCa                                                                <  1:2895270/71‑1 (MQ=255)
  tGTTTAATACCGGGCGCGACCTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCAtt                                                              >  1:3305842/1‑71 (MQ=255)
     ttAATACCGGGCGCGACCTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAg                                                           >  1:94957/1‑71 (MQ=255)
      tAATACCGGGCGCGACCTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGt                                                          >  1:4374833/1‑71 (MQ=255)
      tAATACCGGGCGCGACCTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGt                                                          >  1:3285295/1‑71 (MQ=255)
        aTACCGGGCGCGACCTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGtt                                                         <  1:4281109/70‑1 (MQ=255)
        aTACCGGGCGCGACCTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTc                                                        >  1:206068/1‑71 (MQ=255)
         tACCGGGCGCGACCTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGa                                                            >  1:4438409/1‑66 (MQ=255)
           ccGGGCGCGACCTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGa                                                      >  1:2782653/1‑70 (MQ=255)
           ccGGGCGCGACCTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGa                                                      >  1:2788190/1‑70 (MQ=255)
           ccGGGCGCGACCTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGAt                                                     <  1:117839/71‑1 (MQ=255)
           ccGGGCGCGACCTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGAt                                                     >  1:2158296/1‑71 (MQ=255)
           ccGGGCGCGACCTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGAt                                                     <  1:3404113/71‑1 (MQ=255)
           ccGGGCGCGACCTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGAt                                                     >  1:4070336/1‑71 (MQ=255)
           ccGGGCGCGACCTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGAt                                                     >  1:4139143/1‑71 (MQ=255)
              ggCGCGACCTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGATctg                                                  >  1:2845918/1‑71 (MQ=255)
              ggCGCGACCTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGATc                                                    >  1:4364819/1‑69 (MQ=255)
               gcgcGACCTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGAt                                                     <  1:457459/67‑1 (MQ=255)
                cgcgACCTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGATctgc                                                 >  1:1499247/1‑70 (MQ=255)
                 gcgACCTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGATctgctg                                               >  1:3535733/1‑71 (MQ=255)
                     ccTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGATCTGCTGTTTg                                           >  1:3858881/2‑71 (MQ=255)
                     ccTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGATCTGCTGTTTg                                           >  1:3240178/2‑71 (MQ=255)
                     ccTGGCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGATCTGCTGTTTg                                           >  1:1984867/2‑71 (MQ=255)
                       tGGCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGATCTGCTGtt                                             >  1:3754040/1‑67 (MQ=255)
                        ggCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGATCTGCTGTTTgg                                          >  1:1094897/1‑69 (MQ=255)
                        ggCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCAgg                                                                  <  1:324476/45‑1 (MQ=38)
                         gCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGATCTGCTGt                                              >  1:2279917/1‑64 (MQ=255)
                         gCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGATCTGCTGTTTGGccc                                       >  1:1005310/1‑71 (MQ=255)
                         gCACTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGATCTGCTGTTTGGc                                         <  1:822305/69‑1 (MQ=255)
                          cACTGTTAGGCCGTTTTTACGCTGAAGCTTTGGTGGATTCTGGCATTGAGTTCGATCTGCTGTTTGGCCCt                                      <  1:2045673/71‑1 (MQ=255)
                           aCTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGATCTGCTGttt                                            >  1:2542665/1‑64 (MQ=255)
                           aCTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGATCTGCTGTTTGGCCCTg                                     >  1:2469317/1‑71 (MQ=255)
                           aCTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGCCATTGAGTTCGATCTGCTGTTTGGCCCTg                                     <  1:2698616/71‑1 (MQ=255)
                            cTGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCAtt                                                              <  1:4439147/45‑1 (MQ=38)
                             tGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGATCTGCTGTTTGGcc                                        >  1:3743949/1‑66 (MQ=255)
                             tGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCt                                   >  1:3912772/1‑71 (MQ=255)
                             tGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCt                                   >  1:1626164/1‑71 (MQ=255)
                             tGTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCt                                   >  1:4525788/1‑71 (MQ=255)
                              gTTAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGATCTGCTGt                                              <  1:3533710/59‑1 (MQ=255)
                               ttAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGAt                                                     >  1:1678535/1‑51 (MQ=255)
                               ttAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGATCTGCTGTTTGGCCCTGc                                    >  1:3627518/1‑68 (MQ=255)
                               atAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTgg                                                                  >  1:4006884/2‑38 (MQ=37)
                                tAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCAtt                                                              >  1:408364/1‑41 (MQ=38)
                                tAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTa                                 >  1:2458549/1‑70 (MQ=255)
                                tAGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTa                                 >  1:3367309/1‑70 (MQ=255)
                                 aGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCg                                                       <  1:4322680/47‑1 (MQ=39)
                                 aGGCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTAc                                >  1:1565679/1‑70 (MQ=255)
                                  ggCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACa                               <  1:4333436/70‑1 (MQ=255)
                                   gCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGt                                                          <  1:1034522/42‑1 (MQ=38)
                                   gCCGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACaaa                             <  1:2424188/71‑1 (MQ=255)
                                     cGTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAgg                           >  1:14063/1‑71 (MQ=255)
                                      gTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGAt                                                     >  1:1697544/1‑44 (MQ=38)
                                      gTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAggg                          <  1:1536608/71‑1 (MQ=255)
                                      gTTTTTACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAggg                          <  1:343236/71‑1 (MQ=255)
                                        ttttACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGATCTGCTGTTTGGCCCTGc                                    <  1:941929/59‑1 (MQ=255)
                                         tttACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAg                            >  1:2174700/1‑66 (MQ=255)
                                         tttACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGAt                        <  1:4330933/70‑1 (MQ=255)
                                           tACGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATcc                      <  1:2487395/70‑1 (MQ=255)
                                            aCGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGATCTGCTGTTTg                                           >  1:941734/1‑48 (MQ=39)
                                            aCGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATccc                     >  1:3612353/1‑70 (MQ=255)
                                            aCGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATccc                     <  1:3647606/70‑1 (MQ=255)
                                             cGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGTATCCCGa                   <  1:2673844/71‑1 (MQ=255)
                                             cGCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGa                   <  1:535655/71‑1 (MQ=255)
                                              gCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGa                   <  1:1258216/70‑1 (MQ=255)
                                              gCTGAAGCATTGGTGGATTCTGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGAt                  <  1:2179273/71‑1 (MQ=255)
                                                tGAAGCATTGGTGGATTCTGGCATTGAGTTCGATCTGCTGTTTgg                                          <  1:1583131/45‑1 (MQ=38)
                                                  aaGCATTGGTGGATTCTGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGc               <  1:4045436/70‑1 (MQ=255)
                                                  aaGCATTGGTGGATTCTGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGc               >  1:1197125/1‑70 (MQ=255)
                                                  aaGCATTGGTGGATTCTGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGc               <  1:1900282/70‑1 (MQ=255)
                                                    gcattGGTGGATTCTGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGcgac            <  1:1600689/68‑1 (MQ=255)
                                                     cattGGTGGATTCTGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGcgacc           <  1:732104/69‑1 (MQ=255)
                                                      attGGTGGATTCTGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAgg                           >  1:3771047/2‑54 (MQ=255)
                                                      attGGTGGATTCTGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGcgacca          >  1:4037447/2‑71 (MQ=255)
                                                       ttGGTGGATTCTGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGCgaccac         <  1:3313188/71‑1 (MQ=255)
                                                       ttGGTGGATTCTGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGCgaccac         <  1:2469207/71‑1 (MQ=255)
                                                        tggtggATTCTGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACaaa                             >  1:3942010/1‑50 (MQ=255)
                                                         ggtggATTCTGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGAt                        >  1:2768310/1‑54 (MQ=255)
                                                         ggtggATTCTGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGcgacca          >  1:4353015/1‑68 (MQ=255)
                                                         ggtggATTCTGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGCgaccact        <  1:1173327/70‑2 (MQ=255)
                                                          gtggATTCTGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGCgaccacta       <  1:2385354/70‑3 (MQ=255)
                                                            ggATTCTGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTAc                                <  1:2013871/43‑1 (MQ=255)
                                                            ggATTCTGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGcgacca          >  1:4082331/1‑65 (MQ=255)
                                                            ggATTCTGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGc               <  1:2055809/60‑1 (MQ=255)
                                                             gATTCTGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGa                   <  1:2321605/55‑1 (MQ=255)
                                                              aTTCTGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGCGACCACTACCGcc  >  1:1773308/1‑70 (MQ=255)
                                                               ttCTGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGa                   <  1:1827400/53‑1 (MQ=255)
                                                                  |                                                                  
GCTGTTTAATACCGGGCGCGATCTGGCACTGTTAGGCCGTTTTTACGCTGAAGCGTTGGTGGATTCCGGCATTGAGTTCGATCTGCTGTTTGGCCCTGCTTACAAAGGGATCCCGATTGCCACCACAACCGCT  >  AP009048/3824760‑3824892

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: