Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,920,950	C→T	*80* (TAG→TAA)	atpE ←	F0 sector of membrane‑bound ATP synthase, subunit c

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,920,950	0	C	T	100.0%	157.9 / NA	49	*80* (TAG→TAA)	atpE	F0 sector of membrane‑bound ATP synthase, subunit c
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (22/27);  total (22/27)

CATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTACTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAGAACTGAGTACGCAGCAGAGGAATCAGATCAGGTTGACGCGCTGCGCCTTCC  >  NC_000913/3920823‑3921082                                                                                                                                |                                                                                                                                     cATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGc                                                                                                                           >  2:773921/1‑139 (MQ=255)   tGCAGAACAGAACGAACAGGACAAACGCGATGGTCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGa                                                                                                                         >  1:296052/1‑139 (MQ=255)   tGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAATTAAAAGCAACGCTTATTACGCGACAGCGa                                                                                                                         >  2:599956/1‑139 (MQ=255)   tGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGa                                                                                                                         >  1:120983/1‑139 (MQ=255)   tGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGa                                                                                                                         <  2:120983/139‑1 (MQ=255)   tGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGa                                                                                                                         <  1:825682/139‑1 (MQ=255)     cagagtagaacTAACAGAACAAACGCAATGGCCTCCCCGAGGATTGTTTCGGTAAGATTTACAGCACTATGCCTCTAATGAGTTAGCGTTCTGGTATTGCGCTTTTAATAAAAGCAACGCTTATTACGCGACAGCAAAc                                                                                                                       <  2:69118/133‑1 (MQ=255)             aacgaacAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTa                                                                                                               >  2:99225/1‑139 (MQ=255)                   cAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTAcagacc                                                                                                         <  2:654365/139‑1 (MQ=255)                   cAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTAcagacc                                                                                                         <  1:773921/139‑1 (MQ=255)                      gACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTAcagacccag                                                                                                      <  2:270717/139‑1 (MQ=255)                                aTGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCg                                                                                            <  1:413483/139‑1 (MQ=255)                                              gATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTAcgc                                                                                                                                 <  1:431574/88‑1 (MQ=255)                                              gATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTAcgc                                                                                                                                 >  2:431574/1‑88 (MQ=255)                                                  gttgCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATcc                                                                          >  2:750182/1‑139 (MQ=255)                                                     gCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATccacc                                                                       <  2:467861/139‑1 (MQ=255)                                                               tCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATaa                                                             <  2:296052/139‑1 (MQ=255)                                                                          tGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTAc                                                                                                              <  2:599231/79‑1 (MQ=255)                                                                          tGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTAc                                                                                                              >  1:599231/1‑79 (MQ=255)                                                                          tGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATccac                                                                        <  1:858477/117‑1 (MQ=255)                                                                          tGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATccac                                                                        >  2:858477/1‑117 (MQ=255)                                                                          tGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAGAAc                                                  >  1:223550/1‑139 (MQ=255)                                                                          tGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTCATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCg                                                                                            <  1:842609/97‑1 (MQ=255)                                                                          tGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTCATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCg                                                                                            >  2:842609/1‑97 (MQ=255)                                                                               ctATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGa                                                                                  >  2:624435/1‑102 (MQ=255)                                                                               ctATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGa                                                                                  <  1:624435/102‑1 (MQ=255)                                                                                       ttACCGTTCTGATATTGCTCTTTAAGTAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAGAACTGAGTACgcagca                                     <  1:599956/139‑1 (MQ=255)                                                                                       ttAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATccacc                                                                       >  1:750359/1‑105 (MQ=255)                                                                                       ttAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATccacc                                                                       <  2:750359/105‑1 (MQ=255)                                                                                       ttAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATa                                                        <  1:549620/120‑1 (MQ=255)                                                                                       ttAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCACCCACCAGACCAATAACGATa                                                        >  2:549620/1‑120 (MQ=255)                                                                                           cGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAGAACTGAGTACGCAGCAGAgg                                 <  1:330049/139‑1 (MQ=255)                                                                                             ttCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAGAACTGAGTACGCAGCAGAGGaa                               <  2:686715/139‑1 (MQ=255)                                                                                               cTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTAc                                                                                                              <  2:367900/58‑1 (MQ=255)                                                                                               cTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTAc                                                                                                              >  1:367900/1‑58 (MQ=255)                                                                                                 gATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAGAACTGAGTACg                                          <  1:161627/124‑1 (MQ=255)                                                                                                 gATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAGAACTGAGTACg                                          >  2:161627/1‑124 (MQ=255)                                                                                                        ctctTTAAATAAAAGCAACGCTTATTACGCGACAGCGa                                                                                                                         <  2:389258/38‑1 (MQ=255)                                                                                                        ctctTTAAATAAAAGCAACGCTTATTACGCGACAGCGa                                                                                                                         >  1:389258/1‑38 (MQ=255)                                                                                                        ctctTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATc                                                                                                                    >  1:919523/1‑43 (MQ=255)                                                                                                        ctctTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATc                                                                                                                    <  2:919523/43‑1 (MQ=255)                                                                                                        ctctTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAGAACTGAGTACGCAGCAGAGGAATCAGATCAGGt                    >  1:765814/1‑139 (MQ=255)                                                                                                                     aGCAACGCTTATTACGCGACAGCGAACATCACGTAcagaccca                                                                                                       <  2:341069/43‑1 (MQ=255)                                                                                                                     aGCAACGCTTATTACGCGACAGCGAACATCACGTAcagaccca                                                                                                       >  1:341069/1‑43 (MQ=255)                                                                                                                     aGCAACGCTTATTACGCGACAGCGAACATCACGTAcagacc                                                                                                         <  1:572145/41‑1 (MQ=255)                                                                                                                     aGCAACGCTTATTACGCGACAGCGAACATCACGTAcagacc                                                                                                         >  2:572145/1‑41 (MQ=255)                                                                                                                          cGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAGAACTGAGTACgcagca                                     >  2:152209/1‑104 (MQ=255)                                                                                                                          cGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAGAACTGAGTACgcagca                                     <  1:152209/104‑1 (MQ=255)                                                                                                                          cGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAGAACTGAGTACGCAGCAGAGGAATCAGATCAGGTTGACGCGCTGCGCCTTcc  <  1:750182/139‑1 (MQ=255)                                                                                                                                |                                                                                                                                     CATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTACTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAGAACTGAGTACGCAGCAGAGGAATCAGATCAGGTTGACGCGCTGCGCCTTCC  >  NC_000913/3920823‑3921082

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑