Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	756,643	T→A	14.0%	I246N (ATT→AAT)	sdhA →	succinate dehydrogenase, flavoprotein subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	756,643	0	T	A	14.0%	116.6 / 11.2	43	I246N (ATT→AAT)	sdhA	succinate dehydrogenase, flavoprotein subunit
Reads supporting (aligned to +/- strand):  ref base T (19/18);  new base A (3/3);  total (22/21)
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.34e-01

ACCGTGCTGGCGACTGGCGGAGCAGGGCGTATTTATCAGTCCACCACCAACGCCCACATTAACACCGGCGACGGTGTCGGCATGGCTATCCGTGCCGGCGTACCGGTGCAGGATATGGAAATGTGGCAGTTCCACCCGACCGGCATTGCCGGTGCGGGCGTACTGGTCACCGAAGGTTGCCGTGGTGAAGGCGGTTATCTGCTGAACAAACATGGCGAACGTTTTATGGAGCGTTATGCGCCGAACGCCAAAGACCTGGCGGGCCGTGACGTGGTTGCGCGTTCCATCATGAT  >  NC_000913/756498‑756790                                                                                                                                                  |                                                                                                                                                    aCCGTGCTGGCGACTGGCGGAGCAGGGCGTATTTATCAGTCCACCACCAACGCCCACATTAACACCGGCGACGGTGTCGGCATGGCTATCCGTGCCGGCGTACCGGTGCAGGATATGGAAATGTGGCAGTTCCACCCGACCGGCAATGc                                                                                                                                                  <  2:250425/149‑1 (MQ=255)  ccGTGCTGGCGACTGGCGGAGCAGGGCGTATTTATCAGTCCACCACCAACGCCCACATTAACACCGGCGACGGTGTCGGCATGGCTATCCGTGCCGGCGTACCGGTGCAGGATATGGAAATGTGGCAGTTCCACCCGACCGGCATTGcc                                                                                                                                                 >  2:140609/1‑149 (MQ=255)   cGTGCTGGCGACTGGCGGAGCAGGGCGTATTTATCAGTCCACCACCAACGCCCACATTAACACCGGCGACGGTGTCGGCATGGCTATCCGTGCCGGCGTACCGGTGCAGGATATGGAAATGTGGCAGTTCCACCCGACCGGCATTGCCg                                                                                                                                                <  4:305130/149‑1 (MQ=255)       cTGGCGACTGGCGGAGCAGGGCGTATTTATCAGTCCACCACCAACGCCCACATTAACACCGGCGACGGTGTCGGCATGGCTATCCGTGCCGGCGTACCGGTGCAGGATATGGAAATGTGGCAGTTCCACCCGACCGGCATTGCCGGTGc                                                                                                                                            <  3:8716/149‑1 (MQ=255)            gACTGGCGGAGCAGGGCGTATTTATCAGTCCACCACCAACGCCCACATTAACACCGGCGACGGTGTCGGCATGGCTATCCGTGCCGGCGTACCGGTGCAGGATATGGAAATGTGGCAGTTCCACCCGACCGGCATTGCCGGTGCGGGCg                                                                                                                                       >  4:54041/1‑149 (MQ=255)                 gCGGAGCAGGGCGTATTTATCAGTCCACCACCAACGCCCACATTAACACCGGCGACGGTGTCGGCATGGCTATCCGTGCCGGCGTACCGGTGCAGGATATGGAAATGTGGCAGTTCCACCCGACCGGCATTGCCGGTGCGGGCGTACTg                                                                                                                                  <  1:140609/149‑1 (MQ=255)                      gCAGGGCGTATTTATCAGTCCACCACCAACGCCCACATTAACACCGGCGACGGTGTCGGCATGGCTATCCGTGCCGGCGTACCGGTGCAGGATATGGAAGTGTGGCAGTTCCACTCGACCGGCATTGCCGGTGCGGGCGTACTGGTCAc                                                                                                                             <  4:18130/149‑1 (MQ=255)                          ggCGTATTTATCAGTCCACCACCAACGCCCACATTAACACCGGCGACGGTGTCGGCATGGCTATCCGTGCCGGCGTACCGGTGCAGGATATGGAAATGTGGCAGTTCCACCCGACCGGCATTGCCGGTGCGGGCGTACTGGTCACCGaa                                                                                                                         >  3:283719/1‑149 (MQ=255)                           gCGTATTTATCAGTCCACCACCAACGCCCACATTAACACCGGCGACGGTGTCGGCATGGCTATCCGTGCCGGCGTACCGGTGCAGGATATGGAAATGTGGCAGTTCCACCCGACCGGCATTGCCGGTGCGGGCGTACTGGTCACCGAAg                                                                                                                        >  4:321753/1‑149 (MQ=255)                              tATTTATCAGTCCACCACCAACGCCCACATTAACACCGGCGACGGTGTCGGCATGGCTATCCGTGCCGGCGTACCGGTGCAGGATATGGAAATGTGGCAGTTCCACCCGACCGGCATTGCCGGTGCGGGCGTACTGGTCACCGAAGGtt                                                                                                                     >  1:451856/1‑149 (MQ=255)                                tttATCAGTCCACCACCAACGCCCACATTAACACCGGCGACGGTGTCGGCATGGCTATCCGTGCCGGCGTACCGGTGCAGGATATGGAAATGTGGCAGTTCCACCCGACCGGCATTGCCGGTGAGGGCGTACTGGTCACCGAAGGTTGc                                                                                                                   >  2:155646/1‑149 (MQ=255)                                tttATCAGTCCACCACCAACGCCCACATTAACACCGGCGACGGTGTCGGCATGGCTATCCGTGCCGGCGTACCGGTGCAGGATATGGAAATGTGGCAGTTCCACCCGACCGGCAATGCCGGtcgg                                                                                                                                           >  4:102780/1‑122 (MQ=255)                                tttATCAGTCCACCACCAACGCCCACATTAACACCGGCGACGGTGTCGGCATGGCTATCCGTGCCGGCGTACCGGTGCAGGATATGGAAATGTGGCAGTTCCACCCGACCGGCAATGCCGGtcgg                                                                                                                                           <  3:102780/125‑4 (MQ=255)                                   aTCAGTCCACCACCAACGCCCACATTAACACCGGCGACGGTGTCGGCATGGCTATCCGTGCCGGCGTACCGGTGCAGGATATGGAAATGTGGCAGTTCCACCCGACCGGCATTGCCGGTGCGGGCGTACTGGTCACCGAAGGTTGCCgt                                                                                                                >  2:109652/1‑149 (MQ=255)                                             caccaACGCCCACATTAACACCGGCGACGGTGTCGGCATGGCTATCCGTGCCGGCGTACCGGTGCAGGATATGGAAATGTGGCAGTTCCACCCGACCGGCATTGCCGGTGCGGGCGTACTGGTCACCGAAGGTTGCCGTGGTGAAGGCg                                                                                                      >  4:104251/1‑149 (MQ=255)                                                caACGCCCACATTAACACCGGCGACGGTGTCGGCATGGCTATCCGTGCCGGCGTACCGGTGCAGGATATGGAAATGTGGCAGTTCCACCCGACCGGCATTGCCGGTGCGGGCGTACTGGTCACCGAAGGTTGCCGTGGTGAAGGCGGtt                                                                                                   <  1:246742/149‑1 (MQ=255)                                                caACGCCCACATTAACACCGGCGACGGTGTCGGCATGGCTATCCGTGCCGGCGTACCGGTGCAGGATATGGAAATGTGGCAGTTCCACCCGACCGGCAATGCCGGtcgg                                                                                                                                           >  2:473421/1‑106 (MQ=255)                                                caACGCCCACATTAACACCGGCGACGGTGTCGGCATGGCTATCCGTGCCGGCGTACCGGTGCAGGATATGGAAATGTGGCAGTTCCACCCGACCGGCAATGCCGGtcgg                                                                                                                                           <  1:473421/109‑4 (MQ=255)                                                      ccACATTAACACCGGCGACGGTGTCGGCATGGCTATCCGTGCCGGCGTACCGGTGCAGGATATGGAAATGTGGCAGTTCCACCCGACCGGCATTGCCGGTGCGGGCGTACTGGTCACCGAAGGTTGCCGTGGTGAAGGCGGTTATctgc                                                                                             <  4:265517/149‑1 (MQ=255)                                                        atgccgacacCGTCGCCGGTGTCGGCATGGCTATCCGTGCCGGCGTACCGGTGCAGGATATGGAAATGTGGCAGTTCCACCCGACCGGCATTGCCGGTGCGGGCGTACTGGTCACCGAAGGTTGCCGTGGTGAAGGCg                                                                                                      <  2:102797/132‑1 (MQ=255)                                                        atgccgacacCGTCGCCGGTGTCGGCATGGCTATCCGTGCCGGCGTACCGGTGCAGGATATGGAAATGTGGCAGTTCCACCCGACCGGCATTGCCGGTGCGGGCGTACTGGTCACCGAAGGTTGCCGTGGTGAAGGCg                                                                                                      >  1:102797/7‑138 (MQ=255)                                                            tAACACCGGCGACGGTGTCGGCATGGCTATCCGTGCCGGCGTACCGGTGCAGGATATGGAAATGTGGCAGTTCCACCCGACCGGCATTGCCGGTGCGGGCGTACTGGTCACCGTAGGTTGCCGTGGTGAAGGCGGTTATCTGCTGaaca                                                                                       >  3:325486/1‑149 (MQ=255)                                                              acacCGGCGACGGTGTCGGCATGGCTATCCGTGCCGGCGTACCGGTGCAGGATATGGAAATGTGGCAGTTCCACCCGACCGGCATTGCCGGTGCGGGCGTACTGGTCACCGAAGGTTGCCGTGGTGAAGGCGGTTATCTGCTGaacaaa                                                                                     <  1:109652/149‑1 (MQ=255)                                                                        cGGTGTCGGCATGGCTATCCGTGCCGGCGTACCGGTGCAGGATATGGAAATGTAGCAGTTCCACCCGACCGGCATTGCCGGTGCGGGCGTACTGGTCACCGAAGGTTGCCGTGGTGAAGGCGGTTATCTGCTGAACAAACATGGCGAAc                                                                           <  1:258789/149‑1 (MQ=255)                                                                               ggCATGGCTATCCGTGCCGGCGTACCGGTGCAGGATATGGAAATGTGGCAGTTCCACCCGACCGGCATTGCCGGTGCGGGCGTACTGGTCACCGAAGGTTGCCGTGGTGAAGGCGGTTATCTGCTGAACAAACATGGCGAACGTTTTAt                                                                    <  2:476993/149‑1 (MQ=255)                                                                                gCATGGCTATCCGTGCCGGCGTACCGGTGCAGGATATGGAAATGTGGCAGTTCCACCCGACCGGCATTGCCGGTGCGGGCGTACTGGTCACCGAAGGTTGCCGTGGTGAAGGCGGTTATCTGCTGAACAAACATGGCGAACGTTTTATg                                                                   >  2:5489/1‑149 (MQ=255)                                                                                gCATGGCTATCCGTGCCGGCGTACCGGTGCAGGATATGGAAATGTGGCAGTTCCACCCGACCGGCATTGCCGGTGCGGGCGTACTGGTCACCGAAGGTTGCCGTGGTGAAGGCGGTTATCTGCTGAACAAACATGGCGAACGTTTTATg                                                                   <  3:54041/149‑1 (MQ=255)                                                                                gCATGGCTATCCGTGCCGGCGTACCGGTGCAGGATATGGAAATGTGGCAGTTCCACCCGACCGGCATTGCCGGTGCGGGCGTACTGGTCACCGAAGGTTGCCGTGGTGAAGGCGGTTATCTGCTGAACAAACATGGCGAACGTTTTATg                                                                   >  3:65153/1‑149 (MQ=255)                                                                                  aTGGCTATCCGTGCCGGCGTACCGGTGCAGGATATGGAAATGTGGCAGTTCCACCCGACCGGCATTGCCGGTGCGGGCGTACTGGTCACCGAAGGTTGCCGTGGTGAAGGCGGTTATCTGCTGAACAAACATGGCGAACGTTTTATGGa                                                                 <  1:155646/149‑1 (MQ=255)                                                                                           cGTGCCGGCGTACCGGTGCAGGATATGGAAATGTGGCAGTTCCACCCGACCGGCATTGCCGGTGCGGGCGTACTGGTCACCGAAGGTTGCCGTGGTGAAGGCGGTTATCTGCTGAACAAACATGGCGAACGTTTTATGGAGCGTTATgc                                                        <  2:451856/149‑1 (MQ=255)                                                                                                cGGCGTACCGGTGCAGGATATGGAAATGTGGCAGTTCCACCCGACCGGCATTGCCGGTGCGGGCGTACTGGTCACCGAAGGTTGCCGTGGTGAAGGCGGTTATCTGCTGAACAAACATGGCGAACGTTTTATGGAGCGTTATGCGCCGa                                                   >  1:223770/1‑149 (MQ=255)                                                                                                                gATATGGAAATGTGGCAGTTCCACCCGACCGGCATTGCCGGTGCGGGCGTACTGGTCACCGAAGGTTGCCGTGGTGAAGGCGGTTATCTGCTGAACAAACATGGCGAACGTTTTATGGAGCGTTATGCGCCGAACGCCAAAGACCTGGc                                   <  3:321753/149‑1 (MQ=255)                                                                                                                     ggAAATGTGGCAGTTCCACCCGACCGGCATTGCCGGTGCGGGCGTACTGGTCACCGAAGGTTGCCGTGGTGAAGGCGGTTATCTGCTGAACAAACATGGCGAACGTTTTATGGAGc                                                               >  2:213967/1‑116 (MQ=255)                                                                                                                     ggAAATGTGGCAGTTCCACCCGACCGGCATTGCCGGTGCGGGCGTACTGGTCACCGAAGGTTGCCGTGGTGAAGGCGGTTATCTGCTGAACAAACATGGCGAACGTTTTATGGAGc                                                               <  1:213967/116‑1 (MQ=255)                                                                                                                      gAAATGTGGCAGTTCCACCCGACCGGCATTGCCGGTGCGGGCGTACTGGTCACCGAAGGTTGCCGTGGTGAAGGCGGTTATCTGCTGAACAAACATGGCGAACGTTTTATGGAGCGTTATGCGCCGAACGCCAAAGACCTGGCGGGCCg                             <  3:104251/149‑1 (MQ=255)                                                                                                                              gCAGTTCCACCCGACCGGCATTGCCGGTGCGGGCGTACTGGTCACCGAAGGTTGCCGTGGTGAAGGCGGTTATCTGCTGAACAAACATGGCGAACGTTTTATGGAGCGTTATGCGCCGAACGCCAAAGACCTGGCGGGCCGTGACGTgg                     >  3:71419/1‑149 (MQ=255)                                                                                                                              gCAGTTCCACCCGACCGGCATTGCCGGTGCGGGCGTACTGGTCACCGAAGGTTGCCGTGGTGAAGGCGGTTATCTGCTGAACAAACATGGCGAACGTTTTATGGAGCGTTATGCGCCGAACGCCAAAGACCTGGCGGGCCGTGACGTgg                     >  4:306369/1‑149 (MQ=255)                                                                                                                              gCAGTTCCACCCGACCGGCATTGCCGGTGCGGGCGTACTGGTCACCGAAGGTTGCCGTGGTGAAGGCGGTTATCTGCTGAACAAACATGGCGAACGTTTTATGGAGCGTTATGCGCCGAACGCCAAAGACCTGGCGGGCCGTGACGTgg                     <  1:5489/149‑1 (MQ=255)                                                                                                                                aGTTCCACCCGACCGGCATTGCCGGTGCGGGCGTACTGGTCACCGAAGGTTGCCGTGGTGAAGGCGGTTATCTGCTGAACAAACATGGCGAACGTTTTATGGAGCGTTATGCGCCGAACGCCaaa                                           >  1:202803/1‑125 (MQ=255)                                                                                                                                aGTTCCACCCGACCGGCATTGCCGGTGCGGGCGTACTGGTCACCGAAGGTTGCCGTGGTGAAGGCGGTTATCTGCTGAACAAACATGGCGAACGTTTTATGGAGCGTTATGCGCCGAACGCCaaa                                           <  2:202803/125‑1 (MQ=255)                                                                                                                                      aCCCGACCGGCATTGCCGGTGCGGGCGTACTGGTCACCGAAGGTTGCCGTGGTGAAGGCGGTTATCTGCTGAACAAACATGGCGAACGTTTTATGGAGCGTTATGCGCCGAACGCCAAAGACCTGGCGGGCCGTGACGTGGTTGCGCGt             >  2:128127/1‑149 (MQ=255)                                                                                                                                       ccgcaCCGGCAATGCCGGT‑CGGGCGTACTGGTCACCGAAGGTTGCCGTGGTGAAGGCGGTTATCTGCTGAACAAACATGGCGAACGTTTTATGGAGCGTTATGCGCCGAACGCCAAAGACCTGGCGGGCCGTGACGTGGTTGCGCGTTc           >  3:200006/5‑149 (MQ=255)                                                                                                                                                 aTTGCCGGTGCGGGCGTACTGGTCACCGAAGGTTGCCGTGGTGAAGGCGGTTATCTGCTGAACAAACATGGCGAACGTTTTATGGAGCGTTATGCGCCGAACGCCAAAGACCTGGCGGGCCGTGACGTGGTTGCGCGTTCCATCATGAt  >  2:462990/1‑149 (MQ=255)                                                                                                                                                  ttGCCGGTGCGGGCGTACTGGTCACCGAAGGTTGCCGTGGTGAAGGCGGTTATCTGCTGAACAAACATGGCGAACGttt                                                                       >  4:157334/1‑79 (MQ=255)                                                                                                                                                  ttGCCGGTGCGGGCGTACTGGTCACCGAAGGTTGCCGTGGTGAAGGCGGTTATCTGCTGAACAAACATGGCGAACGttt                                                                       <  3:157334/79‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                    ACCGTGCTGGCGACTGGCGGAGCAGGGCGTATTTATCAGTCCACCACCAACGCCCACATTAACACCGGCGACGGTGTCGGCATGGCTATCCGTGCCGGCGTACCGGTGCAGGATATGGAAATGTGGCAGTTCCACCCGACCGGCATTGCCGGTGCGGGCGTACTGGTCACCGAAGGTTGCCGTGGTGAAGGCGGTTATCTGCTGAACAAACATGGCGAACGTTTTATGGAGCGTTATGCGCCGAACGCCAAAGACCTGGCGGGCCGTGACGTGGTTGCGCGTTCCATCATGAT  >  NC_000913/756498‑756790

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