Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,352,163	G→A	R158C (CGT→TGT)	sapF ←	antimicrobial peptide ABC transporter ATPase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,352,163	0	G	A	100.0%	149.3 / NA	40	R158C (CGT→TGT)	sapF	antimicrobial peptide ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (24/16);  total (24/16)

TAAATATACGAAATGCCCTGTTTTTCCTGTAATTCCAGCATCAGATTAATCAACTGCGAACGCATCGACATATCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACGCTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGTGATCCGGTAGCAGCCCAACCATACGCATCGTTTCAATAATCTGTTTGCGACGCTGTTCCGGTTCCAGGTCGGTGTTCAGGCGCAGTGGAAAATCCAGAATTT  >  NC_000913/1352018‑1352309                                                                                                                                                  |                                                                                                                                                   tAAATATACGAAATGCCCTGTTTTTCCTGTAATTCCAGCATCAGATTAATCAACTGCGAACGCATCGACATATCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGtt                                                                                                                                               >  1:73153/1‑151 (MQ=255)           aaaTGCCCTGTTTTTCCTGTAATTCCAGCATCAGATTAATCAACTGCGAACGCATCGACATATCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCggg                                                                                                                                     >  2:885218/1‑151 (MQ=255)            aaTGCCCTGTTTTTCCTGTAATTCCAGCATCAGATTAATCAACTGCGAACGCATCGACATATCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGATTTTGTCCGGGt                                                                                                                                    <  2:633589/151‑1 (MQ=255)                       tttCCTGTAATTCCAGCATCAGATTAATCAACTGCGAACGCATCGACATATCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCGGGTGCCAACatat                                                                                                                          <  1:257643/150‑1 (MQ=255)                          ccTGTAATTCCAGCATCAGATTAATCAACTGCGAACGCATCGACATATCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCGGGTGCCAACATATGCgg                                                                                                                      >  1:129460/1‑151 (MQ=255)                          ccTGTAATTCCAGCATCAGATTAATCAACTGCGAACGCATCGACATATCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCGGGTGCCAACATATGCgg                                                                                                                      <  2:356657/151‑1 (MQ=255)                               aaTTCCAGCATCAGATTAATCAACTGCGAACGCATCGACATATCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCGGGTGCCAACATATGCGGATAg                                                                                                                  <  1:504073/150‑1 (MQ=255)                               aaTTCCAGCATCAGATTAATCAACTGCGAACGCATCGACATATCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCGGGTGCCAACATATGCGGATAGt                                                                                                                 >  1:125874/1‑151 (MQ=255)                                          cAGATTAATCAACTGCGAACGCATCGACATATCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGTg                                                                                                       >  2:810086/1‑150 (MQ=255)                                               tAATCAACTGCGAACGCATCGACATATCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCGGGTGCCAACATATGCGGTTAGTAACTGACGTGATCCgg                                                                                                 >  2:826277/1‑151 (MQ=255)                                                 aTCAACTGCGAACGCATCGACATATCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGTGATCCGGt                                                                                                >  1:529909/1‑150 (MQ=255)                                                      atGCGAACGCATCGACATATCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGTGATCCGGTagcag                                                                                           >  2:241778/2‑150 (MQ=255)                                                         cGAACGCATCGACATATCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGTGATCCGGTAGCAGCCCa                                                                                       <  1:197774/151‑1 (MQ=255)                                                           aaCGCATCGACATATCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGTGATCCGGTAGCAGCCCAAc                                                                                     >  1:816050/1‑151 (MQ=255)                                                                aTCGACATATCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGTGATCCGGTAGCAGCCCAACCATa                                                                                 >  2:344543/1‑150 (MQ=255)                                                                  cGACATATCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGTGATCCGGTAGCAGCCCAACCATACg                                                                               <  1:384767/150‑1 (MQ=255)                                                                  cGACATATCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGTGATCCGGTAGCAGCCCAACCATACGc                                                                              >  1:521130/1‑151 (MQ=255)                                                                          ccAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGTGATCCGGTAGCAGCCCAACCATACGCATCGTTTc                                                                      >  2:771305/1‑151 (MQ=255)                                                                            aGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGTGATCCGGTAGCAGCCCAACCATACGCATCGTTTCaa                                                                    <  1:137418/151‑1 (MQ=255)                                                                                  gCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGTGATCCGGTAGCAGCCCAACCATACGCATCGTTTCAATAATCt                                                              <  2:528054/151‑1 (MQ=255)                                                                                    gAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGTGATCCGGTAGCAGCCCAACCATACGCATCGTTTCAATAATCTGt                                                            <  1:134054/151‑1 (MQ=255)                                                                                       ggCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCGGGTGCCAACAGATGCGGATAGTAACTGACGTGATCCGGTAGCAGCCCAACCATACGCATCGTTTCAATAATCTGTTTg                                                         <  2:256636/151‑1 (MQ=255)                                                                                           tCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGTGATCCGGTAGCAGCCCAACCATACGCATCGTTTCAATAATCTGTTTGCg                                                       <  2:727795/149‑1 (MQ=255)                                                                                            cATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGTGATCCGGTAGCAGCCCAACCATACGCATCGTTTCAATAATCTGTTTGCGAc                                                     <  2:592842/150‑1 (MQ=255)                                                                                               cGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGTGATCCGGTAGCAGCCCAACCATACGCATCGTTTCAATAATCTGTTTGCGACGCTg                                                 >  2:420126/1‑151 (MQ=255)                                                                                                      aatGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGTGATCCGGTAGCAGCCCAACCATACGCATCGTTTCAATAATCTGTTTGCGACGCTGTTCCgg                                           >  2:190868/1‑150 (MQ=255)                                                                                                       atGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGTGATCCGGTAGCAGCCCAACCATACGCATCGTTTCAATAATCTGTTTGCGACGCTGTTCCGGtt                                         <  2:453992/151‑1 (MQ=255)                                                                                                         gACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGTGATCCGGTAGCAGCCCAACCATACGCATCGTTTCAATAATCTGTTTGCGACGCTGTTCCGGTTcc                                       >  2:890747/1‑151 (MQ=255)                                                                                                         gACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGTGATCCGGTAGCAGCCCAACCATACGCATCGTTTCAATAATCTGTTTGCGACGCTGTTCCGGTTcc                                       >  1:699014/1‑151 (MQ=255)                                                                                                         gACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGTGATCCGGTAGCAGCCCAACCATACGCATCGTTTCAATAATCTGTTTGCGACGCTGTTCCGGTTc                                        <  2:521019/150‑1 (MQ=255)                                                                                                            ttttGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGTGATCCGGTAGCAGCCCAACCATACGCATCGTTTCAATAATCTGTTTGCGACGCTGTTCCGGTTCCAgg                                    >  1:819741/1‑151 (MQ=255)                                                                                                              ttGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGTGATCCGGTAGCAGCCCAACCATACGCATCGTTTCAATAATCTGTTTGCGACGCTGTTCCGGTTCCAGGTc                                  >  1:243975/1‑151 (MQ=255)                                                                                                               tGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGTGATCCGGTAGCAGCCCAACCATACGCATCGTTTCAATAATCTGTTTGCGACGCTGTTCCGGTTCCAGGTCg                                 >  1:738776/1‑151 (MQ=255)                                                                                                                gggCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGTGATCCGGTAGCAGCCCAACCATACGCATCGTTTCAATAATCTGTTTGCGACGCTGTTCCGGTTCCAGGTCgg                                >  1:552121/1‑151 (MQ=255)                                                                                                                            aGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGTGATCCGGTAGCAGCCCAACCATACGCATCGTTTCAATAATCTGTTTGCGACGCTGTTCCGGTTCCAGGTCGGTGTTCAGGCGCa                    <  2:353511/151‑1 (MQ=255)                                                                                                                                  cgcgcCAGACCCAGACACTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGTGATCCGGTAGCAGCCCAACCATACGCATCGTTTCAATAATCTGTTTGCGACGCTGTTCCGGTTCCAGGTCGGTGTTCAGGCGCAGTGGa               >  1:596222/1‑150 (MQ=255)                                                                                                                                     gcCAGACCCAGACACTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGTGATCCGGTAGCAGCCCAACCATACGCATCGTTTCAATAATCTGTTTGCGACGCTGTTCCGGTTCCAGGTCGGTGTTCAGGCGCAGTGGAAAAt           >  1:627643/1‑151 (MQ=255)                                                                                                                                       cagacccagacACTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGTGATCCGGTAGCAGCCCAACCATACGCATCGTTTCAATAATCTGTTTGCGACGCTGTTCCGGTTCCAGGTCGGTGTTCAGGCGCAGTGGAAAATcc         <  2:73136/151‑1 (MQ=255)                                                                                                                                        agacccagacACTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGTGATCCGGTAGCAGCCCAACCATACGCATCGTTTCAATAATCTGTTTGCGACGCTGTTCCGGTTCCAGGTCGGTGTTCAGGCGCAGTGGAAAATcc         >  1:616944/1‑150 (MQ=255)                                                                                                                                              agacACTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGTGATCCGGTAGCAGCCCAACCATACGCATCGTTTCAATAATCTGTTTGCGACGCTGTTCCGGTTCCAGGTCGGTGTTCAGGCGCAGTGGAAAATCCAGAAttt  >  1:487710/1‑151 (MQ=255)                                                                                                                                                  |                                                                                                                                                   TAAATATACGAAATGCCCTGTTTTTCCTGTAATTCCAGCATCAGATTAATCAACTGCGAACGCATCGACATATCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACGCTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGTGATCCGGTAGCAGCCCAACCATACGCATCGTTTCAATAATCTGTTTGCGACGCTGTTCCGGTTCCAGGTCGGTGTTCAGGCGCAGTGGAAAATCCAGAATTT  >  NC_000913/1352018‑1352309

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