Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	1,306,225	A→T	Y5* (TAT→TAA)	cls ←	cardiolipin synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	1,306,225	0	A	T	98.2%	213.2 / ‑3.6	54	Y5* (TAT→TAA)	cls	cardiolipin synthetase
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base T (27/26);  total (27/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.13e-01

CATTAGAATGCGTAAAGTTACGCCTGCAATGAGCAACCAGTATCCCAGAATGGCCAACCAACTCACCAACGTATAAACGGTTGTCATAGATTAGAAATCCTTTTGAAAGCGCATTGTTATGAGTTTACGCATATGGATTCATCT  >  REL606/1306153‑1306296                                                                         |                                                                        cATTAGAATGCGTAAAGTTACGCCTGCAATGAGCAACCAGTATCCCAGAATGGCCAACCAACTCACCAACGTTTaa                                                                      >  1:1113580/1‑76 (MQ=255)      gAATGCGTAAAGTTACGCCTGCAATGAGCAACCAGTATCCCAGAATGGCCAACCAACTCCCCAACGTTTAAACGGt                                                                 <  2:1677573/76‑1 (MQ=255)      gAATGCGTAAAGTTACGCCTGCAATGAGCAACCAGTATCCCAGAATGGCCAACCAACTCACCAACGTTTAAACGGt                                                                 >  2:1706872/1‑76 (MQ=255)      gAATGCGTAAAGTTACGCCTGCAATGAGCAACCAGTATCCCAGAATGGCCAACCAACTCACCAACGTTTAAACGGt                                                                 <  2:1677561/76‑1 (MQ=255)        aTGCGTAAAGTTACGCCTGCAATGAGCAACCAGTATCCCAGAATGGCCAACCAACTCACCAACGTTTAAACGGTTg                                                               <  2:357990/76‑1 (MQ=255)         tGCGTAAAGTTACGCCTGCAATGAGCAACCAGTATCCCAGAATGGCCAACCAACTCACCAACGTATAAACGGTTGt                                                              <  2:606044/76‑1 (MQ=255)           cGTAAAGTTACGCCTGCAATGAGCAACCAGTATCCCAGAATGGCCAACCAACTCACCAACGTTTAAACGGTTGTCa                                                            <  1:809589/76‑1 (MQ=255)            gTAAAGTTACGCCTGCAATGAGCAACCAGTATCCCAGAATGGCCAACCAACTCACCAACGTTTAAACGGTTGTCAt                                                           <  1:554723/76‑1 (MQ=255)            gTAAAGTTACGCCTGCAATGAGCAACCAGTATCCCAGAATGGCCAACCAACTCACCAACGTTTAAACGGTTGTAAt                                                           >  1:1302237/1‑76 (MQ=255)                 gTTACGCCTGCAATGAGCAACCAGTATCCCAGAATGGCCAACCAACTCACCAACGTTTAAACGGTTGTCATAGAtt                                                      <  1:760440/76‑1 (MQ=255)                 gTTACGCCTGCAATGAGCAACCAGTATCCCAGAATGGCCAACCAACTCACCAACGTTTAAACGGTTGTCATAGAtt                                                      <  2:583088/76‑1 (MQ=255)                   tACGCCTGCAATGAGCAACCAGTATCCCAGAATGGCCAACCAACTCACCAACGTTTAAACGGTTGTCATAGATTAg                                                    <  1:939723/76‑1 (MQ=255)                     cGCCTGCAATGAGCAACCAGTATCCCAGAATGGCCAACCAACTCACCAACGTTTAAACGGTTGTCATAGATTAGaa                                                  <  1:571068/76‑1 (MQ=255)                      gCCTGCAATGAGCAACCAGTATCCCAGAATGGCCAACCAACTCACCAACGTTTAAACGGTTGTCATAGATTAGaaa                                                 >  1:157616/1‑76 (MQ=255)                       ccTGCAATGAGCAACCAGTATCCCAGAATGGCCAACCAACTCACCAACGTTTAAACGGTTGTCATAGATTAGAAAt                                                >  1:1300639/1‑76 (MQ=255)                       ccTGCAATGAGCAACCAGTATCCCAGAATGGCCAACCAACTCACCAACGTTTAAACGGTTGTCATAGATCAGAAAt                                                >  2:207462/1‑76 (MQ=255)                        cTGCAATGAGCAACCAGTATCCCAGAATGGCCAACCAACTCACCAACGTTTAAACGGTTGTCATAGATTAGAAATc                                               <  1:1644289/76‑1 (MQ=255)                         tGCAATGAGCAACCAGTATCCCAGAATGGCCAACCAACTCACCAACGTTTAAACGGTTGTCATAGATTAGAAATcc                                              <  1:1405154/76‑1 (MQ=255)                         tGCAATGAGCAACCAGTATCCCAGAATGGCCAACCAACTCACCAACGTTTAAACGGTTGTCATAGATTAGAAATcc                                              >  1:653069/1‑76 (MQ=255)                          gCAATGAGCAACCAGTATCCCAGAATGGCCAACCAACTCACCAACGTTTAAACGGTTGTCATAGATTAGAAATCCt                                             >  2:1420489/1‑76 (MQ=255)                           cAATGAGCAACCAGTATCCCAGAATGGCCAACCAACTCACCAACGTTTAAACGGTTGTCATAGATTAGAAATCCtt                                            >  1:1261629/1‑76 (MQ=255)                           cAATGAGCAACCAGTATCCCAGAATGGCCAACCAACTCACCAACGTTTAAACGGTTGTCATAGATTAGAAATCCtt                                            >  2:733674/1‑76 (MQ=255)                            aaTGAGCAACCAGTATCCCAGAATGGCCAACCAACTCACTAACGTTTAAACGGTTGTCATAGATTAGAAATCCttt                                           >  2:638744/1‑76 (MQ=255)                            aaTGAGCAACCAGTATCCCAGAATGGCCAACCAACTCACCAACGTTTAAACGGTTGTCATAGATTAGAAATCCttt                                           >  1:268577/1‑76 (MQ=255)                            aaTGAGCAACCAGTATCCCAGAATGGCCAACCAACTCACCAACGTTTAAACGGTTGTCATAGATTAGAAATCCttt                                           >  1:1322221/1‑76 (MQ=255)                            aaTGAGCAACCAATAACCCAGAATGGCCAACCAACTCACCAACGTTTAAACGGTTGTCATAGATTAGAAATCCttt                                           <  2:376474/76‑1 (MQ=255)                             ataaGCAACCAGTATCCCAGAATGGCCAACCAACTCACCAACGTTTAAACGGTTGTCATAGATTAGAAATCCtttt                                          <  2:379962/73‑1 (MQ=255)                             aTGAGCAACCAGTATCCCAGAATGGCCAACCAACTCACCAACGTTTAAACGGTTGTCATAGATTAGAAATCCtttt                                          >  1:806859/1‑76 (MQ=255)                              tGAGCAACCAGTATCCCAGAATGGCCAACCAACTCACCAACGTTTAAACGGTTGTCATAGATTAGAAATCCTTTTg                                         <  1:1642551/76‑1 (MQ=255)                               gAGCAACCAGTATCCCAGAATGGCCAACCAACTCACCAACGTTTAAACGGTTGTCATAGATTAGAAATCCTTTTGa                                        <  1:354763/76‑1 (MQ=255)                                aGCAACCAGTATCCCAGAATGGCCAACCAACTCACCAACGTTTAAACGGTTGTCATAGATTAGAAATCCTTTTGaa                                       <  2:20939/76‑1 (MQ=255)                                aGCAACCAGTATCCCAGAATGGCCAACCAACTCACCAACGTTTAAACGGTTGTCATAGATTAGAAATCCTTTTGaa                                       >  1:1304710/1‑76 (MQ=255)                                aGCAACCAGTATCCCAGAATGGCCAACCAACTCACCAACGTTTAAACGGTTGTCATAGATTAGAAATCCTTTTGaa                                       <  1:1246028/76‑1 (MQ=255)                                    aCCAGTATCCCAGAATGGCCAACCAACTCACCAACGTTTAAACGGTTGTCATAGATTAGAAATCCTTTTGAAAgcg                                   <  1:1048651/76‑1 (MQ=255)                                      cAGTATCCCAGAATGGCCAACCAACTCACCAACGTTTAAACGGTTGTCATAGATTAGAAATCCTTTTGAAAGCGCa                                 >  1:673779/1‑76 (MQ=255)                                      cAGTATCCCAGAATGGCCAACCAACTCACCAACGTTTAAACGGTTGTCATAGATTAGAAATCCTTTTGAAAGCGCa                                 >  1:1671772/1‑76 (MQ=255)                                      cAGTATCCCAGAATGGCCAACCAACTCACCAACGTTTAAACGGTTGTCATAGATTAGAAATCCTTTTGAAAGCGCa                                 >  2:1576866/1‑76 (MQ=255)                                        gTATCCCAGAATGGCCAACCAACTCACCAACGTTTAAACGGTTGTCATAGATTAGAAATCCTTTTGAAAGCGCAtt                               >  2:170650/1‑76 (MQ=255)                                          aTCCCAGAATGGCCAACCAACTCACCAACGTTTAAACGGTTGTCATAGATTAGAAATCCTTTTGAAAGCGCATTGt                             >  1:1108348/1‑76 (MQ=255)                                            cccAGAATGGCCAACCAACTCACCAACGTTTAAACGGTTGTCATAGATTAGAAATCCTTTTGAAAGCGCATTGTTa                           <  1:1356886/76‑1 (MQ=255)                                              cAGAATGGCCAACCAACTCACCAACGTTTAAACGGTTGTCATAGATTAGAAATCCTTTTGAAAGCGCATTGTTATg                         <  1:75/76‑1 (MQ=255)                                              cAGAATGGCCAACCAACTCACCAACGTTTAAACGGTTGTCATAGATTAGAAATCCTTTTGAAAGCGCATTGTTATg                         <  2:1710012/76‑1 (MQ=255)                                                 aaTGGCCAACCAACTCACCAACGTTTAAACGGTTGTCATAGATTAGAAATCCTTTTGAAAGCGCATTGTTATGAGt                      >  2:325786/1‑76 (MQ=255)                                                    ggCCAACCAACTCACCAACGTTTAAACGGTTGTCATAGATTAGAAATCCTTTTGAAAGCGCATTGTTATGAGTTTa                   >  1:496661/1‑76 (MQ=255)                                                     gCCAACCAACTCACCAACGTTTAAACGGTTGTCATAGATTAGAAATCCTTTTGAAAGCGCATTGTTATGAGTTTAc                  <  2:1486410/76‑1 (MQ=255)                                                      ccaaccaaCTCACCAACGTTTAAACGGTTGTCATAGATTAGAAATCCTTTTGAAAGCGCATTGTTATGAGTTTACg                 <  1:1007478/76‑1 (MQ=255)                                                       caaccaacTCACCAACGTTTAAACGGTTGTCATAGATTAGAAATCCTTTTGAAAGCGCATTGTTATGAGTTTACGc                >  2:1221187/1‑76 (MQ=255)                                                        aaccaacTCACCAACGTTTAAACGGTTGTCATAGATTAGAAATCCTTTTGAAAGCGCATTGTTATGAGTTTACGCa               >  1:206985/1‑76 (MQ=255)                                                        aaccaacTCACCAACGTTTAAACGGTTGTCATAGATTAGAAATCCTTTTGAAAGCGCATTGTTATGAGTTTACGCa               >  2:886278/1‑76 (MQ=255)                                                           caacTCACCAACGTTTAAACGGTTGTCATAGATTAGAAATCCTTTTGAAAGCGCATTGTTATGAGTTTACGCatat            >  2:1146890/1‑76 (MQ=255)                                                               tCACCAACGTTTAAACGGTTGTCATAGATTAGAAATCCTTTTGAAAGCGCATTGTTATGAGTTTACGCATATGGAt        <  1:1363202/76‑1 (MQ=255)                                                                cACCAACGTTTAAACGGTTGTCATAGATTAGAAATCCTTTTGAAAGCGCATTGTTATGAGTTTACGCATATGGAtt       >  2:1044747/1‑76 (MQ=255)                                                                  ccAACGTTTAAACGGTTGTCATAGATTAGAAATCCTTTTGAAAGCGCATTGTTATGAGTTTACGCATATGGATTCa     <  1:1213826/76‑1 (MQ=255)                                                                     aCGTTTAAACGGTTGTCATAGATTAGAAATCCTTTTGAAAGCGCATTGTTATGAGTTTACGCATATGGATTCATCt  <  1:1478778/76‑1 (MQ=255)                                                                         |                                                                        CATTAGAATGCGTAAAGTTACGCCTGCAATGAGCAACCAGTATCCCAGAATGGCCAACCAACTCACCAACGTATAAACGGTTGTCATAGATTAGAAATCCTTTTGAAAGCGCATTGTTATGAGTTTACGCATATGGATTCATCT  >  REL606/1306153‑1306296

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