Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,917,346	G→A	I10I (ATC→ATT)	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,346	0	G	A	100.0%	142.3 / NA	40	I10I (ATC→ATT)	atpD	ATP synthase F1 complex subunit beta
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (17/23);  total (17/23)

ATTACCATTTTGCACCTCAAGAGCATCGTACACGCGCGGTACGGCATCCTGAGGGAATTCGACGTCAACTACGGCGCCGATTACCTGGACAATCTTTCCAGTAGCCATCTTAAATCCTCTACGAAATAACCTGTTTAAACCGCGGCGGCCCCCGAG  >  NC_000913/3917268‑3917423                                                                               |                                                                              aTTACCATTTTGCACCTCAAGAGCATCGTACACGCGCGGTACGGCATCCTGAGGGAATTCGACGTCAACTACGGCGCCAAt                                                                             >  2:711623/1‑81 (MQ=255)  ttACCATTTTGCACCTCAAGAGCATCGTACACGCGCGGTACGGCATCCTGAGGGAATTCGACGTCAACTACGGCGCCAAtt                                                                            <  2:711568/81‑1 (MQ=255)    aCCATTTTGCACCTCAAGAGCATCGTACACGCGCGGTACGGCATCCTGAGGGAATTCGACGTCAACTACGGCGCCAATTAc                                                                          >  2:711619/1‑81 (MQ=255)      cATTTTGCACCTCAAGAGCATCGTACACGCGCGGTACGGCATCCTGAGGGAATTCGACGTCAACTACGGCGCCAATTACCt                                                                        >  2:711630/1‑81 (MQ=255)      cATTTTGCACCTCAAGAGCATCGTACACGCGCGGTACGGCATCCTGAGGGAATTCGACGTCAACTACGGCGCCAATTACCt                                                                        <  1:711570/81‑1 (MQ=255)       aTTTTGCACCTCAAGAGCATCGTACACGCGCGGTACGGCATCCTGAGGGAATTCGACGTCAACTACGGCGCCAATTACCTg                                                                       <  2:711570/81‑1 (MQ=255)          ttGCACCTCAAGAGCATCGTACACGCGCGGTACGGCATCCTGAGGGAATTCGACGTCAACTACGGCGCCAATTACCTGGa                                                                     >  2:711621/1‑80 (MQ=255)          ttGCACCACAAGAGCATCGTACACTCGCGGTACGGCATCCTGAGGGAATTCGACGTCAACTACGGCGCCAATTACCTGGa                                                                     >  2:711638/1‑80 (MQ=255)            gCACCTCAAGAGCATCGTACACGCGCGGTACGGCATCCTGAGGGAATTCGACGTCAACTACGGCGCCAATTACCTGGACaa                                                                  >  1:711629/1‑81 (MQ=255)               ccTCAAGAGCATCGTACACGCGCGGTACGGCATCCTGAGGGAATTCGACGTCAACTACGGCGCCAATTACCTGGACAATCt                                                               <  1:711572/81‑1 (MQ=255)                       gCATCGTACACGCGCGGTACGGCATCCTGAGGGAATTCGACGTCAACTACGGCGCCAATTACCTGGAAAATCTTTCCAGTa                                                       <  2:711572/81‑1 (MQ=255)                            gTACACGCGCGGTACGGCATCCTGAGGGAATTCGACGTCAACTACGGCGCCAATTACCTGGACAATCTTTCCAGTAGCCAt                                                  >  1:711640/1‑81 (MQ=255)                               cacGCGCGGTACGGCATCCTGAGGGAATTCGACGTCAACTACGGCGCCAATTACCTGGACAATCTTTCCAGTAGCCATCtt                                               <  1:711574/81‑1 (MQ=255)                               cacGCGCGGTACGGCATCCTGAGGGAATTCGACGTCAACTACGGCGCCAATTACCTGGACAATCTTTCCAGTAGCCATCtt                                               <  1:711575/81‑1 (MQ=255)                               cacGCGCGGTACGGCATCCTGAGGGAATTCGACGTCAACTACGGCGCCAATTACCTGGACAATCTTTCCAGTAGCCATCtt                                               <  1:711576/81‑1 (MQ=255)                               cacGCGCGGTACGGCATCCTGAGGGAATTCGACGTCAACTACGGCGCCAATTACCTGGACAATCTTTCCAGTAGCCATCt                                                <  2:711576/80‑1 (MQ=255)                                acGCGCGGTACGGCATCCTGAGGGAATTCGACGTCAACTACGGCGCCAATTACCTGGACAATCTTTCCAGTAGCCATCTTa                                              <  1:711578/81‑1 (MQ=255)                                  gcgcgGTACGGCATCCTGAGGGAATTCGACGTCAACTACGGCGCCAATTACCTGGACAATCTTTCCAGTAGCCATCTTaaa                                            >  2:711631/1‑81 (MQ=255)                                   cgcgGTACGGCATCCTGAGGGAATTCGACGTCAACTACGGCGCCAATTACCTGGACAATCTTTCCAGTAGCCATCTTAAAt                                           <  2:711578/81‑1 (MQ=255)                                        tACGGCATCCTGAGGGAATTCGACGTCAACTACGGCGCCAATTACCTGGACAATCTTTCCAGTAGCCATCTTAAATCctct                                      <  1:711580/81‑1 (MQ=255)                                            gCATCCTGAGGGAATTCGACGTCAACTACGGCGCCAATTACCTGGACAATCTTTCCAGTAGCCATCTTAAATCCTCTACg                                   >  2:711632/1‑80 (MQ=255)                                                 cTGAGGGAATTCGACGTCAACTACGGCGCCAATTACCTGGACAATCTTTCCAGTAGCCATCTTAAATCCTCTACGAAATaa                             <  1:711582/81‑1 (MQ=255)                                                 cTGAGGGAATTCGACGTCAACTACGGCGCCAATTACCTGGACAATCTTTCCAGTAGCCATCTTAAATCCTCTACGAAATaa                             <  2:711580/81‑1 (MQ=255)                                                      ggAATTCGACGTCAACTACGGCGCCAATTACCTGGACAATCTTTCCAGTAGCCATCTTAAATCCTCTACGAAATAACCTGt                        <  2:711582/81‑1 (MQ=255)                                                       gAATTCGACGTCAACTACGGCGCCAATTACCTGGACAATCTTTCCAGTAGCCATCTTAAATCCTCTACGAAATAACCTGtt                       >  1:711626/1‑81 (MQ=255)                                                        aaTTCGACGTCAACTACGGCGCCAATTACCTGGACAATCTTTCCAGTAGCCATCTTAAATCCTCTACGAAATAACCTGttt                      <  2:711583/81‑1 (MQ=255)                                                            cGACGTCAACTACGGCGCCAATTACCTGGACAATCTTTCCAGTAGCCATCTTAAATCCTCTACGAAATAACCTGTTTAAAc                  >  1:711630/1‑81 (MQ=255)                                                             gACGTCAACTACGGCGCCAATTACCTGGACAATCTTTCCAGTAGCCATCTTAAATCCTCTACGAAATAACCTGTTTAAAcc                 <  2:711584/81‑1 (MQ=255)                                                             gACGTCAACTACGGCGCCAATTACCTGGACAATCTTTCCAGTAGCCATCTTAAATCCTCTACGAAATAACCTGTTTAAAcc                 <  1:711586/81‑1 (MQ=255)                                                                gTCAACTACGGCGCCAATTACCTGGACAATCTTTCCAGTAGCCATCTTAAATCCTCTACGAAATAACCTGTTTAAACcgcg              >  1:711627/1‑81 (MQ=255)                                                                 tCAACTACGGCGCCAATTACCTGGACAATCTTTCCAGTAGCCATCTTAAATCCTCTACGAAATAACCTGTTTAAACCgcgg             >  2:711627/1‑81 (MQ=255)                                                                   aaCTACGGCGCCAATTACCTGGACAATCTTTCCAGTAGCCATCTTAAATCCTCTACGAAATAACCTGTTTAAACCgcggc            <  2:711586/80‑1 (MQ=255)                                                                    aCTACGGCGCCAATTACCTGGACAATCTTTCCAGTAGCCATCTTAAATCCTCTACGAAATAACCTGTTTAAACCGcggcgg          >  2:711637/1‑81 (MQ=255)                                                                     cTACGGCGCCAATTACCTGGACAATCTTTCCAGTAGCCATCTTAAATCCTCTACGAAATAACCTGTTTAAACCGCggcggc         <  1:711588/81‑1 (MQ=255)                                                                     cTACGGCGCCAATTACCTGGACAATCTTTCCAGTAGCCATCTTAAATCCTCTACGAAATAACCTGTTTAAACCGCggcggc         <  1:711589/81‑1 (MQ=255)                                                                         ggCGCCAATTACCTGGACAATCTTTCCAGTAGCCATCTTAAATCCTCTACGAAATAACCTGTTTAAACCGCGGCGGccccc     <  2:711589/81‑1 (MQ=255)                                                                         ggCGCCAATTACCTGGACAATCTTTCCAGTAGCCATCTTAAATCCTCTACGAAATAACCTGTTTAAACCGCGGCGGccccc     >  1:711652/1‑81 (MQ=255)                                                                          gcgcCAATTACCTGGACAATCTTTCCAGTAGCCATCTTAAATCCTCTACGAAATAACCTGTTTAAACCGCGGCGGccccc     >  2:711649/1‑80 (MQ=255)                                                                           cgcCAATTACCTGGACAATCTTTCCAGTAGCCATCTTAAATCCTCTACGAAATAACCTGTTTAAACCGCGGCGGccccc     >  1:711636/1‑79 (MQ=255)                                                                            gcCAATTACCTGGACAATCTATCCAGTAGCCATCTTAAATCCTCTGCGAAATAACCTGTTTAAACCGCGGCGGCCCCCgag  <  2:711590/81‑1 (MQ=255)                                                                               |                                                                              ATTACCATTTTGCACCTCAAGAGCATCGTACACGCGCGGTACGGCATCCTGAGGGAATTCGACGTCAACTACGGCGCCGATTACCTGGACAATCTTTCCAGTAGCCATCTTAAATCCTCTACGAAATAACCTGTTTAAACCGCGGCGGCCCCCGAG  >  NC_000913/3917268‑3917423

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