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	95.3%	216.9 / 0.5	65	I338F (ATC→TTC)	clsA	cardiolipin synthase 1
Reads supporting (aligned to +/- strand):  ref base T (1/2);  new base A (34/27);  total (36/29)
Fisher's exact test for biased strand distribution p-value = 5.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

GCGGAAGGATAATACTGACATCCACCCCGCGCTGCGCCGCCGTGCAAATCGCATGAAGTAAATCATCGCTTGGCACAAAGTAGGGCGTGGTCATGATCAAATATTCACGCGCCGAATAAGCCGCAGTCAATAATGCCTGGTGAATGAGATCTTCCGGAAAGCCGGGGCCAGAAGCAATTGTGTGAATGGTGTGA  >  NC_000913/1307538‑1307731                                                                                                 |                                                                                                  gcgGAAGGATAATACTGACATCCACCCCGCGCTGCGCCGCCGTGCAAATCGCATGAAGTAAATCATCGCTTGGCACAAAGTAGGGCGTGGTCATGAACaa                                                                                                <  1:416322/100‑1 (MQ=255)         ataataCTGACATCCACCCCGCGCTGCGCCGCCGTGCAAATCGCATGAAGTAAATCATCGCTTGGCAAAAAGTAGGGCGTGGTCATGAACAAATATTCAc                                                                                        >  2:473204/1‑100 (MQ=255)            ataCTGACATCCACCCCGCGCTGCGCCGCCGTGCAAATCGCATGAAGTAAATCATCGCTTGGCACAAAGTAGGGCGTGGTCATGAACAAATATTCAcgcg                                                                                     >  1:585980/1‑100 (MQ=255)             taCTGACATCCACCCCGCGCTGCGCCGCCGTGCAAATCGCATGAAGTAAATCATCGCTTGGCACAAAGTAGGGCGTGGTCATGAACAAATATTCACgcgc                                                                                    >  1:1355117/1‑100 (MQ=255)              aCTGACATCCACCCCGCGCTGCGCCGCCGTGCAAATCGCATGAAGTAAATCATCGCTTGGCACAAAGTAGGGCGTGGTCATGAACTAATATTCACGCGcc                                                                                   >  1:1325317/1‑100 (MQ=255)               cTGACATCCACCCCGCGCTGCGCCGCCGTGCAAATCGCATGAAGTAAATCATCGCTTGGCACAAAGTAGGGCGTGGTCATGATCAAATATTCACGCGCCg                                                                                  <  2:354186/100‑1 (MQ=255)               cTGACATCCACCCCGCGCTGCGCCGCCGTGCAAATCGCATGAAGTAAATCATCGCTTGGCACAAAGTAGGGCGTGGTCATGAACAAATATTCACGCGCCg                                                                                  >  2:1587700/1‑100 (MQ=255)               cTGACATCCACCCCGCGCTGCGCCGCCGTGCAAATCGCATGAAGTAAATCATCGCTTGGCACAAAGTAGGGCGTGGTCATGAACAAATATTCACGCGCCg                                                                                  >  2:32495/1‑100 (MQ=255)                 gACATCCACCCCTCGCTGCGCCGCCGTTCAAATCGCATGAAGTAAATCATCGCTTGGCACAAAGTAGGGCGTGGTCATGAACAAATATTCACGCGCCg                                                                                  <  1:1842806/98‑1 (MQ=255)                 gACATCCACCCCGCGCTGCGCCGCCGTGCAAATCGCATGAAGTAAATCATCGCTTGGCACAAAGTAGGGCGTGGTCATGAACAAATATTCACGCGCCGaa                                                                                <  2:1372958/100‑1 (MQ=255)                   cATCCACCCCGCGCTGCGCCGCCGTGCAAATCGCATGAAGTAAATCATCGCTTGGCACAAAGTAGGGCGTGGTCATGAACAAATATTCACGCGCCGAATa                                                                              >  2:235955/1‑100 (MQ=255)                    aTCCACCCCGCGCTGCGCCGCCGTGCAAATCGCATGAAGTAAATCATCGCTTGGCACAAAGTAGGGCGTGGTCATGAACAAATATTCACGCGCCGAATaa                                                                             >  2:435359/1‑100 (MQ=255)                      ccACCCCGCGCTGCGCCGCCGTGCAAATCGCATGTAGTAAATCATCGCTTGGCACAAAGTAGGGCGTGGTCATGAACAAATATTCACGCGCCGAATAAGc                                                                           <  1:1027027/100‑1 (MQ=255)                      ccACCCCGCGCTGCGCCGCCGTGCAAATCGCATGAAGTAAATCATCGCTTGGCACAAAGTAGGGCGTGGTCATGAACAAATATTCACGCGCCGAATAAGc                                                                           <  2:1696455/100‑1 (MQ=255)                       cACCGCGCGCTGCGCCGCCGTGCAAATCGCATGAAGTAAATCATCGCTTGGCACAGAGTAGGGCGTGGTCATGAACAAATATTCACGCGCCGAATAAGcc                                                                          >  1:256584/1‑100 (MQ=255)                       cACCCCGCGCTGCGCCGCCGTGCAAATCGCATGAAGTAAATCATCGCTTGGCACAAAGTAGGGCGTGGTCATGAACAAATATTCACGCGCCGAATAAGcc                                                                          >  1:347514/1‑100 (MQ=255)                       cACCCCGCGCTGCGCCGCCGTGCAAATCGCATGAAGTAAATCATCGCTTGGCACAAAGTAGGGCGTGGTCATGAACAAATATTCACGCGCCGAATAAGcc                                                                          >  1:1477520/1‑100 (MQ=255)                         ccccGCGCTGCGCCGCCGTGCAAATCGCATGAAGTAAATCATCGCTTGGCACAAAGTAGGGCGTGGTCATGATCAAATATTCACGCGCCGAATAAGCCGc                                                                        <  1:316902/100‑1 (MQ=255)                         ccccGCGCTGCGCCGCCGTGCAAATCGCATGAAGTAAATCATCGCTTGGCACAAAGTAGGGCGTGGTCATGAACAAATATTCACGCGCCGAATAAGCCGc                                                                        <  2:225917/100‑1 (MQ=255)                         ccccGCGCTGCGCCGCCGTGCAAATCGCATGAAGTAAATCATCGCTTGGCACAAAGTAGGGCGTGGTCATGAACAAATATTCACGCGCCGAATAAGCCGc                                                                        <  2:111645/100‑1 (MQ=255)                         ccccGCGCTGCGCCGCCGTGCAAATCGCATGAAGTAAATCATCGCTTGGCACAAAGTAGGGCGTGGTCATGAACAAATATTCACGCGCCGAATAAGCCGc                                                                        <  2:1608290/100‑1 (MQ=255)                            cgcgCTGCGCCGCCGTGCAAATCGCATGAAGTAAATCATCGCTTGGCACAAAGTAGGGCGTGGTCATGAACAAATATTCACGCGCCGAATAAGCCGCAGt                                                                     <  1:673044/100‑1 (MQ=255)                              cgcTGCGCCGCCGTGCAAATCGCATGAAGTAAATCATCGCTTGGCACAAAGTAGGGCGTGGTCATGAACAAATATTCACGCGCCGAATAAGCCGCAGTCa                                                                   >  2:1093000/1‑100 (MQ=255)                                   cgccgcCGTGCAAATCGCATGAAGTAAATCATCGCTTGGCACAAAGTAGGGCGTGGTCATGAACAAATATTCACGCGCCGAATAAGCCGCAGTCaataat                                                              <  1:993331/100‑1 (MQ=255)                                     ccgccgTGCAAATCGCATGAAGTAAATCATCGCTTGGCACAAAGTAGGGCGTGGTCATGAACAAATATTCACGCGCCGAATAAGCCGCAGTCAATAATGc                                                            <  1:368574/100‑1 (MQ=255)                                          gTGCAAATCGCATGAAGTAAATCATCGCTTGGCACAAAGTAGGGCGTGGTCATGAACAAATATTCACGCGCCGAATAAGCCGCAGTCAATAATGCCTGGt                                                       >  2:445789/1‑100 (MQ=255)                                           tGCAAATCGCATGAAGTAAATCATCGCTTGGCACAAAGTAGGGCGTGGTCATGAACAAATATTCACGCGCCGAATAAGCCGCAGTCAATAATGCCTGGTg                                                      <  2:578539/100‑1 (MQ=255)                                            gCAAATCGCATGAAGTAAATCATCGCTTGGCACAAAGTAGGGCGTGGTCATGAACAGATATTCACGCGCCGAATAAGCCGCAGTCAATAATGCCTGGTGa                                                     >  1:1243315/1‑100 (MQ=255)                                            gCAAATCGCATGAAGTAAATCATCGCTTGGCACAAAGTAGGGCGTGGTCATGAACAAATATTCACGCGCCGAATAAGCCGCAGTCAATAATGCCTGGTGa                                                     >  2:630273/1‑100 (MQ=255)                                            gCAAATCGCATGAAGTAAATCATCGCTTGGCACAAAGTAGGGCGTGGTCATGAACAAATATTCACGCGCCGAATAAGCCGCAGTCAATAATGCCTGGTGa                                                     >  2:466793/1‑100 (MQ=255)                                             cAAATCGCATGAAGTAAATCATCGCTTGGCACAAAGTAGGGCGTGGTCATGAACAAATATTCACGCGCCGAATAAGCCGCAGTCAATAATGCCTGGTGaa                                                    >  1:1284960/1‑100 (MQ=255)                                               aaTCGCATGAAGTAAATCATCGCTTGGCACAAAGTAGGGCGTGGTCATGAACAAATATTCACGCGCCGACTAAGCCGCAGTCAATAATGCCCGGTGAATg                                                  >  2:544380/1‑100 (MQ=255)                                                aTCGCATGAAGTAAATCATCGCTTGGCACAAAGTAGGGCGTGGTCATGAACAAATATTCACGCGCCGAATAAGCCGCAGTCAATAATGCCTGGTGAATga                                                 >  2:134281/1‑100 (MQ=255)                                                    cATGAAGTAAATCATCGCTTGGCACAAAGTAGGGCGTGGTCATGAACAAATATTCACGCGCCGAATAAGCCGCAGTCAATAATGCCTGGTGAATGAGATc                                             <  1:1922957/100‑1 (MQ=255)                                                      tGAAGTAAATCATCGCTTGGCACAAAGTAGGGCGTGGTCATGAACAAATATTCACGCGCCGAATAAGCCGCAGTCAATAATGCCTGGTGAATGAGATCtt                                           >  2:582354/1‑100 (MQ=255)                                                        aaGTAAATCATCGCTTGGCACAAAGTAGGGCGTGGTCATGAACAAATATTCACGCGCCGAATAAGCCGCAGTCAATAATGCCTGGTGAATGAGATCTTcc                                         <  2:1335358/100‑1 (MQ=255)                                                        aaGTAAATCATCGCTTGGCACAAAGTAGGGCGTGGTCATGAACAAATATTCACGCGCCGAATAAGCCGCAGTCAATAATGCCTGGTGAATGAGATCTTcc                                         <  1:1525979/100‑1 (MQ=255)                                                        aaGTAAATCATCGCTTGGCACAAAGTAGGGCGTGGTCATGAACAAATATTCACGCGCCGAATAAGCCGCAGTCAATAATGCCTGGTGAATGAGATCTTcc                                         <  2:1901735/100‑1 (MQ=255)                                                          gTAAATCATCGCTTGGCACAAAGTAGGGCGTGGTCATGAACAAATATTCACGCGCCGAATAAGCCGCAGTCAATAATGCCTGGTGAATGAGATCTTCCgg                                       >  2:1422794/1‑100 (MQ=255)                                                            aaaTCATCGCTTGGCACAAAGTAGGGCGTGGTCATGAACAAATATTCACGCGCCGAATAAGCCGCAGTCAATAATGCCTGGTGAATGAGATCTTCCGGaa                                     >  1:1761111/1‑100 (MQ=255)                                                            aaaTCATCGCTTGGCACAAAGTAGGGCGTGGTCATGAACAAATATTCACGCGCCGAATAAGCCGCAGTCAATAATGCCTGGTGAATGAGATCTTCCGGaa                                     >  2:1038193/1‑100 (MQ=255)                                                                 agcGCTTGGCACAAAGTAGGGCGTGGTCATGAACAAATATTCACGCGCCGAATAAGCCGCAGTCAATAATGCCTGGTGAATGAGATCTTCCGGAAAGCCg                                >  2:1318919/3‑100 (MQ=255)                                                                  tcGCTTGGCACAAAGTAGGGCGTGGTCATGACCAAATATTCACGCGCCGAATAAGCAGCAGTCAATAATGCCTGGTGAATGAGAACTTCCGGAAAGCCgg                               >  1:1160053/1‑100 (MQ=255)                                                                  tcGCTTGGCACAAAGTAGGGCGTGGTCATGAACAAATATTCACGCGCCGAATAAGCCGCAGTCAATAATGCCTGGTGAATGAGATCTTCCGGAAAGCCgg                               <  2:1643287/100‑1 (MQ=255)                                                                   cGCTTGGCACAAAGTAGGGCGTGGTCATGAACAAATATTCACGCGCCGAATAAGCCGCAGTCAATAATGCCTGGTGAATGAGATCTTCCGGAAAGCCggg                              >  2:1764399/1‑100 (MQ=255)                                                                       tGGCACAAAGTAGGGCGTGGTCATGAACAAATATTCACGCGCCGAATAAGCCGCAGTCAATAATGCCTGGTGAATGAGATCTTCCGGAAAGCCGGGGCCa                          <  2:1504226/100‑1 (MQ=255)                                                                       tGGCACAAAGTAGGGCGTGGTCATGAACAAATATTCACGCGCCGAATAAGCCGCAGTCAATAATGCCTGGTGAATGAGATCTTCCGGAAAGCCGGGGCCa                          >  2:1944276/1‑100 (MQ=255)                                                                        ggCACAAAGTAGGGCGTGGTCATGAACAAATATTCACGCGCCGAATAAGCCGCAGTCAATAATGCCTGGTGAATGAGATCTTCCGGAAAGCCGGGGCCAg                         >  1:1730443/1‑100 (MQ=255)                                                                         gCACAAAGTAGGGCGTGGTCATGAACAAATATTCACGCGCCGAATAAGCCGCAGTCAATAATGCCTGGTGAATGAGATCTTCCGGAAAGCCGGGGCCAGa                        >  1:1655367/1‑100 (MQ=255)                                                                         gCACAAAGTAGGGCGTGGTCATGAACAAATATTCACGCGCCGAATAAGCCGCAGTCAATAATGCCTGGTGAATGAGATCTTCCGGAAAGCCGGGGCCAGa                        >  2:1192231/1‑100 (MQ=255)                                                                          cacaAAGTAGGGCGTGGTCATGAACAAATATTCACGCGCCGAATAAGCCGCAGTCAATAATGCCTGGTGAATGAGATCTTCCGGAAAGCCGGGGCCAGaa                       >  2:482649/1‑100 (MQ=255)                                                                          cacaAAGTAGGGCGTGGTCATGAACAAATATTCACGCGCCGAATAAGCCGCAGTCAATAATGCCTGGTGAATGAGATCTTCCGGAAAGCCGGGGCCAGaa                       >  1:46199/1‑100 (MQ=255)                                                                                gTAGGGCGTGGTCATGAACAAATATTCACGCGCCGAATAAGCCGCAGTCAATAATGCCTGGTGAATGAGATCTTCCGGAAAGCCGGGGCCAGAAGCAAtt                 >  2:438568/1‑100 (MQ=255)                                                                                gTAGGGCGTGGTCATGAACAAATATTCACGCGCCGAATAAGCCGCAGTCAATAATGCCTGGTGAATGAGATCTTCCGGAAAGCCGGGGCCAGAAGCAAtt                 <  2:1493296/100‑1 (MQ=255)                                                                                 tAGGGCGTGGTCATGAACAAATATTCACGCGCCGAATAAGCCGCAGTCAATAATGCCTGGTGAATGAGATCTTCCGGAAAGCCGGGGCCAGAAGCAATtg                <  1:1373252/100‑1 (MQ=255)                                                                                   gggCGTGGTCATGAACAAATATTCACGCGCCGAATAAGCCGCAGTCAATAATGCCTGGTGAATGAGATCTTCCGGAAAGCCGGGGCCAGAAGCAATtgtg              <  1:703288/100‑1 (MQ=255)                                                                                     gCGTGGTCATGAACACATATTCACGCGCCGAATAAGCCGCAGGCAATAATGACTGGTGAATGAGATCGTCGGTAAAGCTGGGGGCAGGAGCAATtgtggg            >  1:645945/1‑98 (MQ=255)                                                                                     gCGTGGTCATGAACAAATATTCACGCGCCGAATAAGCCGCAGTCAATAATGCCTGGTGAATGAGATCTTCCGGAAAGCCGGGGCCAGAAGCAATtgtgtg            <  2:211611/100‑1 (MQ=255)                                                                                     gCGTGGTCATGAACAAATATTCACGCGCCGAATAAGCCGCAGTCAATAATGCCTGGTGAATGAGATCTTCCGGAAAGCCGGGGCCAGAAGCAATtgtgtg            <  1:134280/100‑1 (MQ=255)                                                                                      cGTGGTCATGAACAAATATTCACGCGCCGAATAAGCCGCAGTCAATAATGCCTGGTGAATGAGATCTTCCGGAAAGCCGGGGCCAGAAGCAATTGTGTGa           <  1:235953/100‑1 (MQ=255)                                                                                         ggTCATGAACAAATATTCACGCGCCGAATAAGCCGCAGTCAATAATGCCTGGTGAATGAGATCTTCCGGAAAGCCGGGGCCAGAAGCAATTGTGTGAATg        <  2:1325325/100‑1 (MQ=255)                                                                                         ggTCATGAACAAATATTCACGCGCCGAATAAGCCGCAGTCAATAATGCCTGGTGAATGAGATCTTCCGGAAAGCCGGGGCCAGAAGCAATTGTGTGAATg        >  2:927919/1‑100 (MQ=255)                                                                                          gTCATGAACAAATATTCACGCGCCGAATAAGCCGCAGTCAATAATGCCTGGTGAATGAGATCTTCCGGAAAGCCGGGGCCAGAAGCAATTGTGTGAATgg       <  2:850003/100‑1 (MQ=255)                                                                                            cATGAACAAATATTCACGCGCCGAATAAGCCGCAGTCAATAATGCCTGGTGAATGAGATCTTCCGGAAAGCCGGGGCCAGAAGCAATTGTGTGAATGgtg     <  2:1891106/100‑1 (MQ=255)                                                                                               gATCAAATATTCACGCGCCGAATAAGCCGCAGTCAATAATGCCTGGTGAATGAGATCTTCCGGAAAGCCGGGGCCAGAAGCAATTGTGTGAATGGTGTGa  >  1:671437/1‑100 (MQ=255)                                                                                                 |                                                                                                  GCGGAAGGATAATACTGACATCCACCCCGCGCTGCGCCGCCGTGCAAATCGCATGAAGTAAATCATCGCTTGGCACAAAGTAGGGCGTGGTCATGATCAAATATTCACGCGCCGAATAAGCCGCAGTCAATAATGCCTGGTGAATGAGATCTTCCGGAAAGCCGGGGCCAGAAGCAATTGTGTGAATGGTGTGA  >  NC_000913/1307538‑1307731

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