Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	403,036	T→C	D430D (GAT→GAC)	phoA →	alkaline phosphatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	403,036	0	T	C	100.0%	107.4 / NA	30	D430D (GAT→GAC)	phoA	alkaline phosphatase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (12/18);  total (12/18)

AGCTCCGGGCCTCACCCAGGCGCTAAATACCAAAGATGGCGCAGTGATGGTGATGAGTTACGGGAACTCCGAAGAGGATTCACAAGAACATACCGGCAGTCAGTTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGA  >  NC_000913/402958‑403113                                                                               |                                                                              aGCTCCGGGCCTCACCCAGGCGCTAAATACCAAAGATGGCGCAGTGATGGTGATGAGTTACGGGAACTCCGAAGAGGACTc                                                                             <  1:82699/81‑1 (MQ=255) aGCTCCGGGCCTCACCCAGGCGCTAAATACCAAAGATGGCGCAGTGATGGTGATGAGTTACGGGAACTCCGAAGAGGACTc                                                                             <  2:82698/81‑1 (MQ=255)   cTCCGGGCCTCACCCAGGCGCTAAATACCAAAGATGGCGCAGTGATGGTGATGAGTTACGGGAACTCCGAAGAGGACTcac                                                                           <  1:82700/81‑1 (MQ=255)           cTCACCCAGGCGCTAAATACCAAAGATGGCGCAGTGATGGTGATGAGTTACGGGAACTCCGAAGAGGACTCACAAGAACAt                                                                   <  2:82701/81‑1 (MQ=255)                            tACCAAAGATGGCGCAGTGATGGTGATGAGTTACGGGAACTCCGAAGAGGACTCACAAGAACATACCGGCAGTCAGTTGCg                                                  <  1:82702/81‑1 (MQ=255)                            tACCAAAGATGGCGCAGTGATGGTGATGAGTTACGGGAACTCCGAAGAGGACTCACAAGAACATACCGGCAGTCAGTTGCg                                                  >  1:82738/1‑81 (MQ=255)                               cAAAGATGGCGCAGTGATGGTGATGAGTTACGGGAACTCCGAAGAGGACTCACAAGAACATACCGGCAGTCAGTTGCGTAt                                               >  1:82740/1‑81 (MQ=255)                               cAAAGATGGCGAAGTGATGGTGATGAGTTACGGGAACTCCGAAGAGGACTCACAAGAACATACCGGCAGTCAGTTGCGTAt                                               <  1:82703/81‑1 (MQ=255)                                       gcgcAGTGATGGTGATGAGTTACGGGAACTCCGAAGAGGACTCACAAGAACATACCGGCAGTCAGTTGCGTATTgcggcg                                        <  2:82705/80‑1 (MQ=255)                                       gcgcAGTGATGGTGATGAGTTACGGGAACTCCGAAGAGGACTCACAAGAACATACCGGCAGTCAGTTGCGTATTgcggcg                                        <  2:82704/80‑1 (MQ=255)                                       gcgcAGTGATGGTGATGAGTTACGGGAACTCCGAAGAGGACTCACAAGAACATACCGGCAGTCAGTTGCGTATTGCGGCGt                                       <  1:82706/81‑1 (MQ=255)                                         gcAGTGATGGTGATGAGTTACGGGAACTCCGAAGAGGACTCACAAGAACATACCGGCAGTCAGTTGCGTATTGCGGCGTAt                                     >  2:82755/1‑81 (MQ=255)                                         gcAGTGATGGTGATGAGTTACGGGAACTCCGAAGAGGACTCACAAGAACATACCGGCAGTCAGTTGCGTATTGCGGCGTAt                                     >  2:82745/1‑81 (MQ=255)                                         gcAGTGATGGTGATGAGTTACGGGAACTCCGAAGAGGACTCACAAGAACATACCGGCAGTCAGTTGCGTATTGCGGCGTAt                                     <  2:82707/81‑1 (MQ=255)                                            gtgatggtgatgAGTTACGGGAACTCCGAAGAGGACTCACAAGAACATACCGGCAGTCAGTTGCGTATTGCGGCGTATGGc                                  >  1:82743/1‑81 (MQ=255)                                             tgatggtgatgAGTTACGGGAACTCCGAAGAGGACTCACAAGAACATACCGGCAGTCAGTTGCGTATTGCGGCGTATGGcc                                 >  2:82744/1‑81 (MQ=255)                                             tgatggtgatgAGTTACGGGAACTCCGAAGAGGACTCACAAGAACATACCGGCAGTCAGTTGCGTATTGCGGCGTATGGcc                                 <  1:82708/81‑1 (MQ=255)                                                tggtgatgAGTTACGGGAACTCCGAAGAGGACTCACAAGAACATACCGGCAGTCAGTTGCGTATTGCGGCGTATGGCCCGc                              <  1:82709/81‑1 (MQ=255)                                                  gtgatgAGTTACGGGAACTCCGAAGAGGACTCACATGAACATACCGGCAGTCAGTTGCGTATTGCGGCGTATGGCCCGCAt                            >  1:82750/1‑81 (MQ=255)                                                    gatgaGTTACGGGAACTCCGAAGAGGACTCACAAGAACATACCGGCAGTCAGTTGCGTATTGCGGCGTATGGCCCGCATgc                          >  2:82749/1‑81 (MQ=255)                                                       gaGTTACGGGAACTCCGAAGAGGACTCACAAGAACATACCGGCAGTCAGTTGCGTATTGCGGCGTATGGCCCGCATgccgc                       >  2:82742/1‑81 (MQ=255)                                                                   cTCCGAAGAGGACTCACAAGAACATACCGGCGGTCAGTTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTgg           <  2:82710/81‑1 (MQ=255)                                                                    tCCGAAGAGGACTCACAAGAACATACCGGCAGTCAGTTGCGTATTGCGGCGTATGGCCTGCATGCCGCCAATGTTGTTGGa          >  1:82753/1‑81 (MQ=255)                                                                     ccGAAGAGGACTCACAAGAACATACCGGCAGTCAGTTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGAc         >  1:82746/1‑81 (MQ=255)                                                                      cGAAGAGGACTCACAAGAACATACCGGCAGTCAGTTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGAc         <  1:82711/80‑1 (MQ=255)                                                                       gAAGAGGACTCACAAGAACATACCGGCAGTCAGTTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTg       <  1:82712/81‑1 (MQ=255)                                                                       gAAGAGGACTCACAAGAACATACCGGCAGTCAGTTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTg       <  2:82713/81‑1 (MQ=255)                                                                        aaGAGGACTCACAAGAACATACCGGCAGTCAGTTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTga      >  1:82741/1‑81 (MQ=255)                                                                         agagGACTCACAAGAACATACCGGCAGTCAGATGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTgac     <  2:82714/81‑1 (MQ=255)                                                                            ggACTCACAAGAACATACCGGCAGTCAGTTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTgaccga  <  2:82715/81‑1 (MQ=255)                                                                               |                                                                              AGCTCCGGGCCTCACCCAGGCGCTAAATACCAAAGATGGCGCAGTGATGGTGATGAGTTACGGGAACTCCGAAGAGGATTCACAAGAACATACCGGCAGTCAGTTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGA  >  NC_000913/402958‑403113

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