Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 219,287 G→A 5.2% E41K (GAA→AAA)  ivy → inhibitor of vertebrate C‑lysozyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE219,2870GA 5.2% 312.0 / 5.1 97E41K (GAA→AAA) ivyinhibitor of vertebrate C‑lysozyme
Reads supporting (aligned to +/- strand):  ref base G (59/31);  new base A (2/3);  total (63/34)
Fisher's exact test for biased strand distribution p-value = 3.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.44e-01

Only 100 of 101 total aligned reads displayed.

GCCGCTCTGGTCATCGCCACCAGTGCAATGGCGCAGGATGATTTAACCATTAGCAGCCTTGCAAAGGGCGAAACCACCAAAGCTGCATTTAATCAGATGGTACAAGGGCATAAGCTGCCTGCCTGGGTGATGAAAGGCGG  >  minE/219218‑219357
                                                                     |                                                                      
gccgcTCTGGTCATCGCCACCAGTGCAATGGCGCAGGATGATTTAACCATTAGCAGCCTTGCAAAGGGCg                                                                        >  1:3426943/1‑70 (MQ=255)
gccgcTCTGGTCATCGCCACCAGTGCAATGGCGCAGGATGATTTAACCATTAGCAGCCTTGCAAAGGGCg                                                                        >  1:1241034/1‑70 (MQ=255)
 ccgcTCTGGTCATCGCCACCAGTGCAATGGCGCAGGATGATTTAACCATTAGCAGCCTTGCAAAGGGCGaa                                                                      <  1:3863988/71‑1 (MQ=255)
    ctctGGTCATCGCCACCAGTGCAATGGCGCAGGATGATTTAACCATTAGCAGCCTTGCAAAGGGCGAAacc                                                                   <  1:146681/71‑1 (MQ=255)
           cATCGCCACCAGTGCAATGGCGCAGGATGATTTAACCATTAGCAGCCTTGCAAAGGGCGAAACCACCAAAg                                                            >  1:2197563/1‑71 (MQ=255)
            aTCGCCACCAGTGCAATGGCGCAGGATGATTTAACCATTAGCAGCCTTGCAAAGGGCGAAACCACCAAAGc                                                           >  1:2766782/1‑71 (MQ=255)
            aTCGCCACCAGTGCAATGGCGCAGGATGATTTAACCATTAGCAGCCTTGCAAAGGGCAAAACCACCAAAGc                                                           >  1:2258834/1‑71 (MQ=255)
             tCGCCACCAGTGCAATGGCGCAGGATGATTTAACCATTAGCAGCCTTGCAAAGGGCGAAACCACCAAAGc                                                           <  1:3218192/70‑1 (MQ=255)
               gCCACCAGTGCAATGGCGCAGGATGATTTAACCATTAGCAGCCTTGCAAAGGGCGAAACCACCAAAGCTg                                                         >  1:1525711/1‑70 (MQ=255)
                   acaGTGCAATGGCGCAGGATGATTTAACCATTAGCAGCCTTGCAAAGGGCGAAACCACCAAAGCTGCAttt                                                    <  1:3107246/70‑1 (MQ=255)
                    caGTGCAATGGCGCAGGATGATTTAACCATTAGCAGCCTTGCAAAGGGCGAAACCACCAAAGCTGCATTTa                                                   <  1:2921729/71‑1 (MQ=255)
                    caGTGCAATGGCGCAGGATGATTTAACCATTAGCAGCCTTGCAAAGGGCAAAACCACCAAAGCTGCATTTa                                                   <  1:1772975/71‑1 (MQ=255)
                      gTGCAATGGCGCAGGATGATTTAACCATTAGCAGCCTTGCAAAGGGCGAAACCACCAAAGCTGCATTTaa                                                  <  1:855139/70‑1 (MQ=255)
                      gTGCAATGGCGCAGGATGATTTAACCATTAGCAGCCTTGCAAAGGGCGAAACCACCAAAGCTGCATTTAAt                                                 >  1:3439323/1‑71 (MQ=255)
                         cAATGGCGCAGGATGATTTATCCATTAGCAGCCTTGCAAAGGGCGAAACCACCAAAGCTGCATTTAATCAg                                              <  1:2923005/71‑1 (MQ=255)
                         cAATGGCGCAGGATGATTTAACCATTAGCAGCCTTGCAAAGGGCGAAACCACCAAAGCTGCATTTAAt                                                 >  1:2553555/1‑68 (MQ=255)
                         cAATGGCGCAGGATGATTTAACCATTAGCAGCCTTGCAAAGGGCGAAACCACCAAAGCTGCATTTAAt                                                 >  1:1404319/1‑68 (MQ=255)
                         cAATGGCGCAGGATGATTTAACCATTAGCAGCCTTGCAAAGGGCGAAACCACCAAAGCTGCATTTAATCa                                               <  1:2957105/70‑1 (MQ=255)
                         cAATGGCGCACGATGATTTAACCATTAGCAGCCTTGCAAAGGGCGAAACCACCAAAGCTGc                                                        >  1:1253775/1‑61 (MQ=255)
                           aTGGCGCAGGATGATTTAACCATTAGCAGCCTTGCAAAGGGCGAAACCACCAAAGCTGCATTTAATCAGAt                                            >  1:3343001/1‑71 (MQ=255)
                            tGGCGCAGGATGATTTAACCATTAGCAGCCTTGCAAAGGGCGAAACCACCAAAGCTGCATTTa                                                   <  1:1848856/63‑1 (MQ=255)
                            tGGCGCAGGATGATTTAACCATTAGCAGCCTTGCAAAGGGCGAAACCACCAAAGCTGCATTTAATCAGATg                                           >  1:260425/1‑71 (MQ=255)
                             ggCGCAGGATGATTTAACCATTAGCAGCCTTGCAAAGGGCTAAACCACCAAAGCTGCATTTAATCAGATgg                                          >  1:2961580/1‑71 (MQ=255)
                             ggCGCAGGATGATTTAACCATTAGCAGCCTTGCAAAGGGCGAAACCACCAAAGCTGCATTTAATCAGATgg                                          >  1:2916240/1‑71 (MQ=255)
                              gcgcAGGATGATTTAACCATTAGCAGCCTTGCAAAGGGCGAAACCACCAAAGCTGc                                                        >  1:2706765/1‑56 (MQ=255)
                               cgcAGGATGATTTAACCATTAGCAGCCTTGCAAAGGGCGAAACCACCAAAg                                                            >  1:419819/1‑51 (MQ=255)
                               cgcAGGATGATTTAACCATTAGCAGCCTTGCAAAGGGCGAAACCACCAAAGCTGCATTTAATCAGATgg                                          <  1:4052397/69‑1 (MQ=255)
                                gcAGGATGATTTAACCATTAGCAGCCTTGCAAAGGGCGAAACCACCAAAGCTGCATTTAATCAGATGGTa                                        >  1:591345/1‑70 (MQ=255)
                                 cAGGATGATTTAACCATTAGCAGCCTTGCAAAGGGCGAAACCACCAAAGCTGCATTTAATCAGATGGTACa                                      <  1:1113424/71‑1 (MQ=255)
                                  aGGATGATTTAACCATTAGCAGCCTTGCAAAGGGCGAAACCACCAAAGCTGCATTTAATCAGATGGTAc                                       >  1:1705458/1‑69 (MQ=255)
                                    tatgatTTAACCATTAGCAGCCTTGCAAAGGGCGAAACCACCaa                                                              <  1:320448/43‑1 (MQ=255)
                                    gatgatTTAACCATTAGCAGCCTTGCAAAGGGCGAAACCACCAAAGCTGCATTTaa                                                  >  1:2589117/1‑56 (MQ=255)
                                    gatgatTTAACCATTAGCAGCCTTGCAAAGGGCGAAACCACCAAAGCTGCATTTaa                                                  >  1:2584175/1‑56 (MQ=255)
                                    gatgatTTAACCATTAGCAGCCTTGCAAAGGGCGAAACCACCAAAGCTGCATTTaa                                                  >  1:2357068/1‑56 (MQ=255)
                                    gatgatTTAACCATTAGCAGCCTTGCAAAGGGCGAAACCACCAAAGCTGCATTTa                                                   <  1:1144572/55‑1 (MQ=255)
                                    gatgatTTAACCATTAGCAGCCTTGCAAAGGGCGAAACCACCAAAGCTGCATTTAATCAGATGGTAc                                       <  1:1208260/67‑1 (MQ=255)
                                     atgatTTAACCATTAGCAGCCTTGCAAAGGGCGAAACCACCAAAGCTGc                                                        <  1:1066398/49‑1 (MQ=255)
                                      tgatTTAACCATTAGCAGCCTTGCAAAGGTCGAAACCACCAAAGCTGCATTTAATCAGAt                                            <  1:1394363/60‑1 (MQ=255)
                                      tgatTTAACCATTAGCAGCCTTGCAAAGGGCGAAacca                                                                  >  1:451090/1‑38 (MQ=255)
                                      tgatTTAACCATTAGCAGCCTTGCAAAGGGCGAAACCACCAAAGCTGCATTTAAt                                                 >  1:1638491/1‑55 (MQ=255)
                                      tgatTTAACCATTAGCAGCCTTGCAAAGGGCGAAACCACCAAAGCTGCATTTAATCAGATGGTACaa                                     >  1:2352150/1‑67 (MQ=255)
                                      tgatTTAACCATTAGCAGCCTTGCAAAGGGCGAAACCACCAAAGCTGCATTTAATCAGATGGTACaa                                     >  1:3316423/1‑67 (MQ=255)
                                      tgatTTAACCATTAGCAGCCTTGCAAAGGGCGAAACCACCAAAGCTGCATTTAATCAGATGGTACaa                                     >  1:882559/1‑67 (MQ=255)
                                      tgatTTAACCATTAGCAGCCTTGCAAAGGGCGAAACCACCAAAGCTGCATTTAATCAGATGGTACAAggg                                  >  1:2170133/1‑70 (MQ=255)
                                      tgatTTAACCATTAGCAGCCTTGCAAAGGGCGAAACCACCAAAGCTGCATTTAATCAGATGGTACAAggg                                  >  1:1375098/1‑70 (MQ=255)
                                      tgatTTAACCATTAGCAGCCTTGCAAAGGGCGAAACCACCAAAGCTGCATTTAATCAGATGGTACAAggg                                  >  1:1838266/1‑70 (MQ=255)
                                      tgatTTAACCATTAGCAGCCTTGCAAAGGGCGAAACCACCAAAGCTGCATTTAATCAGATGGTACAAggg                                  >  1:528870/1‑70 (MQ=255)
                                      tgatTTAACCATTAGCAGCCTTGCAAAGGGCGAAACCACCAAAGCTGCATTTAATCAGATGGTACAAGGGc                                 <  1:3351086/71‑1 (MQ=255)
                                      tgatTAAACCATTAGCAGCCTTGCAAAGGGCCAAacc                                                                   >  1:3633873/1‑37 (MQ=38)
                                       gatTTAACCATTAGCAGCCTTGCAAAGGGCGAAACCACCAAAGCTGCATTTAATCa                                               <  1:1404202/56‑1 (MQ=255)
                                       gatTTAACCATTAGCAGCCTTGCAAAGGGCGAAACCACCAAAGCTGCATTTAATCAGATGGTACAAGGGCa                                >  1:3887712/1‑71 (MQ=255)
                                        atTTAACCATTAGCAGCCTTGCAAAGGGCGAAACCACCAAAGCTGCATTTAATCAGATGGTACAAGGGCAt                               <  1:535881/71‑1 (MQ=255)
                                         tttAACCATTAGCAGCCTTGCAAAGGGCGAAACCACCAAAGCTGCATTTAATCAGATGGTACAAGGGCAt                               >  1:2705607/1‑70 (MQ=255)
                                          ttAACCATTAGCAGCCTTGCAAAGGGCGAAACCACCAAAGCTGCATTTAATCAGATGGTACAAGGGCAt                               >  1:3975385/1‑69 (MQ=255)
                                             aCCATTAGCAGCCTTGCAAAGGGCGAAACCACCAAAGCTGCATTTAATCAGATGGTACAAGGGCATAAGc                           <  1:3421482/70‑1 (MQ=255)
                                              ccATTAGCAGCCTTGCAAAGGGCGAAACCACCAAAGCTGCATTTAATCAGATGGTACAAGGGCATAAGctg                         >  1:1597184/1‑71 (MQ=255)
                                               cATTAGCAGCCTTGCAAAGGGCGAAACCACCAAAGCTGCATTTAATCAGATGGTACAAGGGCATaa                             <  1:3233345/66‑1 (MQ=255)
                                                 ttAGCAGCCTTGCAAAGGGCGAAACCACCAAAGCTGCATTTAATCAg                                              <  1:2020671/47‑1 (MQ=255)
                                                  tAGCAGCCTTGCAAAGGGCGAAACCACCAAAGCTGCATTTAATCAGATGGTACAAGGGCATAAGctgcctg                     >  1:1233217/1‑71 (MQ=255)
                                                   agcagcCTTGCAAAGGGCGAAACCACCAAAGCTGCATTTAATCAGATGGTACAAGGGCATAAGctgcctgc                    >  1:2102721/1‑71 (MQ=255)
                                                   agcagcCTTGCAAAGGGCGAAACCACCAAAGCTGCATTTAATCAGATGGTACAAGGGCATAAGctgcctg                     >  1:3843011/1‑70 (MQ=255)
                                                    gcagcCTTGCAAAGGGCGAAACCACCAAAGCTGCATTTAATCAGATGGTACAAGGGCATAAGctgcctgc                    <  1:3153756/70‑1 (MQ=255)
                                                      agcCTTGCAAAGGGCGAAACCACCAAAGCTGCATTTAAt                                                 >  1:3336411/1‑39 (MQ=255)
                                                      agcCTTGCAAAGGGCAAAACCACCAAAGCTGCATTTAATCAGATGGTACAAGGGCATAAGCTgcctgcct                  <  1:3171465/70‑1 (MQ=255)
                                                       gcCTTGCAAAGGGCGAAACCACCAAAGCTGCATTTAATCAGATGGTACAAGGGCATAAGCTGCCTGCCTgg                <  1:109592/71‑1 (MQ=255)
                                                       gcCTTGCAAAGGGCGAAACCACCAAAGCTGCATTTAATCAGATGGTACAAGGGCATAAGCTGCCTGCCTgg                <  1:503536/71‑1 (MQ=255)
                                                        ccTTGCAAAGGGCGAAACCACCAAAGCTGCATTTAATCa                                               <  1:2059686/39‑1 (MQ=255)
                                                        ccTTGCAAAGGGCGAAACCACCAAAGCTGCATTTAATCa                                               <  1:2353255/39‑1 (MQ=255)
                                                         cTTGCAAAGGGCGAAACCACCAAAGCTGCATTTAATCAGATGGTACAAGGGCATAAGCTGCCTGCCTGGGt              >  1:713927/1‑71 (MQ=255)
                                                         cTTGCAAAGGGCGAAACCACCAAAGCTGCATTTAATCAGATGGTACAAGGGCATAAGCTGCCTGCCTGGGt              >  1:2619488/1‑71 (MQ=255)
                                                         cTTGCAAAGGGCGAAACCACCAAAGCTGCATTTAATCAGATGGTACAAGGGCATAAGCTGCCTGCCTGGGt              >  1:1724179/1‑71 (MQ=255)
                                                         cTTGCAAAGGGCGAAACCACCAAAGCTGCATTTAATCAGATGGTACAAGGGCATAAGCTGCCTGCCTGGGt              >  1:2088174/1‑71 (MQ=255)
                                                          ttGCAAAGGGCGAAACCACCAAAGCTGCATTTAAt                                                 >  1:3383368/1‑35 (MQ=255)
                                                          ttGCAAAGGGCGAAACCACCAAAGCTGCATTTAAt                                                 >  1:2655535/1‑35 (MQ=255)
                                                          ttGCAAAGGGCGAAACCACCAAAGCTGCATTTAAt                                                 >  1:1066897/1‑35 (MQ=255)
                                                          ttGCAAAGGGCGAAACCACCAAAGCTGCATTTAAt                                                 >  1:1295949/1‑35 (MQ=255)
                                                          ttGCAAAGGGCGAAACCACCAAAGCTGCATTTAATCAGATgg                                          <  1:778787/42‑1 (MQ=255)
                                                          ttGCAAAGGGCGAAACCACCAAAGCTGCATTTAATCAGATGGTACAAGGGCATAAGCTGCCTGCCTGGGt              >  1:1252561/1‑70 (MQ=255)
                                                          ttGCAAAGGGCAAAACCACCAAAGCTGCATTTAATCAGATGGTACAAGGGCATAAGCTgcctgcct                  <  1:3469223/66‑1 (MQ=255)
                                                           tGCAAAGGGCGAAACCACCAAAGCTGCATTTAATCAGATGGTACAAGGGCATAAGc                           >  1:3207416/1‑56 (MQ=255)
                                                           tGCAAAGGGCGAAACCACCAAAGCTGCATTTAATCAGATGGTACAAGGGCATAAGc                           >  1:2520757/1‑56 (MQ=255)
                                                            gCAAAGGGCGAAACCACCAAAGCTGCATTTAATCAg                                              >  1:2533551/1‑36 (MQ=255)
                                                            gCAAAGGGCGAAACCACCAAAGCTGCATTTAATCAGAt                                            >  1:1021744/1‑38 (MQ=255)
                                                            gCAAAGGGCGAAACCACCAAAGCTGCATTTAATCAGATGGTAc                                       <  1:2314067/43‑1 (MQ=255)
                                                            gCAAAGGGCGAAACCACCAAAGCTGCATTTAATCAGATGGTACAAGGGCATAAGc                           >  1:2510266/1‑55 (MQ=255)
                                                            gCAAAGGGCGAAACCACCAAAGCTGCATTTAATCAGATGGTACAAGGGCATAAGc                           >  1:2037755/1‑55 (MQ=255)
                                                            gCAAAGGGCGAAACCACCAAAGCTGCATTTAATCAGATGGTACAAGGGCATAAGCTGCCTGCCTGGGtgat           >  1:2488599/1‑71 (MQ=255)
                                                                aGGGCGAAACCACCAAAGCTGCATTTAATCAGATGGTACaa                                     >  1:3921421/1‑41 (MQ=255)
                                                                aGGGCGAAACCACCAAAGCTGCATTTAATCAGATGGTACAAGGGCATAAGCTGCCTGCCTGGGtgatg          >  1:3373117/1‑68 (MQ=255)
                                                                  ggCAAAACCACCAAAGCTGCATTTAATCAGATGGTACAAGCGCATAAGctg                         >  1:1678763/1‑51 (MQ=255)
                                                                   gCGAAACCACCAAAGCTGCATTTAATCAGATGGTACAAg                                    >  1:3692395/1‑39 (MQ=255)
                                                                   gCGAAACCACCAAAGCTGCATTTAATCAGATGGTACAAGGGCATAAGctgc                        >  1:828525/1‑51 (MQ=255)
                                                                   gCGAAACCACCAAAGCTGCATTTAATCAGATGGTACAAGGGCATAAGctgc                        >  1:1130343/1‑51 (MQ=255)
                                                                   gCGAAACCACCAAAGCTGCATTTAATCAGATGGTACAAGGGCATAAGCTGCCTGCCTGGGTGATGAAAGGc    <  1:452723/71‑1 (MQ=255)
                                                                   gCGAAACCACCAAAGCTGCATTTAATCAGATGGTACAAGGGCATAAGCTGCCTGCCTGGGTGATGAAAGGc    <  1:2042199/71‑1 (MQ=255)
                                                                    cGAAACCACCAAAGCTGCATTTAATCAGATGGTACAAGGGCa                                >  1:2394360/1‑42 (MQ=255)
                                                                    cGAAACCACCAAAGCTGCATTTAATCAGATGGTACAAGGGCATAAGCTGCCTGCCTGGGTGATGAAAGGCg   >  1:3123829/1‑71 (MQ=255)
                                                                    cGAAACCACCAAAGCTGCATTTAATCAGATGGTACAAGGGCATAAGCTGCCTGCCTGGGTGATGAAAGGCg   >  1:1700763/1‑71 (MQ=255)
                                                                     gAAACCACCAAAGCTGCATTTAATCAGATGGTACAAGGGCATAAGCTGCCTGCCTGGGTGATGAAAGGCgg  <  1:3639874/71‑1 (MQ=255)
                                                                     gAAACCACCAAAGCTGCATTTAATCAGATGGTACAAGGGCATAAGCTGCCTGCCTGGGTGATGAAAGGCgg  >  1:906997/1‑71 (MQ=255)
                                                                     |                                                                      
GCCGCTCTGGTCATCGCCACCAGTGCAATGGCGCAGGATGATTTAACCATTAGCAGCCTTGCAAAGGGCGAAACCACCAAAGCTGCATTTAATCAGATGGTACAAGGGCATAAGCTGCCTGCCTGGGTGATGAAAGGCGG  >  minE/219218‑219357

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: