Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	1,305,155	A→C	V362G (GTG→GGG)	cls ←	cardiolipin synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	1,305,155	0	A	C	97.6%	158.5 / ‑3.9	41	V362G (GTG→GGG)	cls	cardiolipin synthetase
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base C (23/17);  total (23/18)
Fisher's exact test for biased strand distribution p-value = 4.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.85e-01

CGTAAAGAATGCGCGACTGGCCCAGCCGACCAGCATCGAGTCA‑TTTTTTCGCGGAAGGATAATACTGACATCCACCCCGCGCTGTGCCGCCGTGCAAATCGCATGAAGTAAATCATCGCTTGGCACAAAGTAGGGCGTTGTCATG  >  REL606/1305082‑1305226                                                                           |                                                                        cGTAAAGAATGCGCGACTGGCCCAGCCGACCAGCATCGAGTCA‑TTTTTTCGCGGAAGGATAATACTGACATCcccc                                                                       >  2:1296622/1‑76 (MQ=255)  gTAAAGAATGCGCGACTGGCCCAGCCGACCAGCATCGAGTCATTTTTTTCGCGGAAGGATAATACTGACATCcccc                                                                       >  2:168417/1‑76 (MQ=255)  gTAAAGAATGCGCGACTGGCCCAGCCGACCAGCATCGAGTCA‑TTTTTTCGCGGAAGGATAATACTGACATcccccc                                                                      >  1:1080161/1‑76 (MQ=255)    aaaGAATGCGCGACTGGCCCAGCCGACCAGCATCGAGTCA‑TTTTTTCGCGGAAGGATAATACTGACATCCCCCCcg                                                                    >  2:1010970/1‑76 (MQ=255)     aaGAATGCGCGACTGGCCCAGCCGACCAGCATCGAGTCA‑TTTTTTCGCGGAAGGATAATACTGACATCCCCCCcgc                                                                   >  2:1873444/1‑76 (MQ=255)        aaTGCGCGACTGGCCCAGCCGACCAGCATCGAGTCA‑TTTTTTCGCGGAAGGATAATACTGACATCCCCCCCGCGCt                                                                >  1:1376766/1‑76 (MQ=255)         aTGCGCGACTGGCCCAGCCGACCAGCATCGAGTCA‑TTTTTTCGCGGAAGGATAATACTGACATCCCCCCCGCGCtg                                                               >  2:699854/1‑76 (MQ=255)         aTGCGCGACTGGCCCAGCCGACCAGCATCGAGTCA‑TTTTTTCGCGGAAGGATAATACTGACATCCCCCCCGCGCtg                                                               <  1:1162605/76‑1 (MQ=255)         aTGCGCGACTGGCCCAGCCGACCAGCATCGAGTCA‑TTTTTTCGCGGAAGGATAATACTGACATCCCCCCCGCGCtg                                                               >  2:461970/1‑76 (MQ=255)          tccgcgACTGGCGCAGCCGACCAGCATCGAGTCA‑TTTTTTCGCGGAAGGATAATACTGACATCCCCCCCGCGCtgt                                                              <  2:953848/74‑1 (MQ=255)          tGCGCGACTGGCCCAGCCGACCAGCATCGAGTCA‑TTTTTTCGCGGAAGGATAATACTGACATCCCCCCCGCGCtgt                                                              >  1:805551/1‑76 (MQ=255)           gcgcGACTGGCCCAGCCGACCAGCATCGAGTCA‑TTTTTTCGCGGAAGGATAATACTGACATCCCCCCCGCGCtgtg                                                             <  1:1488467/76‑1 (MQ=255)               gACTGGCCCAGCCGACCAGCATCGAGTCA‑TTTTTTCGCGGAAGGATAATACTGACATCCCCCCCGCGCTGTgccgc                                                         >  1:1825068/1‑76 (MQ=255)                aCTGGCCCAGCCGACCAGCATCGAGTCA‑TTTTTTCGCGGAAGGATAATACTGACATCCCCCCCGCGCTGTgccgcc                                                        >  1:1661340/1‑76 (MQ=255)                  tGGCCCAGCCGACCAGCATCGAGTCA‑TTTTTTCGCGGAAGGATAATACTGACATCCCCCCCGCGCTGTGCCGCCGt                                                      >  2:301806/1‑76 (MQ=255)                     cccAGCCGACCAGCATCGAGTCA‑TTTTTTCGCGGAAGGATAATACTGACATCCCCCCCGCGCTGTGCCGCCGTGcc                                                   >  1:28716/1‑75 (MQ=255)                     cccAGCCGACCAGCATCGAGTCA‑TTTTTTCGCGGAAGGATAATACTGACATCCCCCCCGCGCTGTGCCGCCGTGCa                                                   >  1:1897110/1‑76 (MQ=255)                       cAGCCGACCAGCATCGAGTCA‑TTTTTTCGCGGAAGGATAATACTGACATCCCCCCCGCGCTGTGCCGCCGTGCaaa                                                 <  1:359935/76‑1 (MQ=255)                        aGCCGACCAGCATCGAGTCA‑TTTTTTCGCGGAAGGATAATACTGACATCCCCCCCGCGCTGTGCCGCCGTGCAAAt                                                >  1:962223/1‑76 (MQ=255)                         gCCGACCAGCATCGAGTCA‑TTTTTTCGCGGAAGGATAATACTGACATCCACCCCGCGCTGTGCCGCCGTGCAAATc                                               <  1:276680/76‑1 (MQ=255)                              ccAGCATCGAGTCA‑TTTTTTCGCGGAAGGATAATACTGACATCCCCCCCGCGCTGTGCCGCCGCGCAAATCGCATg                                          <  1:889120/76‑1 (MQ=255)                               cAGCATCGAGTCA‑TTTTTTCGCGGAAGGATAATACTGACATCCCCCCCGCGCTGTGCCGCCGTGCAAATCGCATGa                                         >  2:629635/1‑76 (MQ=255)                                  gagcGAGTCT‑TTTGTGCGCGGAATGATGATACTGACAGCCCCCCCGCGCTGTGCCGCCGTGCAAATCGCATGAAGt                                      <  1:253144/73‑1 (MQ=255)                                  cATCGAGTCA‑TTTTTTCGCGGAAGGATAATACTGACATCCCCCCCGCGCTGTGCCGCCGTGCAAATCGCATGAAGt                                      <  1:253155/76‑1 (MQ=255)                                   aTCGAGTCA‑TTTTTTCGCGGAAGGATAATACTGACACCCCCCCCGCGCTGTGCCGCCGTGCAAATCGCATGAAGTa                                     <  1:346476/76‑1 (MQ=255)                                    tCGAGTCA‑TTTTTTCGCGGAAGGATAATACTGACATCCCCCCCGCGCTGTGCCGCCGTGCAAATCGCATGAAGTaa                                    >  1:150856/1‑76 (MQ=255)                                     cGAGTCA‑TTTTTTCGCGGAAGGATAATACTGACATCCCCCCCGCGCTGTGCCGCCGTGCAAATCGCATGAAGTaaa                                   <  2:1389062/76‑1 (MQ=255)                                      gAGTCG‑TTTTTTCGCGGAAGGATGATACTGACATCCCCCCCGCGCTGTGCCGCCGTGCAAATCGCATGAAGTAAat                                  <  1:1215782/76‑1 (MQ=255)                                       aGTCA‑TTTTTTCGCGGAAGGATAATACTGACATCCCCCCCGCGCTGTGCCGCCGTGCAAATCGCATGAAGTAAatc                                 <  1:704292/76‑1 (MQ=255)                                           a‑TTTTTTCGCGGAAGGATAATACTGACATCCCCCCCGCGCTGTGCCGCCGTGCAAAt                                                <  2:1698286/57‑1 (MQ=255)                                           a‑TTTTTTCGCGGAAGGATAATACTGACATCCCCCCCGCGCTGTGCCGCCGTGCAAAt                                                >  1:1698286/1‑57 (MQ=255)                                           a‑TTTTTTCGCGGAAGGATAATACTGACATCCCCCCCGCGCTGTGCCGCCGTGCAAATCGCATGAAGTAAATCATCg                             <  1:410814/76‑1 (MQ=255)                                              tttttCGCGGAAGGATAATACTGACATCCCCCCCGCGCTGTGCCGCCGTGCAAATCGCATGAAGTAAATCATCGCt                           >  1:1680789/1‑76 (MQ=255)                                                  tCGCGAAAGGATAATACTGACATCCCCCCCGCGCTGTGCCGCCGTGCAAATCGCATGAAGTAAATCATCGCTTGGc                       <  2:906008/76‑1 (MQ=255)                                                   cgcgGAAGGATAATACTGACATCCCCCCCGCGCTGTGCCGCCGTGCAAATCGCATGAAGTAAATCATCGCTTGGca                      >  2:1277065/1‑76 (MQ=255)                                                   cgcgGAAGGATAATACTGACATCCCCCCCGCGCTGTGCCGCCGTGCAAATCGCATGAAGTAAATCATCGCTTGGca                      <  1:1729639/76‑1 (MQ=255)                                                    gcgGAAGGATAATACTGACATCCCCCCCGCGCTGTGCCGCCGTGCAAATCGCATGAAGTAAATCATCGCTTGGcac                     >  1:651003/1‑76 (MQ=255)                                                     cgGAAGGATAATACTGACATCCCCCCCGCGCTGTGCCGCCGTGCAAATCGCATGAAGTAAATCATCGCTTGGcaca                    >  1:1852584/1‑76 (MQ=255)                                                       gAAGGATAATACTGACATCCCCCCCGCGCTGTGCCGCCGTGCAAATCGCATGAAGTAAATCATCGCTTGGcaccgt                  <  2:101809/76‑4 (MQ=255)                                                       gAAGGACAATACTGACATCCCCCCCGCGCTGTGCCGCCGTGCAAATCGCATGAAGTAAATCATCGCTTGGcaccgt                  <  2:391623/76‑4 (MQ=255)                                                        aaGGATAATACTGACATCCCCCCCGCGCTGTGCCGCCGTGCAAATCGCATGAAGTAAATCATCGCTTGGCACAAAg                 >  2:742937/1‑76 (MQ=255)                                                         aGGATAATACTGACATCCCCCCCGCGCTGTGCCGCCGTGCAAATCGCATGAAGTAAATCATCGCTTGGCACAAAGt                <  1:95212/76‑1 (MQ=255)                                                           gATAATACTGACATCCCCCCCGCGCTGTGCCGCCGTGCAAATCGCATGAAGTAAATCATCGCTTGGCACAAAGTAg              <  1:1126471/76‑1 (MQ=255)                                                              aataCTGACATCCCCCCCGCGCTGTGCCGCCGTGCAAATCGCATGAAGTAAATCATCGCTTGGCACAAAGTAGGGc           >  1:974273/1‑76 (MQ=255)                                                              aataCTGACATCCCCCCCGCGCTGTGCCGCCGTGCAAATCGCATGAAGTAAATCATCGCTTGGCACAAAGTAGGGc           >  1:1002231/1‑76 (MQ=255)                                                               ataCTGACATCCCCCCCGCGCTGTGCCGCCGTGCAAATCGCATGAAGTAAATCATCGCTTGGCACAAAGTAGGGCg          >  1:553639/1‑76 (MQ=255)                                                                   cgACATCCCCCCCGCGCTGTGCCGCCGTGCAAATCGCATGAAGTAAATCATCGCTTGGCACAAAGTAGGGCGTTGt      <  2:75798/75‑1 (MQ=255)                                                                    gACATCCCCCCCGCGCTGTGCCGCCGTGCAAATCGCATGAAGTAAATCATCGCTTGGCACCAAGTAGGGCGTTGTc     >  1:1460516/1‑76 (MQ=255)                                                                       ctCCCCCCCGCGCTGTGCCGCCGTGCAAATCGCATGAAGTAAATCATCGCTTGGCACGAAGTAGGGCGTTGTCATg  <  2:289859/75‑1 (MQ=255)                                                                           |                                                                        CGTAAAGAATGCGCGACTGGCCCAGCCGACCAGCATCGAGTCA‑TTTTTTCGCGGAAGGATAATACTGACATCCACCCCGCGCTGTGCCGCCGTGCAAATCGCATGAAGTAAATCATCGCTTGGCACAAAGTAGGGCGTTGTCATG  >  REL606/1305082‑1305226

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