Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 758,462 G→A intergenic (+57/‑244) sdhB → / → sucA succinate:quinone oxidoreductase, iron‑sulfur cluster binding protein/subunit of E1(0) component of 2‑oxoglutarate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913758,4620GA93.8% 34.3 / ‑6.9 16intergenic (+57/‑244)sdhB/sucAsuccinate:quinone oxidoreductase, iron‑sulfur cluster binding protein/subunit of E1(0) component of 2‑oxoglutarate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (13/2);  minor base T (0/1);  total (13/3)
Fisher's exact test for biased strand distribution p-value = 1.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

TCCTGATTGATAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTCACAGCATCATGAACTGCGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAACCAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTCATTACCCGTAGGCCTGATAAGCGCAGCGCATCAGGCGTAACAAAGAAATGCAGGAAATCTTTAAAAACTGCCCCTGACACTAAGACAGTTTTTAAAGGTTCCTTCGCGAGCCACTACGTAGACAAGAGCTCGCAAGTGAACCCCGGCACGCACATCACTGTGCGTGGTAGTATCCACGGCGAAGTAAGCATAAAAAA  >  NC_000913/758230‑758687
                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                 
tCCTGATTGATAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTCACAGCATCATGAACTGCGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAACCAGTGCCGGATGCGACGTGAACGCCTTATCCgg                                                                                                                                                                                                                 >  1:139199/1‑251 (MQ=255)
                                                 gagggATGCATCCAGCGAACTCCGCTGTCACAGCATCATGAACTGCGTCAGTGGATGCCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAACCCGATGTTGGTGGGACGTAATGAGTAAACCGGAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAACCAGTGCTGGAGGCGTCGTGGACGCCTTATCCGGCCTCCAAGTCATTATCCGTAGGCCTGATAAGCGCAGCGCATCAGGCGt                                                                                                                                                                 <  1:323692/246‑1 (MQ=255)
                                                   gTGATGCATTCAGCGTATTCCGCTGTCACAGCATCATGAACTGCGTCAGTGTCTGTCAGAATGGCCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAACCAGTGCCGGATGCGAAGTGAACGCCTTATCCGGCCTACAAGTCATTACCCGTAGGCCTGATAAGCGCAGCGCATCAGGCGGAAc                                                                                                                                                              >  1:322598/1‑251 (MQ=255)
                                                                        gCTGTCACAGCATCATGAACTGCGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGATGTTGTTGCAACGTAATGAGTAAACCGTAGGCCTAATAAGACGCGCAAGCGTCGCATCAGGCAACCAGTGCCGGATGCGACGTGAACGCCTTATCCGGCCTACAAGTCATTACCCGTAGGCCTGATAAGCGCAGCGCATCAGGCGTAACAAAGAAATGCAGGAAATCtt                                                                                                                                          >  1:20223/1‑250 (MQ=255)
                                                                           gTCACAGCATCCTGAACTGCGTCAGCGTATGTACGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGATACGACGCGCAAGCGTCGCATCAGGCAACCAGTGCCGGATGCGACGTGAACGCCTTATCCGGCCTACAAGTCATTACCCGTAGGCCTGATAAGCGCAGCGCATCAGGCGTAACAAAGAAATGCAGGAAATCTTTaa                                                                                                                                       >  1:112914/1‑250 (MQ=255)
                                                                                 gCATCATGAACTGCGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAACCAGTGCCGGATGCGACGTGAACGCCTTATCCGGCCTACAAGTCATTACCCGTAGGCCTGATAAGCGCAGCGCATCAGGCGTAACAAAGAAATGCAGGAAATCTTTAAAAACTGc                                                                                                                                >  1:249768/1‑251 (MQ=255)
                                                                                  catcatGAACTGCGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGAGCCATCGGCCATATCAAGTCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAACCAGTGCCGGATGCGACGTGAACGCCTTATCCGGCCTACAAGTCATTACCCGTAGGcc                                                                                                                                                                                         >  1:126083/1‑193 (MQ=255)
                                                                                                                      gAACCCGACGCGCTCCATCGGCCATATCAAGTCGATTTTGGTGCAACGTAATGCGTAAACTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAACCAGTTCCGGAGGCGACGTGAACGCCTTATCCGGCCTTCAAGTCATGACCCGTAGGCCTGAGAAGCGCAGCGCATCAGGCGTAACAAAGAAATGCAGGAAATCTTTAAAAACTGCCCCTGACACTAAGACAGTTTTTAAAGGTTCCTTcg                                                                                             <  1:254122/249‑1 (MQ=255)
                                                                                                                                 gcgcCATCGGCCATATCAAGTCGATGTTGTTGCAACGTAATGCGTAAACCGTGGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAACCAGTGCCGGATGCGACGTGAACGCCTTACCCGGCCTACAAGTCATTACCCGTAGGCTTGATAAGCGCAGCGCATCAGGGGTAACAAAGAAATGCAGGAAATCTTTAAAAACTGCCCCTGACACTAAGACAGTTTTTAAAGGTTCCTTCGCGAGCCACTACGt                                                                                <  1:123496/251‑1 (MQ=255)
                                                                                                                                    ccATCGGACATATCAAGTCGATGTTGTTGCAACGTCATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAACCAGTGCCGGATGCGACGTGAACGCCTTATCCGGCCTACAAGTCATTACCCGTAGGCCTGATAAGCGCAGCGCATCAGGCGTAACAAAGAAATGCAGGAAATCTTTAAAAACTGCCCCTGACACTAAGACAGTTTTTAAAGGTTCCTTCGCGAGCCACTACGt                                                                                >  1:281127/1‑248 (MQ=255)
                                                                                                                                           ccATATCAAGTCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAACCAGTGCCGGATGCGACGTGAACGCCTTATCCGGCCTACAAGTCATTACCCGTAGGCCTGATAAGCGCAGCGCATCAGGCGTAACAAAGAAATGCAGGAAATCTTTAAAAACTGCCCCTGACACTAAGACAGt                                                                                                              >  1:125063/1‑211 (MQ=255)
                                                                                                                                                                     gTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAACCAGTGCCGGATGCGACGTGAACGCCTTATCCGGCCTACAAGTCATTACCCGTAGGCCTGATAAGCGCAGCGCATCAGGCGTAACAAAGAAATGCAGGAAATCTTTAAAAACTGCCCCTGACACTAAGACAGTTTTTAAAGGTTCCTTCGCGAGc                                                                                        >  1:23804/1‑207 (MQ=255)
                                                                                                                                                                                aCCGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAACCAGTGCCGGCTGCGACGTGACCGCCTTATCCGGCCTACAAGTCATTACCCGTCGGCCTGATAAGCGCAGCGCATCAGGCGTAACAAAGAAATGCAGGAAATCTTTAAAAACTGCCCCTGACACTAAGACAGTTTTTAAAGGTTCCTTCGCGAGCCACTCCGTAGACAAGAGCTCGCAAGTGAACCCCGGCACGCACATCACTGTGCGTg                                 >  1:242975/1‑251 (MQ=255)
                                                                                                                                                                                                               cATCAGGCAACCAGTGCCGGATGCGACGTGAACGCCTTATCCGGCCTACAAGTCATTACCCGTAGGCCTGATAAGCGCAGCGCATCAGGCGTAACAAAGAAATGCAGGAAATCTTTAAAAACTGCCCCTGACACTAAGACAGTTTTTAAAGGTTCCTTCGCGAGCCACTACGTAGACAAGAGCTCGCAAGTGAACCCCGGCACGCACATCACTGTGCGTGGTAGTATCCACGGCGAAGTCAGCATaaaaa   >  1:253198/1‑250 (MQ=255)
                                                                                                                                                                                                                    ggCAACCAGTGCCGGATGCGACGTGAACGCCTTATCCGGCCTACAAGTCATTACCCGTAGGCCTGATAAGCGCAGCGCATCAGGCGTAACAAAGAAATGCAGGAAATCTTTAAAAACTGCCCCTGACACTAAGACAGTTTTTAAAGGTTCCTTCGCGAGCCACTACGTAGACAAGAGCTCGCAAGTGAACCCCGGCACGCACATCACTGTGCGTGGTAGTATCCACGGCGAAGTAAGCATaaaaaa  >  1:225833/1‑246 (MQ=255)
                                                                                                                                                                                                                               ccGGATGCGACGTGAACGCCTTATCCGGCCTACAAGTCATTACCCGTAGGCCTGATCAGCGCAGCGCATAAGGCGTAACAACGACATGCAGGAAATCTTTcaaa                                                                                                                                     >  1:197952/1‑104 (MQ=37)
                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                 
TCCTGATTGATAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTCACAGCATCATGAACTGCGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAACCAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTCATTACCCGTAGGCCTGATAAGCGCAGCGCATCAGGCGTAACAAAGAAATGCAGGAAATCTTTAAAAACTGCCCCTGACACTAAGACAGTTTTTAAAGGTTCCTTCGCGAGCCACTACGTAGACAAGAGCTCGCAAGTGAACCCCGGCACGCACATCACTGTGCGTGGTAGTATCCACGGCGAAGTAAGCATAAAAAA  >  NC_000913/758230‑758687

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