Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,921,264	T→C	Y263C (TAT→TGT)	atpB ←	ATP synthase Fo complex subunit a

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,921,264	0	T	C	100.0%	129.3 / NA	36	Y263C (TAT→TGT)	atpB	ATP synthase Fo complex subunit a
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (21/15);  total (21/15)

CATGACAGTCTCCAGTTTGTTTCAGTTAAAACGTAGTAGTGTTGGTAAATTAATGTTCTTCAGACGCCATCGACAGATAGACGATCGTCAGAACCATGAAGATGAAGGCTTGCAGCGTAATGATCAGGATGTGGAAAATGGCCCACGGCACATTC  >  NC_000913/3921187‑3921341                                                                              |                                                                              cATGACAGTCTCCAGTTTGTTTCAGTTAAAACGTAGTAGTGTTGGTAAATTAATGTTCTTCAGACGCCATCGAcagacaga                                                                            >  2:765756/1‑81 (MQ=255) cATGACAGTCTCCAGTTTGTTTCAGTTAAAACGTAGTAGTGTTGGTAAATTAATGTTCTTCAGACGCCATCGAcagacaga                                                                            >  2:765734/1‑81 (MQ=255)     aCAGTCTCCAGTTTGTTTCAGTTAAAACGTAGTAGTGTTGGTAAATTAATGTTCTTCAGACGCCATCGACAGACAGACGAt                                                                        >  1:765742/1‑81 (MQ=255)      cAGTCTCCAGTTTGTTTCAGTTAAAACGTAGTAGTGTTGGTAAATTAATGTTCTTCAGACGCCATCGACAGACAGACGATc                                                                       <  2:765702/81‑1 (MQ=255)             cAGTTTGTTTCAGTTAAAACGTAGTAGTGTTGGTAAATTAATGTTCTTCAGACGCCATCGACAGACAGACGAGCGTCa                                                                   <  1:765703/78‑1 (MQ=255)              aGTTTGTTTCAGTTAAAACGTAGTAGTGTTGGTAAATTAATGTTCTTCAGACGCCATCGACAGACAGACGATCGTCAGAAc                                                               <  1:765704/81‑1 (MQ=255)                     ttCAGTTAAAACGTAGTAGTGTTGGTAAATTAATGTTCTTCAGACGCCATCGACAGACAGACGATCGTCAGAACCatg                                                           >  2:765743/1‑78 (MQ=255)                      tCAGTTAAAACGTAGTAGTGTTGGTAAATTGATGTTCTTCAGACGCCATCGACAGACAGACGATCGTCAGAACCatgaag                                                        <  2:765705/80‑1 (MQ=255)                        aGTTAAAACGTAGTAGTGTTGGTAAATTAATGTTCTTCAGACGCCATCGACAGACAGACGATCGTCAGAACCatgaagat                                                      >  1:765746/1‑80 (MQ=255)                           tAAAACGTAGTAGTGTTGGTAAATTAATGTTCTTCAGACGCCATCGACAGACAGACGATCGTCAGAACCatgaagatgaag                                                  >  2:765767/1‑81 (MQ=255)                            aaaaCGTAGTAGTGTTGGTAAATTAATGTTCTTCAGACGCCATCGACAGACAGACGATCGTCAGAACCATGAAGATGAAgg                                                 >  1:765760/1‑81 (MQ=255)                             aaaCGTAGTAGTGTTGGTAAATTAATGTTCTTCAGACGCCATCGACAGACAGACGATCGTCAGAACCATGAAGATGAAgg                                                 <  2:765706/80‑1 (MQ=255)                              aaCGTAGTAGTGTTGGTAAATTAATGTTCTTCAGACGCCATCGACAGACAGACGATCGTCAGAACCATGAAGATGAAGGCt                                               <  1:765707/81‑1 (MQ=255)                               aCGTAGTAGTGTTGGTAAATTAATGTTCTTCAGACGCCATCGACAGACAGACGATCGTCAGAACCATGAAGATGAAGGCtt                                              <  2:765708/81‑1 (MQ=255)                                cGTAGTAGTGTTGGTAAATTAATGTTCTTCAGACGCCATCGACAGACAGACGATCGTCAGAACCATGAAGATGAAGGCTTg                                             <  1:765709/81‑1 (MQ=255)                                 gtagtaGTGTTGGTAAATTAATGTTCTTCAGACGCCATCGACAGACAGACGATCGTCAGAACCATGAAGATGAAGGCTTGc                                            <  2:765710/81‑1 (MQ=255)                                  tagtagTGTTGGTAAATTAATGTTCTTCAGACGCCATCGACAGACAGACGATCGTCAGAACCATGAAGATGAAGGCTTg                                             >  2:765758/1‑79 (MQ=255)                                   agtagaGTTGGTAAATTAATGTTCTTCAGACGCCATCGACAGACAGACGATCGTCAGAACCATGAAGATGAAGGCTTGCAg                                          >  2:765754/1‑81 (MQ=255)                                         gtTGGTAAATTAATGTTCTTCAGACGCCATCGACAGACAGACGATCGTCAGAACCATGAAGATGAAGGCTTGCAGCGTAAt                                    >  1:765749/1‑81 (MQ=255)                                             gTAAATTAATGTTCTTCAGACGCCATCGACAGACAGACGATCGTCAGAACCATGAAGATGAAGGCTTGCAGCGTAATGAt                                 >  2:765772/1‑80 (MQ=255)                                                aaTTAATGTTCTTCAGACGCCATCGACAGACAGACGATCGTCAGAACCATGAAGATGAAGGCTTGCAGCGTAATGATCAg                              >  1:765748/1‑80 (MQ=255)                                                    aaTGTTCTTCAGACGCCATCGACAGACAGACGATCGTCAGAACCATGAAGATGAAGGCTTGCAGCGTAATGATCAGGAtg                          <  2:765711/80‑1 (MQ=255)                                                     aTGTTCTTCAGACGCCATCGACAGACAGACGATCGTCAGAACCATGAAGATGAAGGCTTGCAGCGTAATGATCAGGAtgtg                        <  1:765712/81‑1 (MQ=255)                                                        ttcttcAGACGCCATCGACAGACAGACGATCGTCAGAACCATGAAGATGAAGGCTTGCAGCGTAATGATCAGGATGTGGaa                     <  2:765713/81‑1 (MQ=255)                                                         tcttcAGACGCCATCGACAGACAGACGATCGTCAGAACCATGAAGATGAAGGCTTGCAGCGTAATGATCAGGATGTGGaa                     >  1:765755/1‑80 (MQ=255)                                                           ttcAGACGCCATCGACAGACAGACGATCGTCAGAACCATGAAGATGAAGGCTTGCAGCGTAATGATCAGGATGTGGAAAAt                  >  2:765770/1‑81 (MQ=255)                                                            tcAGACGCCATCGACAGACAGACGATCGTCAGAACCATGAAGATGAAGGCTTGCAGCGTAATGATCAGGATGTGGAAAATg                 >  1:765757/1‑81 (MQ=255)                                                            tcAGACGCCATCGACAGACAGACGATCGTCAGAACCATGAAGATGAAGGCTTGCAGCGTAATGATCAGGATGTGGAAAATg                 >  2:765752/1‑81 (MQ=255)                                                             cAGACGCCATCGACAGACAGACGATCGTCAGAACCATGAAGATGAAGGCTTGCAGCGTAATGATCAGGATGTGGaaaa                   <  1:765714/78‑1 (MQ=255)                                                             cAGACGCCATCGACAGACAGACGATCGTCAGAACCATGAAGATGAAGGCTTGCAGCGTAATGATCAGGATGTGGAAAATgg                >  2:765751/1‑81 (MQ=255)                                                              aGACGCCATCGACAGACAGACGATCGTCAGAACCATGAAGATGAAGGCTTGCAGCGTAATGATCAGGATGTGGAAAATGGc               >  2:765865/1‑81 (MQ=255)                                                                aCGCCATCGACAGACAGACGATCGTCAGAACCATGAAGATGAAGGCTTGCAGCGTAATGATCAGGATGTGGAAAATGGcc              >  2:765775/1‑80 (MQ=255)                                                                 cGCCATCGACAGACAGACGATCGTCAGAACCATGAAGATGAAGGCTTGCAGCGTAATGATCAGGATGTGGAAAATGGCCCa            <  1:765715/81‑1 (MQ=255)                                                                  gCCATCGACAGACAGACGATCGTCAGAACCATGAAGATGAAGGCTTGCAGCGTAATGATCAGGATGTGGAAAATGGCCCAc           >  1:765776/1‑81 (MQ=255)                                                                  gCCATCGACAGACAGACGATCGTCAGAACCATGAAGATGAAGGCTTGCAGCGTAATGATCAGGATGTGGAAAATGGCCCAc           >  1:765781/1‑81 (MQ=255)                                                                     aTCGACAGACAGACGATCGTCAGAACCATGAAGATGAAGGCTTGCAGCGTAATGATCAGGATGTGGAAAATGGCCCACgg         >  1:765783/1‑80 (MQ=255)                                                                     aTCGACAGACAGACGATCGTCAGAACCATGAAGATGAAGGCTTGCAGCGTAATGATCAGGATGTGGAAAATGGCCCACgg         >  1:765762/1‑80 (MQ=255)                                                                     aTCGACAGACAGACGATCGTCAGAACCATGAAGATGAAGGCTTGCAGCGTAATGATCAGGATGTGGAAAATGGCCCACGGc        <  2:765716/81‑1 (MQ=255)                                                                        gacagacaGACGATCGTCAGAACCATGAAGATGAAGGCTTGCAGCGTAATGATCAGGATGTGGAAAATGGCCCACGGcaca     >  2:765788/1‑81 (MQ=255)                                                                         acagacagACGATCGTCAGAACCATGAAGATGAAGGCTTGCAGCGTAATGATCAGGATGTGGAAAATGGCCCACGGCACAt    <  2:765717/81‑1 (MQ=255)                                                                          cagacagaCGATCGTCAGAACCATGAAGATGAAGGCTTGCAGCGTAATGATCAGGATGTGGAAAATGGCCCACGGCACAtt   <  2:765718/81‑1 (MQ=255)                                                                           agacagacGATCGTCAGAACCATGAAGATGAAGGCTTGCAGCGTAATGATCAGGATGTGGAAAATGGCCCACGGCACATTc  <  2:765719/81‑1 (MQ=255)                                                                              |                                                                              CATGACAGTCTCCAGTTTGTTTCAGTTAAAACGTAGTAGTGTTGGTAAATTAATGTTCTTCAGACGCCATCGACAGATAGACGATCGTCAGAACCATGAAGATGAAGGCTTGCAGCGTAATGATCAGGATGTGGAAAATGGCCCACGGCACATTC  >  NC_000913/3921187‑3921341

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