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%	138.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)

AAACCAGCCGCC‑AAAAATACACAGTACAGCCAGCCCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAAT  >  NC_000913/1354372‑1354632                                                                                                                                       |                                                                                                                                aaaaCAGCCGCCAAAAAATACACAGTACAGCCAGCCCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACGATAAAATTTGCGCCAGGAGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATc                                                                                                                             >  1:1208518/1‑139 (MQ=255)  aaCCAGCCGCCAAAAAATACACAGTACAGCCAGCCCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCg                                                                                                                            <  2:1208518/139‑1 (MQ=255)     cAGCCGCC‑AAAAATACACAGTACAGCCAGCCCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAg                                                                                                                        <  2:126731/139‑1 (MQ=255)        ccgcc‑AAAAATACACAGTACAGCCAGCCCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCa                                                                                                                     >  2:370904/1‑139 (MQ=255)           cc‑AAAAATACACAGTACAGCCAGCCCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATAc                                                                                                                  <  2:534729/139‑1 (MQ=255)                            cagccagcCCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCTTATGTTTCa                                                                                                  <  1:1595726/139‑1 (MQ=255)                                  gccCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCg                                                                                                                            >  2:1387445/1‑107 (MQ=255)                                  gccCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCg                                                                                                                            <  1:1387445/107‑1 (MQ=255)                                           aGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATgg                                                                                   <  1:1418189/139‑1 (MQ=255)                                                     gCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCa                                                                                                                     >  1:208341/1‑95 (MQ=255)                                                     gCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCa                                                                                                                     <  2:208341/95‑1 (MQ=255)                                                               gCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCtt                                                                                                          >  2:367417/1‑96 (MQ=255)                                                               gCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCtt                                                                                                          <  1:367417/96‑1 (MQ=255)                                                                   agGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGt                                                           >  1:1050826/1‑139 (MQ=255)                                                                      cGTCACTATAAAATTTACGCCAGGCGGTAGGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCACAATATCAGAAATCACGTTaa                                                        <  2:1275100/139‑1 (MQ=255)                                                                        tCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAg                                                                                                                        <  1:1775935/73‑1 (MQ=255)                                                                        tCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAg                                                                                                                        >  2:1775935/1‑73 (MQ=255)                                                                        tCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGtt                                                          >  1:321573/1‑135 (MQ=255)                                                                        tCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGtt                                                          <  2:321573/135‑1 (MQ=255)                                                                        tCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACa                                                      >  1:187131/1‑139 (MQ=255)                                                                          aCTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACaat                                                    <  2:1508622/139‑1 (MQ=255)                                                                             ataAAATTTGCGCCAGGCGGTACGTATCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAAGCACGTAAACaataat                                                 >  1:1275100/1‑139 (MQ=255)                                                                              taAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATg                                                <  1:370904/139‑1 (MQ=255)                                                                                aaaTTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGAc                                              <  1:1808068/139‑1 (MQ=255)                                                                                  aTTTGCGCCAGGCGCTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAATGCATACCATTCCTTATGTTTCAGATGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCa                                            <  2:1150743/139‑1 (MQ=255)                                                                                      gcgcCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGa                                        >  1:1603764/1‑139 (MQ=255)                                                                                                      aGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATg                                                                                    >  1:429572/1‑79 (MQ=255)                                                                                                      aGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATg                                                                                    <  2:429572/79‑1 (MQ=255)                                                                                                      aGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACt                        >  2:1000573/1‑139 (MQ=255)                                                                                                       gCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACg                                                            >  2:989586/1‑102 (MQ=255)                                                                                                       gCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACg                                                            <  1:989586/102‑1 (MQ=255)                                                                                                       gCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATAACATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTa                                                         <  2:1851347/105‑1 (MQ=255)                                                                                                       gCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATAACATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTa                                                         >  1:1851347/1‑105 (MQ=255)                                                                                                        cGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATc                                                                                                                             <  2:268156/36‑1 (MQ=38)                                                                                                        cGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATc                                                                                                                             >  1:268156/1‑36 (MQ=255)                                                                                                        cGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCTTATGTTTCa                                                                                                  >  1:616999/1‑63 (MQ=255)                                                                                                        cGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCTTATGTTTCa                                                                                                  <  2:616999/63‑1 (MQ=255)                                                                                                         gTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAATGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCAc                                                             <  1:1274169/99‑1 (MQ=255)                                                                                                         gTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCAc                                                             >  2:1274169/1‑99 (MQ=255)                                                                                                          tGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCcgg                    >  1:1377887/1‑139 (MQ=255)                                                                                                           gCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCAcc                                                                               >  1:61202/1‑79 (MQ=255)                                                                                                           gCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCAcc                                                                               <  2:61202/79‑1 (MQ=255)                                                                                                           gCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCcggc                   >  1:1843405/1‑139 (MQ=255)                                                                                                                            cGCTGTATACGTTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAAt  <  2:1843405/139‑1 (MQ=255)                                                                                                                                       |                                                                                                                                AAACCAGCCGCC‑AAAAATACACAGTACAGCCAGCCCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAAT  >  NC_000913/1354372‑1354632

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