Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	1,304,966	C→T	R425Q (CGG→CAG)	cls ←	cardiolipin synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	1,304,966	0	C	T	98.0%	190.4 / ‑3.8	51	R425Q (CGG→CAG)	cls	cardiolipin synthetase
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (26/24);  total (26/25)
Fisher's exact test for biased strand distribution p-value = 4.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.76e-01

AACGGCGGCAAGGTCTGCACCAAAACCTTTATCGTCGATTGCCAGAGTAATCTCGAAATTTAGCCACAGACTCCGCATATCAAGGTTAACGGTGCCAACCAGACTTAGTTCACCATCGACCAGCACGCTCTTGGTATGCAGTAACC  >  REL606/1304893‑1305038                                                                          |                                                                         aaCGGCGGCAAGGTCTGCACCAAAACCTTTATCGTCGATTGCCAGAGTAATCTCGAAATTTAGCCACAGACTCTGc                                                                        <  2:1385598/76‑1 (MQ=255)       ggcAAGGTCTGCACCAAAACCTTTATCGTCGATTGCCAGAGTAATCTCGAAATTTAGCCACAGACTCTGCATATCa                                                                  <  1:797773/76‑1 (MQ=255)        gcAAGGTCTGCACCAAAACCTTTATCGTCGATTGCCAGAGTAATCTCGAAATTTAGCCACAGACTCTGCATATCaa                                                                 <  1:1590022/76‑1 (MQ=255)         cAAGGTCTGCACCAAAACCTTTATCGTCGATTGCCAGAGTAATCTCGAAATTTAGCCACAGACTCTGCATATCAAg                                                                >  2:289465/1‑76 (MQ=255)          aaGGTCTGCACCAAAACCTTTATCGTCGATTGCCAGAGTAATCTCGAAATTTAGCCACAGACTCTGCATATCAAgg                                                               >  1:1150420/1‑76 (MQ=255)          aaGGTCTGCACCAAAACCTTTATCGTCGATTGCCAGAGTAATCTCGAAATTTAGCCACAGACTCTGCATATCAAgg                                                               >  2:1947629/1‑76 (MQ=255)           aGGTCTGCACCAAAACCTTTATCGTCGATTGCCAGAGTAATCTCGAAATTTAGCCACAGACTCTGCATATCAAGGt                                                              <  1:519184/76‑1 (MQ=255)                  cACCAAAACCTTTATCGTCGATTGCCAGAGTAATCTCGAAATTTAGCCACAGACTCTGCATATCAAGGTTAACGGt                                                       <  1:511059/76‑1 (MQ=255)                        aaCCTTTATCGTCGATTGCCAGAGTAATCTCGAAATTTAGCCACAGACTCTGCATATCAAGGTTAACGGTGCCAAc                                                 >  2:1547696/1‑76 (MQ=255)                        aaCCTTTATCGTCGATTGCCAGAGTAATCTCGAAATTTAGCCACAGACTCTGCATATCAAGGTTAACGGTGCCAAc                                                 >  1:1852053/1‑76 (MQ=255)                         aCCTTTATCGTCGATTGCCAGAGTAATCTCGAAATTTAGCCACAGACTCTGCATATCAAGGTTAACGGTGCCAAcc                                                <  2:479023/76‑1 (MQ=255)                          ccTTTATCGTCGATTGCCAGAGTAATCTCGAAATTTAGCCACAGACTCTGCATATCAAGGTTAACGGTGCCAACCa                                               <  1:160884/76‑1 (MQ=255)                          ccTTTATCGTCGATTGCCAGAGTAATCTCGAAATCTAGCCACAGACTCTGCATATCAAGGTTAACGGTGCCAACCa                                               >  2:927747/1‑76 (MQ=255)                           cTTTATCGTCGATTGCCAGAGTAATCTCGAAATTTAGCCACAGACTCTGCATATCAAGGTTAACGGTGCCAACCAg                                              <  1:1713533/76‑1 (MQ=255)                             ttATCGTCGATTGCCAGAGTAATCGCGAAATTTAGCCACAGACTCTGCATATCAAGGTTAACGGTGCCAACCAGAc                                            <  1:925534/76‑1 (MQ=255)                              tATCGTCGATTGCCAGAGTAATCTCGAAATTTAGCCACAGACTCTGCATATCAAGGTTAACGGTGCCAACCAGACt                                           <  1:15087/76‑1 (MQ=255)                              tATCGTCGATTGCCAGAGTAATCTCGAAATTTAGCCACAGACTCTGCATATCAAGGTTAACGGTGCCAACCAGACt                                           >  1:443306/1‑76 (MQ=255)                               aTCGTCGATTGCCAGAGTAATCTCGAAATTTAGCCACAGACTCTGCATATCAAGGTTAACGGTGCCAACCAGACtt                                          >  2:733644/1‑76 (MQ=255)                                 cgtcgATTGCCAGAGTAATCTCGAAATTTAGCCACAGACTCTGCATATCAAGGTTAACGGTGCCAACCAGACTTAg                                        >  1:1618596/1‑76 (MQ=255)                                 cgtcgATTGCCAGAGTAATCTCGAAATTTAGCCACAGACTCTGCATATCAAGGTTAACGGTGCCAACCAGACTTAg                                        <  2:1517997/76‑1 (MQ=255)                                 cgtcgATTGCCAGAGTAATCTCGAAATTTAGCCACAGACTCTGCATATCAAGGTTAACGGTGCCAACCAGACTTAg                                        <  2:1985216/76‑1 (MQ=255)                                    cgATTGCCAGAGTAATCTCGAAATTTAGCCACAGACTCTGCATATCAAGGTTAACGGTGCCAACCAGACTTAGTTc                                     >  1:441300/1‑76 (MQ=255)                                     gATTGCCAGAGTAATCTCGAAATTTAGCCACAGACTCTGCATATCAAGGTTAACGGTGCCAACCAGACTTAGTTCa                                    >  2:1276637/1‑76 (MQ=255)                                        tGCCAGAGTAATCTCGAAATTTAGCCACAGACTCTGCATATCAAGGTTAACGGTGCCAACCAGACTTAGTTCACCa                                 >  2:459782/1‑76 (MQ=255)                                          ccAGAGTAATCTCGAAATTTAGCCACAGACTCTGCATATCAAGGTTAACGGTGCCAACCAGACTTAGTTCACCATc                               >  2:385243/1‑76 (MQ=255)                                           cAGAGTAATCTCGAAATTTAGCCACAGACTCTGCATATCAAGGTTAACGGTGCCAACCAGACTTAGTTCACCATCg                              <  2:891/76‑1 (MQ=255)                                            agagTAATCTCGAAATTTAGCCACAGACTCTGCATATCAAGGTTAACGGTGCCAACCAGACTTAGTTCACCATCGa                             >  1:1017890/1‑76 (MQ=255)                                            agagTAATCTCGAAATTTAGCCACAGACTCTGCATATCAAGGTTAACGGTGCCAACCAGACTTAGTTCACCATCGa                             >  2:579049/1‑76 (MQ=255)                                             gagTAATCTCGAAATTTAGCCACAGACTCTGCATATCAAGGTTAACGGTGCCAACCAGACTTAGTTCACCATCAAc                            >  1:172581/1‑76 (MQ=255)                                               gTAATCTCGAAATTTAGCCACAGACTCTGCATATCAAGGTTAACGGTGCCAACCAGACTTAGTTCACCATCGACCa                          >  2:1821456/1‑76 (MQ=255)                                                 aaTCTCGAAATTTAGCCACAGACTCTGCATATCAAAGTTAACGGTGCCAACCAGACTTAGTTCACCATCGACCAGc                        >  2:1807364/1‑76 (MQ=255)                                                   tctcGAAATTTAGCCACAGACTCTGCATATCAAGGTTAACGGTGCCAACCAGACTTAGTTCACCATCGACCAGCAc                      <  2:1901754/76‑1 (MQ=255)                                                   tctcGAAATTTAGCCACAGACTCCGCATATCAAGGTTAACGGTGCCAACCAGACTTAGTTCACCATCGACCAGCAc                      <  1:1662067/76‑1 (MQ=255)                                                    ctcGAAATTTAGCCACAGACTCTGCATATCAAGGTTAACGGTGCCAACCAGACTTAGTTCACCATCGACCAGCACg                     >  1:382614/1‑76 (MQ=255)                                                    ctcGAAATTTAGCCACAGACTCTGCATATCAAGGTTAACGGTGCCAACCAGACTTAGTTCACCATCGACCAGCACg                     >  1:972510/1‑76 (MQ=255)                                                    ctcGAAATTTAGCCACAGACTCTGCATATCAAGGTTAACGGTGCCAACCAGACTTAGTTCACCATCGACCAGCACg                     <  1:1580026/76‑1 (MQ=255)                                                    ctcGAAATTTAGCCACAGACTCTGCATATCAAGGTTAACGGTGCCAACCAGACTTAGTTCACCATCGACCAGCACg                     <  1:1463367/76‑1 (MQ=255)                                                    ctcGAAATTTAGCCACAGACTCTGCATATCAAGGTTAACGGTGCCAACCAGACTTAGTTCACCATCGACCAGCACg                     >  1:1228923/1‑76 (MQ=255)                                                      cGAAATTTAGCCACAGACTCTGCATATCAAGGTTAACGGTGCCAACCAGACTTAGTTCACCATCGACCAGCACGct                   <  2:940087/76‑1 (MQ=255)                                                         aaTTTAGCCACAGACTCTGCATATCAAGGTTAACGGTGCCAACCAGACTTAGTTCACCATCGACCAGCACGCTCtt                <  2:448779/76‑1 (MQ=255)                                                         aaTTTAGCCACAGACTCTGCATATCAAGGTTAACGGTGCCAACCAGACTTAGTTCACCATCGACCAGCACGCTCtt                >  1:1019234/1‑76 (MQ=255)                                                          aTTTAGCCACAGACTCTGCATATCAAGGTTAACGGTGCCAACCAGACTTAGTTCACCATCGACCAGCACGCTCTTg               >  1:1012080/1‑76 (MQ=255)                                                            ttAGCCACAGACTCTGCATATCAAGGTTAACGGTGCCAACCAGACTTAGTTCACCATCGACCAGCACGCTCTTGGt             >  1:709787/1‑76 (MQ=255)                                                             tAGCCACAGACTCTGCATATCAAGGTTAACGGTGCCAACCAGACTTAGTTCACCATCGACCAGCACGCTCTTGGTa            <  1:260863/76‑1 (MQ=255)                                                              aGCCACAGACTCTGCATATCAAGGTTAACGGTGCCAACCAGACTTAGTTCACCATCGACCAGCACGCTCTTGGTAt           <  1:585523/76‑1 (MQ=255)                                                              aGCCACAGACTCTGCATATCAAGGTTAACGGTGCCAACCAGACTTAGTTCACCATCGACCAGCACGCTCTTGGTAt           <  1:490479/76‑1 (MQ=255)                                                              aGCCACAGACTCTGCATATCAAGGTTAACGGTGCCAACCAGACTTAGTTCACCATCGACCAGCACGCTCTTGGTAt           >  2:799827/1‑76 (MQ=255)                                                              aGCCACAGACTCTGCATATCAAGGTTAACGGTGCCAACCAGACTTAGTTCACCATCGACCAGCACGCTCTTGGCAt           >  2:1488515/1‑76 (MQ=255)                                                                  acaGACTCTGCATATCAAGGTTAACGGTGCCAACCAGACTTAGTTCACCATCGACCAGCACGCTCTTGGTATGCAg       <  1:186530/76‑1 (MQ=255)                                                                    aGACTCTGCATATCAAGGTTAACGGTGCCAACCAGACTTAGTTCACCATCGACCAGCACGCTCTTGGTCTGCAGTa     <  1:1276637/76‑1 (MQ=255)                                                                    aGACTCTGCATATCAAGGTTAACGGTGCCAACCAGACTTAGTTCACCATCGACCAGCACGCTCTTGGTATGCAGTa     <  1:348047/76‑1 (MQ=255)                                                                     gACTCTGCATATCAAGGTTAACGGTGCCGACCAGACTTAGTTCACCATCGACCAGCACGCTCTTGGTATGCAGTaa    <  2:1408252/76‑1 (MQ=255)                                                                     gACTCTGCATATCAAGGTTAACGGTGCCAACCAGACTTAGTTCACCATCGACCAGCACGCTCTTGGTATGCAGTaa    >  1:515371/1‑76 (MQ=255)                                                                       ctctGCATATCAAGGTTAACGGTGCCAACCAGACTTAGTTCACCATCGACCAGCACGCTCTTGGTATGCAGTAAcc  >  1:1024159/1‑76 (MQ=255)                                                                          |                                                                         AACGGCGGCAAGGTCTGCACCAAAACCTTTATCGTCGATTGCCAGAGTAATCTCGAAATTTAGCCACAGACTCCGCATATCAAGGTTAACGGTGCCAACCAGACTTAGTTCACCATCGACCAGCACGCTCTTGGTATGCAGTAACC  >  REL606/1304893‑1305038

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: ‑