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%	246.9 / NA	76	*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 (40/36);  total (40/36)

TTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTACTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAGAACTGAGTACGCAGCAGAGGAATCAGATCAGGTTGACGCGCTGCGCCTTCC  >  NC_000913/3920821‑3921082                                                                                                                                  |                                                                                                                                     ttCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACa                                                                                                                             >  2:730760/1‑139 (MQ=255)     tGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGa                                                                                                                         >  2:1169082/1‑139 (MQ=255)      gCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTCAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTAGGCGACAGCGaa                                                                                                                        >  2:1061560/1‑139 (MQ=255)           acagaacGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACACTTATTACGCGACAGCTAACATCa                                                                                                                   >  1:468745/1‑139 (MQ=255)                     cAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTAcagacc                                                                                                         <  1:1169082/139‑1 (MQ=255)                        gACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTAcagacccag                                                                                                      >  1:918756/1‑139 (MQ=255)                         aCAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTAcgcg                                                                                                                                <  1:158517/112‑1 (MQ=255)                         aCAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTAcgcg                                                                                                                                >  2:158517/1‑112 (MQ=255)                         aCAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTAcagacccaga                                                                                                     >  2:71649/1‑139 (MQ=255)                               gcgATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTAca                                                                                                             >  2:223945/1‑125 (MQ=255)                               gcgATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTAca                                                                                                             <  1:223945/125‑1 (MQ=255)                                  aTGCCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAg                                                                                              <  2:332336/137‑1 (MQ=255)                                  aTGCCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAg                                                                                              >  1:332336/1‑137 (MQ=255)                                       ctggccGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGatcat                                                                                       <  1:71649/139‑1 (MQ=255)                                         ggccGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTAcagacc                                                                                                         <  2:1211009/119‑1 (MQ=255)                                         ggccGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTAcagacc                                                                                                         >  1:1211009/1‑119 (MQ=255)                                            cGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGa                                                                                  <  1:730760/139‑1 (MQ=255)                                               ggATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGa                                                                                                                         <  2:632523/97‑1 (MQ=255)                                               ggATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGa                                                                                                                         >  1:632523/1‑97 (MQ=255)                                               ggATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAg                                                                               >  2:1325870/1‑139 (MQ=255)                                                    gttgCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCg                                                                                     <  1:212741/128‑1 (MQ=255)                                                    gttgCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCg                                                                                     >  2:212741/1‑128 (MQ=255)                                                                 tCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATaa                                                             <  1:1325870/139‑1 (MQ=255)                                                                         caATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCg                                                                                     <  2:328958/107‑1 (MQ=255)                                                                         caATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCg                                                                                     >  1:328958/1‑107 (MQ=255)                                                                          aaTGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCa                                                                                                                   <  1:789541/76‑1 (MQ=255)                                                                          aaTGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCa                                                                                                                   >  2:789541/1‑76 (MQ=255)                                                                                 ctATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAGAACt                                                 >  2:209568/1‑135 (MQ=255)                                                                                 ctATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAGAACt                                                 <  1:209568/135‑1 (MQ=255)                                                                                 ctATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAGAACTGAGt                                             >  2:1296417/1‑139 (MQ=255)                                                                                  tATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAGAACTGAGTa                                            <  1:1004200/139‑1 (MQ=255)                                                                                        gTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAAc                                                                                                                       >  2:1211095/1‑58 (MQ=255)                                                                                        gTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAAc                                                                                                                       <  1:1211095/58‑1 (MQ=255)                                                                                         ttAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTAcagacc                                                                                                         <  2:271417/71‑1 (MQ=255)                                                                                         ttAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTAcagacc                                                                                                         >  1:271417/1‑71 (MQ=255)                                                                                          tAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAGAACTGAGTACGcagcag                                    >  1:319018/1‑139 (MQ=255)                                                                                            aCGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTa                                                                                                               >  2:233812/1‑62 (MQ=255)                                                                                            aCGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTa                                                                                                               <  1:233812/62‑1 (MQ=255)                                                                                            aCGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGc                                                                              >  1:809979/1‑95 (MQ=255)                                                                                            aCGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGc                                                                              <  2:809979/95‑1 (MQ=255)                                                                                            aCGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATaa                                                             >  1:552352/1‑112 (MQ=255)                                                                                            aCGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATaa                                                             <  2:552352/112‑1 (MQ=255)                                                                                              gTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAAc                                                                                                                       <  1:991625/52‑1 (MQ=255)                                                                                              gTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAAc                                                                                                                       >  2:991625/1‑52 (MQ=255)                                                                                                   gATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTAcagacc                                                                                                         <  2:989590/61‑1 (MQ=255)                                                                                                   gATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTAcagacc                                                                                                         >  1:989590/1‑61 (MQ=255)                                                                                                   gATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTAc                                                                                                              >  2:1076804/1‑56 (MQ=255)                                                                                                   gATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTAc                                                                                                              <  1:1076804/56‑1 (MQ=255)                                                                                                       ttGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCa                                                                >  1:1127301/1‑98 (MQ=255)                                                                                                       ttGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCa                                                                <  2:1127301/98‑1 (MQ=255)                                                                                                        tGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAGAACTGAGTACGCAGCAGAGGAat                              >  2:439843/1‑131 (MQ=255)                                                                                                        tGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAGAACTGAGTACGCAGCAGAGGAat                              <  1:439843/131‑1 (MQ=255)                                                                                                         gCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAGAACTGAGTACGCAGCAGAGGAATCAGATCAgg                     >  2:259560/1‑139 (MQ=255)                                                                                                          ctctTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATaa                                                             <  2:1240852/98‑1 (MQ=255)                                                                                                          ctctTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATaa                                                             >  1:1240852/1‑98 (MQ=255)                                                                                                          ctctTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAGAAc                                                  <  2:1272803/109‑1 (MQ=255)                                                                                                          ctctTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAGAAc                                                  >  1:1272803/1‑109 (MQ=255)                                                                                                          ctctTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAGAACTGAGTACGCAGCAGAGGAATCAGATCAGGt                    >  2:1066224/1‑139 (MQ=255)                                                                                                            ctTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCTAGACCTACAGCGATCATCGGGATAGCa                                                                             <  2:802854/80‑1 (MQ=255)                                                                                                            ctTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCa                                                                             >  1:802854/1‑80 (MQ=255)                                                                                                            ctTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAGAACTGAGTACGCAGCAGAGGAATCAGATCAGGTTg                  <  1:1066224/139‑1 (MQ=255)                                                                                                             tttAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAGAACTGAGTACGCAGCAGAGGAATCAGATCAGGTTGa                 <  2:319018/139‑1 (MQ=255)                                                                                                               taaataaaAGCAACGCTTATTACGCGACAGCGAACATCACGTAcagacc                                                                                                         <  1:884455/49‑1 (MQ=255)                                                                                                               taaataaaAGCAACGCTTATTACGCGACAGCGAACATCACGTAcagacc                                                                                                         >  2:884455/1‑49 (MQ=255)                                                                                                               taaataaaAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATaa                                                             <  2:731698/93‑1 (MQ=255)                                                                                                               taaataaaAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATaa                                                             >  1:731698/1‑93 (MQ=255)                                                                                                               taaataaaAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAGa                                                    >  2:660308/1‑102 (MQ=255)                                                                                                               taaataaaAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAGa                                                    <  1:660308/102‑1 (MQ=255)                                                                                                               taaataaaAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAGAACTGAGTACGCAGCAGAGGAATCAGATCAGGTTGAcg               >  1:1170121/1‑139 (MQ=255)                                                                                                                      aaGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCTCCAGACCCATAACGATAAAGAACTGAGTACGCAGCAGAGGAATCAGATCAGGTTGACGCGCTgcg        <  2:918756/139‑1 (MQ=255)                                                                                                                      aaGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAGAACTGAGTACGCAGCAGAGGAATCAGATCAGGTTGACGCGCTgcg        <  1:259560/139‑1 (MQ=255)                                                                                                                       aGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGa                                                                                  <  2:606203/64‑1 (MQ=255)                                                                                                                       aGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGa                                                                                  >  1:606203/1‑64 (MQ=255)                                                                                                                       aGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCa                                                                             >  2:572068/1‑69 (MQ=255)                                                                                                                       aGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCa                                                                             <  1:572068/69‑1 (MQ=255)                                                                                                                            cGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAGAACTGAGCACGCAGCAGAGGAATCAGATCAGGTTGACGCGCTGCGCCTTcc  >  1:432861/1‑139 (MQ=255)                                                                                                                                  |                                                                                                                                     TTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTACTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAGAACTGAGTACGCAGCAGAGGAATCAGATCAGGTTGACGCGCTGCGCCTTCC  >  NC_000913/3920821‑3921082

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