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

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	1,305,211	0	G	.	97.8%	221.5 / ‑1.6	46	coding (1029/1461 nt)	cls	cardiolipin synthetase
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base . (25/20);  total (26/20)
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.66e-01

ATAATACTGACATCCACCCCGCGCTGTGCCGCCGTGCAAATCGCATGAAGTAAATCATCGCTTGGCACAAAGTAGGGCGTTGTCATGATCAAATATTCACGCGCCGAATAAGCCGCAGTCAATAATGCCTGGTGAATGAGA  >  REL606/1305140‑1305280                                                                        |                                                                      ataataCTGACATCCACCCCGCGCTGTGCCGCCGTGCAAATCGCATGAAGTAAATCATCGCTTGGCACAAA‑TAggg                                                                  >  1:391401/1‑76 (MQ=255) ataataCTGACATCCACCCCGCGCTGTGCCGCCGTGCAAATCGCATGAAGTAAATCATCGCTTGGCACAAA‑TAggg                                                                  >  2:525720/1‑76 (MQ=255)   aataCTGANATCCACCCCGCGCTGTGCCGCCGTGCAAATCGCATGAAGTAAATCATCGCTTGGCACAAA‑TAGGGCg                                                                >  1:1754271/1‑76 (MQ=255)         gACATCCACCCCGCGCTGTGCCGCCGTGCAAATCGCATGAAGTAAATCATCGCTTGGCACAAA‑TAGGGCGCTGTCa                                                          >  2:77271/1‑76 (MQ=255)           cATCCACCCCGCGCTGTGCCGCCGTGCAAATCGCATGAAGTAAATCATCGCTTGGCACAAA‑TAGGGCGTTGTCATg                                                        >  1:1464891/1‑76 (MQ=255)            aTCCACCCCGCGCTNTGCCGCCGTTCAAATCGCATGAGGTAAATCATCGCTTGGCACAAA‑TAGGGCGTTGTCATGa                                                       >  1:9044/1‑76 (MQ=255)            aTCCACCCCGCGCTGTGCCGCCGTGCAAATCGCATGAAGTANATCATCGCTTGGCACAAA‑NAGNGCGTTGTCATGa                                                       >  1:1493829/1‑76 (MQ=255)              ccACCCCGCGCTGTGCCGCCGTGCAAATCGCATGAAGTAAATCATCGCTTGGCACAAA‑TAGGGCGTTGTCATGATc                                                     >  2:1722656/1‑76 (MQ=255)                aCCCCGCGCTGTGCCGCCGTGCAAATCGCATGAAGTAAATCATCGCTTGGCACAAA‑TAGGGCGTTGTCATGATCaa                                                   >  1:1615246/1‑76 (MQ=255)                 ccccGCGCTGTGCCGCCGTGCAAATCGCATGAAGTAAATCACCGCTTGGCACAAA‑TAGGGCGTTGTCATGATCaaa                                                  >  2:417552/1‑76 (MQ=255)                  cccGCGCTGTGCCGCCGTGCAAATCGCATGAAGTAAATCATCGCTTGGCACAAA‑TAGGGCGTTGTCATGATCAAat                                                 <  1:127925/76‑1 (MQ=255)                       gcTGTGCCGCCGTGCAAATCGCATGAAGTAAATCATCGCTTGGCACAAA‑TAGGGCGTTGTCATGATCAAATATTCa                                            <  1:799569/76‑1 (MQ=255)                         tgtgCCGCCGTGCAAATCGCATGAAGTAAATCATCGCTTGGCACAAA‑TAGGGCGTTGTCATGATCAAATATTCAcg                                          <  2:735740/76‑1 (MQ=255)                             ccgccgTGCAAATCGCATGAAGTAAATCATCGCTTGGCACAAA‑TAGGGCGTTGTCATGATCAAATATTCACGCGcc                                      <  2:1062775/76‑1 (MQ=255)                                ccgTGCAAATCGCATGAAGTAAATCATCGCTTGGCACAAA‑TAGGGCGTTGTCATGATCAAATATCCACGCGCCGaa                                   <  1:56024/76‑1 (MQ=255)                                  gTGCAAATCGCATGAAGTAAATCATCGCTTGGCACAAA‑TAGGGCGTTGTCATGATCAAATATTCACGCGCCGAATa                                 <  2:747318/76‑1 (MQ=255)                                  gTGCAAATCGCATGAAGTAAATCATCGCTTGGCACAAA‑TAGGGCGTTGTCATGATCAAATATTCACGCGCCGAATa                                 >  1:415263/1‑76 (MQ=255)                                    gCAAATCGCATGAAGTAAATCATCGCTTGGCACAAA‑TAGGGCGTTGTCATGATCAAATATTCACGCGCCGAATAAg                               >  1:705950/1‑76 (MQ=255)                                      aaaTCNCGTTAAGTAAATCATCGCTTGGCCCAAA‑TAGGGCGTTGTCATGATCAAATATTCANGCGCCGAATAAGcc                             >  2:1832126/1‑76 (MQ=255)                                      aaaTCGCATGAAGTAAATGATNGCNTGGCACAGA‑TAGGGCGTTGTCNTGATCGANTNTTCACGCGCCGAATAAGcc                             >  2:258778/1‑76 (MQ=255)                                        aTCGCATGAAGTAAATCATCGCTTGGCACAAA‑TAGGGCGTTGTCATGATCAAATATTCACGCGCCGAATAAGCCGc                           <  2:1059825/76‑1 (MQ=255)                                          cGCATGAAGTAAATCATCGCTTGGCACAAA‑TAGGGCGTTGTCATGATCAAATATTCACGCGCCGAATAAGCCGCAg                         <  1:1377062/76‑1 (MQ=255)                                           gCATGAAGTAAATCATCGCTTGGCACAAA‑TAGGGCGTTGTCATGATCAAATATTCACGCGCCGAATAAGCCGCAGt                        >  2:1354425/1‑76 (MQ=255)                                            cATGAAGTAAATCATCGCTTGGCACAAA‑TAGGGCGTTGTCATGATCAAATATTCACGCGCCGAATAAGCCGCAGTc                       <  1:525721/76‑1 (MQ=255)                                             aTGAAGTAAATCATCGCTTGGCACAGA‑TAGGGCGTTGTCATGATCAAATATTCACGCGCCGAATAAGCCGTAGTCa                      <  1:118969/76‑1 (MQ=255)                                              tGAAGTAAATCATCGCTTGTCACAAA‑TAGGGCGTTGTCATGATCAAATATTCACGCGCCGAATAAGCCGCAGTCaa                     >  1:75539/1‑76 (MQ=255)                                              tGAAGTAAATCATCGCTTGGCACAAA‑TAGGGCGTTGTCATGATCAAATATTCACGCGCCGAATAAGCCGCAGTCaa                     >  1:1360421/1‑76 (MQ=255)                                               gAAGTAAATCATCGCTTGGCACAAA‑TAGGGCGTTGTCATGATCAAATATTCACGCGCCGAATAAGCCGCAGTCaat                    <  1:1175009/76‑1 (MQ=255)                                               gAAGTAAATCATCGCTTGGCACAAA‑TAGGGCGTTGTCATGATCAAATATTCACGCGCCGAATAAGCCGCAGTCaat                    <  1:1243768/76‑1 (MQ=255)                                                aaGTAAATCATCGCTTGGCACAAAGTAGGCACTTGTCATGATCAAATATTCACGCGCCGAATAAGCTGCAGTCaat                    >  2:1305508/1‑76 (MQ=255)                                                aaGTAAATCATCGCTTGGCACAAA‑TAGGGCGTTGTCATGATCAAATATTCACGCGCCGAATAAGCCGCAGTCaata                   >  1:1035327/1‑76 (MQ=255)                                                aaGTAAATCATCGCTTGGCACAAA‑TAGGGCGTTGTCATGATCAAATATTCACGCGCCGAATAAGCCGCAGTCaata                   >  1:1265467/1‑76 (MQ=255)                                                aaGTAAATCATCGCTTGGCACAAA‑TAGGGCGTTGTCATGATCAAATATTCACGCGCCGAATAAGCCGCAGTCaata                   <  2:733620/76‑1 (MQ=255)                                                aaGTAAATCATCGCTTGGCACAAA‑TAGGGCGTTGTCATGATCAAATATTCACGCGCCGAATAAGCCGCAGTCaata                   >  1:183997/1‑76 (MQ=255)                                                 aGTAAATCATCGCTTGGCACAAA‑TAGGGCGTTGTCATGATCAAATATTCACGCGCCGAATAAGCCGCAGTCaataa                  <  1:1415053/76‑1 (MQ=255)                                                  gTAAATCATCGCTTGGCACAAA‑TAGGGCGTTGTCATGATCAAATATTCACGCGCCGAATAAGCCGCAGTCaataat                 <  1:635624/76‑1 (MQ=255)                                                  gTAAATCATCGCGTGGCACAAA‑TAGGGCGTTGTCATGATCAAATATTCACGCGCCGAATAAGCCGCAGTCaataat                 <  2:603839/76‑1 (MQ=255)                                                    aaaTCATCGCTTGGCACAAA‑TAGGGCGTTGTCATGATCAAATATTCACGCGCCGAATAAGCCGCAGTCAATAATGc               <  1:736690/76‑1 (MQ=255)                                                    aaaTCATCGCTTGGCACAAA‑TAGGGCGTTGTCATGATCAAATATTCACGCGCCGAATAAGCCGCAGTCAATAATGc               >  2:1185750/1‑76 (MQ=255)                                                       tcatcGCTTGGCACAAA‑TAGGGCGTTGTCATGATCAAATATTCACGCGCCGAATAAGCCGCAGTCAATAATGCCTg            <  1:370153/76‑1 (MQ=255)                                                        catcGCTTGGCACAAA‑TAGGGCGTTGTCATGATCAAATATTCACGCGCCGAATAAGCCGCAGTCAATAATGCCTgg           >  1:332502/1‑76 (MQ=255)                                                           cGCTTGGCACAAA‑TAGGGCGTTGTCATGATCAAATATTCACGCGCCGAATAAGCCGCAGTCAATAATGCCTGGTGa        >  1:197100/1‑76 (MQ=255)                                                           cGCTTGGCACAAA‑TAGGGCGTTGTCATGATCAAATATTCACGCGCCGAATAAGCCGCAGTCAATAATGCCTGGTGa        >  2:811050/1‑76 (MQ=255)                                                            gCTTGGCACAAA‑TAGGGCGTTGTCATGATCAAATATTCACGCGCCGAATAAGCCGCAGTCAATAATGCCTGGTGaa       <  2:106111/76‑1 (MQ=255)                                                            gCTTGGCACAAA‑TAGGGCGTTGTCATGATCAAATATTCACGCGCCGAATAAGCCGCAGTCAATAATGCCTGGTGaa       >  1:738449/1‑76 (MQ=255)                                                               tGGCACAAA‑TAGGGCGTTGTCATGATCAAATATTCACGCGCCGAATAAGCCGCAGTCAATAATGCCTGGTGAATga    >  1:1398093/1‑76 (MQ=255)                                                                ggCACAAA‑TAGGGCGTTGTCATGATCAAATATTCACGCGCCGAATAAGCCGCAGTCAATAATGCCTGNTGAATgag   <  1:1765055/76‑1 (MQ=255)                                                                 gCACAAA‑TAGGGCGTTGTCATGATCAAATATTCACGCGCCGAATAAGCCGCAGTCAATAATGCCTGGTGAATgaga  <  2:1282016/76‑1 (MQ=255)                                                                        |                                                                      ATAATACTGACATCCACCCCGCGCTGTGCCGCCGTGCAAATCGCATGAAGTAAATCATCGCTTGGCACAAAGTAGGGCGTTGTCATGATCAAATATTCACGCGCCGAATAAGCCGCAGTCAATAATGCCTGGTGAATGAGA  >  REL606/1305140‑1305280

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