Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	753,013	G→T	45.4%	E551* (GAA→TAA)	sdhA →	succinate dehydrogenase, flavoprotein subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	753,013	0	G	T	45.4%	3.8 / 52.7	42	E551* (GAA→TAA)	sdhA	succinate dehydrogenase, flavoprotein subunit
Reads supporting (aligned to +/- strand):  ref base G (7/16);  new base T (8/11);  total (15/27)
Fisher's exact test for biased strand distribution p-value = 5.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.13e-01

ACTTCCAGCGAGTTCAACACCCAGCGCGTTGAGTGCCTGGAACTGGATAACCTGATGGAAACGGCGTATGCAACGGCTGTTTCTGCCAACTTCCGTACCGAAAGCCGTGGCGCGCATAGCCGCTTCGACTTCCCGGATCGTGATGATGAAAACTGGCTGTGCCACTCCCTGTATCTGCCAGAGTCGGAATCCATGACGCGCCGAAGCGTCAACATGGAACCGAAACTGCGCCCGGCATTCCCGCCGAAGATTCGTACTTACTAATGCGGAGACAGGAAAATGAGACTCGAG  >  CP009273/752866‑753156                                                                                                                                                    |                                                                                                                                                acTTCCAGCGAGTTCAACACCCAGCGCGTTGAGTGCCTGGAACTGGATAACCTGATGGAAACGGCGTATGCAACGGCTGTTTCTGCCAACTTCCGTACCGAAAGCCGTGGCGCGCATAGCCGCGTCGACTTCCCGGATCGtgatgctga                                                                                                                                                >  2:754938/1‑149 (MQ=255)  cTTCCAGCGAGTTCAACACCCAGCGCGTTGAGTGCCGGGAACTGGATAACCTGATGGAAACGGCGTATGCAACGGCTGTTTCTGCCAACTTCCGTACCGAAAGCCGTGGCGCGCATAGCCGCTTCGACTTCCCGGATCGTGATGATGaa                                                                                                                                               <  1:24864/149‑1 (MQ=255)             ttCAACACCCAGCGCGTTGAGTGCCTGGAACTGGATAACCTGATGGAAACGGCGTATGCAACGGCTGTTTCTGCCAACTTCCGTACCGAAAGCCGTGGCGCGCATAGCCGCTTCGACTTCCCGGCTCGTGATGATGAAAACTGGCtgtg                                                                                                                                    <  2:647999/149‑1 (MQ=255)              tccacacCCAGCGCGTTGGGTGCCGGGAGCTGGATAACCGGATGGAAACGGCGTATGCAACGGCTGTTTCTGCCAACTTCCGTACCGAAAGCCGTGGCGCGCATAGCCGCTTCGACTTCCCGGATCGTGATGATGAAAACTGGCTGTGc                                                                                                                                   <  1:457365/146‑1 (MQ=255)                aaCACCCAGCGCGTTGAGTGCCTGGGACTGCATAACCTGATGGAAACGGCGTATGCAACGGCTGTTTCTGCCAACTTCCGTACCGAAAGCCGTGGCGCGCATAGCCGCTTCGACTTCCCGGATCGTGATGATGAAAACTGGCTGTGCCa                                                                                                                                 >  2:449438/1‑149 (MQ=255)                            gTTGAGTGCCTGGAACTGGATAACCTGATGGAAACGGCGTATGCAACGGCTGTTTCTGCCAACGTCCGGACCGAAAGCCGTGGCGCGCCTAGCCGCGTCGACTTCCGGGGACGTGATGATGAAAACGGGGTGAGCGACTGACTGTAGCt                                                                                                                     >  1:741020/1‑149 (MQ=255)                              tGAGGGCCTGGAATCGGAGAACCTGATGGAAACGGCGTATGCAACGGCTGTTTCTGCCAAGTTCCGTACCGAAAGCCGTGGCGCGCATAGCCGCTTCGACTTCCCGGATCGTGATGATGAAAACTGGCTGTGCCACGCCCTGTATCTGc                                                                                                                   <  1:71592/149‑1 (MQ=255)                                agctccTGCAACTGGATAACCCAATGGAAACGGCGTATCCAACGGCTGTTTCTGCCAACTTCCGTACCGAAAGCCGTGGCGCGCATAGCCGCTTCGACTCCCCGGATCGTGATGATTAAAACTGGCTGTGCCACTCCCTGTATCTGCCa                                                                                                                 <  1:439564/145‑1 (MQ=255)                                agcgcccgtaaCTCGTTAACCTTATGGAAACGGCGTATGCAACGGCTGTTTCTGCACACGTTATGACCGAAAGCCGTGGCGCGCATAGCTGCTTGGACTTCCCGGATCGTGATGATGAAAACTGGCTGTGCCACTCCCTGTATCTGCCa                                                                                                                 <  2:418077/140‑1 (MQ=255)                                agcgCCGGGAACTGGAGAACCTGATGGAAACGGCGTATGCAACGGCTGTTTCTGCCAACTTCCGTACCGAAAGCCGTGGCGCGCATAGCCGCTTCGACTTCCCGGATCGTGATGATGAAAACTGGCTGTGCCACTCCCTGTATCTGCCa                                                                                                                 <  2:235475/146‑1 (MQ=255)                                aGTGCCTGGCACTGGATAACCTGATGGAAACGGCGTATGCAACGGCTGTTTCTGCCAACTTCCGTACCGAAAGCCGTGGCGCGCATAGCCGCTTCGACTTCCCGGATCGTGATGATTAAAACTGGCTGTGCCACTCCCTGTATCTGCCa                                                                                                                 <  1:833224/149‑1 (MQ=255)                                aGTCCCGGGAACTGGATAACCTGATGGAAACGGCGTATGCAACGGCTGTTTCTGCCAACTTCCGTACCGAAAGCCGTGGCGCGCATAGCCGCTTCGACTTCCCGGATCGTGATGATGAAAACTGGCTGTGCCACTCCCTGTATCTGCCa                                                                                                                 <  2:131276/149‑1 (MQ=255)                                     cTGGAACTGGATAACCTGATGGAAACGGCGTATGCAACGGCTGTTTCTGCCAACTTCCGTACCGAAAGCCGTGGCGCGCATAGCCGCTTCGACTTCCCGGATCGTGATGATTAAAACTGGCTGTGCCACTCCCTGTATCTGCCAGAGTc                                                                                                            <  1:658033/149‑1 (MQ=255)                                     cTGGAACTGGATAACCTGATGGAAACGGCGTATGCAACGGCTGTTTCTGCCAACTTCCGTACCGAAAGCCGTGGCGCGCATAGCCGCTTCGACTTCCCGGATCGTGATGATTAAAACTGGCTGTGCCACTCCCTGTATCTGCCAGAGTc                                                                                                            <  1:195555/149‑1 (MQ=255)                                         aaCTGGATAACCTGATGGAAACGGCGTATGCAACGGCTGTTTCTGCCAACTTCCGTACCGAAAGCCGTGGCGCGCATAGCCGCTTCGACTTCCCGGATCGTGATGATTAAAACTGGCTGTGCCACTCCCTGTATCTGCCAGAGTCGGaa                                                                                                        >  2:665270/1‑149 (MQ=255)                                         aaCTGGATAACCTGATGGAAACGGCGTATGCAACGGCTGTTTCTGCCAACTTCCGTACCGAAAGCCGTGGCGCGCATAGCCGCTTCGACTTCCCGGATCGTGATGATTAAAACTGGCTGTGCCACTCCCTGTATCTGCCAGAGTCGGaa                                                                                                        >  2:294388/1‑149 (MQ=255)                                           cTGGATAACCTGATGGAAACGGCGTATGCAACGGCTGTTTCTGCCAACTTCCGTACCGAAAGCCGTGGCGCGCATAGCCGCTTCGACTTCCCGGATCGTGATGATTAAAACTGGCTGTGCCACTCCCTGTATCTGCCAGAGTCGGAATc                                                                                                      <  2:558541/149‑1 (MQ=255)                                                tAACCTGATGGAAACGGCGTATGCAACCGCTGTTTCTGCCAACTTCCGTACCGAACGCCGTGGCGCGCATAGCCGCTTCGACTTCCCGGATCGTGATGATGAAAACTGGCTGTGCCACTCCCTGTATCTGCCAGAGTCGGAATCCATGa                                                                                                 <  2:824842/149‑1 (MQ=255)                                                    cTGATGGAAACGGCGTATGCAACGGCAGTTTCTGCCAACTTCCGTACCGAAAGCCGTGGCGCGCATAGCCGCTTCGACTTCCCGGATCGTGATGATGAAAACTGGCTGTGCCACTCCCTGTATCTGCCAGAGTCGGAATCCATGAcgcg                                                                                             >  2:465066/1‑149 (MQ=255)                                                     tGATGGAAACGGCGTATGCAACGGCTGTTTCTGCCAACTTCCGTACCGAAAGCCGTGGCGCGCATAGCCGCTTCGACTTCCCGGATCGTGATGATTAAAACTGGCTGTGCCACTCCCTGTATCTGCCAGAGTCGGAATCCATGACgcgc                                                                                            <  1:209236/149‑1 (MQ=255)                                                      gATGGAAACGGCGTATGCAACGGCTGTTTCTGCCAACTTCCGTACCGAAAGCCGTGGCGCGCATAGCCGCTTCGACTTCCCGGATCGTGATGATGAAAACTGGCTGTGCCACTCCCTGTATCTGCCAGAGTCGGAATCCATGACGCGcc                                                                                           >  2:44255/1‑149 (MQ=255)                                                       aTGGAAACGGCGTATGCAACGGCTGTTTCTGCCAACTTCCGTACCGAAAGCCGTGGCGCGCATAGCCGCTTCGACTTCCCGGATCGTGATGATTAAAACTGGCTGTGCCACTCCCTGTATCTGCCCGAGCCGGAATCCATGACGCGCCg                                                                                          >  2:446297/1‑149 (MQ=255)                                                         ggAAACGGCGTATGCAACGGCTGTTTCTGCCAACTTCCGTACCGAAAGCCGTGGCGCGCATAGCCGCTTCGACTTCCCGGATCGTGATGATTAAAACTGGCTGTGCCACTCCCTGTATCTGCCAGAGTCGGAATCCATGACGCGCCGaa                                                                                        >  2:344588/1‑149 (MQ=255)                                                         ggAAACGGCGTATGCAACGGCTGTTTCTGCCAACTTCCGTACCGAAAGCCGTGGCGCGCATAGCCGCTTCGACTTCCCGGATCGTGATGATGAAAACTGGCTGTGCCACTCCCTGTATCTGCCAGAGTCGGAATCCATGACGCGCCGaa                                                                                        <  2:472472/149‑1 (MQ=255)                                                           aaaCGGCGTATGCAACGGCTGTTTCTGCCAACTTCCGTACCGAAAGCCGTGGCGCGCATAGCCGCTTAGACTTCCCGGATCGTGATGATGAAAACTGGCTGTGCCACTCCCTGTATCTGCCAGAGTCGGAATCCATGACGCGCCGAAGc                                                                                      >  2:525938/1‑149 (MQ=255)                                                             aCGGCGTATGCAACGGCTGTTTATGCCAACTTCCGTACCGAAAGCCGTGGCGCGCATCGCCGCTTCGACTTCCCGGATCGGGATGATGAAAACTGGCTGTGCCACTCCCTGTATCTGCCAGAGTCGGAATCCCTGACGCGCCGAAGCGt                                                                                    <  2:549564/149‑1 (MQ=255)                                                             aCGGCGTATGCAACGGCTGTTGCTCCCAACCTCCGTACCGAAAGCCGTGGCGCGCATAGCCGCTTCGACTTCCCGGATCGTGATGATTAAAACTGGCTGTGCCACTCCCTGTATCTGCCAGAGTCGGAATCCATGACGCGCCGAAGCGt                                                                                    <  2:311505/149‑1 (MQ=255)                                                                         aCGGCTGTTTCTGCCAACTTCCGTACCGAAAGCCGTGGCGCGCATAGCCGCTTCGACGTCCCGGAGCATGATTATTAAAACTGGCTGTGCCACTCCCTGTATCTGCCAGAGTCGGAATCCGTGACGCGCCGAAGCGTCAACATGGAAcc                                                                        <  2:650114/149‑1 (MQ=255)                                                                                tttCTGCCAACTTCCGTACCGAAAGCCGTGGCGCGCATAGCCGCTTCGACTTCCCGGATCGTGATGATTAAAACTGGCTGTGCCACTCCCTGTATCTGCCAGAGTCGGAATCCATGACGCGCCGAAGCGTCAACATGGAACCGAAACTg                                                                 >  2:293869/1‑149 (MQ=255)                                                                                      ccAACTTCCGTACCGAAAGCCGTGGCGCGCATAGCCGCTTCGACTTCCCGGATCGTGATGATTAAAACTGGCTGTGCCACTCCCTGTATCTGCCAGAGTCGGCATACATTACGCGCCGAAGGGTCAACATGGAACCGAAACTGCGCCCg                                                           >  2:162502/1‑149 (MQ=255)                                                                                         aCTTCCGTACCGAAAGCCGTGGCGCGCATAGCCGCTTCGACTTCCCGGATCGTGATGATGAAAACTGGCTGTGCCACTCCCTGTATCTGCCAGAGTCGGAATCCATGACGCGCCGAAGCGTCAACATGGAACCGAAACTGCGCCCGGCa                                                        >  2:705252/1‑149 (MQ=255)                                                                                                tACCGAAAGCCGTGGCGCGCATAGCCGCTTCGACTTCCCGGATCGTGATGATTAAAACTGGCTGTGCCACTCCCTGTATCTGCCAGAGTCGGAATCCATGACGCGCCGAAGCGTCAACATGGAACCGAAACTGCGCCCGGCATTCccgc                                                 <  1:294388/149‑1 (MQ=255)                                                                                                    gAAAGCCGTGGCGCGCATAGCCGCTTCGACTTCCCGGATCGTGATGATTAAAACTGGCTGTGCCACTCCCTGTATCTGCCAGAGTCGGAATCCATGACGCGCCGAAGCGTCAACATGGAACCGAAACTGCGCCCGGCATTCCCGCCGaa                                             <  2:347328/149‑1 (MQ=255)                                                                                                                    aTAGCCGCTTCGACTTCCCGGATCGTGATGATTAAAACTGGCTGTGCCACTCCCTGTATCTGCCAGAGTCGGAATCCATGACGCGCCGAAGCGTCAACATGGAACCGAAACTGCGCCCGGCATTCCCGCCGAAGATTCGTACTTACTaa                             >  2:860968/1‑149 (MQ=255)                                                                                                                        ccGCTTCGACTTCCCGGATCGTGATGATGAAAACTGGCTGTGCCACTCCCTGTATCTGCCAGAGTCGGAATCCATGACGCGCCGAAGCGTCAACATGGAACCGAAAATGCGTACCGCATGCCCGCCGAAGATTCGTACTTACTAATGCg                         <  2:425434/149‑1 (MQ=255)                                                                                                                              cGACTTCCCGGATCGTGATGATGAAAACTGGCTGTGCCACTCCCTGTATCTGCCAGAGTCGGAATGCATGACGCGCCGAAGCGTCAACCTGGCACCGAAACTGCGCCCGGCATTCCCCCCGAAGATTCGTACTTACTAAAGCGGCGAcc                   >  1:162622/1‑148 (MQ=255)                                                                                                                               gACTTCCCGGATCGTGATGATTAAAACTGGCTGTGCCACTCCCTGTATCTGCCAGAGTCGGAATCCATGACGCGCCGAAGCGTCAACATGGAACCGAAACTGCGCCCGGCATTCCCGCCGAAGATTCGTACTTACTAATGCGGAGACAg                  >  1:264658/1‑149 (MQ=255)                                                                                                                                aCTTCCCGGATCGTGATGATTAAAACTGGCTGTGCCACGCCCTGTATCTCCCAGAGTCGGAATCCATGACGCGCCGAAGCGTCAACATGGAACCGAAACTGCGCCCGGCATTCCCGCCGAAGATTCGTACTTACTAATGCGGAGACAgg                 <  2:547276/149‑1 (MQ=255)                                                                                                                                aCTTCCCGGATCGTGATGATGAAAACTGGCTGTGCCACTCCCTGTATCTGCCAGAGTCGGAATCCATGACGCGCCGAAGCGTCAACATGGAACCGAAACTGCGCCCGGCATTCCCGCCGAAGATTCGTACTTACTAATGCGGAGACAgg                 <  2:607790/149‑1 (MQ=255)                                                                                                                                aCTTCCCGGATCGTGATGATGAAAACGGGCTGTGCCACTCCCTGTATCTGCCAGAGTCGGAATCCATGACGCGCCGAAGCGTCAACATGGAACCGAAACTGCGCCCGGCATTCCCGCCGAAGATTCGTACTTACTAATGCGGAGACAgg                 <  2:291825/149‑1 (MQ=255)                                                                                                                                aCTCCCCGGACCGCGATGATGAAAACTGGCTGTGCCACTCCATGTATCTGCCAGAGTCGGAATCCATGACGCGCCGAAGCGTCAACATGGAACCGAAACTGCGCCCGGCATTCCCGCCGAAGATTCGTACTTACTAATGCGGAGACAgg                 <  1:259323/149‑1 (MQ=255)                                                                                                                                      cGGATCGTGATGATGAAAACTGGCTGTGCCACTCCCTGTATCTGCCAGAGTCGGAATCCATGACGCGCCGAAGCGTCAACATGGAACCGAAACTGCGCCCGGCATTCCCGCCGAAGATTCGTACTTACTAATGCGGAGACAGGAAAATg           <  1:465066/149‑1 (MQ=255)                                                                                                                                         aTCGTGATGATGAAAACTGGCTGTGCCACTCCCTGTATCTGCCAGAGTCGGAATCCATGACGCGCCGAAGCGTCAACATGGAACCGAAACTGCGCCCGGCATTCCCGCCGAAGATTCGTACTTACTAATGCGGAGACAGGAAAATgaga        <  1:754938/149‑1 (MQ=255)                                                                                                                                               atgatgaAAACTGGCTGTGCCACTCCCTGTATCTGCCAGAGTCGGAATCCATGACGCGCCGAAGCGTCAACATGGAACCGAAACTGCGCCCGGCATTCCCGCCGAAGATTCGTACTTACTAATGCGGAGACAGGAAAATGAGACTCGAg  >  2:543079/1‑149 (MQ=255)                                                                                                                                               atgatgaAAACTGGCTGTGCCACTCCCTGTATCTGCCAGAGTCGGAATCCATGACGCGCCGAAGCGTCAACATGGAACCGAAACTGCGCCCGGCATTCCCGCCGAAGATTCGTACTTACTAATGCGGAGACAGGAAAATGAGACTCGAg  >  1:554397/1‑149 (MQ=255)                                                                                                                                               atgatgaAAACTGGCTGTGCCACTCCCTGTATCTGCCAGAGTCGGAATCCATGACGCGCCGAAGCGTCAACATGGAACCGAAACTGCGCCCGGCATTCCCGCCGAAGATTCGTACTGACTAAGGCGGCGACAGGCAAATGAGACTCGAg  >  2:788554/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                ACTTCCAGCGAGTTCAACACCCAGCGCGTTGAGTGCCTGGAACTGGATAACCTGATGGAAACGGCGTATGCAACGGCTGTTTCTGCCAACTTCCGTACCGAAAGCCGTGGCGCGCATAGCCGCTTCGACTTCCCGGATCGTGATGATGAAAACTGGCTGTGCCACTCCCTGTATCTGCCAGAGTCGGAATCCATGACGCGCCGAAGCGTCAACATGGAACCGAAACTGCGCCCGGCATTCCCGCCGAAGATTCGTACTTACTAATGCGGAGACAGGAAAATGAGACTCGAG  >  CP009273/752866‑753156

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑