Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,337,996	C→T	G156S (GGC→AGC)	mlaD ←	intermembrane phospholipid transport system, substrate binding protein MlaD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,337,996	0	C	T	100.0%	102.9 / NA	30	G156S (GGC→AGC)	mlaD	intermembrane phospholipid transport system, substrate binding protein MlaD
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (18/12);  total (18/12)

GTACCCACAGGTTCAGTGGTTTCATTATTACCTGGCGCAGCAGCTGGCGCATCGCCACTATTCTTATTGTCATCGCCTTTACTACCGTAAAGGAACTGACCAATGAGATCTTCCAGCACCATCGCAGACTTAGTGTCCTGAATTGTATCGCCAT  >  NC_000913/3337920‑3338073                                                                             |                                                                              gTACCCACAGGTTCAGTGGTTTCATTATTACCTGGCGCAGCAGCTGGCGCATCGCCACTATTCTTATTGTCATCGCtttt                                                                            <  2:586730/80‑1 (MQ=255) gTACCCACAGGTTCAGTGGTTTCATTATTACCTGGCGCAGCAGCTGGCGCATCGCCACTATTCTTATTGTCATCGCTTTta                                                                           >  2:586772/1‑81 (MQ=255)      cacaGGTTCAGTGGTTTCATTATTACCTGGCGCAGCAGCTGGCGCATCGCCACTATTCTTATTGTCATCGCTTTTACTAcc                                                                      <  2:586731/81‑1 (MQ=255)       acaGGTTCAGTGGTTTCATTATTACCTGGCGCAGCAGCTGGCGCATCGCCACTATTCTTATTGTCATCGCTTTTACTACCg                                                                     <  1:586733/81‑1 (MQ=255)       acaGGTTCAGTGGTTTCATTATTACCTGGCGCAGCAGCTGGCGCATCGCCACTATTCTTATTGTCATCGCTTTTACTACCg                                                                     <  2:586734/81‑1 (MQ=255)       acaGGTTCAGTGGTTTCATTATTACCTGGCGCAGCAGCTGGCGCATCGCCACTATTCTTATTGTCATCGCTTTTACTACCg                                                                     <  1:586732/81‑1 (MQ=255)             tCAGTGGTTTCATTATTACCTGGCGCAGCAGCTGGCGCATCGCCACTATTCTTATTGTCATCGCTTTTACTACCGTAAAg                                                                >  2:586773/1‑80 (MQ=255)              cAGTGGTTTCATTATTACCTGGCGCAGCAGCTGGCGCATCGCCACTATTCTTATTGTCATCGCTTTTACTACCGTAAAgg                                                               >  1:586768/1‑80 (MQ=255)               aGTGGTTTCATTATTACCTGGCGCAGCAGCTGGCGCATCGCCACTATTCTTATTGTCATCGCTTTTACTACCGTAAAGGaa                                                             >  2:586764/1‑81 (MQ=255)                        attattACCTGGCGCAGCAGCTGGCGCATCGCCACTATTCTTATTGTCATCGCTTTTACTACCGTAAAGGAACTGACCAAt                                                    >  1:586794/1‑81 (MQ=255)                         ttattaCCTGGCGCAGCAGCTGGCGCATCGCCACTATTCTTATTGTCATCGCTTTTACTACCGTAAAGGAACTGACCAATg                                                   <  1:586735/81‑1 (MQ=255)                         ttattaCCTGGCGCAGCAGCTGGCGCATCGCCACTATTCTTATTGTCATCGCTTTTACTACCGTAAAGGAACTGACCAATg                                                   <  1:586736/81‑1 (MQ=255)                           attaCCTGGCGCAGCAGCTGGCGCATCGCCACTATTCTTATTGTCATCGCTTTTACTACCGTAAAGGAACTGACCAAt                                                    >  2:586774/1‑78 (MQ=255)                            ttaCCTGGCGCAGCAGCTGGCGCATCGCCACTATTCTTATTGTCATCGCTTTTACTACCGTAAAGGAACTGACCAAt                                                    >  2:586780/1‑77 (MQ=255)                             taCCTGGCGCAGCAGCTGGCGCATCGCCACTATTCTTATTGTCATCGCTCTTACTACCGTAAAGGAACTGACCAAt                                                    >  2:586761/1‑76 (MQ=255)                               ccTGGCGCAGCAGCTGGCGCATCGCCACTATTCTTATTGTCATCGCTTTTACTACCGTAAAGGAACTGACCAATGAGATc                                              <  1:586737/80‑1 (MQ=255)                                  ggCGCAGCAGCTGGCGCATCGCCACTATTCTTATTGTCATCGCTTTTACTACCGTAAAGGAACTGACCAATGAGATCTTc                                           >  1:586792/1‑80 (MQ=255)                                  ggCGCAGCAGCTGGCGCATCGCCACTATTCTTATTGTCATCGCTTTTACTACCGTAAAGGAACTGACCAATGAGATCTTc                                           >  1:586791/1‑80 (MQ=255)                                  ggCGCAGCAGCTGGCGCATCGCCACTATTCTTATTGTCATCGCTTTTACTACCGTAAAGGAACTGACCAATGAGATCTTc                                           >  1:586760/1‑80 (MQ=255)                                   gcgcAGCAGCTGGCGCATCGCCACTATTCTTATTGTCATCGCTTTTACTACCGTAAAGGAACTGACCAATGAGATCTTcc                                          <  1:586738/80‑1 (MQ=255)                                     gcagcaGCTGGCGCATCGCCACTATTCTTATTGTCATCGCTTTTACTACCGTAAAGGAACTGACCAATGAGATCTTCCAg                                        >  2:586779/1‑80 (MQ=255)                                         cagcTGGCGCATCGCCACTATTCTTATTGTCATCGCTTTTACTACCGTAAAGGAACTGACCAATGAGATCTTCCAGCACCa                                   >  1:586782/1‑81 (MQ=255)                                          agcTGGCGCATCGCCACTATTCTTATTGTCATCGCTTTTACTACCGTAAAGGAACTGACCAATGAGATCTTCCAGCACCAt                                  <  2:586739/81‑1 (MQ=255)                                           gcTGGCGCATCGCCACTATTCTTATTGTCATCGCTTTTACTACCGTAAAGGAACTGACCAATGAGATCTTCCAGCACCATc                                 <  2:586740/81‑1 (MQ=255)                                              ggCGCATCGCCACTATTCTTATTGTCATCGCTTTTACTACCGTAAAGGAACTGACCAATGAGATCTTCCAGCACCATCg                                >  2:586783/1‑79 (MQ=255)                                                  cATCGCCACTATTCTTATTGTCATCGCTTTTACTACCGTAAAGGAACTGACCAATGAGATCTTCCAGCACCATCGCAGACt                          <  1:586741/81‑1 (MQ=255)                                                  cATCGCCACTATTCTTATTGTCATCGCTTTTACTACCGTAAAGGAACTGACCAATGAGATCTTCCAGCACCATCGCAGACt                          >  1:586786/1‑81 (MQ=255)                                                   aTCGCCACTATTCTTATTGTCATCGCTTTTACTACCGTAAAGGAACTGACCAATGAGATCTTCCAGCACCATCGCAGACt                          >  1:586787/1‑80 (MQ=255)                                                     cGCCACTATTCTTATTGTCATCGCTTTTACTACCGTAAAGGAACTGACCAATGAGATCTTCCAGCACCATCGCAGACTTAg                       >  2:586784/1‑81 (MQ=255)                                                             ttCTTATTGTCATCGCTTTTACTACCGTAAAGGAACTGACCAATGAGATCTTCCAGCACCATCGCAGACTTAGTGTCCTg                <  2:586742/80‑1 (MQ=255)                                                                          cGCTTTTACTACCGTAAAGGAACTGACCAATGAGATCTTCCAGCACCATCGCAGACTTAGTGTCCTGAATTGTATCGCCAt  >  2:586788/1‑81 (MQ=255)                                                                             |                                                                              GTACCCACAGGTTCAGTGGTTTCATTATTACCTGGCGCAGCAGCTGGCGCATCGCCACTATTCTTATTGTCATCGCCTTTACTACCGTAAAGGAACTGACCAATGAGATCTTCCAGCACCATCGCAGACTTAGTGTCCTGAATTGTATCGCCAT  >  NC_000913/3337920‑3338073

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: ‑