Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,917,313	A→G	D21D (GAT→GAC)	atpD ←	ATP synthase F1 complex subunit beta

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,917,313	0	A	G	100.0%	154.0 / NA	45	D21D (GAT→GAC)	atpD	ATP synthase F1 complex subunit beta
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (21/24);  total (21/24)

ACCGCGACGCAGACCGTCGGAGGAACCCATTGCGATGGTACGTACGATACCGCCGCCGAGCTGCTGCTGAACTTCCAGCACCAGACGCTCATTACCATTTTGCACCTCAAGAGCATCGTACACGCGCGGTACGGCATCCTGAGGGAATTCGACGTCAACTACGGCGCCGATTACCTGGACAATCTTTCCAGTAGCCATCTTAAATCCTCTACGAAATAACCTGTTTAAACCGCGGCGGCCCCCGAGACGATCTCGGTGAGTTCCTGAG  >  NC_000913/3917178‑3917445                                                                                                                                        |                                                                                                                                     aCCGCGACGCAGACCGTCGGAGGAACCCATTGCGATGGTACGTACGATACCGCCGCCGAGCTGCTGCTGAACTTCCAGCACCAGACGCTCATTACCATTTTGCACCTCAAGAGCATCGTACACGCGCGGTACGGCGTCCt                                                                                                                                  >  1:211012/1‑140 (MQ=255)     cgACGCAGACCGTCGGAGGAACCCATTGCGATGGTACGTACGATACCGCCGCCGAGCTGCTGCTGAACTTCCAGCACCAGACGCTCATTACCATTTTGCACCTCAAGAGCATCGTACACGCGCGGTACGGCGTCCTGAgg                                                                                                                              >  1:27708/1‑140 (MQ=255)       aCGCAGACCGTCGGAGGAACCCATTGCGATGGTACGTACGATACCGCCGCCGAGCTGCTGCTGAACTTCCAGCACCAGACGCTCATTACCATTTTGCACCTCAAGAGCATCGTACACGCGCGGTACGGCGTCCTGAGGGa                                                                                                                            >  2:165725/1‑140 (MQ=255)        cGCAGACCGTCGGAGGAACCCATTGCGATGGTACGTACGATACCGCCGCCGAGCTGCTGCTGAACTTCCAGCACCAGACGCTCATTACCATTTTGCACCTCAAGAGCATCGTACACGCGCGGTACGGCGTCCTGAGGGaa                                                                                                                           >  1:672493/1‑140 (MQ=255)           aGACCGTCGGAGGAACCCATTGCGATGGTACGTACGATACCGCCGCCGAGCTGCTGCTGAACTTCCAGCACCAGACGCTCATTACCATTTTGCACCTCAAGAGCATCGTACACGCGCGGTACGGCGTCCTGAGGGAATTc                                                                                                                        <  1:846066/140‑1 (MQ=255)             aCCGTCGGAGGAACCCATTGCGATGGTACGTACGATACCGCCGCCGAGCTGCTGCTGAACTTCCAGCACCAGACGCTCATTACCATTTTGCACCTCAAGAGCATCGTACACGCGCGGTACGGCGTCCTGAGGGAATTCg                                                                                                                       <  2:445304/139‑1 (MQ=255)             aCCGTCGGAGGAACCCATTGCGATGGTACGTACGATACCGCCGCCGAGCTGCTGCTGAACTTCCAGCACCAGACGCTCATTACCATTTTGCACCTCAAGAGCATCGTACACGCGCGGTACGGCGTCCTGAGGGAATTCg                                                                                                                       >  1:445304/1‑139 (MQ=255)                 tCGGAGGAACCCATTGCGATGGTACGTACGATACCGCCGCCGAGCTGCTGCTGAACTTCCAGCACCAGACGCTCATTACCATTTTGCACCTCAAGAGCATCGTACACGCGCGGTACGGCGTCCTGAGGGAATTCGACGTc                                                                                                                  >  2:446197/1‑140 (MQ=255)                      ggaACCCATTGCGATGGTACGTACGATACCGCCGCCGAGCTGCTGCTGAACTTCCAGCACCAGACGCTCATTACCATTTTGCACCTCAAGAGCATCGTACACGCGCGGTACGGCGTCCTGAGGGAATTCGACGTCAACTa                                                                                                             <  2:710026/140‑1 (MQ=255)                              ttGCGATGGTACGTACGATACCGCCGCCGAGCTGCTGCTGAACTTCCAGCACCAGACGCTCATTACCATTTTGCACCTCAAGAGCATCGTACACGCGCGGTACGGCGTCCTGAGGGAATTCGACGTCAACTACGGCGCCg                                                                                                     <  1:53315/140‑1 (MQ=255)                              ttGCGATGGTACGTACGATACCGCCGCCGAGCTGCTGCTGAACTTCCAGCACCAGACGCTCATTACCATTTTGCACCTCAAGAGCATCGTACACGCGCGGTACGGCGTCCTGAGGGAATTCGACGTCAACTACGGCGCCg                                                                                                     <  2:672493/140‑1 (MQ=255)                              ttGCGATGCTACGTACAGTACCTCCGCCGACCTGCGGCTTAACTTCCATCACCAGACGCTCATTACCATTTTGCACCTCAAGAGCATCGTACACGCGCGGTACGGCGTCCTGAGGGAATTCGACGTCAACTACGGCGCCg                                                                                                     <  2:407180/140‑1 (MQ=255)                               tGCGATGGTACGTACGATACCGCCGCCGAGCTGCTGCTGAACTTCCAGCACCAGACGCTCATTACCATTTTGCACCTCAAGAGCATCGTACGCGCGCGGTTCGGCGGCCTGAGGGAAATCGAGGTCAACTAGGGCGCCGa                                                                                                    >  2:413477/1‑140 (MQ=255)                                         cgtacgATACCGCCGCCGAGCTGCTGCTGAACTTCCAGCACCAGACGCTCATTACCATTTTGCACCTCAAGAGCATCGTACACGCGCGGTACGGCGTCCTGAGGGAATTCGACGTCAACTACGGCGCCGATTACCTGGAc                                                                                          >  2:195403/1‑140 (MQ=255)                                            acgATACCGCCGCCGATATGCTCCTGAACTTGCTGCACCTGACGCTCATTACCATTTTGCACCTCAAGAGCATCGTACACGCGCGGTACGGCGTCCTGAGGGAATTCGACGTCAACTACGGCGCCGATTACCTGGACAAt                                                                                       <  2:389177/140‑1 (MQ=255)                                                tACCGCCGCCGAGCTGCTGCTGAACTTCCAGCACCAGACGCTCATTACCATTTTGCACCTCAAGAGCATCGTACACGCGCGGTACGGCGTCCTGAGGGAATTCGACGTCAACTACGGCGCCGATTACCTGGACAATCttt                                                                                   >  2:609605/1‑140 (MQ=255)                                                 aCCGCCGCCGAGCTGCTGCTGAACTTCCAGCACCAGACGCTCATTACCATTTTGCACCTCAAGAGCATCGTACACGCGCGGTACGGCGTCCTGAGGGAATTCGACGTCAACTACGGCGCCGATTACCTGGACAATCTTTc                                                                                  <  1:830341/140‑1 (MQ=255)                                                   cgccgccgAGCTGCTGCTGAACTTCCAGCACCAGACGCTCATTACCATTTTGCACCTCAAGAGCATCGTACACGCGCGGTACGGCGTCCTGAGGGAATTCGACGTCAACTACGGCGCCGATTACCTGGACAATCTTTCCa                                                                                >  1:212994/1‑140 (MQ=255)                                                     ccgccgAGCTGCTGCTGAACTTCCAGCACCAGACGCTCATTACCATTTTGCACCTCAAGAGCATCGTACACGCGCGGTACGGCGTCCTGAGGGAATTCGACGTCAACTACGGCGCCGATTACCTGGACAATCTTTCCAGt                                                                              >  2:735471/1‑140 (MQ=255)                                                        ccgAGCTGCTGCTGAACTTCCAGCACCAGACGCTCATTACCATTTTGCACCTCAAGAGCATCGTACACGCGCGGTACGGCGTCCTGAGGGAATTCGACGTCAACTACGGCGCCGATTACCTGGACAATCTTTCCAGTAGc                                                                           <  1:571799/140‑1 (MQ=255)                                                         cgAGCTGCTGCTGAACTTCCAGCACCAGACGCTCATTACCATTTTGCACCTCAAGAGCATCGTACACGCGCGGTACGGCGTCCTGAGGGAATTCGACGTCAACTACGGCGCCGATTACCTGGACAATCTTTCCAGTAGcc                                                                          <  1:735471/140‑1 (MQ=255)                                                         cgAGCTGCTGCTGAACTTCCAGCACCAGACGCTCATTACCATTTTGCACCTCAAGAGCATCGTACACGCGCGGTACGGCGTCCTGAGGGAATTCGACGTCAACTACGGCGCCGATTACCTGGACAATCTTTCCAGTAGcc                                                                          <  2:211012/140‑1 (MQ=255)                                                         cgAGCTGCTGCTGAACTTCCAGCACCAGACGCTCATTACCATTTTGCACCTCAAGAGCATCGTACACGCGCGGTACGGCGTCCTGAGGGAATTCGACGTCAACTACGGCGCCGATTACCTGGACAATCTTTCCAGTAGcc                                                                          <  1:195403/140‑1 (MQ=255)                                                                   ctgAACTTCCAGCACCAGACGCTCATTACCATTTTGCACCTCAAGAGCATCGTACACGCGCGGTACGGCGTCCTGAGGGAATTCGACGTCAACTACGGCGCCGATTACCTGGACAATCTTTCCAGTAGCCATCTTAAATc                                                                <  1:687081/140‑1 (MQ=255)                                                                        cTTCCAGCACCAGACGCTCATTACCATTTTGCACCTCAAGAGCATCGTACACGCGCGGTACGGCGTCCTGAGGGAATTCGACGTCAACTACGGCGCCGATTACCTGGACAATCTTTCCAGTAGc                                                                           <  2:799446/124‑1 (MQ=255)                                                                        cTTCCAGCACCAGACGCTCATTACCATTTTGCACCTCAAGAGCATCGTACACGCGCGGTACGGCGTCCTGAGGGAATTCGACGTCAACTACGGCGCCGATTACCTGGACAATCTTTCCAGTAGc                                                                           >  1:799446/1‑124 (MQ=255)                                                                          tCCAGCACCAGACGCTCATTACCATTTTGCACCTCAAGAGCATCGTACACGCGCGGTACGGCGTCCTGAGGGAATTCGACGTCAACTACGGCGCCGATTACCTGGACAATCTTTCCAGTAGCCATCTTAAATc                                                                >  1:620029/1‑133 (MQ=255)                                                                          tCCAGCACCAGACGCTCATTACCATTTTGCACCTCAAGAGCATCGTACACGCGCGGTACGGCGTCCTGAGGGAATTCGACGTCAACTACGGCGCCGATTACCTGGACAATCTTTCCAGTAGCCATCTTAAATc                                                                <  2:620029/133‑1 (MQ=255)                                                                            cAGCACCAGACGCTCATTACCATTTTGCACCTCAAGAGCATCGTACACGCGCGGTACGGCGTCCTGAGGGAATTCGACGTCAACTACGGCGCCGATTACCTGGACAATCTTTCCAGTAGCCATCTTAAATCCTCTACGaa                                                       <  2:212994/140‑1 (MQ=255)                                                                             aGCACCAGACGCTCATTACCATTTTGCACCTCAAGAGCATCGTACACGCGCGGTACGGCGTCCTGAGGGAATTCGACGTCAACTACGGCGCCGATTACCTGGACAATCTTTCCAGTAGCCATCTTAAATCCTCTACGaaa                                                      <  1:444255/140‑1 (MQ=255)                                                                              gCACCAGACGCTCATTACCATTTTGCACCTCAAGAGCATCGTACACGCGCGGTACGGCGTCCTGAGGGAATTCGACGTCAACTACGGCGCCGATTACCTGGACAATCTTTCCAGTAGCCATCTTAAATCCTCTACGAAAt                                                     <  1:609605/140‑1 (MQ=255)                                                                              gCACCAGACGCTCATTACCATTTTGCACCTCAAGAGCATCGTACACGCGCGGTACGGCGTCCTGAGGGAATTCGACGTCAACTACGGCGCCGATTACCTGGACAATCTTTCCAGTAGCCATCTTAAATCCTCTACGAAAt                                                     <  2:622596/140‑1 (MQ=255)                                                                                aCCAGACGCTCATTACCATTTTGCACCTCAAGAGCATCGTACACGCGCGGTACGGCGTCCTGAGGGAATTCGACGTCAACTACGGCGCCGATTACCTGGACAATCTTTCCAGTAGCCATCTTAAATCCTCTACGAAATaa                                                   <  2:569284/140‑1 (MQ=255)                                                                                    gACGCTCATTACCATTTTGCACCTCAAGAGCATCGTACACGCGCGGTACGGCGTCCTGAGGGAATTCGACGTCAACTACGGCGCCGATTACCTGGACAATCTTTCCAGTAGCCATCTTAAATCCTCTACGAAATAACCTg                                               >  1:597108/1‑140 (MQ=255)                                                                                              aCCATTTTGCACCTCAAGAGCATCGTACACGCGCGGTACGGCGTCCTGAGGGAATTCGACGTCAACTACGGCGCCGATTACCTGGACAATCTTTCCAGTAGCCATCTTAAATCCTCTACGAAATAACCTGTTTAAACcgc                                     >  1:719902/1‑140 (MQ=255)                                                                                                  ttttGCACCTCAAGAGCATCGTACACGCGCGGTACGGCGTCCTGAGGGAATTCGACGTCAACTACGGCGCCGATTACCTGGACAATCTTTCCAGTAGCCATCTTAAATCCTCTACGAAATAACCTGTTTAAACCgcggcg                                 <  2:597108/140‑1 (MQ=255)                                                                                                     tGCACCTCAAGAGCATCGTACACGCGCGGTACGGCGTCCTGAGGGAATTCGACGTCAACTACGGCGCCGATTACCTGGACAATCTTTCCAGTAGCCATCTTAAATCCTCTACGAAATAACCTGTTTAAACCGCGGCGGcc                              >  1:49152/1‑140 (MQ=255)                                                                                                       cACCTCAAGAGCATCGTACACGCGCGGTACGGCGTCCTGAGGGAATTCGACGTCAACTACGGCGCCGATTACCTGGACAATCTTTCCAGTAGCCATCTTAAATCCTCTACGAAATAACCTGTTTAAACCGCGGCGGcccc                            >  2:586195/1‑140 (MQ=255)                                                                                                         ccTCAAGAGCATCGTACACGCGCGGTACGGCGTCCTGAGGGAATTCGACGTCAACTACGGCGCCGATTACCTGGACAATCTTTCCAGTAGCCATCTTAAATCCTCTACGAAATAACCTGTTTAAACCGcggcgg                                <  2:660912/134‑1 (MQ=255)                                                                                                         ccTCAAGAGCATCGTACACGCGCGGTACGGCGTCCTGAGGGAATTCGACGTCAACTACGGCGCCGATTACCTGGACAATCTTTCCAGTAGCCATCTTAAATCCTCTACGAAATAACCTGTTTAAACCGcggcgg                                >  1:660912/1‑134 (MQ=255)                                                                                                                         cacGCGCGGTACGGCGTCCTGAGGGAATTCGACGTCAACTACGGCGCCGATTACCTGGACAATCTTTCCAGTAGCCATCTTAAATCCTCTACGAAATAACCTGTTTAAACCGCGGCGGCCCCCGAGACGATCTCGGTGAg          >  1:21407/1‑140 (MQ=255)                                                                                                                            gcgcgGTACGGCGTCCTGAGGGAATTCGACGTCAACTACGGCGCCGATTACCTGGACAATCTTTCCAGTAGCCATCTTAAATCCTCTACGAAATAACCTTTTTAAACCGCGGCGGCCCCCGAGACGATCTCGGTGAGTTc       <  2:605816/140‑1 (MQ=255)                                                                                                                               cgGTACGGCGTCCTGAGGGAATTCGACGTCAACTACGGCGCCGATTACTTGGACAATCTTTCCAGTAGCCATCTTAAATCCTCTACGAAATAACCTGTTTAAACCGCGGCGGCCCCCGAGACGATCTCGGTGAGTTCCTg    >  1:207318/1‑140 (MQ=255)                                                                                                                               cgGTACGGCGTCCTGAGGGAATTCGACGTCAACTACGGCGCCGATTACCTGGACAATCTTTCCAGTAGCCATCTTAAATCCTCTACGAAATAACCTGTTTAAACCGCGGCGGCCCCCGAGACGATCTCGGTGAGTTCCTg    >  1:797979/1‑140 (MQ=255)                                                                                                                                 gTACGGCGTCCTGAGGGAATTCTACTTCAACTACGGCACCGATTACCTGGACAATCTTTCCAGTAGCCATCTTAAATCCTCTACGAAATAACCTGTTTAAACCGCGGCGGCCCCCGAGACGATCTCGGTGAGTTCCTGAg  <  1:413477/140‑1 (MQ=255)                                                                                                                                        |                                                                                                                                     ACCGCGACGCAGACCGTCGGAGGAACCCATTGCGATGGTACGTACGATACCGCCGCCGAGCTGCTGCTGAACTTCCAGCACCAGACGCTCATTACCATTTTGCACCTCAAGAGCATCGTACACGCGCGGTACGGCATCCTGAGGGAATTCGACGTCAACTACGGCGCCGATTACCTGGACAATCTTTCCAGTAGCCATCTTAAATCCTCTACGAAATAACCTGTTTAAACCGCGGCGGCCCCCGAGACGATCTCGGTGAGTTCCTGAG  >  NC_000913/3917178‑3917445

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