Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	1,305,683	Δ1 bp	coding (557/1461 nt)	cls ←	cardiolipin synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	1,305,683	0	G	.	96.2%	112.1 / ‑1.5	26	coding (557/1461 nt)	cls	cardiolipin synthetase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base . (13/12);  minor base A (1/0);  total (14/12)
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.73e-01

TGACCTTTAAGGCTTCGACCACTTCAATACCGGCATTACGCATTAGCTCGGGCCACGGGCTGCGGAAAAAAGCCACACTCCCGGCGGAGTCGAGCATCAATCGGCAATGAATGCCGCGTCGCGCAGCGGCCATTAATGATTCAG  >  REL606/1305612‑1305755                                                                        |                                                                         tGACCTTTAAGGCTTCGACCACTTCAATACCGGCATTACGCATTAGCTCGGGCCACGGGCTGCGGAAAAAA‑Ccaca                                                                     <  1:451706/76‑1 (MQ=255)      tttAAGGCTTCGACCACTTCAATACCGGCATTACGCATTAGCTCGGGCCACGGGCTGCGGAAAAAA‑CCACACTccc                                                                >  1:1262656/1‑76 (MQ=255)        tAAGGCTTCGACCACTTCAATACCGGCATTACGCATTAGCTCGGGCTACGGGCTGCGGAAAAAA‑CCACACTCCcgg                                                              >  2:1446670/1‑76 (MQ=255)          aGGCTTCGACCACTTCAATACCGGCATTACGCATTAGCTCGGGCCACGGGCTGCGGAAAAAA‑CCACACTCCcggcg                                                            >  2:205889/1‑76 (MQ=255)             cTTCGACCACTTCAATACCGGCATTACGCATTAGCTCGGGCCACGGGCTGCGGAAAAGA‑CCACACGCCCGGCGGAg                                                         <  1:1618694/76‑1 (MQ=255)               tCGACCACTTCNATNCCGGCATTACGTATTAGCTCGGGCCACGGGCTGCGGAAAAAA‑CGACACTCCAGGCGGAGTc                                                       <  2:257514/76‑1 (MQ=255)                cGACCACTTCAATACCGGCATTACGCATTAGCTCGGGCCACGGGCTGCGGAAAAAA‑CCACACTCCCGGCGGAGTCg                                                      <  2:817688/76‑1 (MQ=255)                 gACCACTTCAATACCGGCATTACGCATTAGCTCGGGCCACGGGCTGCGGAAAAAA‑CCACACTCCCGGCGGAGTCGa                                                     <  2:939012/76‑1 (MQ=255)                  aCCACTTCAATACCGGCATTACGCATTAGCTCGGGCCACGGGCTGCGGAAAAAA‑CCACACTCCCGGCGGAGTCGAg                                                    <  1:1210266/76‑1 (MQ=255)                   ccACTTCAATACTGGCATTACGCGTTAGCTCGGGCCACGGGCTGTGGAAAAAA‑CCACANTCCCGGCGGAGTCGAGc                                                   <  1:211974/76‑1 (MQ=255)                   ccACTTCAATACCGGCATTACGCATTAGCTCGGGCCACGGGCTGCGGAAAAAA‑CCACACTCCCGGCGGAGTCGAGc                                                   <  1:1066249/76‑1 (MQ=255)                         cAATACCGGCATTACGCATTAGCTCGGGCCACGGGCTGCGGAAAAAA‑CCACACTCCCGGCGGAGTCGAGCATCAAt                                             >  1:769208/1‑76 (MQ=255)                                     tACGCATTAGCTCGGGCCACGGGCTGCGGAAAAAA‑CCACACTCCCGGCGGAGTCGAGCATCAATCGGCaatgaatg                                 <  2:177340/76‑1 (MQ=255)                                       cGCATTAGCTCGGGCCACGGGCTGCGGAAAAAA‑CCACACTCCCGGCGGAGTCGAGCATCAATCGGCAATGAATGcc                               >  1:1550721/1‑76 (MQ=255)                                        gCATTAGCTCGGGCCACGGGCTGCGGAAAAAA‑CCACACTCCCGGCGGANACGAGCATCAATCGGCAATGAATGCcg                              <  2:4701/76‑1 (MQ=255)                                        gCATTAGCTCGGGCCACGGGCTGCGGAAAAAA‑CCACACTCCCGGCGGAGTCGAGCATCAATCGGCAATGAATGCcg                              >  1:1557734/1‑76 (MQ=255)                                             aGCTCGGGCCACGGGCTGCGGAAAAAA‑CCACACTCCCGGCGGAGTCGAGCATCAATCGGCAATGAATGCCGCGTcg                         <  2:204296/76‑1 (MQ=255)                                               cTCGGGCCACGGGCTGCGGAAAAAA‑CCACACTCCCGGCGGAGTCGAGCATCAATCGGCAATGAATGCCGCGTcgcg                       <  2:362674/76‑1 (MQ=255)                                                 cGGGCCACGGGCTGCGGAAAAAA‑CCACACTCCCGGCGGAGTCGAGCATCAATCGGCAATGAATGCCGCGTCGCGCa                     >  2:658753/1‑76 (MQ=255)                                                  gggCCACGAGCTGCGGAAAAAA‑CCACGGACCTGGCAGCGTCGCGCAGCAAGCGGCATTGAGTGTCGCGTCGCGCAg                    <  1:1657788/76‑1 (MQ=255)                                                     ccACGGGCTGCGGAAAAAA‑CCACACTCCCGGCGGAGTCGAGCATCAATCGGCAATGAATGCCGCGTCGCGCAGCgg                 >  1:75433/1‑76 (MQ=255)                                                        cGGGCTGCGGAAAAAA‑CCACACTCCCGGCGGAGTCGAGCATCAATCGGGAATGAATGCCGCGTCGCGCAGCGGCCa              >  1:164185/1‑76 (MQ=255)                                                        cGGGCTGCGCAAAAAA‑CGGCACTCCCGGCGGAGTCGAGCATCAATCGGCAATGAATGCCGCNTCNCGCAGCGGCCa              >  2:256882/1‑76 (MQ=255)                                                          ggCTGCGGAAAAAA‑CCACACTCCCGGCGGAGTCGAGCATCAATCGGCAATGAATGCCGCGTCGCGCAGCGGCCCtt            >  1:1295785/1‑76 (MQ=255)                                                          ggCTGCGGAAAAAA‑CCACACTCCCGGCGGAGTCGAGCATCAATCGGCAATGAATGCCGCGTCGCGCAGCGGCCAtt            <  1:1217110/76‑1 (MQ=255)                                                              gCGGAAAAAA‑CCACACTCCCGGCGGAGTCGAGCATCAATCGGCAATGAATGCCGCGTCGCGCAGCGGCCATTAATg        >  2:188727/1‑76 (MQ=255)                                                              gCGGAAAAAA‑CCACACTCCCGGCGGAGTCGAGCATCAATCGGCAATGAATGCCGCGTCGCGCAGCGGCCATTAATg        >  2:805376/1‑76 (MQ=255)                                                                 gganaaaACCACACTCCCGGCGGAGTCGAGCATCAATCGGCAATGAATGCCGCGTCGCGCAGCGGCCATTAATGAt      >  1:1543063/3‑76 (MQ=255)                                                                     aaaaCCACACTCCCGGCGGAGTCGAGCATCAATCGGCAATGAATGCCGCGTCGCGCAGCGGCCATTAATGATTCAg  <  1:910905/76‑1 (MQ=255)                                                                     aaaaCCACACTCCCGGCGGAGTCGAGCATCAATCGGCAATGAATGCCGCGTCGAGCAGCGGCCATTAATGATTCAg  >  2:141204/1‑76 (MQ=255)                                                                        |                                                                         TGACCTTTAAGGCTTCGACCACTTCAATACCGGCATTACGCATTAGCTCGGGCCACGGGCTGCGGAAAAAAGCCACACTCCCGGCGGAGTCGAGCATCAATCGGCAATGAATGCCGCGTCGCGCAGCGGCCATTAATGATTCAG  >  REL606/1305612‑1305755

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