Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,308,318	G→C	100%	R110G (CGT→GGT)	clsA ←	cardiolipin synthase 1

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,308,318	0	G	C	96.9%	225.7 / ‑0.2	64	R110G (CGT→GGT)	clsA	cardiolipin synthase 1
Reads supporting (aligned to +/- strand):  ref base G (2/0);  new base C (35/27);  total (37/27)
Fisher's exact test for biased strand distribution p-value = 5.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.65e-01

GCGCGAGCTGGATGTCGCGGATTAACGCCTGCATCACATCATCTGACTCGGTCATCAGTTGTAGCTGATTCCCTTTGACCCCAGCGATCCCCTGACGACGCTCGCAAAGCTTGAATAATGGCGCAGCGACACTGCTATTTTCTTCGGCGAAGATATGCTTACAGGCTTTAAGGTCGTTAAGCCATTTTGCG  >  NC_000913/1308222‑1308412                                                                                                 |                                                                                               gcgcgAGCTGGATGTCGCGGATTAACGCCTGCATCACATCATCTGACTCGGTCATCAGTTGTAGCTGATTCCCTTTGACCCCAGCGATCCCCTGaccacg                                                                                             <  2:1741016/100‑1 (MQ=255)   gcgAGCTGGATGTCGCGGATTAACGCCTGCATCACATCATCTGACTCGGTCATCAGTTGTAGCTGATTCCCTTTGACCCCAGCGATCCCCTGACCACGCt                                                                                           <  1:1019498/100‑1 (MQ=255)      aGCTGGATGTCGCGGATTAACGCCTGCATCACATCATCTGACTCGGTCATCAGTTGTAGCTGATTCCCTTTGACCCCAGCGATCCCCTGACCACGCTCGc                                                                                        <  1:1675553/100‑1 (MQ=255)       gCTGGATGTCGCGGATTAACGCCTGCATCACATCATCTGACTCGGTCATCAGTTGTAGCTGATTCCCTTTGACCCCAGCGATCCCCTGACCACGCTCGCa                                                                                       >  1:545965/1‑100 (MQ=255)        cTGGATGTCGCGGATTAACGCCTGCATCACATCATCTGACTCGGTCATCAGTTGTAGCTGATTCCCTTTGACCCCAGCGATCCCCTGACCACGCTCGCaa                                                                                      >  2:1960960/1‑100 (MQ=255)           gATGTCGCGGATTAACGCCTGCATCACATCATCTGACTCGGTCATCAGTTGTAGCTGATTCCCTTTGACCCCAGCGATCCCCTGACCACGCTCGCAAAGc                                                                                   <  1:190965/100‑1 (MQ=255)              gTCGCGGATTAACGCCTGCATCACATCATCTGACTCGGTCATCAGTTGTAGCTGATTCCCTTTGACCCCAGCGATCCCCTGACCACGCTCGCAAAGCTTg                                                                                <  2:326762/100‑1 (MQ=255)              gTCGCGGATTAACGCCTGCATCACATCATCTGACTCGGTCATCAGTTGTAGCTGATTCCCTTTGACCCCAGCGATCCCCTGACCACGCTCGCAAAGCTTg                                                                                >  1:332710/1‑100 (MQ=255)               tCGCGGATTAACGCCTGCATCACATCATCTGACTCGGTCATCAGTTGTAGCTGATTCCCTTTGACCCCAGCGATCCCCTGACCACGCTCGCAAAGCTTGa                                                                               <  2:389320/100‑1 (MQ=255)               tCGCGGATTAACGCCTGCATCACATCATCTGACTCGGTCATCAGTTGTAGCTGATTCCCTTTGACCCCAGCGATCCCCTGACCACGCTCGCAAAGCTTGa                                                                               <  2:139475/100‑1 (MQ=255)               tCGCGGATTAACGCCTGCATCACATCATCTGACTCGGTCATCAGTTGTAGCTGATTCCCTTTGACCCCAGCGATCCCCTGACCACGCTCGCAAAGCTTGa                                                                               >  1:1958985/1‑100 (MQ=255)                cgcgGATTAACGCCTGCATCACATCATCTGACTCGGTCATCAGTTGTAGCTGATTCCCTTTGACCCCAGCGATCCCCTGACCACGCTCGCAAAGCTTGaa                                                                              >  1:1813632/1‑100 (MQ=255)                   ggATTAACGCCTGCATCACATCATCTGACTCGGTCATCAGTTGTAGCTGATTCCCTTTGACCCCAGCGATCCCCTGACCACGCTCGCAAAGCTTGaataa                                                                           >  2:2053135/1‑100 (MQ=255)                      ttAACGCCTGCATCACATCATCTGACTCGGTCATCAGTTGTAGCTGATTCCCTTTGACCCCAGCGATCCCCTGACCACGCTCGCAAAGCTTGAATAATgg                                                                        <  2:684749/100‑1 (MQ=255)                        aaCGCCTGCATCACATCATCTGACTCGGTCATCAGTTGTAGCTGATTCCCTTTGACCCCAGCGATCCCCTGACCACGCTCGCAAAGCTTGAATAATGgcg                                                                      >  2:925057/1‑100 (MQ=255)                        aaCGCCTGCATCACATCATCTGACTCGGTCATCAGTTGTAGCTGATTCCCTTTGACCCCAGCGATCCCCTGACCACGCTCGCAAAGCTTGAATAATGgcg                                                                      >  2:924854/1‑100 (MQ=255)                            ccTGCATCACATCATCTGACTCGGTCATCAGTTGTAGCTGATTCCCTTTGACCCCAGCGATCCCCTGACCACGCTCGCAAAGCTTGAATAATGGCGCAGc                                                                  >  1:834835/1‑100 (MQ=255)                            ccTGCATCACATCATCTGACTCGGTCATCAGTTGTAGCTGATTCCCTTTGACCCCAGCGATCCCCTGACCACGCTCGCAAAGCTTGAATAATGGCGCAGc                                                                  <  1:288381/100‑1 (MQ=255)                               gCATCACATCATCTGACTCGGTCATCAGTTGTAGCTGATTCCCTTTGACCCCAGCGATCCCCTGACCACGCTCGCAAAGCTTGAATAATGGCGCAGCGac                                                               >  2:740909/1‑100 (MQ=255)                               gCATCACATCATCTGACTCGGTCATCAGTTGTAGCTGATTCCCTTTGACCCCAGCGATCCCCTGACCACGCTCGCAAAGCTTGAATAATGGCGCAGCGac                                                               >  2:1143619/1‑100 (MQ=255)                                   cacatcaTCTGACTCGGTCATCAGTTGTAGCTGATTCCCTTTGACCCCAGCGATCCCCTGACCACGCTCGCAAAGCTTGAATAATGGCGCAGCGACACTg                                                           >  1:134448/1‑100 (MQ=255)                                     catcatCTGACTCGGTCATCAGTTGTAGCTGATTCCCTTTGACCCCAGCGATCCCCTGACCACGCTCGCAAAGCTTGAATAATGGCGCAGCGACACTGCt                                                         >  2:702543/1‑100 (MQ=255)                                        catcTGACTCGGTCATCAGTTGTAGCTGATTCCCTTTGACCCCAGCGATCCCCTGACCACGCTCGCAAAGCTTGAATAATGGCGCAGCGACACTGCTAtt                                                      >  2:527128/1‑100 (MQ=255)                                            tGACTCGGTCATCAGTTGTAGCTGATTCCCTTTGACCCCAGCGATCCCCTGACCACGCTCGCAAAGCTTGAATAATGGCGCAGCGACACTGCTATTttct                                                  <  1:426667/100‑1 (MQ=255)                                             gACTCGGTCATCAGTTGTAGCTGATTCCCTTTGACCCCAGCGATCCCCTGACCACGCTCGCAAAGCTTGAATAATGGCGCAGCGACACTGCTATTttctt                                                 <  2:1021019/100‑1 (MQ=255)                                               cTCGGTCATCAGTTGTAGCTGATTCCCTTTGACCCCAGCGATCCCCTGACCACGCTCGCAAAGCTTGAATAATGGCGCAGCGACACTGCTATTTTCTTCg                                               >  1:193834/1‑100 (MQ=255)                                               cTCGGTCATCAGTTGTAGCTGATTCCCTTTGACCCCAGCGATCCCCTGACCACGCTCGCAAAGCTTGAATAATGGCGCAGCGACACTGCTATTTTCTTCg                                               >  1:1853437/1‑100 (MQ=255)                                                  ggTCATCAGTTGTAGCTGATTCCCTTTGACCCCAGCGATCCCCTGACCACGCTCGCAAAGCTTGAATAATGGCGCAGCGACACTGCTATTTTCTTCGGCg                                            >  1:247276/1‑100 (MQ=255)                                                  ggTCATCAGTTGTAGCTGATTCCCTTTGACCCCAGCGATCCCCTGACCACGCTCGCAAAGCTTGAATAATGGCGCAGCGACACTGCTATTTTCTTCGGCg                                            <  2:1853423/100‑1 (MQ=255)                                                   gTCATCAGTTGTAGCTGATTCCCTTTGACCCCAGCGATCCCCTGACCACGCTCGCAAAGCTTGAATAATGGCGCAGCGACACTGCTATTTTCTTCGGCGa                                           >  1:76136/1‑100 (MQ=255)                                                     catcaGTTGTAGCTGATTCCCTTTGACCCCAGCGATCCCCTGACCACGCTCGCAAAGCTTGAATAATGGCGCAGCGACACTGCTATTTTCTTCGGCGAAg                                         >  2:650160/1‑100 (MQ=255)                                                      atcaGTTGTAGCTGATTCCCTTTGACCCCAGCGATCCCCTGACCACGCTCGCAAAGCTTGAATAATGGCGCAGCGACACTGCTATTTTCTTCGGCGAAGa                                        <  1:2008664/100‑1 (MQ=255)                                                       tcaGTTGTAGCTGATTCCCTTTGACCCCAGCGATCCCCTGACCACGCTCGCAAAGCTTGAATAATGGCGCAGCGACACTGCTATTTTCTTCGGCGAAGat                                       >  2:86196/1‑100 (MQ=255)                                                            tGTAGCTGATTCCCTTTGACCCCAGCGATCCCCTGACCACGCTCGCAAAGCTTGAATAATGGCGCAGCGACCCTGCTATTTTCTTCGGCGAAGATATGc                                   >  2:85991/1‑99 (MQ=255)                                                            tGTAGCTGATTCCCTTTGACCCCAGCGATCCCCTGACCACGCTCGCAAAGCTTGAATAATGGCGCAGCGACACTGCTATTTTCTTCGGCGAAGATATGCt                                  <  1:1673308/100‑1 (MQ=255)                                                            tGTAGCTGATTCCCTTTGACCCCAGCGATCCCCTGACCACGCTCGCAAAGCTTGAATAATGGCGCAGCGACACTGCTATTTTCTTCGGCGAAGATATGCt                                  <  2:693451/100‑1 (MQ=255)                                                             gTATCTGATTCCCTTTGACCCCAGCGATCCCCTGACCACGCTCGCAAAGCTTGAATAATGGCGCAGCGACACTGCTATTTTCTTCGGCGAAGATATGCtt                                 >  2:283708/1‑100 (MQ=255)                                                               aGCTGATTCCCTTTGACCCCAGCGATCCCCTGACGACGCTCGCAAAGCTTGAATAATGGCGCAGCGACACTGCTATTTTCTTCGGCGAAGATATGCTTAc                               >  1:425939/1‑100 (MQ=255)                                                                   gATTCCCTTTGACCCCAGCGATCCCCTGACCACGCTCGCAAAGCTTGAATAATGGCGCAGCGACACTGCTATTTTCTTCGGCGAAGATATGCTTACAGGc                           <  1:1950639/100‑1 (MQ=255)                                                                   gATTCCCTTTGACCCCAGCGATCCCCTGACCACGCTCGCAAAGCTTGAATAATGGCGCAGCGACACTGCTATTTTCTTCGGCGAAGATATGCTTACAGGc                           <  1:603982/100‑1 (MQ=255)                                                                   gATTCCCTTTGACCCCAGCGATCCCCTGACCACGCTCGCAAAGCTTGAATAATGGCGCAGCGACACTGCTATTTTCTTCGGCGAAGATATGCTTACAGGc                           <  2:1531664/100‑1 (MQ=255)                                                                    aTTCCCTTTGACCCCAGCGATCCCCTGACCACGCTCGCAAAGCTTGAATAATGGCGCAGCGACACTGCTATTTTCTTCGGCGAAGATATGCTTACAGGCt                          <  2:869708/100‑1 (MQ=255)                                                                     ttCCCTTTGACCCCAGCGATCCCCTGACCACGCTCGCAAAGCTTGAATAATGGCGCAGCGACACTGCTATTTTCTTCGGCGAAGATATGCTTACAGGCtt                         <  1:747754/100‑1 (MQ=255)                                                                     ttCCCTTTGACCCCAGCGATCCCCTGACCACGCTCGCAAAGCTTGAATAATGGCGCAGCGACACTGCTATTTTCTTCGGCGAAGATATGCTTACAGGCtt                         >  2:1639892/1‑100 (MQ=255)                                                                       cccTTTGACCCCAGCGATCCCCTGACCACGCTCGCAAAGCTTGAATAATGGCGCAGCGACACTGCTATTTTCTTCGGCGAAGATATGCTTACAGGCTTTa                       <  2:276372/100‑1 (MQ=255)                                                                             gACCCCAGCGATCCCCTGACCACGCTCGCAAAGCTTGAATAATGGCGCAGCGACACTGCTATTTTCTTCGGCGCAGATATGCTTACAGGCTTTAAGGTCg                 >  1:676937/1‑100 (MQ=255)                                                                              aCCCCAGCGATCCCCTGACCACGCTCGCAAAGCTTGAATAATGGCGCAGCGACACTGCTATTTTCTTCGGCGAAGATATGCTTACAGGCTTTAAGGt                   <  2:1958969/97‑1 (MQ=255)                                                                              aCCCCAGCGATCCCCTGACCACGCTCGCAAAGCTTGAATAATGGCGCAGCGACACTGCTATTTTCTTCGGCGAAGATATGCTTACAGGCTTTAAGGTCGt                <  2:103363/100‑1 (MQ=255)                                                                               ccccAGCGATCCCCTGACCACGCTCGCAAAGCTTGAATAATGGCGCAGCGACACTGCTATTTTCTTCGGCGAAGATATGCTTACAGGCTTTAAGGTCGtt               >  2:982430/1‑100 (MQ=255)                                                                               ccccAGCGATCCCCTGACCACGCTCGCAAAGCTTGAATAATGGCGCAGCGACACTGCTATTTTCTTCGGCGAAGATATGCTTACAGGCTTTAAGGTCGtt               >  1:875625/1‑100 (MQ=255)                                                                               ccccAGCGATCCCCTGACCACGCTCGCAAAGCTTGAATAATGGCGCAGCGACACTGCTATTTTCTTCGGCGAAGATATGCTTACAGGCTTTAAGGTCGtt               >  2:492661/1‑100 (MQ=255)                                                                               ccccAGCGATCCCCTGACCACGCTCGCAAAGCTTGAATAATGGCGCAGCGACACTGCTATTTTCTTCGGCGAAGATATGCTTACAGGCTTTAAGGTCGtt               >  1:678064/1‑100 (MQ=255)                                                                               ccccAGCGATCCCCTGACCACGCTCGCAAAGCTTGAATAATGGCGCAGCGACACTGCTATTTTCTTCGGCGAAGATATGCTTACAGGCTTTAAGGTCGtt               >  1:1981403/1‑100 (MQ=255)                                                                                 ccAGCGATCCCCTGACCACGCTCGCAAAGCTTGAATAATGGCGCAGCGACACTGCTATTTTCTTCGGCGAAGATATGCTTACAGGCTTTagatcggaaga             >  1:224230/1‑90 (MQ=255)                                                                                  cAGCGATCCCCTGACCACGCTCGCAAAGCTTGAATAATGGCGCAGCGACACTGCTATTTTCTTCGGCGAAGATATGCTTACAGGCTTTAAGGTCGTTAAg            <  1:216727/100‑1 (MQ=255)                                                                                  cAGCGATCCCCTGACCACGCTCGCAAAGCTTGAATAATGGCGCAGCGACACTGCTATTTTCTTCGGCGAAGATATGCTTACAGGCTTTAAGGTCGTTAAg            >  1:1976663/1‑100 (MQ=255)                                                                                  cAGCGATCCCCTGACCACGCTCGCAAAGCTTGAATAATGGCGCAGCGACACTGCTATTTTCTTCGGCGAAGATATGCTTACAGGCTTTAAGGTCGTTAAg            <  1:1861202/100‑1 (MQ=255)                                                                                   aGCGATCCCCTGACCACGCTCGCAAAGCTTGAATAATGGCGCAGCGACACTGCTATTTTCTTCGGCGAAGATATGCTTACAGGCTTTAAGGTCGTTAAGc           >  2:766989/1‑100 (MQ=255)                                                                                   aGCGATCCCCTGACCACGCTCGCAAAGCTTGAATAAAGGCGCAGTGACACTGCTATTTTCTTCGGCGAAGATATGCTTACAGGCTTTAAGGTCGTTAAGc           <  1:1973199/100‑1 (MQ=255)                                                                                    gCGATCCCCTGACCACGCTCGCAAAGCTTGAATAATGGCGCAGCGACACTGCTATTTTCTTCGGCGAAGATATGCTTACAGGCTTTAAGGTCGTTAAGcn          <  1:2048113/100‑2 (MQ=255)                                                                                    gCGATCCCCTGACCACGCTCGCAAAGCTTGAATAATGGCGCAGCGACACTGCTATTTTCTTCGGCGAAGATATGCTTACAGGCTTTAAGGTCGTTAAGcc          <  1:1604744/100‑1 (MQ=255)                                                                                      gATCCCCTGACCACGCTCGCAAAGCTTGAATAATGGCGCAGCGACACTGCTATTTTCTTCGGCGAAGATATGCTTACAGGCTTTAAGGTCGTTAAGCCAt        >  2:647240/1‑100 (MQ=255)                                                                                          cccTGACCACGCTCGCAAAGCTTGAATAATGGCGCAGCGACACTGCTATTTTCTTCGGCGTAGATATACTTACAGGCTTTAAGGTCGTTAAGCCATTTTg    >  2:1455852/1‑100 (MQ=255)                                                                                          cccTGACCACGCTCGCAAAGCTTGAATAATGGCGCAGCGACACTGCTATTTTCTTCGGCGAAGATATGCTTACAGGCTTTAAGGTCGTTAAGCCATTTTg    >  1:1106234/1‑100 (MQ=255)                                                                                            cTGACGACGCTCGCAAAGCTTGAATAATGGCGCAGCGAAACTGCTATTTTCTTCGGCGAAGATATGCTTACAGGCTTTAAGGTCGTTAAGCCATTTTGCg  >  2:135955/1‑100 (MQ=255)                                                                                                 |                                                                                               GCGCGAGCTGGATGTCGCGGATTAACGCCTGCATCACATCATCTGACTCGGTCATCAGTTGTAGCTGATTCCCTTTGACCCCAGCGATCCCCTGACGACGCTCGCAAAGCTTGAATAATGGCGCAGCGACACTGCTATTTTCTTCGGCGAAGATATGCTTACAGGCTTTAAGGTCGTTAAGCCATTTTGCG  >  NC_000913/1308222‑1308412

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: ‑