Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	1,304,826	G→A	Q472* (CAA→TAA)	cls ←	cardiolipin synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	1,304,826	0	G	A	97.8%	169.0 / ‑3.6	46	Q472* (CAA→TAA)	cls	cardiolipin synthetase
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (24/21);  total (25/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.20e-01

CCTGTTTAACCTCTGTTGGCGACGTTTTACAGCAACGGGCTGAAGAAGTAAAACAGTCGCTCGGCGACACGTTGCCATAATGGACGTTTTAGCCATAAACGGGCATCGAGCAGACGTGAACGCGAAATATAATCGTCCTGAACG  >  REL606/1304753‑1304896                                                                          |                                                                       ccTGTTTAACCTCTGTTGGCGACGTTTTACAGCAACGGGCTGAAGAAGTAAAACAGTCGCTCGGCGACACGTTAcc                                                                      <  2:1607864/76‑1 (MQ=255) ccTGTTTAACCTCTGTTGGCGACGTTTTACAGCAACGGGCTGAAGAAGTAAAACAGTCGCTCGGCGACACGTTAcc                                                                      >  2:432986/1‑76 (MQ=255)   tGTTTAACCTCTGTTGGCGACGTTTTACAGCAACGGGCTGAAGAAGTAAAACAGTCGCTCGGCGACACGTTACCAt                                                                    >  1:154389/1‑76 (MQ=255)    gTTTAACCTCTGTTGGCGACGTTTTACAGCAACGGGCTGAAGAAGTAAAACAGTCGCTCGGCGACACGTTACCATa                                                                   >  1:1388252/1‑76 (MQ=255)     tttAACCTCTGTTGGCGACGTTTTACAGCAACGGGCTGAAGAAGTAAAACAGTCGCTCGGCGACACGTTACCATaa                                                                  <  1:1084243/76‑1 (MQ=255)     tttAACCTCTGTTGGCGACGTTTTACAGCAACGGGCTGAAGAAGTAAAACAGTCGCTCGGCGACACGTTACCATaa                                                                  >  2:1435399/1‑76 (MQ=255)      ttAACCTCTGTTGGCGACGTTTTACAGCAACGGGCTGAAGAAGTAAAACAGTCGCTCGGCGACACGTTACCATAAt                                                                 <  1:290537/76‑1 (MQ=255)      ttAACCTCTGTTGGCGACGTTTTACAGCAACGGGCTGAAGAAGTAAAACAGTCGCTCGGCGACACGTTACCATAAt                                                                 <  2:199722/76‑1 (MQ=255)       tAACCTCTGTTGGCGACGTTTTACAGCAACGGGCTGAAGAAGTAAAACAGTCGCTCGGCGACACGTTACCATAATg                                                                >  2:609946/1‑76 (MQ=255)        aaCCTCTGTTGGCGACGTTTTACAGCAACGGGCTGAAGAAGTAAAACAGTCGCTCGGCGACACGTTACCATAATgg                                                               <  1:1319521/76‑1 (MQ=255)          ccTCTGTTGGCGACGTTTTACAGCAACGGGCTGAAGAAGTAAAACAGTCGCTCGGCGACATGTTACCATAATGGAc                                                             >  2:6827/1‑76 (MQ=255)             ctGTTGGCGACGTTTTACAGCAACGGGCTGAAGAAGTAAAACAGTCGCTCGGCGACACGTTACCATAATGGACGtt                                                          >  2:1095971/1‑76 (MQ=255)                 tGGCGACGTTTTACAGCAACGGGCTGAAGAAGTAAAACAGTCGCTCGGCGACACGTTACCATAATGGACGTTTTAg                                                      >  2:1318119/1‑76 (MQ=255)                         ttttACACCACCGGGCTGAAGAAGAAAAACAGTCGCTCGGCGACACGTTACCATAATGGACGTTTTAGCCATAAAc                                              <  1:594398/76‑1 (MQ=255)                          tttACAGCAACGGGCTGAAGAAGTAAAACAGTCGCTCGGCGACACGTTACCATAATGGACGTTTTAGCCATAAACg                                             >  1:149446/1‑76 (MQ=255)                           ttACAGCAACGGGCTGAAGAAGTAAAACAGTCGCTCGGCGACACGTTACCATAATGGACGTTTTAGCCATAAACgg                                            >  1:639428/1‑76 (MQ=255)                            tACAGCAACGGGCTGAAGAAGTAAAACAGTCGCTCGGCGACACGTTACCATAATGGACGTTTTAGCCATAAACggg                                           <  1:1426892/76‑1 (MQ=255)                             aCAGCATCGGGCGAAAAAAGAAAAACAGTCGCTCGGCGGCACGTTACCATAATGGACGTTTTAGCCATAAACGGGc                                          <  2:302454/76‑1 (MQ=255)                              cAGCAACGGGCTGAAGAAGTAAAACAGTCGCTCGGCGACACGTTACCATAATGGACGTTTTAGCCATAAACGGGCa                                         >  2:511002/1‑76 (MQ=255)                               aGCAACGGGCTGAAGAAGTAAAACAGTCGCTCGGCGACACGTTACCATAATGGACGTTTTAGCCATAAACGGGCAt                                        <  2:833918/76‑1 (MQ=255)                                 cAACGGGCTGAAGAAGTAAAACAGTCGCTCGGCGACACGTTACCATAATGGACGTTTTAGCCATAAACGGGCATCg                                      <  1:1655096/76‑1 (MQ=255)                                  aaCGGGCTGAAGAAGTAAAACAGTCGCTCGGCGACACGTTACCATAATGGACGTTTTAGCCATAAACGGGCATCGa                                     >  1:819678/1‑76 (MQ=255)                                    cGGGCTGAAGAAGTAAAACAGTCGCTCGGCGACACGTTACCATAATGGACGTTTTAGCCATAAACGGGCATCGAGc                                   <  1:118392/76‑1 (MQ=255)                                     gggCTGAAGAAGTAAAACAGTCGCTCGGCGACACGTTACCATAATGGACGTTTTAGCCATAAACGGGCATCGAGCa                                  >  1:1400674/1‑76 (MQ=255)                                      ggCTGAAGAAGTAAAACAGTCGCTCGGCGACACGTTACCATAATGGACGTTTTAGCCATAAACGGGCATCGAGCac                                 >  1:479716/1‑75 (MQ=255)                                        cTGAAGAAGTAAAACAGTCGCTCGGCGACACGTTACCATAATGGACGTTTTAGCCATAAACGGGCATCCAGCAGAc                               >  1:925317/1‑76 (MQ=255)                                          gaagaaGTAAAACAGTCGCTCGGCGACACGTTACCATAATGGACGTTTTAGCCATAAACGGGCATCGAGCAGACgc                             >  1:953448/1‑75 (MQ=255)                                          gaagaaGTAAAACAGTCGCTCGGCGACACGTTACCATAATGGACGTTTTAGCCATAAACGGGCATCGAGCAGACGt                             >  1:614822/1‑76 (MQ=255)                                           aaaaagTAAAACAGTCGCTCGGCGACACGTTACCATAATGGACGTTTTAGCCATAAACGGGCATCGAGCAGACGTg                            <  1:1099434/73‑1 (MQ=255)                                            agaagTAAAACAGTCGCTCGGCGACACGTTACCATAATGGACGTTTTAGCCATAAACGGGCATCGAGCAGACGTGa                           >  1:1593348/1‑76 (MQ=255)                                               agTAAAACAGTCGCTCGGCGACACGTTACCATAATGGACGTTTTAGCCATAAACGGGCATCGAGCAGACGTGAAcg                        >  2:1399889/1‑76 (MQ=255)                                                 tAAAACAGTCGCTCGGCGACACGTTACCATAATGGACGTTTTAGCCATAAACGGGCATCGAGCAGACGTGAAcgcg                      >  1:582448/1‑76 (MQ=255)                                                 tAAAACAGTCGCTCGGCGACACGTTACCATAATGGACGTTTTAGCCATAAACGGGCATCGAGCAGACGTGAAcgcg                      <  2:523980/76‑1 (MQ=255)                                                  aaaaCAGTCGCTCGGCGACACGTTGCCATAATGGACGTTTTAGCCATAAACGGGCATCGAGCAGACGTGAACGCGa                     >  1:946158/1‑76 (MQ=255)                                                   aaaCAGTCGCTCGGCGACACGTTACCATAATGGACGTTTTAGCCATAAACGGGCATCGAGCAGACGTGACAGCGaa                    >  2:274618/1‑76 (MQ=255)                                                   aaaCAGTCGCTCGGCGACACGTTACCATAATGGACGTTTTAGCCATAAACGGGCATCGAGCAGACGTGAACGCGaa                    <  2:1429514/76‑1 (MQ=255)                                                      cAGTCGCTCGGCGACACGTTACCATAATGGACGTTTTAGCCATAAACGGGCATCGAGCAGACGTGAACGCGAAata                 <  2:1079956/76‑1 (MQ=255)                                                      cAGTCGCTCGGCGACACGTTACCATAATGGACGTTTTAGCCATAAACGGGCATCGAGCAGACGTGAACGCGAAata                 >  1:249603/1‑76 (MQ=255)                                                       ggTCGCTCGGCGACACGTTACCATAATGGACGTTTTAGCCATAAACGGGCATCGAGCAGACGTGAACGCGAAatat                <  1:237391/75‑1 (MQ=255)                                                        gTCGCTCGGCGACACGTTACCATAATGGACGTTTTAGCCATAAACGGGCATCGAGCAGACGTGAACGCGAAAtata               <  1:593030/76‑1 (MQ=255)                                                          cGCTCGGCGACACGTTACCATAATGGACGTTTTAGCCATAAACGGGCATCGAGCAGACGTGAACGCGAAATATAAt             <  2:1193804/76‑1 (MQ=255)                                                          cGCTCGGCGACACGTTACCATAATGGACGTTTTAGCCATAAACGGGCATCGAGCAGACGTGAACGCGAAATATAAt             >  1:807483/1‑76 (MQ=255)                                                          cGCTCGGCGACACGTTACCATAATGGACGTTTTAGCCATAAACGGGCATCGAGCAGACGTGAACGCGAAATATAAt             <  2:691134/76‑1 (MQ=255)                                                          cGCTCGGCGACACGTTACCATAATGGACGCTTTAGCCATAAACGGGCATCGAGCAGACGTGAACGCGAAATATAAt             >  1:492371/1‑76 (MQ=255)                                                            cTCGGCGACACGTTACCATAATGGACGTTTTAGCCATAAACGGGCATCGAGCAGACGTGAACGCGAAATATAATCg           <  1:1099532/76‑1 (MQ=255)                                                                gCGACACGTTACCATAATGGACGTTTTAGCCATAAACGGGCATCGAGCAGACGTGAACGCGAAATATAATCGTCCt       <  2:255739/76‑1 (MQ=255)                                                                    cacGTTACCATAATGGGCGTTTTAGCCATAAACGGGCATCGAGCAGACGTGAACGCGAAATATAATCGTCCTGAAc   <  1:744297/76‑1 (MQ=255)                                                                     acGTTACCATAATGGACGTTTTAGCCATAAACGGGCATCGAGCAGACGTGAACGCGAAATATAATCGTCCTGAAcg  >  1:480899/1‑76 (MQ=255)                                                                          |                                                                       CCTGTTTAACCTCTGTTGGCGACGTTTTACAGCAACGGGCTGAAGAAGTAAAACAGTCGCTCGGCGACACGTTGCCATAATGGACGTTTTAGCCATAAACGGGCATCGAGCAGACGTGAACGCGAAATATAATCGTCCTGAACG  >  REL606/1304753‑1304896

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