Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	1,305,551	A→G	I230T (ATC→ACC)	cls ←	cardiolipin synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	1,305,551	0	A	G	98.3%	234.4 / ‑3.7	59	I230T (ATC→ACC)	cls	cardiolipin synthetase
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (24/34);  total (25/34)
Fisher's exact test for biased strand distribution p-value = 4.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.10e-01

ACGCCCGCATCTTGTTTGAAGTAGCGAGGGTCGACCATATTCATGCTGCCGGTGTACGCGATGTAATTATCGATCATGATCATCTTGCGATGTTGGCGCAGGTCCATACGGCGTAAAAACACACGCATCAGATTGACCTTTAA  >  REL606/1305479‑1305621                                                                         |                                                                       aCGCCCGCATCTTGTTTGAAGTAGCGAGGGTCGACCATATTCATGCTGCCGGTGTACGCGATGTAATTATCggtc                                                                      <  1:1076235/75‑1 (MQ=255)   gCCCGCATCTTGTTTGAAGTAGCGAGGGTCGACCATATTCATGCTGCCGGTGTACGCGATGTAATTATCggtcat                                                                    <  2:65130/75‑1 (MQ=255)   gCCCGCATCTTGTTTGAAGTAGCGAGGGTCGACCATATTCATGCTGCCGGTGTACGCGATGTAATTATCggtcat                                                                    >  2:1632077/1‑75 (MQ=255)    cccGCATCTTGTTTGAAGTAGCGAGGTTCGACCATATTCATGCTGCCGGTGTACGCGATGTAATTATCggtcatg                                                                   <  1:152682/75‑1 (MQ=255)    cccGCATCTTGTTTGAAGTAGCGAGGGTCGACCATATTCATGCTGCCGGTGTACGCGATGTAATTATCggtcatg                                                                   >  1:92523/1‑75 (MQ=255)    cccGCATCTTGTTTGAAGTAGCGAGGGTCGACCATATTCATGCTGCCGGTGTACGCGATGTAATTATCggtcatg                                                                   <  2:510395/75‑1 (MQ=255)     ccGCATCTTGTTTGAAGTAGCGAGGGTCGACCATATTCATGCTGCCGGTGTACGCGATGTAATTATCggtcatga                                                                  >  1:930799/1‑75 (MQ=255)     ccGCATCTTGTTTGAAGTAGCGAGGGTCGACCATATTCATGCTGCCGGTGTACGCGATATAATTATCggtcatga                                                                  <  1:1305108/75‑1 (MQ=255)      cGCATCTTGTTTGAAGTAGCGAGGGTCGACCATATTCATGCTGCCGGTGTACGCGATGTAATTATCggtcatgat                                                                 <  1:759288/75‑1 (MQ=255)       gCATCTTGTTTGAAGTAGCGAGGGTCGACCATATTCATGCTGCCGGTGTACGCGATGTAATTATCggtcatgatc                                                                >  2:338383/1‑75 (MQ=255)       gCATCTTGTTTGAAGTAGCGAGGGTCGACCATATTCATGCTGCCGGTGTACGCGATGTAATTATCggtcatgatc                                                                >  2:1003108/1‑75 (MQ=255)        cATCTTGTTTGAAGTAGCGAGGGTCGACCATATTCATGCTGCCGGTGTACGCGATGTAATTATCggtcatgatca                                                               <  2:759254/75‑1 (MQ=255)        cATCTTGTTTGAAGTAGCGAGGGTCGACCATATTCATGCTGCCGGTGTACGCGATGTAATTATCggtcatgatca                                                               <  2:1295809/75‑1 (MQ=255)          tCTTGTTTGAAGTAGCGAGGGTCGACCATATTCATGCTGCCGGTGTACGCGATGTAATTATCGGTCATGatcatc                                                             <  2:1029816/75‑1 (MQ=255)               tttGAAGTAGCGAGGGTCGACCATATTCATGCTGCCGGTGTACGCGATGTAATTATCGGTCATGATCATCTTGCg                                                        >  1:508147/1‑75 (MQ=255)                  gAAGTAGCGAGGGTCGACCATATTCATGCTGCCGGTGTACGCGATGTAATTATCGGTCATGATCATCTTGCGATg                                                     >  1:1208628/1‑75 (MQ=255)                   aaGTAGCGAGGGTCGACCATATTCATGCTGCCGGTGTACGCGATGTAATTATCGGTCATGATCATCTTGCGATGt                                                    >  1:771463/1‑75 (MQ=255)                     gTAGCGAGGGTCGACCATATTCATGCTGCCGGTGTACGCGATGTAATTATCGGTCATGATCATCTTGCGATGTTg                                                  >  1:1686851/1‑75 (MQ=255)                     gTAGCGAGGGTCGACCATATTCATGCTGCCGGTGTACGCGATGTAATTATCGGTCATGATCATCTTGCGATGTTg                                                  >  1:2030297/1‑75 (MQ=255)                       aGCGAGGGTCGACCATATTCATGCTGCCGGTGTACGCGATGTAATTATCGGTCATGATCATCTTGCGATGTTGgc                                                >  2:542054/1‑75 (MQ=255)                       aGCGAGGGTCGACCATATTCATGCTGCCGGTGTACGCGATGTAATTATCGGTCATGATCATCTTGCGATGTTGgc                                                >  2:661010/1‑75 (MQ=255)                        gCGAGGGTCGACCATATTCATGCTGCCGGTGTACGCGATGTAATTATCGGTCATGATCATCTTGCGATGTTGgcg                                               <  2:312836/75‑1 (MQ=255)                        gCGAGGGTCGACCATATTCATGCTGCCGGTGTACGCGATGTAATTATCGGTCATGATCATCTTGCGATGTTGgcg                                               >  2:769957/1‑75 (MQ=255)                         cGAGGGTCGACCATATTCATGCTGCCGGTGTACGCGATGTAATTATCGGTCATGATCATCTTGCGATGTTGgcgc                                              <  2:1326528/75‑1 (MQ=255)                          gAGGGTCGACCATATTCATGCTGCCGGTGTACGCGATGTAATTATCGGTCATGATCATCTTGCGATGTTGGCGCa                                             <  1:920631/75‑1 (MQ=255)                           aGGGTCGACCATATTCATGCTGCCGGTGTACGCGATGTAATTATCGGTCATGATCATCTTGCGATGTTGGCGCAg                                            >  2:167328/1‑75 (MQ=255)                            gggTCGACCATATTCATGCTGCCGGTGTACGCGATGTAATTATCGGTCATGATCATCTTGCGATGTTGGCGCAgg                                           <  1:1390486/75‑1 (MQ=255)                              gTCGACCATATTCATGCTGCCGGTGTACGCGATGTAATTATCGGTCATGATCATCTTGCGATGTTGGCGCAGGTc                                         <  1:1918330/75‑1 (MQ=255)                                cGACCATATTCATGCTGCCGGTGTACGCGATGTAATTATCGGTCATGATCATCTTGCGATGTTGGCGCAGGTCCa                                       <  2:1539218/75‑1 (MQ=255)                                cGACCATATTCATGCTGCCGGTGTACGCGATGTAATTATCGGTCATGATCATCTTGCGATGTTGGCGCAGGTCCa                                       <  1:526796/75‑1 (MQ=255)                                    cATATTCATGCTGCCGGTGTACGCGATGTAATTATCGGTCATGATCATCTTGCGATGTTGGCGCAGGTCCATACg                                   <  2:1130053/75‑1 (MQ=255)                                      tatTCATGCTGCCGGTGTACGCGATGTAATTATCGGTCATGATCATCTTGCGATGTTGGCGCAGGTCCATACGGc                                 >  1:1934833/1‑75 (MQ=255)                                      tatTCATGCTGCCGGTGTACGCGATGTAATTATCGGTCATGATCATCTTGCGATGTTGGCGCAGGTCCATACGGc                                 <  2:573101/75‑1 (MQ=255)                                        ttCTTGCTGCCGGTGTACGCGATGTAAATATCAGTCATGATCATCTTGCGATGTTGTCGCATGTCCATCCGGCGt                               <  2:1459362/75‑1 (MQ=255)                                        ttCATGCTGCCGGTGTACGCGATGTAATTATCGGTCATGATCATCTTGCGATGTTGGCGCAGGTCCATACGGCGt                               <  2:1459353/75‑1 (MQ=255)                                         tCATGCTGCCGGTGTACGCGATGTAATTATCGGTCATGATCATCTTGCGATGTTGGCGCAGGTCCATACGGCGTa                              >  1:1541912/1‑75 (MQ=255)                                           aTGCTGCCGGTGTACGCGATGTAATTATCGGTCATGATCATCTTGCGATGTTGGCGCAGGTCCATACGGCGTaaa                            <  2:1330495/75‑1 (MQ=255)                                           aTGCTGCCGGTGTACGCGATGTAATTATCGGTCATGATCATCTTGCGATGTTGGCGCAGGTCCATACGGCGTaaa                            <  1:1661657/75‑1 (MQ=255)                                            tgctgcCGGTGTACGCGATGTAATTATCGGTCATGATCATCTTGCGATGTTGGCGCAGGTCCATACGGCGTaaaa                           >  1:725029/1‑75 (MQ=255)                                              ctgcCGGTGTACGCGATGTAATTATCGGTCATGATCATCTTGCGATGTTGGCGCAGGTCCATACGGCGTAAAAac                         <  1:2047967/75‑1 (MQ=255)                                               tgcCGGTGTACGCGATGTAATTATCGGTCATGATCATCTTGCGATGTTGGCGCAGGTCCATACGGCGTAAAAaca                        >  2:1055921/1‑75 (MQ=255)                                                gcCGGTGTACGCGATGTAATTATCGGTCATGATCATCTTGCGATGTTGGCGCAGGTCCATACGGCGTAAAAacac                       <  1:1544215/75‑1 (MQ=255)                                                gcCGGTGTACGCGATGTAATTATCGGTCATGATCATCTTGCGATGTTGGCGCAGGTCCATACGGCGTAAAAacac                       <  1:1378842/75‑1 (MQ=255)                                                gcCGGTGTACGCGATGTAATTATCGGTCATGATCATCTTGCGATGTTGGCGCAGGTCCATACGGCGTAAAAacac                       >  2:1477373/1‑75 (MQ=255)                                                  cGGTGTACGCGATGTAATTATCGGTCATGATCATCTTGCGATGTTGGCGCAGGTCCATACGGCGTAAAAacacac                     <  2:596851/75‑1 (MQ=255)                                                  cGGTGTACGCGATGTAATTATCGGTCATGATCATCTTGCGATGTTGGCGCAGGTCCATACGGCGTAAAAacacac                     <  2:1536232/75‑1 (MQ=255)                                                  cGGTGTACGCGATGTAATTATCGGTCATGATCATCTTGCGATGTTGGCGCAGGTCCATACGGCGTAAAAacacac                     <  1:661010/75‑1 (MQ=255)                                                     tgtACGCGATGTAATTATCGGTCATGATCATCTTGCGATGTTGGCGCAGGTCCATACGGCGTAAAAACACACGCa                  <  2:376241/75‑1 (MQ=255)                                                     tgtACGCGATGTAATTATCGGTCATGATCATCTTGCGATGTTGGCGCAGGTCCATACGGCGTAAAAACACACGCa                  <  1:1632077/75‑1 (MQ=255)                                                     tgtACGCGATGTAATTATCGGTCATGATCATCTTGCGATGTTGGCGCAGGTCCATACGGCGTAAAAACACACGCa                  >  2:753148/1‑75 (MQ=255)                                                      gtACGCGATGTAATTATCGGTCATGATCATCTTGCGATGTTGGCGCAGGTCCATACGGCGTAAAAACACACGCAt                 >  2:1183153/1‑75 (MQ=255)                                                      gtACGCGATGTAATTATCGGTCATGATCATCTTGCGATGTTGGCGCAGGTCCATACGGCGTAAAAACACACGCAt                 >  1:1999971/1‑75 (MQ=255)                                                      gtACGCGATGTAATTATCGGTCATGATCATCTTGCGATGTTGGCGCAGGTCCATACGGCGTAAAAACACACGCAt                 <  1:219927/75‑1 (MQ=255)                                                      gtACGCGATGTAATTATCGATCATGATCATCTTGCGATGTTGGCGCAGGTCCATACGGCGTAAAAACACACGCAt                 >  1:454508/1‑75 (MQ=255)                                                       tACGCGATGTAATTATCGGTCATGATCATCTTGCGATGTTGGCGCAGGTCCATACGGCGTAAAAACACACGCATc                <  1:828363/75‑1 (MQ=255)                                                       tACGCGATGTAATTATCGGTCATGATCATCTTGCGATGTTGGCGCAGGTCCATACGGCGTAAAAACACACGCATc                <  1:21742/75‑1 (MQ=255)                                                       tACGCGATGTAATTATCGGTCATGATCATCTTGCGATGTTGGCGCAGGTCCATACGGCGTAAAAACACACGCATc                >  1:1317049/1‑75 (MQ=255)                                                        aCGCGATGTAATTATCGGTCATGATCATCTTGCGATGTTGGCGCAGGTCCATACGGCGTAAAAACACACGCATCa               <  1:1320031/75‑1 (MQ=255)                                                        aCGCGATGTAATTATCGGTCATGATCATCTTGCGATGTTGGCGCAGGTCCATACGGCGTAAAAACACACGCATCa               <  1:1465667/75‑1 (MQ=255)                                                        aCGCGATGTAATTATCGGTCATGATCATCTTGCGATGTTGGCGCAGGTCCATACGGCGTAAAAACACACGCATCa               <  2:930799/75‑1 (MQ=255)                                                        aCGCGATGTAATTATCGGTCATGATCATCTTGCGATGTTGGCGCAGGTCCATACGGCGTAAAAACACACGCATCa               >  2:1893667/1‑75 (MQ=255)                                                         cgcgATGTAATTATCGGTCATGATCATCTTGCGATGTTGGCGCAGGTCCATACGGCGTAAAAACACACGCATCAg              >  1:796640/1‑75 (MQ=255)                                                          gcgATGTAATTATCGGTCATGATCATCTTGCGATGTTGGCGCAGGTCCATACGGCGTAAAAACACACGCACCAGa             >  1:950072/1‑75 (MQ=255)                                                           cgATGTAATTATCGGTCATGATCATCTTGCGATGTTGGCGCAGGTCCATACGGCGTAAAAACACACGCATCAGAt            <  1:408339/75‑1 (MQ=255)                                                            gATGTAATTATCGGTCATGATCATCTTGCGATGTTGGCGCAGGTCCATACGGCGTAAAAACACACGCATCAGAtt           <  1:1829701/75‑1 (MQ=255)                                                            gATGTAATTATCGGTCATGATCATCTTGCGATGTTGGCGCAGGTCCATACGGCGTAAAAACACACGCATCAGAtt           <  1:361191/75‑1 (MQ=255)                                                                tAATTATCGGTCATGATCATCTTGCGATGTTGGCGCAGGTCCATACGGCGTAAAAACACACGCATCAGATTGAcc       <  1:374157/75‑1 (MQ=255)                                                                  aTTATCGGTCATGATCATCTTGCGATGTTGGCGCAGGTCCATACGGCGTAAAAACACACGCATCAGATTGACCtt     <  2:855586/75‑1 (MQ=255)                                                                   ttATCGGTCATGATCATCTTGCGATGTTGGCGCAGGTCCATACGGCGTAAAAACACACGCATCAGATTGACCttt    >  1:1500802/1‑75 (MQ=255)                                                                   ttATCGGTCATGATCATCTTGCGATGTTGGCGCAGGTCCATACGGCGTAAAAACACACGCATCAGATTGACCttt    >  1:138646/1‑75 (MQ=255)                                                                   ttATCGGTCATGATCATCTTGCGATGTTGGCGCAGGTCCATACGGCGTAAAAACACACGCATCAGATTGACCttt    >  2:720208/1‑75 (MQ=255)                                                                     aTCGGTCATGATCATCTTGCGATGTTGGCGCAGGTCCATACGGCGTAAAAACACACGCATCAGATTGACCTTTaa  <  2:184098/75‑1 (MQ=255)                                                                     aTCGGTCATGATCATCTTGCGATGTTGGCGCAGGTCCATACGGCGTAAAAACACACGCATCAGATTGACCTTTaa  <  2:638725/75‑1 (MQ=255)                                                                     aTCGATCTTGATCATCTTGCGATGTTGGCGCAGGTCCATACGGCGTAAAAACACACGCATCAGATTGACCTTTaa  <  1:148388/75‑1 (MQ=255)                                                                         |                                                                       ACGCCCGCATCTTGTTTGAAGTAGCGAGGGTCGACCATATTCATGCTGCCGGTGTACGCGATGTAATTATCGATCATGATCATCTTGCGATGTTGGCGCAGGTCCATACGGCGTAAAAACACACGCATCAGATTGACCTTTAA  >  REL606/1305479‑1305621

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