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

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	1,305,475	0	C	.	91.7%	93.2 / 4.0	24	coding (765/1461 nt)	cls	cardiolipin synthetase
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base . (14/8);  total (15/9)
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.92e-01

CAGTCGCAGGAATAAATAATCCCCATCGCGGTGGCGATGGGGCCTTCCATACGCGCCATCAGATCAATCCATTGCCCTACGCCCGCATCTTGTTTGAAGTAGCGAGGGTCGACCATATTCATGCTGCCGGTGTACGCGATGTAATT  >  REL606/1305401‑1305546                                                                           |                                                                        cAGTCGCAGGAATAAATAATCCCCATCGCGGTGGCGATGGGGCCTTCCATACGCGCCATCAGATCAATCCATTGcc                                                                        >  1:809224/1‑76 (MQ=255)  aGTCGCAGGAATAAATAACCCCCATCGCGGGGGCGATGGGGCCTTCCATACGCGCCATCAGATCAATCCATTGcct                                                                       <  2:61348/76‑2 (MQ=255)       cAGGAATAAATAATCCCCATCGCGGCGGCGATGGGGCCTTCCATACGCGCCATCAGATCAATCCATTG‑CCTACGcc                                                                 >  1:1094878/1‑76 (MQ=255)            ataaataaTCCCCATCGCGGTGGCGATGGGGCCTTCCATACGCGCCATCAGATCAATCCATTG‑CCTACGCCCGCAt                                                            >  2:926960/1‑76 (MQ=255)              aaataaTCCCCATCGCGGTGGCGATGGGGCCTTCCATACGCGCCATCAGATCAATCCATTGCCCTACGCCCGCATc                                                           >  2:476897/1‑76 (MQ=255)              aaataaTCCCCATCGCGGTGGCGATGGGGCCTTCCATACGCGCCATCAGATCAATCCATTG‑CCTACGCCCGCATCt                                                          >  2:678441/1‑76 (MQ=255)                   atCCCCATCGCGGTGGCGATGGGGCCTTCCATACGCGCCATCAGATCAATCCATTG‑CCTACGCCCGCATCTTGttt                                                     >  2:1037044/1‑76 (MQ=255)                   acccccACCGGGGGGGCGAGGGGGCCTTCCATACGCGCCATCAGATCAATCCATTG‑CCTACGCCCGCATCTTGttt                                                     <  2:1531493/74‑1 (MQ=255)                    tCCCCATCGCCGAGGCGATGGGGCCTTCCATACGCCCCATCAGATCAATCCATTG‑CCTACGCCCGCGTCACGTTTg                                                    >  2:1277424/1‑76 (MQ=255)                     ccccATCGCGGTGGCGATGGGGCCTTCCATACGCGCCATCAGATCAATCCATTG‑CCTACGCCCGCATCTTGTTTGa                                                   >  2:766167/1‑76 (MQ=255)                     ccccATCGCGGTGGCGATGGGGCCTTCCATACGCGCCATCAGATCAATCCATTG‑CCTACGCCCGCATCTTGTTTGa                                                   >  2:1373346/1‑76 (MQ=255)                      cccATCGCGGTGGCGATTGGGCCTTCCATACGCGCCATCAGATCAATCCCTTG‑CCTACGCCCGCATCTTGTTTGaa                                                  >  2:512983/1‑76 (MQ=255)                           cccgGTGGTGATTGCCCCTTCCATACGCGCCATCAGATCAATCCATTGCCCTACGCCCGCATCTTGTTTGAAGTAg                                              <  2:1011095/74‑1 (MQ=255)                                    gATGGGGCCTTCCATACGCGCCATCAGATCAATCCATTG‑CCTACGCCCGCATCTTGTTTGAAGTAGCGAGGGTCGa                                    >  1:492239/1‑76 (MQ=255)                                       ggggCCTTCCATACGCGCCATCAGATCAATCCATTG‑CCTACGCCCGCATCTTGTTTGAAGTAGCGAGGGGCGACCa                                 >  1:495552/1‑76 (MQ=255)                                         ggCCTTCCATACGCGCCATCAGATCAATCCATTG‑CCTACGCCCGCATCTTGTTTGAAGTAGCGAGGGTCGACCata                               <  1:14537/76‑1 (MQ=255)                                           ccTTCCATACGCGCCATCAGATCAATCCATTG‑CCTACGCCCGCATCTTGTTTGAAGTAGCGAGGGTCGACCATAtt                             <  2:888967/76‑1 (MQ=255)                                           ccTTCCATACGCGCCATCAGATCAATCCATTG‑CCTACGCCCGCATCTTGTTTGAAGTAGCGAGGGTCGACCATAtt                             <  2:800122/76‑1 (MQ=255)                                            cTTCCATACGCGCCATCAGATCAATCCATTG‑CCTACGCCCGCATCTTGTTTGAAGTAGCGAGGGTCGACCATATTc                            >  1:19643/1‑76 (MQ=255)                                               ccATACGCGCCATCAGATCAATCCATTG‑CCTACGCCCGCATCTTGTTTGAAGTAGCGAGGGTCGACCATATTCAtg                         >  1:42026/1‑76 (MQ=255)                                                 aTACGCGCCATCAGATCAATCCATTG‑CCTACGCCCGCATCTTGTTTGAAGTAGCGAGGGTCGACCATATTCAtgct                       >  1:378662/1‑76 (MQ=255)                                                   aCGCGCCATCAGATCAATCCATTG‑CCTACGCCCGCATCTTGTTTGAAGTAGCGAGGGTCGACCATATTCAtgctgc                     >  1:433978/1‑76 (MQ=255)                                                      cgcCATCAGATCAATCCATTG‑CCTACGCCCGCATCTTGTTTGAAGTAGCGAGGGTCGACCATATTCATGCTGCCgg                  <  2:523234/76‑1 (MQ=255)                                                         cATCAGATCAATCCATTG‑CCTACGCCCGCATCTTGTTTGAAGTAGCGAGGGTCGACCATATTCATGCTGCCGgtgt               >  1:931432/1‑76 (MQ=255)                                                         cATCAGATCAATCCATTG‑CCTACGCCCGCATCTTGTTTGAAGTAGCGAGGGTCGACCATATTCATGCTGCCGgtgt               >  1:643438/1‑76 (MQ=255)                                                          aTCAGATCAATCCATTG‑CCTACGCCCGCATCTTGTTTGAAGTAGCGAGGGTCGACCATATTCATGCTGCCGGTGTa              >  1:554992/1‑76 (MQ=255)                                                                   aTCCATTG‑CCTACGCCCGCATCTTGTTTGAAGTAGCGAGGGTCGACCATATTCATGCTGCCGGTGTACGCGATGTa     <  2:1390642/76‑1 (MQ=255)                                                                    tCCATTG‑CCTACGCCCGCATCTTGTTTGAAGTAGCGAGGGTCGACCATATTCATGCTGCCGGTGTACGCGATGTaa    <  2:1575779/76‑1 (MQ=255)                                                                      cATTG‑CCTACGCCCGCATCTTGTTTGAAGTAGCGAGGGTCGACCATATTCATGCTGCCGGTGTACGCGATGTAAtt  <  2:221767/76‑1 (MQ=255)                                                                           |                                                                        CAGTCGCAGGAATAAATAATCCCCATCGCGGTGGCGATGGGGCCTTCCATACGCGCCATCAGATCAATCCATTGCCCTACGCCCGCATCTTGTTTGAAGTAGCGAGGGTCGACCATATTCATGCTGCCGGTGTACGCGATGTAATT  >  REL606/1305401‑1305546

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