Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 3,490,530 (T)6→7 intergenic (‑106/‑84) glpE ← / → glpD thiosulfate sulfurtransferase/sn‑glycerol‑3‑phosphate dehydrogenase, aerobic, FAD/NAD(P)‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6063,490,5241.T85.7% 108.1 / 13.2 35intergenic (‑100/‑90)glpE/glpDthiosulfate sulfurtransferase/sn‑glycerol‑3‑phosphate dehydrogenase, aerobic, FAD/NAD(P)‑binding
Reads supporting (aligned to +/- strand):  ref base . (3/2);  new base T (15/15);  total (18/17)
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.63e-01

GGGCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAA  >  REL606/3490451‑3490593
                                                                          |                                                                       
gggCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGtt                                                                      >  2:494392/1‑76 (MQ=255)
       tttACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGtttttttca                                                              >  2:865897/1‑73 (MQ=255)
             gCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTTCAATGTTa                                                         <  2:867067/76‑1 (MQ=255)
               acaGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCt                                                      >  1:808146/1‑76 (MQ=255)
                  gTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTTCAATGTTACCTaa                                                    >  1:57915/1‑76 (MQ=255)
                  gTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTTCAATGCTACCTaa                                                    <  1:791280/76‑1 (MQ=255)
                        gTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTTCAATGTTACCTAAAGcgcg                                              >  1:440532/1‑76 (MQ=255)
                            aaGCGGCAGATAAACGCCATAATGTTATACATATCACTCTCAAATGTTTTTTTGAATGTTACCTAAAGCGCGATTc                                          >  1:864141/1‑76 (MQ=255)
                               cGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTTCAATGTTACCTAAAGCGCGATTCttt                                       >  2:946613/1‑76 (MQ=255)
                               cGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTTCAATGTTACCTAAAGCGCGATTCttt                                       <  1:647011/76‑1 (MQ=255)
                                  cAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTCTTTTTNAATGTTNCCTAAAGCGCGATTCTCTGcn                                    >  1:998708/1‑75 (MQ=255)
                                     aTAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAat                                 >  2:895531/1‑76 (MQ=255)
                                      tAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAata                                <  2:666372/76‑1 (MQ=255)
                                         actccATAATGTTATACATATCACTCTTAAATGTTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGt                             <  1:557256/73‑1 (MQ=255)
                                            ccATAATGTTATACATATCACGCTAAAATGTTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCg                          <  2:1005969/76‑1 (MQ=255)
                                              aTAATGTTATACATATCACTCTAAAATGTTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGAt                        <  1:214661/76‑1 (MQ=255)
                                              aTAATGTTATACATATCACTCTAAAATGTTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGAt                        >  2:1016246/1‑76 (MQ=255)
                                               tAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCGNAAGGGCGATTCTTTGCTAATATGTTCGATaa                      <  1:93720/76‑1 (MQ=255)
                                                 aTGTTATACATATCACTCTAAAATGTTTTTTTCGATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAAc                     <  1:776102/76‑1 (MQ=255)
                                                 aTGTTATACATATCACTCTAAAATGTTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAAc                     <  2:817520/76‑1 (MQ=255)
                                                 aTGTTATACATATCACTCTAAAATGTTCTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAAc                     <  1:494392/76‑1 (MQ=255)
                                                   gTTATACATATCACTCTAAAATGTTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGa                   >  2:571762/1‑76 (MQ=255)
                                                    ttATACATATCACTCTAAAATGTTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGaa                  >  1:878639/1‑76 (MQ=255)
                                                      ataCATATCACTCTAAAATGTTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATCACGAACa                >  2:766560/1‑76 (MQ=255)
                                                       taCATATCACTCTAAAATGTTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACAt               >  1:316704/1‑76 (MQ=255)
                                                        aCATATCACTCTAAAATGTTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACAtt              <  2:72392/76‑1 (MQ=255)
                                                         cATATCACTCTAAAATGTTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACAttt             >  2:9032/1‑76 (MQ=255)
                                                          atatCACTCTAAAATGTTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTa            >  1:891421/1‑76 (MQ=255)
                                                           tatCACTCTAAAATGTTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTAt           >  1:790328/1‑76 (MQ=255)
                                                            atCACTCTAAAATGTTTTTTTCAATGCTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATg          <  1:569685/76‑1 (MQ=255)
                                                                 tctAAAATGTTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAGTATGTTCGATAACGAACATTTATGAGCtt     <  1:1166825/76‑1 (MQ=255)
                                                                  ctAAAATGTTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCttt    <  2:360688/76‑1 (MQ=255)
                                                                  ctAAAATGTTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCttt    >  2:329197/1‑76 (MQ=255)
                                                                  ctAAAATGTTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCttt    >  1:898504/1‑76 (MQ=255)
                                                                   tAAAATGTTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTa   <  1:1119040/76‑1 (MQ=255)
                                                                   tAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTaa  >  2:500515/1‑76 (MQ=255)
                                                                    aaaaTGTTTTTTTCAATGTTACCTAAAGAGCGATTCTTTGCTAATATGTTCGATAACTAACATTTATGAGCTTTaa  <  2:1081511/76‑1 (MQ=255)
                                                                          |                                                                       
GGGCAAATTTACCGCACAGTTTACGTCGAAGCGGCAGATAAACGCCATAATGTTATACATATCACTCTAAAATGTTTTTTCAATGTTACCTAAAGCGCGATTCTTTGCTAATATGTTCGATAACGAACATTTATGAGCTTTAA  >  REL606/3490451‑3490593

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