Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 3,948,122 C→T intergenic (‑316/+370) aslA ← / ← hemY acrylsulfatase‑like enzyme/predicted protoheme IX synthesis protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6063,948,1220CT98.3% 213.4 / NA 60intergenic (‑316/+370)aslA/hemYacrylsulfatase‑like enzyme/predicted protoheme IX synthesis protein
Reads supporting (aligned to +/- strand):  ref base C (1/0);  major base T (38/19);  minor base G (1/0);  total (41/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACGGCAATTAATATATATGTTATTAATTGAATGAATTTTATTATTCATTATATATATGTGTAGAATCGTGCGCAGGAGAAATATTCACTCAGGAAGTTATTACTCAGGAAGCAAAGAGGATTACAGAATTATCTCATAACAAGTG  >  REL606/3948049‑3948194
                                                                         |                                                                        
aaCGGCAATTAATATATATGTTATTAATTGAATGAATTTTATTATTCATTATATATATGTGTAGAATCGTGCGTAg                                                                        <  2:804420/76‑1 (MQ=255)
  cGGCAATTAATATATATGTTATTAATTGAATGAATTTTATTATTCATTATATATATGTGTAGAATCGTGCGTAGga                                                                      <  1:819333/76‑1 (MQ=255)
     catttAATATATATGTTATTAATTGANTGAATTTTATTATTCATTATATATATGTGTAGAATCGTGCGTAGGAGaa                                                                   >  2:932489/4‑76 (MQ=255)
     cAATTAATATATATGTTATTAATTGAATGAATTTTATTATTCATTATATATATGTGTAGAATCGTGCGTAGGAGaa                                                                   <  1:1045343/76‑1 (MQ=255)
      aaTTAATATATATGTTATTAATTGAATGAATTTTATTATTCATTATATATATGTGTAGAATCGTGCGTAGGAGaaa                                                                  >  1:495758/1‑76 (MQ=255)
      aaTTAATATATATGTTATTAATTGAATGAATTTTATTATTCATTATATATATGTGTAGAATCGTGCGTAGGAGaaa                                                                  >  1:281246/1‑76 (MQ=255)
      aaTTAATATATATGTTATTAATTGAATGAATTTTATTATTCATTATATATATGTGTAGAATCGTGCGTAGGAGaaa                                                                  >  2:647375/1‑76 (MQ=255)
      aaTTAATATATATGTTATTAATTGAATGAATTTTAGTATTCATTATCTATATGTATAGAATCGTGCGTAGgagata                                                                  >  1:119759/1‑74 (MQ=255)
        ttAATATATATGTTATTAATTGAATGAATNTTATTATGCATTGTANATATGTGTAGTAGCGTGCGTAGGAGNAnta                                                                <  2:912369/76‑1 (MQ=255)
           atatatatGTTATTAATTGAATGAATTTTATTATTCATTATATATATGTGTAGAATCGTGCGTAGGAGAAATATTc                                                             >  1:672233/1‑76 (MQ=255)
           atatatatGTTATTAATTGAATGAATTTTATTATTCATTATATATATGTGTAGAATCGTGCGTAGGAGAAATATTc                                                             <  2:137707/76‑1 (MQ=255)
             atatatGTTATTAATTGAATGAATTTTACTATTCATTATATATATGTGTAGAATCGTGCGCGGGAGAAATATTCAc                                                           >  2:1051263/1‑76 (MQ=255)
              tatatGTTATTAATTGAATGAATTTTATTATTCATTATATATATGTGTAGAATCGTGCGAAGGAGAAATATTCACt                                                          >  2:441543/1‑76 (MQ=255)
               atatGTTATTAATTGAATGAATTTTATTATTCATTATATATATGTGTAGAATCGTGCGTAGGAGAAATATTCACTc                                                         >  1:620918/1‑76 (MQ=255)
                tatGTTATTAATTGAATGAATTTTATTATTCATTATATATATGTGTAGAATCGTGCGTAGGAGAAATATTCACTCa                                                        <  2:790731/76‑1 (MQ=255)
                 atGTTATTAATTGAATGAATTTTATTNTTCATTCTATATATGTGTAGAATCGTGCGTAGGAGAAATATTCACCCAg                                                       >  2:933676/1‑76 (MQ=255)
                 atGTTATTAATTGAATGAATTTTATTNTTCATTCTATATATGTGTAGAATCGTGCGTAGGAGAAATATTCACACAg                                                       >  2:933683/1‑76 (MQ=255)
                   gTTATTAATTGAATGAATTTTATTATTCATTATATATATGTGTAGAATCGTGCGTAGGAGAAATATTCACTCAGGa                                                     <  1:573372/76‑1 (MQ=255)
                    ttattaGTTGAATGAATTTTATTATTCATTATATATATGTGTAGAATCGTGCGTAGGAGAAATATTCACTCAGGaa                                                    >  1:158053/1‑76 (MQ=255)
                    ttattaATTGAATGAATTTTATTATTCATTATATATATGTGTAGAATCGTGCGTAGGAGAAATATTCACTCAGGaa                                                    >  1:580640/1‑76 (MQ=255)
                    ttattaAATGAATGAATTTTATTATTCATTATATATATGTGTAGAATCGTGCGTAGGAGAAATATTCACTCAGGaa                                                    >  1:709898/1‑76 (MQ=255)
                        taATTGAATGAATTTTATTATTCATTATATATATGTGTAGAATCGTGCGTAGGAGAAATATTCACTCAGGAAGtta                                                >  2:400124/1‑76 (MQ=255)
                           ttGAATGAATTTTATTATTCATTATATATATGTGTAGAATCGTGCGTAGGAGAAATATTCACTCAGGAAGttatta                                             >  1:230457/1‑76 (MQ=255)
                             gaatgaatTTTATTATTCATTATATATATGTGTAGAATCGTGCGTAGGAGAAATATTCACTCAGGAAGTTATTACt                                           >  1:477930/1‑76 (MQ=255)
                               atgaatTTTATTATTCATTATATATATGTGTAGAATCGTGCGTAGGAGAAATATTCACTCAGGAAGTTATTACTCa                                         <  2:580640/76‑1 (MQ=255)
                                  aatTTTATTATTCATTATATATATGTGTAGAATCGTGCGTAGGAGAAATATTCACTCAGGAAGTTATTACTCAGGa                                      <  1:665571/76‑1 (MQ=255)
                                   atTTTATTATTCATTATATATATGTGTAGAATCGTGCGTAGGAGAAATATTCACTCAGGAAGTTATTACTCAGGaa                                     >  2:64036/1‑76 (MQ=255)
                                       tattatTCATTATATATATGTGTAGAATCGTGCGTAGGAGAAATATTCACTCAGGAAGTTATTACTCAGGAAGCaa                                 >  1:487621/1‑76 (MQ=255)
                                       tattatTCATTATATATATGTGTAGAATCGTGCGTAGGAGAAATATTCACTCAGGAAGTTATTACTCAGGAAGCaa                                 >  1:487633/1‑76 (MQ=255)
                                       tattatTCATTATATATATGTGTAGAATCGTGCGTAGGAGAAATATTCACTCAGGAAGTTATTACTCAGGAAGCaa                                 >  1:324236/1‑76 (MQ=255)
                                         ttattCATTATATATATGTGTAGAATCGTGCGTAGGAGAAATATTCACTCAGGAAGTTATTACTCAGGAAGCAAag                               <  1:704192/76‑1 (MQ=255)
                                          tattCATTATATATATGTGTAGAATCGTGCGTAGGAGAAATCTTCACNCAGGAANTTATTACTCAGGAAGCAAagn                              >  1:907539/1‑75 (MQ=255)
                                          tattCATTATATATATGTGTAGAATCGTGCGTAGGAGAAATATTCACTCAGNAAGTTATTACTCAGGAAGCAAaga                              <  2:188065/76‑1 (MQ=255)
                                            ttCATTATATATATGTGTAGAATCGTGCGTAGGAGAAATATTCACTCAGGAAGTTATTACTCAGGAAGCAAAGAgg                            >  2:356750/1‑76 (MQ=255)
                                            ttCATTATATATATGTGTAGAATCGTGCGTAGGAGAAATATTCACTCAGGAAGTTATTACTCAGGAAGCAAAGAgg                            >  2:356769/1‑76 (MQ=255)
                                               aTTATATATATGTGTAGAATCGTGCGTAGGAGAAATATTCACTCAGGAAGTTATTACTCAGGAAGCAAAGAGGAtt                         >  2:1074324/1‑76 (MQ=255)
                                                ttATCTATATGTGTAGAATCGTGAGGAGTAGATNTATTCACTCAGGAAGTTATTACTNAGGAAGCAAAGAGGATTa                        >  2:109846/1‑76 (MQ=255)
                                                ttATATATATGTGTAGAATCGTGCGTAGGAGAAATATTCACTCAGGTAGTTATTACTCAGGAAGCAAAGAGGATTa                        <  1:1054562/76‑1 (MQ=255)
                                                ttATATATATGTGTAGAATCGTGCGTAGGAGAAATATTCACTCAGGAAGTTATTACTCAGGAAGCAAAGAGGATTa                        >  1:646802/1‑76 (MQ=255)
                                                ttATATATATGTGTAGAATCGTGCGTAGGAGAAATATTCACTCAGGAAGTTATTACTCAGGAAGCAAAGAGGATTa                        >  1:31054/1‑76 (MQ=255)
                                                 tatatataTGTGTAGAATCGTGCGTAGGAGAAATATTCGCTCAGGAAGTTATTACTCAGGAAGCAAAGAGGATTAc                       >  1:791324/1‑76 (MQ=255)
                                                   tatatatGTGTAGAATCGTGCGTAGGAGAAATATTCACTCAGGAAGTTATTACTCAGGAAGCAAAGAGGATTACAg                     >  2:1009676/1‑76 (MQ=255)
                                                     tatatGTGTAGAATCGTNCGTAGGAGAAATATTCACTCAGGAAGTTATTACNGAGGCAGCAGCGAGGATTACAGaa                   >  2:892971/1‑76 (MQ=255)
                                                       tatGTGCAGAATCGTGCGTAGGAGAAATATTCACTCAGGAAGTTATTACTCAGGAAGCAAAGAGGATTACAGAAtt                 >  1:165244/1‑76 (MQ=255)
                                                         tgtgTAGAATCGTGCGTAGGAGAAATATTCACTCAGGAAGTTATTACTCAGGAAGCAAAGAGGATTACAGAATTAt               >  2:782858/1‑76 (MQ=255)
                                                          gtgtAGAATCGTGCGTAGGAGAAATATTCACTCAGGAGGTTATTACTCAGGAAGCAAAGAGGATTACAGAATTAtc              >  1:578817/1‑76 (MQ=255)
                                                           tgtAGAATCGTGCGTAGGAGAAATATTCACTCAGGAAGTTATTACTCAGGAAGCAAAGAGGATTACAGAATTAtct             <  1:802888/76‑1 (MQ=255)
                                                           tgtAGAATCGTGCGTAGGAGAAATATTCACTCAGGAAGTTATTACTCAGGAACCAAAGAGGATTACAGAATTAtct             <  1:841955/76‑1 (MQ=255)
                                                             tAGAATCGTGCGTAGGAGAAATATTCACTCAGGAAGTTATTACTCAGGAAGCAAAGAGGATTACAGAATTATCTCa           >  1:648209/1‑76 (MQ=255)
                                                              aGAATCGTGCGTAGGAGAAATATTCACTCAGGAAGTTATTACTCAGGAAGCAAAGAGGATTACAGAATTATCTCAt          <  1:676033/76‑1 (MQ=255)
                                                              aGAATCGTGCGTAGGAGAAATATTCACTCAGGAAGTTATTACTCAGGAAGCAAAGAGGATTACAGAATTATCTCAt          >  1:414898/1‑76 (MQ=255)
                                                               gAATCGTGCGTAGGAGAAATATTCACTCAGGAAGTTATTACTCAGGAAGCAAAGAGGATTACAGAATTATCTCATa         >  1:700342/1‑76 (MQ=255)
                                                               gAATCGTGCGTAGGAGAAATATTCACTCAGGAAGTTATTACTCAGGAAGCAAAGAGGATTACAGAATTATCTCATa         >  1:700321/1‑76 (MQ=255)
                                                               gAATCGTGCGTAGGAGAAATATTCACTCAGGAAGTTATTACTCAGGAAGCAAAGAGGATTACAGAATTATCTCATa         <  2:594045/76‑1 (MQ=255)
                                                                aaTCGTGCGTAGGAGAAATATTCACTCAGGAAGTTATTACTCAGGAAGCAAAGAGGATTACAGAATTATCTCATaa        >  1:476980/1‑76 (MQ=255)
                                                                   cGTGCGTAGGAGAAATATTCACTCAGGAAGTTATTACTCAGGAAGCAAAGAGGATTACAGAATTATCTCATAACaa     >  2:591783/1‑76 (MQ=255)
                                                                   cGTGCGTAGGAGAAATATTCACTCAGGAAGTTATTACTCAGGAAGCAAAGAGGATTACAGAATTATCTCATAACaa     >  1:676748/1‑76 (MQ=255)
                                                                     tGCGTAGGAGAAATATTCACTCAGGAAGTTATTACTCAGGAAGCAAAGAGGATTACAGAATTATCTCATAACAAgt   <  1:585269/76‑1 (MQ=255)
                                                                     tGCGTAGGAGAAATATTCACTCAGGAAGTTATTACTCAGGAAGCAAAGAGGATTACAGAATTATCTCATAACAAgt   <  1:351109/76‑1 (MQ=255)
                                                                     tGCGTAGGAGAAATATTCACTCAGGAAGTTATTACTCAGGAAGCAAAGAGGATTACAGAATTATCTCATAACAAgt   <  1:884405/76‑1 (MQ=255)
                                                                      gcgtAGGAGAAATATTCACTCAGGAAGTTATTACTCAGGAAGCAAAGAGGATTACAGAATTATCTCATAACAAgtg  <  2:457874/76‑1 (MQ=255)
                                                                         |                                                                        
AACGGCAATTAATATATATGTTATTAATTGAATGAATTTTATTATTCATTATATATATGTGTAGAATCGTGCGCAGGAGAAATATTCACTCAGGAAGTTATTACTCAGGAAGCAAAGAGGATTACAGAATTATCTCATAACAAGTG  >  REL606/3948049‑3948194

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