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%	240.7 / NA	73	*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 (39/34);  total (39/34)

CGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTACTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAGAACTGAGTACGCAGCAGAGGAATCAGATCAGGTTGACGCGCTGCG  >  NC_000913/3920816‑3921076                                                                                                                                       |                                                                                                                               cGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTAcg                                                                                                                            <  2:653023/139‑1 (MQ=255)         aTGCAGAACAGAACGAACAGTACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCg                                                                                                                    <  1:382721/139‑1 (MQ=255)          tGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTCGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGa                                                                                                                   >  1:194511/1‑139 (MQ=255)          tGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGa                                                                                                                   >  2:1387002/1‑139 (MQ=255)          tGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGa                                                                                                                   <  2:194511/139‑1 (MQ=255)          tGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGa                                                                                                                   >  2:233728/1‑139 (MQ=255)                         acAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGTGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTAcgcg                                                                                                                          <  1:70639/117‑1 (MQ=255)                         acAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTAcgcg                                                                                                                          >  2:70639/1‑117 (MQ=255)                           aGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTAcagaccc                                                                                                  >  1:453627/1‑139 (MQ=255)                               cAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTAcagacccagac                                                                                              >  1:70610/1‑139 (MQ=255)                                       aTGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACCGACCCAGACCAAAAGCg                                                                                      >  1:1303422/1‑139 (MQ=255)                                       aTGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCg                                                                                      >  2:1209348/1‑139 (MQ=255)                                        tggcctggccGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGa                                                                                     <  2:70610/139‑1 (MQ=255)                                                   aGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATa                                                                          <  1:1387002/139‑1 (MQ=255)                                                    ggATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAGCGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAg                                                                         >  1:1052726/1‑139 (MQ=255)                                                      aTTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCa                                                                       <  2:453627/139‑1 (MQ=255)                                                            gCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATccacc                                                                 <  2:706188/139‑1 (MQ=255)                                                             cGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATccacca                                                                <  2:1303422/139‑1 (MQ=255)                                                             cGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATccacca                                                                <  2:846351/139‑1 (MQ=255)                                                              gTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACACGACAGCGAACATCACGTACAGACCCAGACCTACAGCGatcat                                                                                 <  1:1419926/121‑1 (MQ=255)                                                              gTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACACGACAGCGAACATCACGTACAGACCCAGACCTACAGCGatcat                                                                                 >  2:1419926/1‑121 (MQ=255)                                                               ttAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGa                                                              <  1:837191/139‑1 (MQ=255)                                                                       cacaGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAAc                                                      >  1:358800/1‑139 (MQ=255)                                                                            cacaATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGc                                                                        >  2:771549/1‑116 (MQ=255)                                                                            cacaATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGc                                                                        <  1:771549/116‑1 (MQ=255)                                                                            cacaATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCa                                                          >  1:1044655/1‑130 (MQ=255)                                                                            cacaATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCa                                                          <  2:1044655/130‑1 (MQ=255)                                                                            cacaATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTAAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATaa                                                       <  2:549631/133‑1 (MQ=255)                                                                            cacaATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTAAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATaa                                                       >  1:549631/1‑133 (MQ=255)                                                                              caATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCTCGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAg                                               >  1:1055643/1‑139 (MQ=255)                                                                               aaTGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAGa                                              <  2:1052726/139‑1 (MQ=255)                                                                                 tGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTAcagacc                                                                                                   <  2:682107/84‑1 (MQ=255)                                                                                 tGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTAcagacc                                                                                                   >  1:682107/1‑84 (MQ=255)                                                                                 tGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAg                                                                                        >  2:118370/1‑95 (MQ=255)                                                                                 tGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAg                                                                                        <  1:118370/95‑1 (MQ=255)                                                                                 tGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCAt                                                         <  1:26881/126‑1 (MQ=255)                                                                                 tGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCAt                                                         >  2:26881/1‑126 (MQ=255)                                                                                            aGTTAACGTTCTGATATTGCCCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAGAACTGAGTACgcgg                                 >  1:865499/1‑137 (MQ=255)                                                                                             gTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAGAACTGAGTACgcagc                                >  1:968861/1‑139 (MQ=255)                                                                                             gTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAGAACTGAGTACgcagc                                >  1:428439/1‑139 (MQ=255)                                                                                              ttAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACGTACAGCGATCATCGGGa                                                                            <  2:921024/94‑1 (MQ=255)                                                                                              ttAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGa                                                                            >  1:921024/1‑94 (MQ=255)                                                                                              ttAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAGa                                              >  1:1204607/1‑124 (MQ=255)                                                                                              ttAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAGa                                              <  2:1204607/124‑1 (MQ=255)                                                                                                aaCGTTCTGATATTGCTCTTTAAATAAAATCAACGCTTATTACGCGACACCGAACATCACGTACAGACAAAGGCCTACAGCGATCAGCGCGATATCATCCACCAGACTCATGACGAGAAACATCTGAGCACGCAGCaga                             >  1:304307/1‑139 (MQ=255)                                                                                                aaCGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAGAACTGAGTACGCAGCaga                             >  2:320961/1‑139 (MQ=255)                                                                                                   gTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTAcagacc                                                                                                   <  1:1188988/66‑1 (MQ=255)                                                                                                   gTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTAcagacc                                                                                                   >  2:1188988/1‑66 (MQ=255)                                                                                                      cTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCg                                                                                                                    >  1:1468728/1‑46 (MQ=255)                                                                                                      cTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCg                                                                                                                    <  2:1468728/46‑1 (MQ=255)                                                                                                      cTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAGAACTGAGTACGCAGCaga                             >  1:255612/1‑133 (MQ=255)                                                                                                      cTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAGAACTGAGTACGCAGCaga                             <  2:255612/133‑1 (MQ=255)                                                                                                      cTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAGAACTGAGTACGCAGCaga                             <  2:930355/133‑1 (MQ=255)                                                                                                      cTGATATTGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAGAACTGAGTACGCAGCaga                             >  1:930355/1‑133 (MQ=255)                                                                                                             tGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCg                                                                                                                    >  1:548167/1‑39 (MQ=255)                                                                                                             tGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCg                                                                                                                    <  2:548167/39‑1 (MQ=255)                                                                                                             tGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAGaa                                             <  2:1243598/110‑1 (MQ=255)                                                                                                             tGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAGaa                                             >  1:1243598/1‑110 (MQ=255)                                                                                                             tGCTCTTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAGAACTGAGTACGCAGCAGAGGAatcagatcag                >  1:474833/1‑139 (MQ=255)                                                                                                               ctctTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTAcagacc                                                                                                   >  2:782177/1‑54 (MQ=255)                                                                                                               ctctTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTAcagacc                                                                                                   <  1:782177/54‑1 (MQ=255)                                                                                                               ctctTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAGAACTGAGTACGCAGCAGAGGAATCAGATCAGGt              >  2:1284831/1‑139 (MQ=255)                                                                                                               ctctTTAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAGAACTGAGTACGCAGCAGAGGAATCAGATCAGGt              >  2:139022/1‑139 (MQ=255)                                                                                                                   ttAAATAAAAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAGAACTGAGTACGCAGCAGAGGAATCAGATCAGGTTGAc          >  2:202475/1‑139 (MQ=255)                                                                                                                    taaataaaAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAGAACTGAGTACgcagca                               <  2:1153283/117‑1 (MQ=255)                                                                                                                    taaataaaAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAGAACTGAGTACgcagca                               >  1:1153283/1‑117 (MQ=255)                                                                                                                      aataaaAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAGAACTGAGTACGCAGCaga                             >  2:1259359/1‑117 (MQ=255)                                                                                                                      aataaaAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAGAACTGAGTACGCAGCaga                             <  1:1259359/117‑1 (MQ=255)                                                                                                                      aataaaAGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAGAACTGAGTACGCAGCAGAGGAATCAGATCAGGTTGAcgcg       <  2:358800/139‑1 (MQ=255)                                                                                                                           aaGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAGAACTGAGTACGCAGCAGAGGAATCAGATCAGGTTGACGCGCTgcg  <  2:474833/139‑1 (MQ=255)                                                                                                                           aaGCAACGCTTATTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAGAACTGAGTACGCAGCAGAGGAATCAGATCAGGTTGACGCGCTgcg  <  2:1055643/139‑1 (MQ=255)                                                                                                                            aGCAACGCTTATTACGCGACAGCGAACATCACGTAcagacc                                                                                                   >  2:689734/1‑41 (MQ=255)                                                                                                                            aGCAACGCTTATTACGCGACAGCGAACATCACGTAcagacc                                                                                                   <  1:689734/41‑1 (MQ=255)                                                                                                                                       |                                                                                                                               CGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCAACGCTTACTACGCGACAGCGAACATCACGTACAGACCCAGACCTACAGCGATCATCGGGATAGCATCCACCAGACCCATAACGATAAAGAACTGAGTACGCAGCAGAGGAATCAGATCAGGTTGACGCGCTGCG  >  NC_000913/3920816‑3921076

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