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%	169.5 / NA	52	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 (24/28);  total (24/28)

CGCAAACCAGCCGCCAAAAATACACAGTACAGCCAGCCCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCC  >  NC_000913/1354369‑1354634                                                                                                                                         |                                                                                                                                  cgcAAACCAGCCGCCAAAAATACACAGTACAGCCAGCCCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTa                                                                                                                                 >  1:71080/1‑139 (MQ=255)      aCCAGCCGCCAAAAATACACAGTACAGCCAGCCCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTa                                                                                                                            <  2:661969/139‑1 (MQ=255)      aCCAGCCGCCAAAAATACACAGTACAGCCAGCCCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTa                                                                                                                            >  1:207640/1‑139 (MQ=255)                aaaaaTACACAGTACAGCCAGCCCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATACCa                                                                                                                  <  2:643947/139‑1 (MQ=255)                              cagccagcCCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCTTATGTTTCa                                                                                                    <  1:1252871/139‑1 (MQ=255)                                 ccagccCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCTTATGTTTCa                                                                                                    >  2:313319/1‑136 (MQ=255)                                 ccagccCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCTTATGTTTCa                                                                                                    <  1:313319/136‑1 (MQ=255)                                      cccGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTa                                                                                            >  1:1030753/1‑139 (MQ=255)                                           gcAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCAt                                                                                       <  2:323289/139‑1 (MQ=255)                                                       gCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAg                                                                                                                          <  2:802136/92‑1 (MQ=255)                                                       gCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAg                                                                                                                          >  1:802136/1‑92 (MQ=255)                                                       gCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCa                                                                               <  2:907123/135‑1 (MQ=255)                                                       gCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCa                                                                               >  1:907123/1‑135 (MQ=255)                                                                       gCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTa                                                           <  2:598352/139‑1 (MQ=255)                                                                        cGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTaa                                                          <  2:591661/139‑1 (MQ=255)                                                                          tCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCa                                                                                   >  1:365210/1‑112 (MQ=255)                                                                          tCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCa                                                                                   <  2:365210/112‑1 (MQ=255)                                                                                 aaaaTTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCg                                                                                                                              >  2:881892/1‑62 (MQ=255)                                                                                 aaaaTTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCg                                                                                                                              <  1:881892/62‑1 (MQ=255)                                                                                 aaaaTTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGa                                                 >  2:1367822/1‑139 (MQ=255)                                                                                    aTTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCa                                              <  1:526501/139‑1 (MQ=255)                                                                                        gcgcCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCTTATGTTTCAGAggg                                                                                               <  2:1488665/86‑1 (MQ=255)                                                                                        gcgcCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCTTATGTTTCAGAggg                                                                                               >  1:1488665/1‑86 (MQ=255)                                                                                        gcgcCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCaa                                                                              >  1:29009/1‑103 (MQ=255)                                                                                        gcgcCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCaa                                                                              <  2:29009/103‑1 (MQ=255)                                                                                         cgcCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAg                                         <  2:207640/139‑1 (MQ=255)                                                                                                        aGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCaa                                                                              >  1:75412/1‑87 (MQ=255)                                                                                                        aGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCaa                                                                              <  2:75412/87‑1 (MQ=255)                                                                                                        aGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACt                          >  1:351516/1‑139 (MQ=255)                                                                                                             gCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAg                                                                                                                          >  1:454563/1‑38 (MQ=255)                                                                                                             gCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAg                                                                                                                          <  2:454563/38‑1 (MQ=255)                                                                                                             gCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCAt                                                                                                                      <  1:1013434/42‑1 (MQ=255)                                                                                                             gCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCAt                                                                                                                      >  2:1013434/1‑42 (MQ=255)                                                                                                             gCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATAc                                                                                                                    >  1:556985/1‑44 (MQ=255)                                                                                                             gCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATAc                                                                                                                    <  2:556985/44‑1 (MQ=255)                                                                                                             gCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCTTAt                                                                                                          <  2:1421071/54‑1 (MQ=255)                                                                                                             gCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCTTAt                                                                                                          >  1:1421071/1‑54 (MQ=255)                                                                                                             gCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCa                                                                                        >  2:644071/1‑72 (MQ=255)                                                                                                             gCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCa                                                                                        <  1:644071/72‑1 (MQ=255)                                                                                                             gCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACa                                                        >  1:778734/1‑104 (MQ=255)                                                                                                             gCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACa                                                        <  2:778734/104‑1 (MQ=255)                                                                                                             gCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGt                                            >  1:1464606/1‑116 (MQ=255)                                                                                                             gCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGt                                            <  2:1464606/116‑1 (MQ=255)                                                                                                             gCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCAt                              <  1:760856/130‑1 (MQ=255)                                                                                                             gCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCAt                              >  2:760856/1‑130 (MQ=255)                                                                                                               ccGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACTCAAAatat                                                                         >  2:1185030/1‑85 (MQ=255)                                                                                                               ccGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAatat                                                                         <  1:1185030/85‑1 (MQ=255)                                                                                                                       cgcTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAAt                                                                  >  1:1238929/1‑84 (MQ=255)                                                                                                                       cgcTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAAt                                                                  <  2:1238929/84‑1 (MQ=255)                                                                                                                                cTGTATACGTTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATcc  <  2:1030753/139‑1 (MQ=255)                                                                                                                                   tataCGTTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTa                                                           <  1:749372/79‑1 (MQ=255)                                                                                                                                   tataCGTTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTa                                                           >  2:749372/1‑79 (MQ=255)                                                                                                                                         |                                                                                                                                  CGCAAACCAGCCGCCAAAAATACACAGTACAGCCAGCCCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCC  >  NC_000913/1354369‑1354634

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