Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,354,505	C→T	S5N (AGC→AAC)  *322* (TAG→TAA)	sapC ← sapB ←	antimicrobial peptide transport ABC transporter permease antimicrobial peptide transport ABC transporter permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,354,505	0	C	T	100.0%	137.7 / NA	44	S5N (AGC→AAC)  *322* (TAG→TAA)	sapC sapB	antimicrobial peptide transport ABC transporter permease antimicrobial peptide transport ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (21/23);  total (21/23)

CGCAAACCAGCCGCCAAAAATACACAGTACAGCCAGCCCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCT  >  NC_000913/1354369‑1354625                                                                                                                                         |                                                                                                                         cgcAAACCAGCCGCCAAAAATACACAGTACAGCCAGCCCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTa                                                                                                                        >  2:716235/1‑139 (MQ=255)   cAAACCAGCCGCCAAAAATACACAGTACAGCCAGCCCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATc                                                                                                                      <  1:959566/139‑1 (MQ=255)    aaaCCAGCCGCCAAAAATACACAGTACAGCCAGCCCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCg                                                                                                                     >  2:1602102/1‑139 (MQ=255)              ccAAAAATACACAGTACAGCCAGCCCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATAc                                                                                                           <  1:716235/139‑1 (MQ=255)                      acacAGTACAGCCAGCCCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCtt                                                                                                   <  1:586735/139‑1 (MQ=255)                              cagccagcCCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCTTATGTTTCa                                                                                           <  2:340323/139‑1 (MQ=255)                                    gccCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATAc                                                                                                           >  2:807516/1‑117 (MQ=255)                                    gccCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATAc                                                                                                           <  1:807516/117‑1 (MQ=255)                                    gccCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGt                                                                                     <  1:1602102/139‑1 (MQ=255)                                          cgcAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGtt                                                                                    <  2:157628/134‑1 (MQ=255)                                          cgcAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGtt                                                                                    >  1:157628/1‑134 (MQ=255)                                               ccGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGcc                                                                          >  1:935397/1‑138 (MQ=255)                                                      ggCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCaaaa                                                                   <  1:118769/139‑1 (MQ=255)                                                         cGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAatat                                                                >  2:853759/1‑139 (MQ=255)                                                         cGACCATCGCAGAGCAGTAGCTATAAAATTTGCGCCAGGCGGTACGCCGCGTTCCCGGCGGGCGATTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAatat                                                                >  2:1101148/1‑139 (MQ=255)                                                           aCCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCtt                                                                                                   >  1:1263096/1‑102 (MQ=255)                                                           aCCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCtt                                                                                                   <  2:1263096/102‑1 (MQ=255)                                                                cGCAGATGCGTCACTATAAAACTTCCGCCAGGCGGTACGAAGCGTGCCCGGCGGGCGGTTGTGGCAGTAGACGTTATCGTACAGCATGCTACAGCTTATGGTGCTGGGACATCGGCATGACACCCAAAATATCaaagat                                                         <  1:1101148/139‑6 (MQ=255)                                                                          tCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTAt                                                                                                                       <  1:294079/67‑1 (MQ=255)                                                                          tCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTAt                                                                                                                       >  2:294079/1‑67 (MQ=255)                                                                          tCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATAc                                                                                                           <  1:1623278/79‑1 (MQ=255)                                                                          tCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATAc                                                                                                           >  2:1623278/1‑79 (MQ=255)                                                                          tCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCAcc                                                                        <  2:1329833/114‑1 (MQ=255)                                                                          tCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCAcc                                                                        >  1:1329833/1‑114 (MQ=255)                                                                                    aTTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCa                                     <  1:853759/139‑1 (MQ=255)                                                                                      ttGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGt                                   >  2:261501/1‑139 (MQ=255)                                                                                       tGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTaa                                                 >  1:1247560/1‑124 (MQ=255)                                                                                       tGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTaa                                                 <  2:1247560/124‑1 (MQ=255)                                                                                                    aCGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCAt                     <  1:29225/139‑1 (MQ=255)                                                                                                        aGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCTTAt                                                                                                 >  1:1468028/1‑59 (MQ=255)                                                                                                        aGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCTTAt                                                                                                 <  2:1468028/59‑1 (MQ=255)                                                                                                         gCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCTTATGTTTCAGAggg                                                                                      <  2:1226979/69‑1 (MQ=255)                                                                                                         gCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCTTATGTTTCAGAggg                                                                                      >  1:1226979/1‑69 (MQ=255)                                                                                                          cGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCAc                                                      <  2:310395/100‑1 (MQ=255)                                                                                                          cGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCAc                                                      >  1:310395/1‑100 (MQ=255)                                                                                                             gCCGGGCGGGCGCTTTTCGCAGTATACGTTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCAGACACCATCACTCcggc            >  2:1133517/1‑139 (MQ=255)                                                                                                             gCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTa                                                                                                                   <  2:455502/36‑1 (MQ=255)                                                                                                             gCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTa                                                                                                                   >  1:455502/1‑36 (MQ=255)                                                                                                             gCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTa                                                  <  2:1478983/101‑1 (MQ=255)                                                                                                             gCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTa                                                  >  1:1478983/1‑101 (MQ=255)                                                                                                                     ggcgCTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATgg    <  1:261501/139‑1 (MQ=255)                                                                                                                       cgcTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCt  >  1:1048459/1‑139 (MQ=255)                                                                                                                            ttCGCTGTATACGTTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGaa                                                           <  2:329557/77‑1 (MQ=255)                                                                                                                            ttCGCTGTATACGTTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGaa                                                           >  1:329557/1‑77 (MQ=255)                                                                                                                                         |                                                                                                                         CGCAAACCAGCCGCCAAAAATACACAGTACAGCCAGCCCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCT  >  NC_000913/1354369‑1354625

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