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

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	1,305,545	0	T	.	97.8%	217.7 / ‑0.8	46	coding (695/1461 nt)	cls	cardiolipin synthetase
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base . (21/24);  total (22/24)
Fisher's exact test for biased strand distribution p-value = 4.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.49e-01

ATTGCCCTACGCCCGCATCTTGTTTGAAGTAGCGAGGGTCGACCATATTCATGCTGCCGGTGTACGCGATGTAATTATCGATCATG‑‑‑ATCATCTTGCGATGTTGGCGCAGGTCCATACGGCGTAAAAACACACGCATCAGATTGACC  >  REL606/1305471‑1305616                                                                           |                                                                           aTTGCCCTACGCCCGCATCTTGTTTGAAGTAGCGAGGGTCGACCATATTCATGCTGCCGGTGTACGCGATGTAAta                                                                           <  2:1148111/76‑2 (MQ=255)   tGCCCTATGCCCGCATGTTGTTTGAAGTAGCGAGGGTCGTCCATATTCATGCTGCCGGTGTACGCGATGTAAtatc                                                                         <  1:29724/76‑4 (MQ=255)   tGCCCTACGTCCGTATCTTGTTTGAAGTAGCNCGGGACGACCATATTCATGCTGCCGGTGTACTNGATGTAAtatc                                                                         >  1:243191/1‑73 (MQ=255)       cTACGCCCGCATCTTGTTTGAAGTAGCGAGGGTCGACCATATTCATGCTGCCGGTGTACGCGATGTAA‑TATCgatc                                                                    <  1:1862403/76‑1 (MQ=255)       cTACGCCCGCATCTTGTTTGAAGTAGCGAGGGTCAACCATATTCATGCTGCCGGTGTACGCGATGTAA‑TATCgatc                                                                    >  2:1165418/1‑76 (MQ=255)        tACGCCCGCATCTTGTTTGAAGTAGCGAGGGTCGACCATATTCATGCTGCCGGTGTACGCGATGTAA‑TATCgatca                                                                   >  1:556163/1‑76 (MQ=255)         aCGCCCGCATCTTGTTTGAAGTAGCGAGGGTCGACCATATTCATGCTGCCGGTGTACGCGATGTAA‑TATCgatcat                                                                  >  2:1449845/1‑76 (MQ=255)         aCGCCCGCATCTTGTTTGAAGTAGCGAGGGTCGACCATATTCATGCTGCCGGTGTACGCGATGTAA‑TATCgatcat                                                                  >  1:1738841/1‑76 (MQ=255)           gCCCGCATCTTGTTTGAAGTAGCGAGGGTCGACCATATTCATGCTGCCGGTGTACGCGATGTAA‑TATCgatcatg‑‑‑a                                                             <  1:1023794/76‑1 (MQ=255)           gCCCGCATCTTATTTGAAGTAGCGAGGNTCGACCANNATCATGCTGCCTGTGTACGCAATGTAA‑TATCgatcatg‑‑‑a                                                             >  1:2007733/1‑76 (MQ=255)             ccGCATCTTGTTTGAAGTAGCGAGGGTCGACCATATTCATGCTGCCGGTGTACGCGATGTAA‑TATCgatcatg‑‑‑atc                                                           <  1:336537/76‑1 (MQ=255)                  tCTTGTTTGAAGTAGCGAGGGTCGACCATATTCATGCTGCCGGTGTACGCGATGTAA‑TATCGATCATG‑‑‑ATCATCtt                                                      <  1:278193/76‑1 (MQ=255)                  tCTTGTTTGAAGTAGCGAGGGTCAACCATATTCATGCTGCCGGTGTACNCGATGTAA‑TATCGATCATG‑‑‑ATCATCtt                                                      <  2:8702/76‑1 (MQ=255)                     tGTTTGAAGTAGCGAGGGTCGACCATATTCATGCTGCCGGTGTACGCGATGTAA‑TATCGATCATG‑‑‑ATCATCTTGCg                                                   <  1:855449/76‑1 (MQ=255)                     tGTTTGAAGTAGCGAGGGTCGACCATATTCATGCTGCCGGTGTACGCGATGTAA‑TATCGATCATG‑‑‑ATCATCTTGCg                                                   >  2:595439/1‑76 (MQ=255)                     tGTTTGAAGTAGCGAAGGTCGACCATATTCATGCTGCCGGTGTACGCGATGTAA‑TATCGATCATG‑‑‑ATCATCTTGCg                                                   <  1:1694115/76‑1 (MQ=255)                      gTTTGAAGTAGCGAGGGTCGACCATATTCATGCTGCCGGTGTACGCGATGTAA‑TATCGATCATG‑‑‑ATCATCTTGCGa                                                  <  2:1789732/76‑1 (MQ=255)                       tttGAAGTAGCGAGGGTGGACCATATTCATGCTGCCGGTGTACGCGATGTAA‑TATAGATCATG‑‑‑ATCATCTTGCGAt                                                 >  1:2016377/1‑76 (MQ=255)                       tttGAAGTAGCGAGGGTCGACCATATTCATGCTGCCGGTGTACGCGATGTAA‑TATCGATCATG‑‑‑ATCATCTTGCGAt                                                 >  1:1284874/1‑76 (MQ=255)                         tGAAGTAGCGAGGGTCGACCATATTCATGCTGCCGGTGTACGCGATGTAA‑TATCGATCATG‑‑‑ATCATCTTGCGATGt                                               >  1:1130895/1‑76 (MQ=255)                          gAAGTAGCGAGGGTCGACCATATTCATGCTGCCGGTGTACGCGATGTAA‑TATCGATCATG‑‑‑ATCATCTTGCGATGtt                                              <  1:1496083/76‑1 (MQ=255)                          gAAGTAGCGAGGGTCGACCATATTCATGCTGCCGGTGTACGCGATGTAA‑TATCGATCATG‑‑‑ATCATCTTGCGATGtt                                              >  1:1177704/1‑76 (MQ=255)                             gTAGCGAGGGTCGACCATATTCATGCTGCCGGTGTACGCGATGTAA‑TATCGATCATG‑‑‑ATCATCTTGCGATGTTGgc                                           <  1:874306/76‑1 (MQ=255)                             gTAGCGAGGGTCGACCATATTCATGCTGCCGGTGTACGCGATGTAA‑TATCGATCATG‑‑‑ATCATCTTGCGATGTTGgc                                           >  1:769455/1‑76 (MQ=255)                               aGCGAGGGTCGACCNTATTCATGCTGCCGGTGTACGCGATGTGA‑TATCGATCATG‑‑‑ATCATCTTCCGATGTTGgcgc                                         <  1:1756015/76‑1 (MQ=255)                               aGCGAGGGTCGACCATATTCATGCTGCCGGTGTACGCGATGTAA‑TATCGATCATG‑‑‑ATCATTTTGCGATGTTGgcgc                                         >  2:206588/1‑76 (MQ=255)                               aGCCAGGGGCGACCATATTCATGCTGCCGGTGTACGCGATGTAATTATCGATCATG‑‑‑ATCATCTTGCGATGTTGgcg                                          >  1:1901071/1‑76 (MQ=255)                                 cGAGGGTCGACCATATTCATGCTGCCGGTGTACGCGATGTAA‑TATCGATCATG‑‑‑ATCATCTTGCGATGTTGGCGCAg                                       >  2:730548/1‑76 (MQ=255)                                   aGGGTCGACCATATTCATGCTGCCGGTGTACGCGATGTAA‑TATCGATCATG‑‑‑ATCATCTTGCGATGTTGGCGCAGGt                                     <  1:1251697/76‑1 (MQ=255)                                     ggTCGACCATATTCATGCTGCCGGTGTACGCGATGTAA‑TATCGATCATG‑‑‑ATCATCTTGCGATGTTGGCGCAGGTcc                                   <  2:881446/76‑1 (MQ=255)                                       tCGACCATATTCATGCTGCCGGTGTACGCGATGTAA‑TATCGATCATG‑‑‑ATCATCTTGCGATGTTGGCGCAGGTCCAt                                 <  1:1380399/76‑1 (MQ=255)                                        cGACCATATTCATGCTGCCGGTGTACGCGATGTAA‑TATCGATCATG‑‑‑ATCATCTTGCGATGTTGGCGCAGGTCCATa                                >  1:860590/1‑76 (MQ=255)                                        cGACCATATTAATGCTGCCGGTGTACGCGATGTAA‑TATCGATCATGATCATCATCTTGCGAGGTTGGCGCAGGTcc                                   >  1:168394/1‑76 (MQ=255)                                         gACCATATTCATGCTGCCGGTGTACGCGATGTAA‑TATCGATCATG‑‑‑ATCATCTTGCGATGTTGGCGCAGGTCCATAc                               <  1:922617/76‑1 (MQ=255)                                              tatTCATGCTGCCGGTGTACGCGATGTAA‑TATCGATCATG‑‑‑ATCATCTTGCGATGTTGGCGCAGGTCCATACGGCGt                          >  1:984155/1‑76 (MQ=255)                                              tatTCATGCTGCCGGTGTACGCGATGTAA‑TATCGATCATG‑‑‑ATCATCTTGCGATGTTGGCGCAGGTCCATACGGCGt                          >  2:1347039/1‑76 (MQ=255)                                               atTCATGCTGCCGGTGTACGCGATGTAA‑TATCGATCATG‑‑‑ATCATCTTGCGATGTTGGCGCAGGTCCATACGGCTTa                         <  2:676026/76‑1 (MQ=255)                                                 tCATGCTGCCGGTGTACGCGATGTAA‑TATCGATCATG‑‑‑ATCATCTTGCGATGTTGGCGCAGGTCCATACGGCGTaaa                       >  2:1230033/1‑76 (MQ=255)                                                    tgctgcCGGTGTACGCGATGTAA‑TATCGATCATG‑‑‑ATCATCTTGCGATGTTGGCGCAGGTCCATACGGCGTAAAAac                    <  2:1474166/76‑1 (MQ=255)                                                       tgcCGGTGTACGCGATGTAA‑TATCGATCATG‑‑‑ATCATCTTGCGATGTTGGCGCAGGTTCATACGGCGTAAAAacaca                 <  1:1725693/76‑1 (MQ=255)                                                          cGGTGTACGCGATGTAA‑TATCGATCATG‑‑‑ATCATCTTGCGATGTTGGCGCAGGTCCATACGGCGTAAAAACACACGc              <  1:1505449/76‑1 (MQ=255)                                                          cGGTGTACGCGATGTAA‑TATCGATCATG‑‑‑ATCATCTTGCGATGTTGGCGCAGGTCCATACGGCGTAAAAACACACGc              <  1:976140/76‑1 (MQ=255)                                                            gtggACGCGGTGTAA‑TATCGATCATG‑‑‑ATCATCTTGCGATGTTGGCGCAGGTCCATACGGCGTAAAAACACACGCAt            <  2:129727/76‑1 (MQ=255)                                                             tgtACGCGATGTAA‑TATCGATCATG‑‑‑ATCATCTTGCGATGTTGGCGCAGGTCCATACGGCGTAAAAACACACGCATc           <  1:1022246/76‑1 (MQ=255)                                                                aCGCGATGTAA‑TATCGATCATG‑‑‑ATCATCTTGCGATGTTGGCGCAGGTCCATACGGCGTAAAAACACACGCATCAGa        >  1:1109082/1‑76 (MQ=255)                                                                  gcgATGTAA‑TATCGATCATG‑‑‑ATCATCTTGCGATGTTGGCGCAGGTCCATACGGCGTAAAAACACACGCATCAGAtt      >  1:1631432/1‑76 (MQ=255)                                                                   cgATGTAA‑TATCGATCATG‑‑‑ATCATCTTGCGATGTTGGCGCAGGTCCATACGGCGTAAAAACACACGCATCAGATTg     >  2:1343890/1‑76 (MQ=255)                                                                   cgATGTAA‑TATCGATCATG‑‑‑ATCATCTTGCGATGTTGGCGCAGGTCCATACGGCGTAAAAACACACGCATCAGATTg     <  1:1392012/76‑1 (MQ=255)                                                                   cgATGTAA‑NATCGATCATG‑‑‑ATCATCTTGCGATGGTAGCGCAGGTCCATACGGCGTAAAAACACACGCATCAGATTg     <  2:345078/76‑1 (MQ=255)                                                                     aTGTAA‑TATCGATCATG‑‑‑ATCATCTTGCGATGTTGGCGCAGGTCCATACGGCGTAAAAACACACGCATCAGATTGAc   <  1:914232/76‑1 (MQ=255)                                                                      tGTAA‑TATCGATCATG‑‑‑ATCATCTTGCGATGTTGGCGCAGGTCCATACGGCGTAAAAACACACGCATCAGATTGAcc  >  1:1376404/1‑76 (MQ=255)                                                                           |                                                                           ATTGCCCTACGCCCGCATCTTGTTTGAAGTAGCGAGGGTCGACCATATTCATGCTGCCGGTGTACGCGATGTAATTATCGATCATG‑‑‑ATCATCTTGCGATGTTGGCGCAGGTCCATACGGCGTAAAAACACACGCATCAGATTGACC  >  REL606/1305471‑1305616

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