Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	1,305,551	A→T	I230N (ATC→AAC)	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	T	97.4%	152.2 / ‑3.4	39	I230N (ATC→AAC)	cls	cardiolipin synthetase
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base T (20/18);  total (21/18)
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.88e-01

TACGCCCGCATCTTGTTTGAAGTAGCGAGGGTCGACCATATTCATGCTGCCGGTGTACGCGATGTAATTATCGATCATGATCATCTTGCGATGTTGGCGCAGGTCCATACGGCGTAAAAACACACGCATCAGATTGACCTTTA  >  REL606/1305478‑1305620                                                                          |                                                                      tACGCCCGCATCTTGTTTGAAGTAGCGAGGGTCGACCATATTCATGCTGCCGGTGTACGCGATGTAATTATCgttc                                                                     <  1:1008219/76‑1 (MQ=255)  aCGCCCGCATCTTGTTTGAAGTAGCGAGGGTCGACCATATTCATGCTGCCGGTGTACGCGATGTAATTATCgttca                                                                    <  1:1231279/76‑1 (MQ=255)     cccGCATCTTGTTTGAAGTAGCGAGGGTCGACCATATTCATGCTGCCGGTGTACGCGATGTAATTATCgttcatga                                                                 <  2:225510/76‑1 (MQ=255)        gCATCTTGTTTGAAGTAGCGAGGGTCGACCATATTCATGCTGCCGGTGTACGCGATGTAATTATCgttcatgatca                                                              <  2:1477903/76‑1 (MQ=255)            cTTGTTTGAAGTAGCGAGGGTCGACCATATTCATGCTGCCGGTGTACGCGATGTAATTATCGTTCATGCTCATCtt                                                          >  1:218045/1‑76 (MQ=255)            cTTGTTTGAAGTAGCGAGGGTCGACCATATTCATGCTGCCGGTGTACGCGATGTAATTATCGTTCATGATCATCtt                                                          >  1:1503084/1‑76 (MQ=255)            cTTGTTTGAAGTAGCGAGGGTCGACCATATTCATGCTGCCGGTGTACGCGATGTAATTATCGTTCATGATCATCtt                                                          <  1:660632/76‑1 (MQ=255)             ttGTTTGAAGTAGCGAGGGTCGACCATATTCATGCTGCCGGTGTACGCGATGTAATTATCGTTCATGATCATCTTg                                                         <  2:1315482/76‑1 (MQ=255)             ttGTTTGAAGTAGCGAGGGTCGACCATATTCATGCTGCCGGTGTACGCGATGTAATTATCGTTCATGATCATCTTg                                                         <  2:327702/76‑1 (MQ=255)                   gAAGTAGCGAGGGTCGACCATATTCATGCTGCCGGTGTACGCGATGTAATTATCGTTCATGATCATCTTGCGATGt                                                   >  2:1021934/1‑76 (MQ=255)                       tAGCGAGGGTCGACCATATTCATGCTGCCGGTGTACGCGATGTAATTATCGTTCATGATCATCTTGCGATGTTGgc                                               <  1:349170/76‑1 (MQ=255)                        aGCGAGGGTCGACCATATTCATGCTGCCGGTGTACGCGATGTAATTATCGTTCATGATCATCTTGCGATGTTGgcg                                              >  1:977988/1‑76 (MQ=255)                         gCGAGGGTCGACCATATTCATGCTGCCGGTGTACGCGATGTAATTATCGTTCATGATCATCTTGCGATGTTGgcgc                                             <  1:439008/76‑1 (MQ=255)                         gCGAGGGTCGACCATATTCATGCTGCCGGTGTACGCGATGTAATTATCGTTCATGATCATCTTGCGATGTTGgcgc                                             >  1:126940/1‑76 (MQ=255)                          cGAGGGTCGACCATATTCATGCTGCCGGTGTACGCGATGTAATTATCGTTCATGATCATCTTGCGATGTTGGCGCa                                            >  1:1240309/1‑76 (MQ=255)                            aGGGTCGACCATATTCATGCTGCCGGTGTACGCGATGTAATTATCGTTCATGATCATCTTGCGATGTTGGCGCAgg                                          <  2:418408/76‑1 (MQ=255)                                 cGACCATATTCATGCTGCCGGTGTACGCGATGTAATTATCGTTCATGATCATCTTGCGATGTTGGCGCAGGTCCAt                                     <  1:899392/76‑1 (MQ=255)                                  gACCATATTCATGCTGCCGGTGTACGCGATGTAATTATCGTTCATGATCATCTTGCGATGTTGGCGCAGGTCCATa                                    >  2:1252253/1‑76 (MQ=255)                                    ccATATTCATGCTGCCGGTGTACGCGATGTAATTATCGTTCATGATCATCTTGNGATGTTGGCGCAGGTCCATACg                                  <  2:221789/76‑1 (MQ=255)                                        atTCATGCTGTCGGTGTACGCGATGTAATTATCGTTCATGATCATCTTGCGATGTTGGCGCAGGTCCATACGGCGt                              <  1:755169/76‑1 (MQ=255)                                         ttCATGCTGCCGGTGTACGCGATGTAATTATCGTTCATGATCATCTTGCGATGTTGGCGCAGGTCCATACGGCGTa                             <  1:503102/76‑1 (MQ=255)                                            aTGCTGCCGGTGTACGCGATGTAATTATCGTTCATGATCATCTTGCGATGTTGGCGCAGGTCCATACGGCGTaaaa                          >  2:1619831/1‑76 (MQ=255)                                             tgctgcCGGTGTACGCGATGTAATTATCGTTCATGATCATCTTGCGATGTTGGCGCAGGTCCATACGGCGTaaaaa                         <  1:1481976/76‑1 (MQ=255)                                                tgcCGGTGTACGCGATGTAATTATCGTTCATGATCATCTTGCGATGTTGGCGCAGGTCCATACGGCGTAAAAacac                      >  2:1285729/1‑76 (MQ=255)                                                 gcCGGTGTACGCGATGTAATTATCGTTCATGATCATCTTGCGATGTTGGCGCAGGTCCATACGGCGTAAAAacaca                     >  2:166293/1‑76 (MQ=255)                                                  ccGGTGTACGCGATGTAATTATCGTTCATGATCATCTTGCGATGTTGGCGCAGGTCCATACGGCGTAAAAacacac                    <  2:154762/76‑1 (MQ=255)                                                  ccGGTGTACGCGATGTAATTATCGTTCATGATCATCTTGCGATGTTGGCGCAGGTCCATACGGCGTAAAAacacac                    >  1:285724/1‑76 (MQ=255)                                                  ccGGTGTACGCGATGTAATTATCGTTCATGATCATCTTGCGATGTTGGCGCAGGTCCATACGGCGTAAAAacacac                    <  1:658030/76‑1 (MQ=255)                                                    ggTGTACGCGATGTAATTATCGTTCATGATCATCTTGCGATGTTGGCGCAGGTCCATACGGCGTAAAAACACACGc                  >  1:623242/1‑76 (MQ=255)                                                    ggTGTACGCGATGTAATTATCGTTCATGATCATCTTGCGATGTTGGCGCAGGTCCATACGGCGTAAAAACACACGc                  >  1:30625/1‑76 (MQ=255)                                                      tgtACGCGATGTAATTATCGTTCATGATCATCTTGCGATGTTGGCGCAGGTCCATACGGCGTAAAAACACACGCAt                >  2:476705/1‑76 (MQ=255)                                                       gtACGCGATGTAATTATCGATCATGATCATCTTGCGATGTTGGCGCAGGTCCATACGGCGTAAAAACACACGCATc               >  1:1463015/1‑76 (MQ=255)                                                        tACGCGATGTAATTATCGTTCATGATCATCTTGCGATGTTGGCGCAGGTCCATACGGCGTAAAAACACACGCATCa              <  2:1180666/76‑1 (MQ=255)                                                            cgATGTAATTATCGTTCATGATCATCTTGCGATGTTGGCGCAGGTCCATACGGCGTAAAAACACACGCATCAGAtt          <  1:321857/76‑1 (MQ=255)                                                             gATGTAATTATCGTTCATGATCATCTTGCGATGTTGGCGCAGGTCCATACGGCGTAAAAACACACGCATCAGATTg         >  1:830383/1‑76 (MQ=255)                                                                gTAATTATCGTTCATGATCATCTTGCGATGTTGGCGCAGGTCCATACGGCGTAAAAACACACGCATCAGATTGAcc      >  1:1414004/1‑76 (MQ=255)                                                                 tAATTATCGTTCATGATCATCTTGCGATGTTGGCGCAGGTCCATACGGCGTAAAAACACACGCATCAGATTGACCt     >  1:777378/1‑76 (MQ=255)                                                                 tAATTATCGTTCATGATCATCTTGCGATGTTGGCGCAGGTCCATACGGCGTAAAAACACACGCATCAGATTGACCt     >  1:1616024/1‑76 (MQ=255)                                                                  aaTTATCGTTCATGATCATCTTGCGATGTTGGCGCAGGTCCATACGGCGTAAAAACACACGCATCAGATTGACCtt    >  1:930137/1‑76 (MQ=255)                                                                   aTTATCGTTCATGATCATCTTGCGATGTTGGCGCAGGTCCATACGGCGTAAAAACACACGCATCAGATTGACCttt   <  2:144766/76‑1 (MQ=255)                                                                   aTTATCGTTCATGATCATCTNNCGATGTTGGCGCAGNTNNATACGGCGTAAAAACACACGCATCAGATTGACCttt   <  1:1188887/76‑1 (MQ=255)                                                                    ttATCGTTCATGATCATCTTGCGATGTTGGCGCAGGTCCATACGGCGTAAAAACACACGCATCAGATTGACCTTTa  >  1:579161/1‑76 (MQ=255)                                                                    ttATCGTTCATGATCATCTTGCGATGTTGGCGCAGGTCCATACGGCGTAAAAACACACGCATCAGATTGACCTTTa  <  2:285724/76‑1 (MQ=255)                                                                          |                                                                      TACGCCCGCATCTTGTTTGAAGTAGCGAGGGTCGACCATATTCATGCTGCCGGTGTACGCGATGTAATTATCGATCATGATCATCTTGCGATGTTGGCGCAGGTCCATACGGCGTAAAAACACACGCATCAGATTGACCTTTA  >  REL606/1305478‑1305620

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