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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0090483,825,0510GA100.0% 249.7 / NA 81E135E (GAG→GAApyrEorotate phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (34/47);  total (34/47)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GGTAGCGCGTTACAAGGACGCGTAATGCTGGTAGATGATGTGATCACCGCCGGAACGGCGATTCGCGAGTCGATGGAGATTATTCAGGCCAATGGCGCGACGCTTGCTGGCGTGTTGATTTCGCTCGATCGTC  >  AP009048/3824983‑3825115
                                                                    |                                                                
ggtAGCGCGTTACAAGGACGCGTAATGCTGGTAGATGATGTGATCACCGCCGGAACGGCGATTCGCGAATc                                                                >  1:2165954/1‑71 (MQ=255)
 gtAGCGCGTTACAAGGACGCGTAATGCTGGTAGATGATGTGATCACCGCCGGAACGGCGATTCGCGAATc                                                                <  1:2494876/70‑1 (MQ=255)
 gtAGCGCGTTACAAGGACGCGTAATGCTGGTAGATGATGTGATCACCGCCGGAACGGCGATTCGCGAATCg                                                               <  1:4512522/71‑1 (MQ=255)
  tAGCGCGTTACAAGGACGCGTAATGCTGGTAGATGATGTGATCACCGCCGGAACGGCGATTCGCGAATCGa                                                              <  1:445446/71‑1 (MQ=255)
  tAGCGCGTTACAAGGACGCGTAATGCTGGTAGATGATGTGATCACCGCCGGAACGGCGATTCGCGAATCGa                                                              >  1:984052/1‑71 (MQ=255)
    gcgcGTTACAAGGACGCGTAATGCTGGTAGATGATGTGATCACCGCCGGAACGGCGATTCGCGAATCGAt                                                             <  1:2491985/70‑1 (MQ=255)
    gcgcGTTACAAGGACGCGTAATGCTGGTAGATGATGTGATCACCGCCGGAACGGCGATTCGCGAATCGAt                                                             >  1:4236615/1‑70 (MQ=255)
    gcgcGTTACAAGGACGCGTAATGCTGGTAGATGATGTGATCACCGCCGGAACGGCGATTCGCGAATCGATg                                                            <  1:2626659/71‑1 (MQ=255)
       cgTTACAAGGACGCGTAATGCTGGTAGATGATGTGATCACCGCCGGAACGGCGATTCGCGAATCGATg                                                            >  1:2295082/1‑68 (MQ=255)
       cgTTACAAGGACGCGTAATGCTGGTAGATGATGTGATCACCGCCGGAACGGCGATTCGCGAATCGATGgag                                                         <  1:632462/71‑1 (MQ=255)
             aaGGACGCGTAATGCTGGTAGATGATGTGATCACCGCCGGAACGGCGATTCGCGAATCGATGGAGattatt                                                   >  1:151699/1‑71 (MQ=255)
              aGGACGCGTAATGCTGGTAGATGATGTGATCACCGCCGGAACGGCGATTCGCGAATCGATGGAGattatt                                                   <  1:525086/70‑1 (MQ=255)
              aGGACGCGTAATGCTGGTAGATGATGTGATCACCGCCGGAACGGCGATTCGCGAATCGATGGAGattatt                                                   <  1:4097897/70‑1 (MQ=255)
              aGGACGCGTAATGCTGGTAGATGATGTGATCACCGCCGGAACGGCGATTCGCGAATCGATGGAGATTATTc                                                  <  1:1325170/71‑1 (MQ=255)
               ggACGCGTAATGCTGGTAGATGATGTGATCACCGCCGGAACGGCGATTCGCGAATCGATGGAGattatt                                                   <  1:771205/69‑1 (MQ=255)
                  cgcgTAATGCTGGTAGATGATGTGATCACCGCCGGAACGGCGATTCGCGAATCGATGGAGATTATTCAgg                                               <  1:2115368/70‑1 (MQ=255)
                     gTAATGCTGGTAGATGATGTGATCACCGCCGGAACGGCGATTCGCGAATCGATGGAGATTATTCAGGct                                             >  1:1640609/1‑68 (MQ=255)
                     gTAATGCTGGTAGATGATGTGATCACCGCCGGAACGGCGATTCGCGAATCGATGGAGATTATTCAGGCTaa                                           >  1:451702/1‑71 (MQ=255)
                     gTAATGCTGGTAGATGATGTGATCACCGCCGGAACGGCGATTCGCGAATCGATGGAGATTATTCAGGCTaa                                           <  1:291911/71‑1 (MQ=255)
                      tAATGCTGGTAGATGATGTGATCACCGCCGGAACGGCGATTCGCGAATCGATGGAGATTATTCAGGct                                             <  1:81738/68‑2 (MQ=255)
                      tAATGCTGGTAGATGATGTGATCACCGCCGGAACGGCGATTCGCGAATCGATGGAGATTATTCAGGCTaa                                           <  1:2912019/70‑1 (MQ=255)
                       aaTGCTGGTAGATGATGTGATCACCGCCGGAACGGCGATTCGCGAATCGATGGAGATTATTCAGGCTAATg                                         >  1:406701/1‑71 (MQ=255)
                        aTGCTGGTAGATGATGTGATCACCGCCGGAACGGCGATTCGCGAATCGATGgaga                                                        <  1:4013524/55‑1 (MQ=255)
                         tGCTGGTAGATGATGTGATCACCGCCGGAACGGCGATTCGCGAATc                                                                <  1:4019314/46‑1 (MQ=255)
                         tGCTGGTAGATGATGTGATCACCGCCGGAACGGCGATTCGCGAATc                                                                <  1:3280093/46‑1 (MQ=255)
                         tGCTGGTAGATGATGTGATCACCGCCGGAACGGCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGgc                                       <  1:1275995/71‑1 (MQ=255)
                          gCTGGTAGATGATGTGATCACCGCCGGAACGGCGATTCGCGAATCGATGGAGattat                                                    <  1:3334781/57‑1 (MQ=255)
                           cTGGTAGATGATGTGATCACCGCCGGAACGGCGATTCGCGAATCGATGGAGattat                                                    <  1:1166505/56‑1 (MQ=255)
                           cTGGTAGATGATGTGATCACCGCCGGAACGGCGATTCGCGAATCGATGGAGatta                                                     <  1:1052266/55‑1 (MQ=255)
                             ggTAGATGATGTGATCACCGCCGGAACGGCGATTCGCGAATCGATGGAGATTATTCAg                                                <  1:4355147/58‑1 (MQ=255)
                              gTAGATGATGTGATCACCGCCGGAACGGCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGAc                                  <  1:3705377/71‑1 (MQ=255)
                              gTAGATGATGTGATCACCGCCGGAACGGCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGAc                                  >  1:2804411/1‑71 (MQ=255)
                               tAGATGATGTGATCACCGCCGGAACGGCGATTCGCGAATCGATGGAGATTATTCAGGc                                              >  1:4231823/1‑58 (MQ=255)
                               tAGATGATGTGATCACCGCCGGAACGGCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGAc                                  >  1:1504258/1‑70 (MQ=255)
                                 gatgatGTGATCACCGCCGGAACGGCGATTCGCGAATCGATGGAGattatt                                                   <  1:4116019/51‑1 (MQ=255)
                                   tgatgTGATCACCGCCGGAACGGCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGcgt                              >  1:1320742/1‑68 (MQ=255)
                                   tgatgTGATCACCGCCGGAACGGCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTg                             <  1:549396/71‑1 (MQ=255)
                                      tgtgATCACCGCCGGAACGGCGATTCGCGAATCGATGgag                                                         <  1:1147421/40‑1 (MQ=255)
                                      tgtgATCACCGCCGGAACGGCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGcgcg                                    <  1:2105981/61‑1 (MQ=255)
                                      tgtgATCACCGCCGGAACGGCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTg                          >  1:3300849/1‑71 (MQ=255)
                                      tgtgATCACCGCCGGAACGGCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTg                          <  1:2417389/71‑1 (MQ=255)
                                       gtgATCACCGCCGGAACGGCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTg                             <  1:885950/67‑1 (MQ=255)
                                        tgATCACCGCCGGAACGGCGATTCGCGAATCGATGGAGattatt                                                   <  1:637518/44‑1 (MQ=255)
                                        tgATCACCGCCGGAACGGCGATTCGCGAATCGATGGAGATTATTCAGGCTaa                                           <  1:3738126/52‑1 (MQ=39)
                                        tgATCACCGCCGGAACGGCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGATGGc                        >  1:4149700/1‑71 (MQ=255)
                                           tCACCGCCGGAACGGCGATTCGCGAATCGATGGAGatta                                                     >  1:3489574/1‑39 (MQ=255)
                                           tCACCGCCGGAACGGCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCgta                     >  1:1649644/1‑70 (MQ=255)
                                           tCACCGCCGGAACGGCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGATGGCgta                     <  1:1124544/71‑2 (MQ=255)
                                            cACCGCCGGAACGGCGATTCGCTAATCGATGGAGattatt                                                   <  1:3819686/40‑1 (MQ=38)
                                            cACCGCCGGAACGGCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCgtat                    >  1:1586172/1‑69 (MQ=255)
                                             aCCGCCGGAACGGCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTAtt                   >  1:3431158/1‑71 (MQ=255)
                                              ccgccgGAACGGCGATTCGCGAATCGATGGAGatta                                                     <  1:2559948/36‑1 (MQ=255)
                                               cgccgGAACGGCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGa                 >  1:212467/1‑71 (MQ=255)
                                                 ccgGAACGGCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGCCGCTTGCTgg                         >  1:2632190/1‑61 (MQ=255)
                                                   ggAACGGCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTc             <  1:2888825/71‑1 (MQ=255)
                                                    gAACGGCGATTCGCGAATCGATGGAGATTATTCAGGct                                             <  1:876001/38‑2 (MQ=38)
                                                     aaCGGCGATTCGCGAATCGATGGAGATTATTCAGGct                                             >  1:3459322/1‑36 (MQ=38)
                                                     aaCGGCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGc           <  1:1641416/71‑1 (MQ=255)
                                                     aaCGGCGATTCGCGAATCGATGGAGATTATTCAGGATAATGGCGCGa                                   >  1:2647512/1‑47 (MQ=37)
                                                      aCGGCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCt                           >  1:1367202/1‑54 (MQ=39)
                                                      aCGGCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCt          >  1:3202809/1‑71 (MQ=255)
                                                       cGGCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCtt         <  1:4295501/71‑2 (MQ=255)
                                                       cGGCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCtt         >  1:1341393/1‑70 (MQ=255)
                                                       cGGCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCtt         <  1:4522355/71‑2 (MQ=255)
                                                        ggCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGc           <  1:162912/68‑1 (MQ=255)
                                                        ggCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCtt         >  1:1754405/1‑69 (MQ=255)
                                                         gCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCt          >  1:646567/1‑68 (MQ=255)
                                                         gCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGa       >  1:2788483/1‑71 (MQ=255)
                                                         gCGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGa       >  1:3304369/1‑71 (MQ=255)
                                                          cGATTCGCGAATCGATGGAGATTATTCAGGCTAATgg                                        <  1:1714816/37‑1 (MQ=37)
                                                          cGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCt          >  1:3030462/1‑67 (MQ=255)
                                                          cGATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGa       >  1:4403973/1‑70 (MQ=255)
                                                           gATTCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATc     >  1:2647793/1‑71 (MQ=255)
                                                             ttCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTg                             <  1:1324806/45‑1 (MQ=38)
                                                             ttCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGt   <  1:517736/71‑1 (MQ=255)
                                                             ttCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGt   <  1:2687081/71‑1 (MQ=255)
                                                             ttCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGt   <  1:3742535/71‑1 (MQ=255)
                                                              tCGCGAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGt   >  1:4349746/1‑70 (MQ=255)
                                                               cgcgAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTc             >  1:866565/1‑59 (MQ=255)
                                                               cgcgAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCg    <  1:1977238/68‑1 (MQ=255)
                                                               cgcgAATCGATGGAGATTATTCAGGCTAATGGCGCGACGCTTGCTGGCGTATTGATTTCGCTTGATCGTc  <  1:2946924/70‑1 (MQ=255)
                                                                    |                                                                
GGTAGCGCGTTACAAGGACGCGTAATGCTGGTAGATGATGTGATCACCGCCGGAACGGCGATTCGCGAGTCGATGGAGATTATTCAGGCCAATGGCGCGACGCTTGCTGGCGTGTTGATTTCGCTCGATCGTC  >  AP009048/3824983‑3825115

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: