Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	1,305,909	G→A	Q111* (CAG→TAG)	cls ←	cardiolipin synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	1,305,909	0	G	A	98.0%	174.8 / ‑3.2	49	Q111* (CAG→TAG)	cls	cardiolipin synthetase
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (23/25);  total (24/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 = 3.96e-01

ATTAACGCCTGCATCACATCATCTGACTCGGTCATCAGTTGTAGCTGATTCCCTTTGACCCCAGCGATCCCCTGACGACGCTCGCAAAGCTTGAATAATGGCGCAGCGACACTGCTATTTTCTTCGGCGAAGATATGCTTACAGGC  >  REL606/1305836‑1305981                                                                          |                                                                         aTTAACGCCTGCATCACATCATCTGACTCGGTCATCAGTTGTAGCTGATTCCCTTTGACCCCAGCGATCCCCTaac                                                                        >  1:664779/1‑76 (MQ=255)      cGCCTGCATCACATCATCTGACTCGGTCATCAGTTGTAGCTGATTCCCTTTGACCCCAGCGATCCCCTAACGACGc                                                                   >  1:1379278/1‑76 (MQ=255)       gCCTGCATCACATCATCTGACTCGGTCATCAGTTGTAGCTGATTCCCTTTGACCCCAGCGATCCCCTAACGACGCt                                                                  <  2:154/76‑1 (MQ=255)        ccTGCATCACATCATCTGACTCGGTCATCAGTTGTAGCTGATTCCCTTTGACCCCAGCGATCCCCTAACGACGCTc                                                                 <  2:141081/76‑1 (MQ=255)           gCATCACATCATCTGACTCGGTCATCANTTGTAGCTGATTCCCTTTGACCCCAGCGATCCCCTAACGACGCTCGCa                                                              >  1:1606458/1‑76 (MQ=255)           gCATCACATCATCTGACTCGGTCATCAGTTGTAGCTGATTCCCTTTGACCCCAGCGATCCCCTAACGACGCTCGCa                                                              >  2:875757/1‑76 (MQ=255)            catcacatcaTCTGACTCGGTCATCAGTTGTAGCTGATTCCCTTTGACCCCAGCGATCCCCTAACGACGCTCGCaa                                                             >  2:194048/1‑76 (MQ=255)              tcacatcaTCTGNCTCGGTCATCAGTTGTAGCTGATTCCCTTTGACCCCAGCGATCCCCTAACGACGCTCGCAAAg                                                           >  1:19898/1‑76 (MQ=255)              tcacatcaTCTGACTCGGTCATCAGTTGTAGCTGATTCCCTTTGACCCCAGCGATCCCCTAACGACGCTCGCAAAg                                                           <  1:826651/76‑1 (MQ=255)               cacatcaTCTGACTCGGTCATCAGTTGTAGCTGATTCCCTTTGACCCCAGCGATCCCCTAACGACGCTCGCAAAGc                                                          <  1:225980/76‑1 (MQ=255)                 catcctATGACTCGGTCATCAGTTGTAGCTGATTCCCTTTGACCCCAGGGATCCCCTAACGACGCTCGCAAAGCtt                                                        <  1:1589032/76‑1 (MQ=255)                   tcatcTGACTCGGTCATTAGTTGTAGCTGATTCCCTTTGACCCCAGCGATCCCCTAACGACGCTCGCAAAGCTTGa                                                      <  2:332898/76‑1 (MQ=255)                    catcTGACTCGGTCATCAGTTGTAGCTGATTCCCTTTGACCCCAGCGATCCCCTAACGACGCTCGCAAAGCTTGaa                                                     >  1:1247191/1‑76 (MQ=255)                    catcTGACTCGGTCATCAGTTGTAGCTGATTCCCTTTGACCCCAGCGATCCCCTAACGACGCTCGCAAAGCTTGaa                                                     <  2:682899/76‑1 (MQ=255)                          aCTCGGTCATCAGTTGTAGCTGATTCCCTTTGACCCCAGCGATCCCCTAACGACGCTCGCAAAGCTTGAATAATgg                                               <  1:548712/76‑1 (MQ=255)                           cTCGGTCATCAGTTGTAGCTGATTCCCTTTGACCCCAGCGATCCCCTAACGACGCTCGCAAAGCTTGAATAATGgc                                              >  1:786147/1‑76 (MQ=255)                             cGGTCATCAGTTGTAGCTGATTCCCTTTGACCCCAGCGATCCCCTAACGACGCTCGCAAAGCTTGAATAATGgcgc                                            >  2:1569827/1‑76 (MQ=255)                              ggTCATCAGTTGTAGCTGATTCCCTTTGACCCCAGCGATCCCCTAACGACGCTCGCAAAGCTTGAATAATGGCGCa                                           <  2:160390/76‑1 (MQ=255)                              ggTCATCAGTTGTAGCTGATTCCCTTTGACCCCAGCGATCCCCTAACGACGCTCGCAAAGCTTGAATAATGGCGCa                                           >  1:476143/1‑76 (MQ=255)                              ggTCATCAGTTGTAGCTGATTCCCTTTGACCCCAGCGATCCCCTAACGACGCTCGCAAAGCTTGAATAATGGCGCa                                           >  1:1055096/1‑76 (MQ=255)                               gTCATCAGTTGTAGCTGATTCCCTTTGACCCCAGCGATCCCCTAACGACGCTCGCAAAGCTTGAATAATGGCGCAg                                          >  2:1231766/1‑76 (MQ=255)                                tcatcaGTTGTAGCTGATTCCCTTTGACCCCAGCGATCCCCTAACGACGCTCGCAAAGCTTGAATAATGGCGCAGc                                         >  2:751049/1‑76 (MQ=255)                                    caGTTGTAGCTGATTCCCTTTGACCCCAGCGATCCCCTAACGACGCTCGCAAAGCTTGAATAATGGCGCAGCGaca                                     >  1:486354/1‑76 (MQ=255)                                      gTTGTAGCTGATTCCCTTTGACCCCAGCGATCCCCTAACGACGCTCGCAAAGCTTGAATAATGGCGCAGCGACACt                                   >  1:560474/1‑76 (MQ=255)                                      gTTGTAGCTGATTCCCTTTGACCCCAGCGATCCCCTAACGACGCTCGCAAAGCTTGAATAATGGCGCAGCGACACt                                   <  1:934371/76‑1 (MQ=255)                                       ttGTAGCTGATTCCNTTTGACCCCAGCGGTCCCCNAACGATGCNCGCAAAGCTTGAATAATGGCGCAGCGACACTg                                  <  2:1345058/76‑1 (MQ=255)                                       ttGTAGCTGATTCCCTTTGACCCCAGCGATCCCCTAACGACGCTCGCAAAGCTTGAATAATGGCGCAGCGACACTg                                  <  1:598766/76‑1 (MQ=255)                                          tAGCTGATTCCCTTTGACCCCAGCGATCCCCTAACGACGCTCGCAAAGCTTGAATAATGGCGCAGCGACACTGCTa                               <  1:693706/76‑1 (MQ=255)                                              tGATTCCCTTTGACCCCAGCGATCCCCTAACGACGCTCGCAAAGCTTGAATAATGGCGCAGCGACACTGCTAtttt                           <  1:1061438/76‑1 (MQ=255)                                               gATTCCCTTTGACCCCAGCGATCCCCTAACGACGCTCGCAAAGCTTGAATAATGGCGCAGCGACACTGCTATTttc                          >  2:401756/1‑76 (MQ=255)                                               gATTCCCTTTGACCCCAGCGATCCCCTAACGACGCTCGCAAAGCTTGAATAATGGCGCAGCGACACTGCTATTttc                          >  2:401754/1‑76 (MQ=255)                                                aTTCCCTTTGACCCCAGCGATCCCCTAACGACGCTCGCAAAGCTTGAATAATGGCGCAGCGACACTGCTATTttct                         >  2:1140417/1‑76 (MQ=255)                                                aTTCCCTTTGACCCCAGCGATCCCCTAACGACGCTCGCAAAGCTTGAATAATGGCGCAGCGACACTGCTATTttct                         >  2:1042664/1‑76 (MQ=255)                                                  tCCCTTTGACCCCAGCGATCCCCTAACGACGCTCGCAAAGCTTGAATAATGGCGCAGCGACACTGCTATTttcttc                       <  1:243024/76‑1 (MQ=255)                                                   cccTTTGACCCCAGCGATCCCCTAACGACGCTCGCAAAGCTTGAATAATGGCGCAGCGACACTGCTATTTTCTTCg                      <  2:973616/76‑1 (MQ=255)                                                      tttGACCCCAGCGATCCCCTGACGACGCTCGCAAAGCTTGAATAATGGCGCAGCGACACTGCTATTTTCTTCGGCg                   >  2:1048673/1‑76 (MQ=255)                                                        tGACCCCAGCGATCCCCTAACGACGCTCGCAAAGCTTGAATAATGGCGCAGCGACACTGCTATTTTCTTCGGCGaa                 <  1:396287/76‑1 (MQ=255)                                                         gACCCCAGCGATCCCCTAACGACGCTCGCAAAGCTTGAATAATGGCGCAGCGACACTGCTATTTTCTTCGGCGAAg                <  2:932588/76‑1 (MQ=255)                                                            cccAGCGATCCCCTAACGACGCTCGCAAAGCTTGAATAATGGCGCAGCNACACTGCTATTTTCTTCGGCGAAGata             <  1:1606323/76‑1 (MQ=255)                                                            cccAGCGATCCCCTAACGACGCTCGCAAAGCTTGAATAATGGCGCAGCGACACTGCTATTTTCTTCGGCGAAGata             <  1:850274/76‑1 (MQ=255)                                                              cAGCGATCCCCTAACGACGCTCGCAAAGCTTGAATAATGGCGCAGCGACACTGCTATTTTCTTCGGCGAAGATATg           <  1:332156/76‑1 (MQ=255)                                                              cAGCGATCCCCTAACGACGCTCGCAAAGCTTGAATAATGGCGCAGCGACACTGCTATTTTCTTCGGCGAAGATATg           <  1:715398/76‑1 (MQ=255)                                                              cAGCGATCCCCTAACGACGCTCGCAAAGCTTGAATAATGGCGCAGCGACACTGCTATTTTCTTCGGCGAAGATATg           >  2:985559/1‑76 (MQ=255)                                                                 cGATCCCCTAACGACGCTCGCAAAGCTTGAATAATGGCGCAGCGACACTGCTATTTTCTTCGGCGAAGATATGCtt        <  2:996501/76‑1 (MQ=255)                                                                  gATCCCCTAACGACGCTCGCTANNCTTGAATAATNNAGCAGCGACACTGCGATTTTCTTCGGCGAAGATATGCTTa       <  1:1450433/76‑1 (MQ=255)                                                                  gATCCCCTAACGACGCTCGCAAAGCTTGAATAATGGCGCAGCGACACTGCTATTTTCTTCGGCGAAGATATGCTTa       <  2:666514/76‑1 (MQ=255)                                                                  gATCCCCTAACGACGCTCGCAAAGCTTGAATAATGGCGCAGCGACACTGCTATTTTCTTCGGCGAAGATATGCTTa       >  2:1442815/1‑76 (MQ=255)                                                                     ccccTAACGACGCTCGCAAAGCTTGAATAATGGCGCAGCGACACTGCTATTTTGTTCGGCGAAGATATGCTTACAg    >  2:56429/1‑76 (MQ=255)                                                                     ccccTAACGACGCTCGCAAAGCTTGAATAATGGCGCAGCGACACTGCTATTTTCTTCGGCGAAGATATGCTTACAg    >  2:704809/1‑76 (MQ=255)                                                                       ccTAACGACGCTCGCAAAGCTTGAATAATGGCGCAGCGACACTGCTATTTTCTTCGGCGAAGATATGCTTACAGGc  >  1:773613/1‑76 (MQ=255)                                                                          |                                                                         ATTAACGCCTGCATCACATCATCTGACTCGGTCATCAGTTGTAGCTGATTCCCTTTGACCCCAGCGATCCCCTGACGACGCTCGCAAAGCTTGAATAATGGCGCAGCGACACTGCTATTTTCTTCGGCGAAGATATGCTTACAGGC  >  REL606/1305836‑1305981

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