Predicted mutation
evidence seq id position mutation annotation gene description
RA 1 991,761 C→T intergenic (‑394/+467) unknown ← / ← unknown hypothetical protein/Thiol:disulfide oxidoreductase associated with MetSO reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*1991,7610CT97.2% 88.7 / ‑5.9 36intergenic (‑394/+467)unknown/unknownhypothetical protein/Thiol:disulfide oxidoreductase associated with MetSO reductase
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (18/17);  total (19/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 = 1.00e+00

CGACTTGCTGACCGCTAAATAATAGCCAACTGGCTAAGCAGTTTAAATCCTTATAACTATCAAATAGTTCGATTTCTTTGGATGGTTTTCAAGGCACTATAACCATATGTACAATTACTGTCAAAAATGATTTTAAGGTTATCGGCACAAGTGCACTGGTAAGTAAGTAAAACTACAATACAGAGATCTGTAAACAGATCGCTTATCAAAATTTGATAAATTGTGATAGCTAGAGGAGCATTTATTAAAACAGCGTATTTATGAACAATGTGACTAAAAAGTTATTCTACGCAAGGTTGGTTTATTTACTTAAGCCGATTTTGCAAAATATTAACTGTTGAAGATAATATGATAATGAGTTATCTAACTATATGATGATATACCAAAATGGTTTATTTTATTTTGATGTAGAAATTGAATTTGATATAGTGCGTAGA  >  1/991523‑991959
                                                                                                                                                                                                                                              |                                                                                                                                                                                                      
cGACTTGCTGACCGCTAAATAATAGCCAACTGGCTAAGCAGTTTAAATCCTTATAACTATCAAATAGTTCGATTTCTTTGGATGGTTTTCAAGGCACTATAACTATATGTACAATTACTGTCAAAAATGATTTTAAGGTTATCGGCACAAGTGCACTGGTAAGTAAGTAAAACTACAATACAGAGATCTGTAAACAGATCGCTTATCAAAATTTGATAAATTGTGATAGCTAGAGGAGTATTTATTAAAAc                                                                                                                                                                                            <  1:137036/251‑1 (MQ=255)
                               ggCTAAGCAGTTTAAATCCTTATAACTATCAAATAGTTCGATTTCTTTGGATGGTTTTCAAGGCACTATAACTATATGTACAATTACTGTCAAAAATGATTTTAAGGTTATCGGCACAAGTGCACTGGTAAGTAAGTAAAACTACAATACAGAGATCTGTAAACAGATCGCTTATCAAAATTTGATAAATTGTGATAGCTAGAGGAGTATTTATTAAAAc                                                                                                                                                                                            <  1:593446/220‑1 (MQ=255)
                                              aTCCTTATAACTATCAAATAGTTCGATTTCTTTGGATGGTTTTCAAGGCACTATAACTATATGTACAATTACTGTCAAAAATGATTTTAAGGTTATCGGCACAAGTGCACTGGTAAGTAAGTAAAACTACAATACAGAGATCTGTAAACAGATCGCTTATCAAAATTTGATAAATTGTGATAGCTAGAGGAGTATTTATTAAAACAGCGTATTTATGAACAATGTGACTAAAAAGTTATTCTACGTAAggt                                                                                                                                              <  1:493623/251‑1 (MQ=255)
                                              aTCCTTATAACTATCAAATAGTTCGATTTCTTTGGATGGTTTTCAAGGCACTATAACTATATGTACAATTACTGTCAAAAATGATTTTAAGGTTATCGGCACAAGTGCACTGGTAAGTAAGTAAAACTACAATACAGAGATCTGTAAACAGATCGCTTATCAAAATTTGATAAATTGTGATAGCTAGAGGAGTATTTATTAAAACAGCGTATTTATGAACAATGTGACTAAAAAGTTATTCTACGTAAggt                                                                                                                                              <  1:644796/251‑1 (MQ=255)
                                                ccTTATAACTATCAAATAGTTCGATTTCTTTGGATGGTTTTCAAGGCACTATAACTGTATGTACAATTACTGTCAAAAATGATTTTAAGGTTATCGGCACAAGTGCACTGGTAAGTAAGTAAAACTACAATACAGAGATCTGTAAACAGATCGCTTATCAAAATTTGATAAATTGTGATAGCTAGAGGAGTATTTATTAAAACAGCGTATTTATGAACAATGTGACTAAAAAGTTATTCTACGTAAggt                                                                                                                                              >  1:207292/1‑249 (MQ=255)
                                                ccTTATAACTATCAAATAGTTCGATTTCTTTGGATGGTTTTCAAGGCACTATAACTATATGTACAATTACTGTCAAAAATGATTTTAAGGTTATCGGCACAAGTGCACTGGTAAGTAAGTAAAACTACAATACAGAGATCTGTAAACAGATCGCTTATCAAAATTTGATAAATTGTGATAGCTAGAGGAGtatttatt                                                                                                                                                                                                 >  1:754214/1‑198 (MQ=255)
                                                ccTTATAACTATCAAATAGTTCGATTTCTTTGGATGGTTTTCAAGGCACTATAACTATATGTACAATTACTGTCAAAAATGATTTTAAGGTTATCGGCACAAGTGCACTGGTAAGTAAGTAAAACTACAATACAGAGATCTGTAAACAGATCGCTTATCAAAATTTGATAAATTGTGATAGCTAGAGGAGTATTTATTAAAAc                                                                                                                                                                                            >  1:145881/1‑203 (MQ=255)
                                                        cTATCAAATAGTTCGATTTCTTTGGATGGTTTTCAAGGCACTATAACTATATGTACAATTACTGTCAAAAATGATTTTAAGGTTATCGGCACAAGTGCACTGGTAAGTAAGTAAAACTACAATACAGAGATCTGTAAACAGATCGCTTATCAAAATTTGATAAATTGTGATAGCTAGAGGAGTATTTATTAAAACAGCGTATTTATGAACAATGTGACTAAAAAGTTATTCTACGTAAGGTTGGtttattt                                                                                                                                    <  1:202339/251‑1 (MQ=255)
                                                          aTCAAATAGTTCGATTTCTTTGGATGGTTTTCAAGGCACTATAACTATATGTACAATTACTGTCAAAAATGATTTTAAGGTTATCGGCACAAGTGCACTGGTAAGTAAGTAAAACTACAATACAGAGATCTGTAAACAGATCGCTTATCAAAATTTGATAAATTGTGATAGCTAGAGGAGtatttatt                                                                                                                                                                                                 >  1:42690/1‑188 (MQ=255)
                                                                                    gTTTTCAAGGCACTATAACTATATGTACAATTACTGTCAAAAATGATTTTAAGGTTATCGGCACAAGTGCACTGGTAAGTAAGTAAAACTACAATACAGAGATCTGTAAACAGATCGCTTATCAAAATTTGATAAATTGTGATAGCTAGAGGAGTATTTATTAAAACAGCGTATTTATGAACAATGTGACTAAAAAGTTATTCTACGTAAggt                                                                                                                                              >  1:670643/1‑213 (MQ=255)
                                                                                     ttttCAAGGCACTATAACTATATGTACAATTACTGTCAAAAATGATTTTAAGGTTATCGGCACAAGTGCACTGGTAAGTAAGTAAAACTACAATACAGAGATCTGTAAACAGATCGCTTATCAAAATTTGATAAATTGTGATAGCTAGAGGAGTATTTATTAAAACAGCGTATTTATGAACAATGTGACTAAAAAGTTATTCTACGTAAGGTTGGTTTATTTACTTAAGCCGATTTTGCAAAATATTAAc                                                                                                        <  1:462281/250‑1 (MQ=255)
                                                                                     ttttCAAGGCACTATAACTATATGTACAATTACTGTCAAAAATGATTTTAAGGTTATCGGCACAAGTGCACTGGTAAGTAAGTAAAACTACAATACAGAGATCTGTAAACAGATCGCTTATCAAAATTTGATAAATTGTGATAGCTAGAGGAGTATTTATTAAAACAGCGTATTTATGAACAATGTGACTAAAAAGTTATTCTACGTAAGGTTGGTTTATTTACTTAAGCCGATTTTGCAAAATATTAACt                                                                                                       <  1:727808/251‑1 (MQ=255)
                                                                                                cTATAACTATATGTACAATTACTGTCAAAAATGATTTTAAGGTTATCGGCACAAGTGCACTGGTAAGTAAGTAAAACTACAATACAGAGATCTGTAAACAGATCGCTTATCAAAATTTGATAAATTGTGATAGCTAGAGGAGTATTTATTAAAAc                                                                                                                                                                                            <  1:219575/155‑1 (MQ=255)
                                                                                                cTATAACTATATGTACAATTACTGTCAAAAATGATTTTAAGGTTATCGGCACAAGTGCACTGGTAAGTAAGTAAAACTACAATACAGAGATCTGTAAACAGATCGCTTATCAAAATTTGATAAATTGTGATAGCTAGAGGAGTATTTATTAAAACAGCGTATTTATGAAc                                                                                                                                                                             >  1:443755/1‑170 (MQ=255)
                                                                                                  ataACTATATGTACAATTACTGTCAAAAATGATTTTAAGGTTATCGGCACAAGTGCACTGGTAAGTAAGTAAAACTACAATACAGAGATCTGTAAACAGATCGCTTATCAAAATTTGATAAATTGTGATAGCTAGAGGAGTATTTATTAAAAc                                                                                                                                                                                            >  1:26492/1‑153 (MQ=255)
                                                                                                   taACTATATGTACAATTACTGTCAAAAATGATTTTAAGGTTATCGGCACAAGTGCACTGGTAAGTAAGTAAAACTACAATACAGAGATCTGTAAACAGATCGCTTATCAAAATTTGATAAATTGTGATAGCTAGAGGAGTATTTATTAAAACAGCGTATTTATGAACAATGTGACTAAAAAGTTATTCTACGTAAGGTTGGTTTATTTACTTAAGCCGATTTTGCAAAATATTAACTGt                                                                                                     <  1:834970/239‑1 (MQ=255)
                                                                                                    aactatatGTACAATTACTGTCAAAAATGATTTTAAGGTTATCGGCACAAGTGCACTGGTAAGTAAGTAAAACTACAATACAGAGATCTGTAAACAGATCGCTTATCAAAATTTGATAAATTGTGATAGCTAGAGGAGTATTTATTAAAACAGCGTATTTATGAAc                                                                                                                                                                             >  1:531882/5‑166 (MQ=255)
                                                                                                                    aCTGTCAAAAATGATTTTAAGGTTATCGGCACAAGTGCACTGGTAAGTAAGTAAAACTACAATACAGAGATCTGTAAACAGATCGCTTATCAAAATTTGATAAATTGTGATAGCTAGAGGAGTATTTATTAAAACAGCGTATTTATGAACAATGTGACTAAAAAGTTATTCTACGTAAGGTTGGTTTATTTACTTAAGCCGATTTTGCAAAATATTAACt                                                                                                       >  1:91544/1‑220 (MQ=255)
                                                                                                                                 aTTTTAAGGTTATCGGCACAAGTGCACTGGTAAGTAAGTAAAACTACAATACAGAGATCTGTAAACAGATCGCTTATCAAAATTTGATAAATTGTGATAGCTAGAGGAGTAttt                                                                                                                                                                                                    >  1:730565/1‑114 (MQ=255)
                                                                                                                                         gTTATCGGCACAAGTGCACTGGTAAGTAAGTAAAACTACAATACAGAGATCTGTAAACAGATCGCTTATCAAAATTTGATAAATTGTGATAGCTAGAGGAGTATTTATTAAAACAGCGTATTTATGAACAATGTGACTAAAAAGTTATTCTACGTAAGGTTGGTTTATTTACTTAAGCCGATTTTGCAAAATATTAACt                                                                                                       <  1:808402/199‑1 (MQ=255)
                                                                                                                                                      gTGCACTGGTAAGTAAGTAAAACTACAATACAGAGATCTGTAAACAGATCGCTTATCAAAATTTGATAAATTGTGATAGCTAGAGGAGTATTTATTAAAACAGCGTATTTATGAACAATGTGACTAAAAAGTTATTCTACGTAAGGTTGGTTTATTTACTTAAGCCGATTTTGCAAAATATTAACTGTTGAAGat                                                                                              <  1:703871/195‑1 (MQ=255)
                                                                                                                                                                      gtaaAACTACAATACAGAGATCTGTAAACAGATCGCTTATCAAAATTTGATAAATTGTGATAGCTAGAGGAGTATTTATTAAAACAGCGTATTTATGAACAATGTGACTAAAAAGTTATTCTACGTAAGGTTGGTTTATTTACTTAAGCCGATTTTGCAAAATATTAACTGTTGAAGATAATATGAt                                                                                      >  1:389162/1‑187 (MQ=255)
                                                                                                                                                                                 aTACAGAGATCTGTAAACAGATCGCTTATCAAAATTTGATAAATTGTGATAGCTAGAGGAGTATTTATTAAAACAGCGTATTTATGAACAATGTGACTAAAAAGTTATTCTACGTAAGGTTGGTTTATTTACTTAAGCCGATTTTGCAAAATATTAACTGTTGAAGCTAATATGATAATGAGTTATCTAACTATATGATGATATACCAAAATGGTTTATTTTATTTTGATGTAGAAATTGAATTTGAtata           >  1:758282/1‑251 (MQ=255)
                                                                                                                                                                                 aTACAGAGATCTGTAAACAGATCGCTTATCAAAATTTGATAAATTGTGATAGCTAGAGGAGTATTTATTAAAACAGCGTATTTATGAACAATGTGACTAAAAAGTTATTCTACGTAAGGTTGGTTTATTTACTTAAGCCGATTTTGCAAAATATTAACTGTTGAAGATAATATGATAATGAGt                                                                               >  1:618154/1‑183 (MQ=255)
                                                                                                                                                                                   aCAGAGATCTGTAAACAGATCGCTTATCAAAATTTGATAAATTGTGATAGCTAGAGGAGTATTTATTAAAACAGCGTATTTATGAACAATGTGACTAAAAAGTTATTc                                                                                                                                                        <  1:128696/108‑1 (MQ=255)
                                                                                                                                                                                       agaTCTGTAAACAGATCGCTTATCAAAATTTGATAAATTGTGATAGCTAGAGGAGTATTTATTAAAACAGCGTATTTATGAACAATGTGACTAAAAAGTTATTCTACGTAAGGTTGGTTTATTTACTTAAGCCGATTTTGCAAAATATTAACTGTTGAAGATAATATGATAATGAGTTATCTAACTATATGATGATATACCAAAATGGTTTATTTTATTTTGATGTAGAAATTGAATTTGATATAGTGc       <  1:808041/249‑1 (MQ=255)
                                                                                                                                                                                        gaTCTGTAAACAGATCGCTTATCAAAATTTGATAAATTGTGATAGCTAGAGGAGtatttatt                                                                                                                                                                                                 >  1:139395/1‑62 (MQ=255)
                                                                                                                                                                                         aTCTGTAAACAGATCGCTTATCAAAATTTGATAAATTGTGATAGCTAGAGGAGTATTTATTAAAACAGCGTATTTATGAACAATGTGACTAAAAAGTTATTCTACGTAAGGTTGGTTTATTTACTTAAGCCGATTTTGCAAAATATTAACTGTTGAAGATAATATGATAATGAGTTATCTAACTATATGATGATATACCAAAATGGTTTATTTTATTTTGATGTAGAAATTGAATTTGATATAGTGCGt     <  1:873730/249‑1 (MQ=255)
                                                                                                                                                                                          tCTGTAAACAGATCGCTTATCAAAATTTGATAAATTGTGATAGCTAGAGGAGTATTTATTAAAACAGCGTATTTATGAACAATGTGACTAAAAAGTTATTCTACGTAAggt                                                                                                                                              >  1:583202/1‑111 (MQ=255)
                                                                                                                                                                                          tCTGTAAACAGATCGCTTATCAAAATTTGATAAATTGTGATAGCTAGAGGAGTATTTATTAAAACAGCGTATTTATGAACAATGTGACTAAAAAGTTATTCTACGTAAGGTTGGTTTATTTACTTAAGCCGATTTTGCAAAATATTAACTg                                                                                                      <  1:607989/151‑1 (MQ=255)
                                                                                                                                                                                          tCTGTAAACAGATCGCTTATCAAAATTTGATAAATTGTGATAGCTAGAGGAGTATTTATTAAAACAGCGTATTTATGAACAATGTGACTAAAAAGTTATTCTACGTAAGGTTGGTTTATTTACTTAAGCCGATTTTGCAAAATATTAACTGTTGAAGATAATATGATaa                                                                                    >  1:287519/1‑169 (MQ=255)
                                                                                                                                                                                          tCTGTAAACAGATCGCTTATCAAAATTTGATAAATTGTGATAGCTAGAGGAGTATTTATTAAAACAGCGTATTTATGAACAATGTGACTAAAAAGTTATTCTACGTAAGGTTGGTTTATTTACTTAAGCCGATTTTGCAAAATATTAACTGTTGAAGATAATATGATAATGAGTTATCTAACTATatgat                                                               >  1:741156/1‑190 (MQ=255)
                                                                                                                                                                                          tCTGTAAACAGATCGCTTATCAAAATTTGATAAATTGTGATAGCTAGAGGAGTATTTATTAAAACAGCGTATTTATGAACAATGTGACTAAAAAGTTATTCTACGTAAGGTTGGTTTATTTACTTAAGCCGATTTTGCAAAATATTAACTGTTGAAGATAATATGATAATGAGTTATCTAACTATATGATGATATACCAAAATGGTTTATTTTATTTTGATGTAGAAATTGAATTTGat              <  1:538900/239‑1 (MQ=255)
                                                                                                                                                                                          tCTGTAAACAGATCGCTTATCAAAATTTGATAAATTGTGATAGCTAGAGGAGTATTTATTAAAACAGCGTATTTATGAACAATGTGACTAAAAAGTTATTCTACGTAAGGTTGGTTTATTTACTTAAGCCGATTTTGCAAAATATTAACTGTTGAAGATAATATGATAATGAGTTATCTAACTATATGATGATATACCAAAATGGTTTATTTTATTTTGATGTAGAAATTGAATTTGATATAGTGCGTAGa  >  1:595933/1‑251 (MQ=255)
                                                                                                                                                                                          tCTGTAAACAGATCGCTTATCAAAATTTGATAAATTGTGATAGCTAGAGGAGCATTTATTAAAACAGCGTATTTATGAACAATGTGACTAAAAAGTTATTCTACGTAAGGTTGGTTTATTTACTTAAGCCGATTTTGCAAAATATTAACTGTTGAAGATAATATGATAATGAGTTATCTAAc                                                                       >  1:870771/1‑182 (MQ=255)
                                                                                                                                                                                                                      gATAAATTGTGATAGCTAGAGGAGTATTTATTAAAACAGCGTATTTa                                                                                                                                                                                  >  1:348602/1‑47 (MQ=255)
                                                                                                                                                                                                                          aaTTGTGATAGCTAGAGGAGTATTTATTAAAACAGCGTATTTATGAACAATGTGACTAAAAAGTTATTCTACGTAAGGTTGGTTTATTTACTTAAGCCGATTTTGCAAAATATTAACTGTTGAAGATAATATGATAATGAGt                                                                               >  1:656502/1‑142 (MQ=255)
                                                                                                                                                                                                                                   aGCTAGAGGAGTATTTATTAAAACAGCGTATTTATGAACAATGTGACTAAAAAGTTATTCTACGTAAGGTTGGTTTATTTACTTAAGCCGATTTTGCAAAATATTAACTGTTGAAGATAATATGATAATGAGTTATCTAAc                                                                       <  1:317342/141‑1 (MQ=255)
                                                                                                                                                                                                                                   aGCTAGAGGAGTATTTATTAAAACAGCGTATTTATGAACAATGTGACTAAAAAGTTATTCTACGTAAGGTTGGTTTATTTACTTAAGCCGATTTTGCAAAATATTAACTGTTGAAGATAATATGATAATGAGTTATCTAACTATATGATGatat                                                          >  1:331543/1‑154 (MQ=255)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                      
CGACTTGCTGACCGCTAAATAATAGCCAACTGGCTAAGCAGTTTAAATCCTTATAACTATCAAATAGTTCGATTTCTTTGGATGGTTTTCAAGGCACTATAACCATATGTACAATTACTGTCAAAAATGATTTTAAGGTTATCGGCACAAGTGCACTGGTAAGTAAGTAAAACTACAATACAGAGATCTGTAAACAGATCGCTTATCAAAATTTGATAAATTGTGATAGCTAGAGGAGCATTTATTAAAACAGCGTATTTATGAACAATGTGACTAAAAAGTTATTCTACGCAAGGTTGGTTTATTTACTTAAGCCGATTTTGCAAAATATTAACTGTTGAAGATAATATGATAATGAGTTATCTAACTATATGATGATATACCAAAATGGTTTATTTTATTTTGATGTAGAAATTGAATTTGATATAGTGCGTAGA  >  1/991523‑991959

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑