Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,307,634	T→A	I338F (ATC→TTC)	clsA ←	cardiolipin synthase 1

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,307,634	0	T	A	97.8%	154.4 / NA	45	I338F (ATC→TTC)	clsA	cardiolipin synthase 1
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base A (27/17);  total (27/18)
Fisher's exact test for biased strand distribution p-value = 4.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCGGAAGGATAATACTGACATCCACCCCGCGCTGCGCCGCCGTGCAAATCGCATGAAGTAAATCATCGCTTGGCACAAAGTAGGGCGTGGTCATGATCAAATATTCACGCGCCGAATAAGCCGCAGTCAATAATGCCTGGTGAATGAGATCTTCCGGAAAGCCGGGGCCAGAAGCAATTGTGTGAATGGTGT  >  NC_000913/1307537‑1307729                                                                                                  |                                                                                                cgcgGAAGGATAATACTGACATCCACCCCGCGCTGCGCCGCCGTGCAAATCGCATGAAGTAAATCATCGCTTGGCACAAAGTAGGGCGTGGTCATGAACa                                                                                               <  2:1460794/100‑1 (MQ=255)   cgGAAGGATAATACTGACATCCACCCCGCGCTGCGCCGCCGTGCAAATCGCATGAAGTAAATCATCGCTTGGCACAAAGTAGGGCGTGGTCATGAACaaa                                                                                             <  2:1256424/100‑1 (MQ=255)   cgGAAGGATAATACTGACATCCACCCCGCGCTGCGCCGCCGTGCAAATCGCATGAAGTAAATCATCGCTTGGCACAAAGTAGGGCGTGGTCATGAACaaa                                                                                             <  1:283597/100‑1 (MQ=255)    ggAAGGATAATACTGACATCCACCCCGCGCTGCGCCGCCGTGCAAATCGCATGAAGTAAATCATCGCTTGGCACAAAGTAGGGCGTGGTCATGAACAAat                                                                                            >  2:1280137/1‑100 (MQ=255)              taCTGACATCCACCCCGCGCTGCGCCGCCGTGCAAATCGCATGAAGTAAATCATCGCTTGGCACAAAGTAGGGCGTGGTCATGAACAAATATTCACgcgc                                                                                  >  2:543057/1‑100 (MQ=255)              taCTGACATCCACCCCGCGCTGCGCCGCCGTGCAAATCGCATGAAGTAAATCATCGCTTGGCACAAAGTAGGGCGTGGTCATGAACAAATATTCACgcgc                                                                                  >  2:456267/1‑100 (MQ=255)                cTGACATCCACCCCGCGCTGCGCCGCCGTGCAAATCGCATGAAGTAAATCATCGCTTGGCACAAAGTAGGGCGTGGTCATGAACAAATATTCACGCGCCg                                                                                >  2:1381165/1‑100 (MQ=255)                   aCATCCACCCCGCGCTGCGCCGCCGTGCAAATCGCATGAAGTAAATCATCGCTTGGCACAAAGTAGGGCGTGGTCATGAACAAATATTCACGCGCCGAAt                                                                             >  1:1267735/1‑100 (MQ=255)                    cATCCACCCCGCGCTGCGCCGCCGTGCAAATCGCATGAAGTAAATCATCGCTTGGCACAAAGTAGGGCGTGGTCATGAACAAATATTCACGCGCCGAATa                                                                            >  2:1488808/1‑100 (MQ=255)                     aTCCACCCCGCGCTGCGCCGCCGTGCAAATCGCATGAAGTAAATCATCGCTTGGCACAAAGTAGGGCGTGGTCATGAACAAATATTCACGCGCCGAATaa                                                                           <  1:1509344/100‑1 (MQ=255)                      tCCACCCCGCGCTGCGCCGCCGTGCAAATCGCATGAAGTAAATCATAGCTTGGCACAAAGTAGGGCGTGGTCATGAAGAAATATTCACGCGCCGAATa                                                                            >  1:171305/1‑98 (MQ=255)                       ccACCCCGCGCTGCGCCGCCGTGCAAATCGCATGAAGTAAATCATCGCTTGGCACAAAGTAGGGCGTGGTCATGAACAAATATTCACGCGCCGAATAAGc                                                                         >  2:194710/1‑100 (MQ=255)                       ccACCCCGCGCTGCGCCGCCGTGCAAATCGCATGAAGTAAATCATCGCTTGGCACAAAGTAGGGCGTGGTCATGAACAAATATTCACGCGCCGAATAAGc                                                                         >  2:738384/1‑100 (MQ=255)                         aCCCCGCGCTGCGCCGCCGTGCAAATCGCATGAAGTAAATCATCGCTTGGCACAAAGTAGGGCGTGGTCATGAACAAATATTCACGCGCCGAATAAGCCg                                                                       >  1:800666/1‑100 (MQ=255)                            ccGCGCTGCGCCGCCGTGCAAATCGCATGAAGTAAATCATCGCTTGGCACAAAGTAGGGCGTGGTCATGAACAAATATTCACGCGCCGAATAAGCCGCAg                                                                    <  1:1244410/100‑1 (MQ=255)                            ccGCGCTGCGCCGCCGTGCAAATCGCATGAAGTAAATCATCGCTTGGCACAAAGTAGGGCGTGGTCATGAACAAATATTCACGCGCCGAATAAGCCGCAg                                                                    <  1:552032/100‑1 (MQ=255)                            ccGCGCTGCGCCGCCGTGCAAATCGCATGAAGTAAATCATCGCTTGGCACAAAGTAGGGCGTGGTCATGAACAAATATTCACGCGCCGAATAAGCCGCAg                                                                    >  2:1017915/1‑100 (MQ=255)                             cgcgCTGCGCCGCCGTGCAAATCGCATGAAGTAAATCATCGCTTGGCACAAAGTAGGGCGTGGTCATGAACAAATATTCACGCGCCGAATAAGCCGCAGt                                                                   <  2:580821/100‑1 (MQ=255)                              gcgcTGCGCCGCCGTGCAAATCGCATGAAGTAAATCATCGCTTGGCACAAAGTAGGGCGTGGTCATGAACAAATATTCACGCGCCGAATAAGCCGCAGTc                                                                  <  2:1564161/100‑1 (MQ=255)                              gcgcTGCGCCGCCGTGCAAATCGCATGAAGTAAATCATCGCTTGGCACAAAGTAGGGCGTGGTCATGAACAAATATTCACGCGCCGAATAAGCCGCAGTc                                                                  <  2:1564128/100‑1 (MQ=255)                                gcTGCGCCGCCGTGCAAATCGCGTGAAGTAAATCATCGCTTGGCACAAAGTATGGCGTAGTCAAGAACAAATATTCACGCGCCGAATAAGCCGCAGACaa                                                                >  2:751445/1‑100 (MQ=255)                                 cTGCGCCGCCGTGCAAATCGCATGAAGTAAATCATCGCTTGGCACAAAGTAGGGCGTGGTCATGAACAAATATTCACGCGCCGAATAAGCCGCAGTCaat                                                               >  2:265758/1‑100 (MQ=255)                                         ccgTGCAAATCGCATGAAGTAAATCATCGCTTGGCACAAAGTAGGGCGTGGTCATGAACAAATATTCACGCGCCGAATAAGCCGCAGTCAATAATGCCTg                                                       <  1:1504818/100‑1 (MQ=255)                                          cgTGCAAATCGCATGAAGTAAATCATCGCTTGGCACAAAGTAGGGCGTGGTCATGAACAAATATTCACGCGCCGAATAAGCCGCAGTCAATAATGCCTgg                                                      >  2:728713/1‑100 (MQ=255)                                            tGCAAATCGCATGAAGTAAATCATCGCTTGGCACAAAGTAGGGCGTGGTCATGAACAAATATTCACGCGCCGAATAAGCCGCAGTCAATAATGCCTGGTg                                                    >  1:586697/1‑100 (MQ=255)                                            tGCAAATCGCATGAAGTAAATCATCGCTTGGCACAAAGTAGGGCGTGGTCATGAACAAATATTCACGCGCCGAATAAGCCGCAGTCAATAATGCCTGGTg                                                    <  2:171305/100‑1 (MQ=255)                                            tGCAAATCGCATGAAGTAAATCATCGCTTGGCACAAAGTAGGGCGTGGTCATGAACAAATATTCACGCGCCGAATAAGCCGCAGTCAATAATGCCTGGTg                                                    <  1:472108/100‑1 (MQ=255)                                             gCAAATCGCATGAAGTAAATCATCGCTTGGCACAAAGTAGGGCGTGGTCATGAACAAATATTCACGCGCCGAATAAGCCGCAGTCAATAATGCCTGGTGa                                                   >  2:646847/1‑100 (MQ=255)                                                   cGCATGAAGTAAATCATCGCTTGGCACAAAGTAGGGCGTGGTCATGAACAAATATTCACGCGCCGAATAAGCCGCAGTCAATAATGCCTGGTGAATgaga                                             >  1:568973/1‑100 (MQ=255)                                                    gCATGAAGTAAATCATCGCTTGGCACAAAGTAGGGCGTGGTCATGAACAAATATTCACGCGCCGAATAAGCCGCAGTCAATAATGCCTGGTGAATGAGAt                                            <  1:751909/100‑1 (MQ=255)                                                          aGTAAATCATCGCTTGGCACAAAGTAGGGCGTGGTCATGAACAAATATTCACGCGCCGAATAAGCCGCAGTCAATAATGCCTGGTGAATGAGATCTTcc                                       >  1:656384/1‑99 (MQ=255)                                                                   tcGCTTGGCACAAAGTAGGGCGTGGTCATGAACAAATATTCACGCGCCGAATAAGCCGCAGTCAATAATGCCTGGTGAATGAGATCTTCCGGAAAGCCgg                             <  1:1488817/100‑1 (MQ=255)                                                                         ggCACAAAGTAGGGCGTGGTCATGAACAAATATTCACGCGCCGAATAAGCCGCAGTCAATAATGCCTGGTGAATGAGATCTTCCGGAAAGCCGGGGCCAg                       <  2:1285005/100‑1 (MQ=255)                                                                           cacaAAGTAGGGCGTGGTCATGAACAAATATTCACGCGCCGAATAAGCCGCAGTCAATAATGCCTGGTGAATGAGATCTTCCGGAAAGCCGGGGCCAGaa                     >  2:1493065/1‑100 (MQ=255)                                                                             caAAGTAGGGCGTGGTCATGAACAAATATTCACGCGCCGAATAAGCCGCAGTCAATAATGCCTGGTGAATGAGATCTTCCGGAAAGCCGGGGCCAGAAGc                   >  1:352368/1‑100 (MQ=255)                                                                             caAAGTAGGGCGTGGTCATGAACAAATATTCACGCGCCGAATAAGCCGCAGTCAATAATGCCTGGTGAATGAGATCTTCCGGAAAGCCGGGGCCAGAAGc                   >  1:1549681/1‑100 (MQ=255)                                                                                   aGGGCGTGGTCATGAACAAATATTCACGCGCCGAATAAGCCGCAGTCAATAATGCCTGGTGAATGAGATCTTCCGGAAAGCCGGGGCCAGAAGCAATtg              >  2:1115735/1‑99 (MQ=255)                                                                                    gggCGTGGTCATGAACAAATATTCACGCGCCGAATAAGCCGCAGTCAATAATGCCTGGTGAATGAGATCTTCCGGAAAGCCGGGGCCAGAAGCAATtgtg            >  2:505627/1‑100 (MQ=255)                                                                                         tGGTCATGAACAAATATTCACGCGCCGAATAAGCCGCAGTCAATAATGCCTGGTGAATGAGATCTTCCGGAAAGCCGGGGCCAGAAGCAATTGTGTGAAt       <  1:265758/100‑1 (MQ=255)                                                                                          ggTCATGAACAAATATTCACGCGCCGAATAAGCCGCAGTCAATAATGCCTGGTGAATGAGATCTTCCGGAAAGCCGGGGCCAGAAGCAATTGTGTGAATg      <  2:445905/100‑1 (MQ=255)                                                                                          ggTCATGAACAAATATTCACGCGCCGAATAAGCCGCAGTCAATAATGCCTGGTGAATGAGATCTTCCGGAAAGCCGGGGCCAGAAGCAATTGTGTGAATg      >  2:1245204/1‑100 (MQ=255)                                                                                             cATGATCAAATATTCACGCGCCGAATAAGCCGCAGTCAATAATGCCTGGTGAATGAGATCTTCCGGAAAGCCGGGGCCAGAAGCAATTGTGTGAATGgtg   <  2:266502/100‑1 (MQ=255)                                                                                             cATGAACAAATATTCACGCGCCGAATAAGCCGCAGTCAATAATGCCTGGTGAATGAGATCTTCCGGAAAGCCGGGGCCAGAAGCAATTGTGTGAATGgtg   >  1:1277793/1‑100 (MQ=255)                                                                                             cATGAACAAATATTCACGCGCCGAATAAGCCGCAGTCAATAATGCCTGGTGAATGAGATCTTCCGGAAAGCCGGGGCCAGAAGCAATTGTGTGAATGgtg   >  1:131186/1‑100 (MQ=255)                                                                                              aTGAACAAATATTCACGCGCCGAATAAGCCGCAGTCAATAATGCCTGGTGAATGAGATCTTCCGGAAAGCCGGGGCCAGAAGCAATTGTGTGAATGgtgt  >  2:163587/1‑100 (MQ=255)                                                                                                  |                                                                                                CGCGGAAGGATAATACTGACATCCACCCCGCGCTGCGCCGCCGTGCAAATCGCATGAAGTAAATCATCGCTTGGCACAAAGTAGGGCGTGGTCATGATCAAATATTCACGCGCCGAATAAGCCGCAGTCAATAATGCCTGGTGAATGAGATCTTCCGGAAAGCCGGGGCCAGAAGCAATTGTGTGAATGGTGT  >  NC_000913/1307537‑1307729

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