Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 464,884:1 +T 24.1% intergenic (+24/‑277) sdhB → / → sucA succinate dehydrogenase, FeS subunit/2‑oxoglutarate decarboxylase, thiamin‑requiring

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported464,8841.T24.1% 122.0 / 65.4 58intergenic (+24/‑277)sdhB/sucAsuccinate dehydrogenase, FeS subunit/2‑oxoglutarate decarboxylase, thiamin‑requiring
Reads supporting (aligned to +/- strand):  ref base . (31/13);  new base T (9/5);  total (40/18)
Fisher's exact test for biased strand distribution p-value = 7.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGCCATCGGCCATATCAAGTCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAACCAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTCATTACCCGT  >  Exported/464814‑464954
                                                                       |                                                                       
gcgcCATCGGCCATATCAAGTCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCa                                                                          <  1:5465554/71‑1 (MQ=255)
  gcCATCGGCCATATCAAGTCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAg                                                                      >  1:624876/1‑71 (MQ=255)
     aTCGGCCATATCAAGTCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGt                                                                   >  1:1090347/1‑71 (MQ=255)
     aTCGGCCATATCAAGTCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGt                                                                   >  1:35157/1‑71 (MQ=255)
        ggCCATATCAAGTCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCg                                                                 >  1:4829995/1‑70 (MQ=255)
        ggCCATATCAAGTCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGc                                                                <  1:1246320/71‑1 (MQ=255)
            atatCAAGTCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGCATCa                                                            <  1:3049843/71‑1 (MQ=255)
             tatCAAGTCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAg                                                           >  1:5031890/1‑71 (MQ=255)
               tCAAGTCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCaa                                                                        >  1:1826767/1‑57 (MQ=255)
               tCAAGTCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCaa                                                                        >  1:324804/1‑57 (MQ=255)
               tCAAGTCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAg                                                                      >  1:367953/1‑58 (MQ=255)
               tCAAGTCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAg                                                                      >  1:5551587/1‑58 (MQ=255)
               tCAAGTCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAg                                                                      >  1:440743/1‑58 (MQ=255)
               tCAAGTCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAg                                                                      >  1:3356972/1‑58 (MQ=255)
               tCAAGTCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGCAt                                                              >  1:4153199/1‑66 (MQ=255)
               tCAAGTCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGc                                                         <  1:5116807/71‑1 (MQ=255)
               tCAAGTCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGc                                                         <  1:4598662/71‑1 (MQ=255)
                cAAGTCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGc                                                         >  1:2353200/1‑70 (MQ=255)
                  aGTCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAAc                                                      >  1:1794064/1‑71 (MQ=255)
                  aGTCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAAc                                                      >  1:799971/1‑71 (MQ=255)
                      gATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAACCAg                                                   >  1:4378387/1‑70 (MQ=255)
                       aTGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAAc                                                      >  1:2290297/1‑66 (MQ=255)
                       aTGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAAc                                                      >  1:5328296/1‑66 (MQ=255)
                       aTGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAAc                                                      >  1:1624072/1‑66 (MQ=255)
                        tgttgtTGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAACCAGTGc                                                >  1:5400891/1‑71 (MQ=255)
                         gttgttGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAACCAGTGcc                                               >  1:2855310/1‑71 (MQ=255)
                          ttgttgCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAACCAGTg                                                 <  1:1834866/68‑1 (MQ=255)
                            gttgCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCg                                                                 >  1:1339669/1‑50 (MQ=255)
                            gttgCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCg                                                                 >  1:777925/1‑50 (MQ=255)
                               tacacttcgcccaTAAACCGTAGGCCTGATAAGACGCGCa                                                                          <  1:192668‑M1/27‑1 (MQ=255)
                               tacacttcgcccaTAAACCGTAGGCCTGATAAGACGCGCATAGCGTc                                                                  >  1:1920950‑M1/14‑47 (MQ=255)
                               tacacttcgcccaTAAACCGTAGGCCTGATAAGACGCGCATAGCGTc                                                                  >  1:4022236‑M1/14‑47 (MQ=255)
                               tacacttcgcccaTAAACCGTAGGCCTGATAAGACGCGCATAGCGTc                                                                  >  1:373375‑M1/14‑47 (MQ=255)
                               tacacttcgcccaTAAACCGTAGGCCTGATAAGACGCGCATAGCGTc                                                                  >  1:5000250‑M1/14‑47 (MQ=255)
                               gCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGc                                                                >  1:2531986/1‑48 (MQ=255)
                               gCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGc                                                                >  1:3950061/1‑48 (MQ=255)
                               gCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGc                                                                >  1:1913771/1‑48 (MQ=255)
                               gCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGc                                                                >  1:1694434/1‑48 (MQ=255)
                                 cacttcgcccgTAAGCCGTAGGCCTGATAAGACGCGCAAAGCGTCGCAt                                                              <  1:3447958‑M1/38‑1 (MQ=255)
                                 cacttcgcccaTAAACCGTAGGCCTGATAAGACGCGCATAGCGTCGc                                                                <  1:5248528‑M1/36‑1 (MQ=255)
                                 cacttcgcccaTAAACCGTAGGCCTGATAAGACGCGCATAGCGTCGc                                                                <  1:2213374‑M1/36‑1 (MQ=255)
                                 cacttcgcccaTAAACCGTAGGCCTGATAAGACGCGCATAGCGTCGc                                                                <  1:5409598‑M1/36‑1 (MQ=255)
                                 cacttcgcccaTAAACCGTAGGCCTGATAAGACGCGCATAGCGTCGCAt                                                              <  1:4947273‑M1/38‑1 (MQ=255)
                                   cttcgcccaTAAACCGTAGGACTGATAAGACGCGCATAGCGTCGCATCAGGCa                                                        >  1:2507995‑M1/10‑53 (MQ=255)
                                     tcgcccaTAAACCGTAGGCCTGATAAGACGCGCATAGCGTCGCAt                                                              <  1:309600‑M1/38‑1 (MQ=255)
                                     tAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAACCAGTGCCGGATGCGGCGTg                                   <  1:3840306/71‑1 (MQ=255)
                                     tAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAACCAGTGCCGGATGCGGCGTg                                   <  1:3085767/71‑1 (MQ=255)
                                      cgcccgTAAGCCGTAGGCCTGATAAGACGCGCAAAGCGTCGCAt                                                              <  1:276795‑M1/38‑1 (MQ=255)
                                       aTGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAACCAGTGCCGGATGCGGCGTGaa                                 <  1:2752619/71‑1 (MQ=255)
                                       aTGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAACCAGTGCCGGATGCGGCGTGaa                                 >  1:4074436/1‑71 (MQ=255)
                                           aTAAACCGTAGGCCTGATAAGACGCGCATAGCGTCGc                                                                >  1:98135‑M1/2‑37 (MQ=255)
                                           aTAAACCGTAGGCCTGATAAGACGCGCATAGCGTCGc                                                                >  1:2218124‑M1/2‑37 (MQ=255)
                                           aTAAACCGTAGGCCTGATAAGACGCGCATAGCGTCGc                                                                >  1:951124‑M1/2‑37 (MQ=255)
                                           aTAAACCGTAGGCCTGATAAGACGCGCATAGCGTCGc                                                                >  1:3048109‑M1/2‑37 (MQ=255)
                                                ccGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAACCAGTGCCGGATGCGGCGTGAACGCCTTATc                        <  1:2056247/71‑1 (MQ=255)
                                                           aTAAGACGCGCAAGCGTCGCATCAGGCAACCAGTGCCGGATGCGGCGTGAACGCGTTATCCGGCCTACAAg             >  1:5401865/1‑71 (MQ=255)
                                                           aTAAGACGCGCAAGCGTCGCATCAGGCAACCAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACa               <  1:4838075/69‑1 (MQ=255)
                                                           aTAAGACGCGCAAGCGTCGCATCAGGCAACCAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAg             >  1:4453097/1‑71 (MQ=255)
                                                           aTAAGACGCGCAAGCGTCGCATCAGGCAACCAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAg             >  1:351885/1‑71 (MQ=255)
                                                               gACGCGCAAGCGTCGCATCAGGCAACCAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTCAt         <  1:4442509/71‑1 (MQ=255)
                                                                aCGCGCAAGCGTCGCATCAGGCAACCAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTCAt         >  1:1394081/1‑70 (MQ=255)
                                                                aCGCGCAAGCGTCGCATCAGGCAACCAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTCAt         >  1:884615/1‑70 (MQ=255)
                                                                      aaGCGTCGCATCAGGCAACCAGTGCCGGATGCGGCGTGAACGCCTTATCCgg                     >  1:4186831/1‑52 (MQ=255)
                                                                      aaGCGTCGCATCAGGCAACCAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTCATTACCCGt  <  1:3415362/71‑1 (MQ=255)
                                                                       |                                                                       
GCGCCATCGGCCATATCAAGTCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAACCAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTCATTACCCGT  >  Exported/464814‑464954

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: