Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,921,515	C→T	T179T (ACG→ACA)	atpB ←	F0 sector of membrane‑bound ATP synthase, subunit a

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,921,515	0	C	T	100.0%	127.5 / NA	41	T179T (ACG→ACA)	atpB	F0 sector of membrane‑bound ATP synthase, subunit a
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (25/16);  total (25/16)

AATCAGCTCACCGGCATACATGTTACCGAACAGTCGCAAACCGAGTGAAACTGGTTTGGACAGCAGGCTTACCCCTTCAAGGATTAAGTTGACAGGAATGAACGCCCAGTGATTGAACGGCTGCAGCGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCATTTTGATGCTGTAGAACAGAATCAGGATAAATACGCCCAGTGCCATAGACAGCGTTACGTTCACGTCCGCAGACGGAACCACACGCAGTGCAGGCAGA  >  NC_000913/3921389‑3921647                                                                                                                               |                                                                                                                                     aaTCAGCTCACCGGCATACATGTTACCGAACAGTCGCAAACCGAGTGAAACTGGTTTGGACAGCAGGCTTACCCCTTCAAGGATTAAGTTGACAGGAATGAACGCCCAGTGATTGAACGGCTGCAGTGTCAACTCTTTc                                                                                                                          >  1:245560/1‑139 (MQ=255)                         aCCGAACAGTCGCAAACCGAGTGAAACTGGTTTGGACAGCAGGCTTACCCCTTCAAGGATTAAGTTGACAGGAATGAACGCCCAGTGATTGAACGGCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCttt                                                                                                      >  2:784473/1‑135 (MQ=255)                         aCCGAACAGTCGCAAACCGAGTGAAACTGGTTTGGACAGCAGGCTTACCCCTTCAAGGATTAAGTTGACAGGAATGAACGCCCAGTGATTGAACGGCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCttt                                                                                                      <  1:784473/135‑1 (MQ=255)                            gAACAGTCGCAAACCGAGTGAAACTGGTTTGGACAGCAGGCTTACCCCTTCAAGGATTAAGTTGACAGGAATGAACGCCCAGTGATTGAACGGCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCATTTTg                                                                                               >  2:693358/1‑139 (MQ=255)                             aaCAGTCGCAAACCGAGTGAAACTGGTTTGGACAGCAGGCTTACCCCTTCAAGGATTAAGTTGACAGGAATGAACGCCCAGTGATTGAACGGCTGCAGTGTCAACTCTTTCGTGAAGccgccg                                                                                                              >  1:474937/1‑123 (MQ=255)                             aaCAGTCGCAAACCGAGTGAAACTGGTTTGGACAGCAGGCTTACCCCTTCAAGGATTAAGTTGACAGGAATGAACGCCCAGTGATTGAACGGCTGCAGTGTCAACTCTTTCGTGAAGccgccg                                                                                                              <  2:474937/123‑1 (MQ=255)                             aaCAGTCGCAAACCGAGTGAAACTGGTTTGGACAGCAGGCTTACCCCTTCAAGGATTAAGTTGACAGGAATGAACGCCCAGTGATTGAACGGCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCATTTTGa                                                                                              >  1:835396/1‑139 (MQ=255)                                aGTCGCAAACCGAGTGAAACTGGTTTGGACAGCAGGCTTACCCCTTCAAGGATTAAGTTGACAGGAATGAACGCCCAGTGATTGAACGGCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCATTTTGATGc                                                                                           >  1:964003/1‑139 (MQ=255)                                 gTCGCAAACCGAGTGAAACTGGTTTGGACAGCAGGCTTACCCCTTCAAGGATTAAGTTGACAGGAATGAACGCCCAGTGATTGAACGGCTGCAGTGTCaa                                                                                                                                 >  2:895010/1‑100 (MQ=255)                                 gTCGCAAACCGAGTGAAACTGGTTTGGACAGCAGGCTTACCCCTTCAAGGATTAAGTTGACAGGAATGAACGCCCAGTGATTGAACGGCTGCAGTGTCaa                                                                                                                                 <  1:895010/100‑1 (MQ=255)                                  tCGCAAACCGAGTGAAACTGGTTTGGACAGCAGGCTTACCCCTTCAAGGATTAAGTTGACAGGAATGAACGCCCAGTGATTGAACGGCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCATTTTGATGCTg                                                                                         >  2:288071/1‑139 (MQ=255)                                       aaCCGAGTGAAACTGGTTTGGACAGCAGGCTTACCCCTTCAAGGATTAAGTTGACAGGAATGAACGCCCAGTGATTGAACGGCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCATTTTGATGCTGTAGaa                                                                                    <  2:835396/139‑1 (MQ=255)                                         ccGAGTGAAACTGGTTTGGACAGCAGGCTTACCCCTTCAAGGATTAAGTTGACAGGAATGAACGCCCAGTGATTGAACGGCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCATTTTGATGCTGTAGAACa                                                                                  <  2:709776/139‑1 (MQ=255)                                           gAGTGAAACTGGTTTGGACAGCAGGCTTACCCCTTCAAGGATTAAGTTGACAGGAATGAACGCCCAGTGATTGAACGGCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCATTTTGATGCTGTAGAACAGa                                                                                <  2:964003/139‑1 (MQ=255)                                                            aCAGCAGGCTTACCCCTTCAAGGATTAAGTTGACAGGAATGAACGCCCAGTGATTGAACGGCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCATTTTGATGCTGTAGAACAGAATCAGGATAAATACGcc                                                               <  1:498948/139‑1 (MQ=255)                                                            aCAGCAGGCTTACCCCTTCAAGGATTAAGTTGACAGGAATGAACGCCCAGTGATTGAACGGCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCATTTTGATGCTGTAGAACAGAATCAGGATAAATACGcc                                                               >  2:498948/1‑139 (MQ=255)                                                              agcagGCTTACCCCTTCAAGGATTAAGTTGACAGGAATGAACGCCCAGTGATTGAACGGCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCATTTTGATGCTGTAGAACAGAATCAGGATAAATACGCCCa                                                             <  2:857011/139‑1 (MQ=255)                                                                      tACCCCTTCAAGGATTAAGTTGACAGGAATGAACGCCCAGTGATTGAACGGCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCt                                                                                                        >  1:298838/1‑88 (MQ=255)                                                                      tACCCCTTCAAGGATTAAGTTGACAGGAATGAACGCCCAGTGATTGAACGGCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCt                                                                                                        <  2:298838/88‑1 (MQ=255)                                                                             tCAAGGATTAAGTTGACAGGAATGAACGCCCAGTGATTGAACGGCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCATTTTGATGCTGTAGAACAGAATCAGGATAAATACGCCCAGTGCCATAGACAGCg                                              >  2:357814/1‑139 (MQ=255)                                                                                   aTTAAGTTGACAGGAATGAACGCCCAGTGATTGAACGGCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCATTTTGATGCTGTAGAACAGAATCAGGATAAATACGCCCa                                                             >  2:89169/1‑118 (MQ=255)                                                                                   aTTAAGTTGACAGGAATGAACGCCCAGTGATTGAACGGCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCATTTTGATGCTGTAGAACAGAATCAGGATAAATACGCCCa                                                             <  1:89169/118‑1 (MQ=255)                                                                                   aTTAAGTTGACAGGAATGAACGCCCAGTGATTGAACGGCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCATTTTGATGCTGTAGAACAGAATCAGGATAAATACGCCCAGTGCCATAGACAGCGTTACGt                                        >  1:961549/1‑139 (MQ=255)                                                                                            aCAGGAATGAACGCCCAGTGATTGAACGGCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCATTTTGATGCTGTAGAACAGAATCAGGATAAATACGCCCAGTGCCATAGACAGCGTTACGTTCACGTCCg                               >  2:951092/1‑139 (MQ=255)                                                                                                  aTGAACGCCCAGTGATTGAACGGCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCATTTTGATGCTGTAGAACAGAATCAGGATAAATACGCCCAGTGCCATAGACAGCGTTACGTTCACGTCCGCAGACg                         <  1:951092/139‑1 (MQ=255)                                                                                                          ccAGTGATTGAACGGCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCATTTTGATGCTGTAGAACAGAAt                                                                              >  2:642673/1‑78 (MQ=255)                                                                                                          ccAGTGATTGAACGGCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCATTTTGATGCTGTAGAACAGAAt                                                                              <  1:642673/78‑1 (MQ=255)                                                                                                          ccAGTGATTGAACGGCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCATTTTGATGCTGTAGAACAGAATCAGGATAAATACGcc                                                               >  1:223236/1‑93 (MQ=255)                                                                                                          ccAGTGATTGAACGGCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCATTTTGATGCTGTAGAACAGAATCAGGATAAATACGcc                                                               <  2:223236/93‑1 (MQ=255)                                                                                                          ccAGTGATTGAACGGCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCATTTTGATGCTGTAGAACAGAATCAGGATAAATACGCCCAGTGCCATAGACAGCGTTACGTTCACGTCCGCAGACGGAACcaca                 >  1:544782/1‑139 (MQ=255)                                                                                                          ccAGTGATTGAACGGCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCATTTTGATGCTGTAGAACAGAATCAGGATAAATACGCCCAGTGCCATAGACAGCGTTACGTTCACGTCCGCAGACGGAACcaca                 >  2:19778/1‑139 (MQ=255)                                                                                                           cAGTGATTGAACGGCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCATTTTGa                                                                                              >  2:305923/1‑61 (MQ=255)                                                                                                           cAGTGATTGAACGGCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCATTTTGa                                                                                              <  1:305923/61‑1 (MQ=255)                                                                                                            aGTGATTGAACGGCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCATTTTGATGCTGTAGAACAGAATCAGGATAAATACGCCCAGTGCCATAGACa                                                 <  2:960878/105‑1 (MQ=255)                                                                                                            aGTGATTGAACGGCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCATTTTGATGCTGTAGAACAGAATCAGGATAAATACGCCCAGTGCCATAGACa                                                 >  1:960878/1‑105 (MQ=255)                                                                                                                      cGGCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCATTTTGATGCTGTAGAACAGAATCAGGATAAATACGCCCAGTGCCATAGACAGCGTTACGTTCACGTCCGCAGACGGAACCACACGCAGTgcaggc     >  2:430949/1‑139 (MQ=255)                                                                                                                       ggCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCATTTTGATGCTGTAGAACAGAATCAGGATAAATACGCCCAGTGCCATAGACAGCGTTACGTTCACGTCCGCAGACGGAACCACACGCAGTgcaggca    <  1:547248/139‑1 (MQ=255)                                                                                                                        gCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCATTTTGATGCTGTAGAACAGAATCAGGATAAATACGCCCAGTGCCATAGACAGCGTTACGTTCACGTCCGCAGACGGAACCACACGCAGTgcaggcag   >  2:364795/1‑139 (MQ=255)                                                                                                                        gCTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCATTTTGATGCTGTAGAACAGAATCAGGATAAATACGCCCAGTGCCATAGACAGCGTTACGTTCACGTCCGCAGACGGAACCACACGCAGTgcaggcag   >  2:244882/1‑139 (MQ=255)                                                                                                                         cTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCATTTTGATGCTGTAGAACAGAATCAGGATAAATACGCCCAGTGCCATAGACAGCGTTACGTTCACGTCCGCAGACGGAACCACACGCAGTGCAGGCAGa  >  2:611470/1‑139 (MQ=255)                                                                                                                         cTGCAGTGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCATTTTGATGCTGTAGAACAGAATCAGGATAAATACGCCCAGTGCCATAGACAGCGTTACGTTCACGTCCGCAGACGGAACCACACGCAGTGCAGGCAGa  >  1:909889/1‑139 (MQ=255)                                                                                                                               |                                                                                                                                     AATCAGCTCACCGGCATACATGTTACCGAACAGTCGCAAACCGAGTGAAACTGGTTTGGACAGCAGGCTTACCCCTTCAAGGATTAAGTTGACAGGAATGAACGCCCAGTGATTGAACGGCTGCAGCGTCAACTCTTTCGTGAAGCCGCCGATGCCTTTCATTTTGATGCTGTAGAACAGAATCAGGATAAATACGCCCAGTGCCATAGACAGCGTTACGTTCACGTCCGCAGACGGAACCACACGCAGTGCAGGCAGA  >  NC_000913/3921389‑3921647

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