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%	194.4 / NA	59	*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 (30/29);  total (30/29)

CAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTACTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAGAACTGAGTACGCAGCAGAGGAATCAGATCAGGTTGACGCGCTGCGCCTTC  >  NC_000913/3920827‑3921081                                                                                                                            |                                                                                                                                    cagaacagaaCGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAAc                                                                                                                      >  1:1275649/1‑139 (MQ=255)   gaacagaacGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACAt                                                                                                                    >  1:990813/1‑139 (MQ=255)                   aCAAACGCGATGGCCTGGACGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTAcagacccaga                                                                                                    >  2:838328/1‑139 (MQ=255)                       aCGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTAcg                                                                                                                                 >  2:1172676/1‑106 (MQ=255)                       aCGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTAcg                                                                                                                                 <  1:1172676/106‑1 (MQ=255)                            aTGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCg                                                                                           <  1:375816/139‑1 (MQ=255)                                 ctggccGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGatcat                                                                                      <  2:1275649/139‑1 (MQ=255)                                               ttgCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTAcagacccag                                                                                                     >  2:1140277/1‑110 (MQ=255)                                               ttgCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTAcagacccag                                                                                                     <  1:1140277/110‑1 (MQ=255)                                                tgCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATccac                                                                       <  2:990813/139‑1 (MQ=255)                                                tgCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATccac                                                                       <  1:838328/139‑1 (MQ=255)                                                 gCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATccacc                                                                      <  1:458566/139‑1 (MQ=255)                                                  cGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATccacca                                                                     <  1:1036010/139‑1 (MQ=255)                                                                 cacaATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATaa                                                      >  1:1032164/1‑139 (MQ=255)                                                                    aaTGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGAcc                                                                 >  1:628490/1‑125 (MQ=255)                                                                    aaTGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGAcc                                                                 <  2:628490/125‑1 (MQ=255)                                                                           ctATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTAc                                                                                                             >  2:133218/1‑74 (MQ=255)                                                                           ctATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTAc                                                                                                             <  1:133218/74‑1 (MQ=255)                                                                              tttAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGTAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGGTAGCATCCACCAGACCCATAACGATAAAGAACt                                                >  1:1009737/1‑132 (MQ=255)                                                                              tttAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATccacc                                                                      <  1:613640/110‑1 (MQ=255)                                                                              tttAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATccacc                                                                      >  2:613640/1‑110 (MQ=255)                                                                              tttAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAGAACt                                                <  2:1009737/132‑1 (MQ=255)                                                                                  gTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCg                                                                                                                         <  2:269663/55‑1 (MQ=255)                                                                                  gTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCg                                                                                                                         >  1:269663/1‑55 (MQ=255)                                                                                      gcGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCa                                                                            <  1:292751/95‑1 (MQ=255)                                                                                      gcGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCa                                                                            >  2:292751/2‑96 (MQ=255)                                                                                       cGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTa                                                                                                              >  1:596720/1‑61 (MQ=255)                                                                                       cGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTa                                                                                                              <  2:596720/61‑1 (MQ=255)                                                                                        gTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAGAACTGAGTACGCAGCAGAGGa                               >  1:1271005/1‑139 (MQ=255)                                                                                         ttCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAGAACTGAGTACGCAGCAGAGGaa                              >  1:572683/1‑139 (MQ=255)                                                                                           cTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAGAACTGAGTACGCAGCAGAGGAatc                            <  2:572683/139‑1 (MQ=255)                                                                                            tGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAGAACTGAGTACGCAGCAGAGGAatca                           >  1:780337/1‑139 (MQ=255)                                                                                                   gCTCTTTAAATAAAAGCAACGCTTATTACGCGACAg                                                                                                                           <  1:369148/36‑1 (MQ=255)                                                                                                   gCTCTTTAAATAAAAGCAACGCTTATTACGCGACAg                                                                                                                           >  2:369148/1‑36 (MQ=255)                                                                                                   gCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTAcagaccc                                                                                                       <  1:800790/56‑1 (MQ=255)                                                                                                   gCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTAcagaccc                                                                                                       >  2:800790/1‑56 (MQ=255)                                                                                                   gCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAGAACTGAGTACGCAGCAGAGGAat                             >  2:1101619/1‑130 (MQ=255)                                                                                                   gCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAGAACTGAGTACGCAGCAGAGGAat                             <  1:1101619/130‑1 (MQ=255)                                                                                                    ctctTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTAcagacccagac                                                                                                   >  1:327989/1‑59 (MQ=255)                                                                                                    ctctTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTAcagacccagac                                                                                                   <  2:327989/59‑1 (MQ=255)                                                                                                    ctctTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATccac                                                                       <  1:262484/87‑1 (MQ=255)                                                                                                    ctctTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATccac                                                                       >  2:262484/1‑87 (MQ=255)                                                                                                    ctctTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCAt                                                              >  1:307181/1‑96 (MQ=255)                                                                                                    ctctTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCAt                                                              <  2:307181/96‑1 (MQ=255)                                                                                                    ctctTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATa                                                       >  2:1189158/1‑103 (MQ=255)                                                                                                    ctctTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATa                                                       <  1:1189158/103‑1 (MQ=255)                                                                                                      ctTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATc                                                                                                                   >  2:1165333/1‑41 (MQ=255)                                                                                                      ctTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATc                                                                                                                   <  1:1165333/41‑1 (MQ=255)                                                                                                        ttAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCa                                                                            <  1:198996/78‑1 (MQ=255)                                                                                                        ttAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCa                                                                            >  2:198996/1‑78 (MQ=255)                                                                                                               aaaGCAACGCTTATTACGCGACAGCGAACATCACGTa                                                                                                              <  2:1057977/37‑1 (MQ=255)                                                                                                               aaaGCAACGCTTATTACGCGACAGCGAACATCACGTa                                                                                                              >  1:1057977/1‑37 (MQ=255)                                                                                                                aaGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGa                                                                                 <  2:1128179/65‑1 (MQ=255)                                                                                                                aaGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGa                                                                                 >  1:1128179/1‑65 (MQ=255)                                                                                                                 aGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAGAACTGAGTACGCAGCAGAGGAATCAGATCAGGTTGACGCGCTgcgc      >  1:674043/1‑139 (MQ=255)                                                                                                                 aGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAGAACTGAGTACGCAGCAGAGGAATCAGATCAGGTTGACGCGCTgcgc      >  1:1131609/1‑139 (MQ=255)                                                                                                                     aCGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACg                                                          <  2:1393135/83‑1 (MQ=255)                                                                                                                     aCGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACg                                                          >  1:1393135/1‑83 (MQ=255)                                                                                                                     aCGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAGAACTGAGTACGCAGCAGAGGAATCAGATCAGGTTGACGCGCTGCGCCTTc  <  2:1032164/139‑1 (MQ=255)                                                                                                                            |                                                                                                                                    CAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTACTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAGAACTGAGTACGCAGCAGAGGAATCAGATCAGGTTGACGCGCTGCGCCTTC  >  NC_000913/3920827‑3921081

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