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_0009131,354,5050CT100.0% 123.9 / NA 39S5N (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 (20/19);  total (20/19)

CGCAAACCAGCCGCCAAAAATACACAGTACAGCCAGCCCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAAT  >  NC_000913/1354369‑1354621
                                                                                                                                        |                                                                                                                    
cgcAAACCAGCCGCCAAAAATACACAGTACAGCCAGCCCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTa                                                                                                                    >  1:944727/1‑139 (MQ=255)
        aGCCGCCAAAAATACACAGTACAGCCAGCCCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAgg                                                                                                            <  2:399079/139‑1 (MQ=255)
           cgccAAAAATACACAGTACAGCCAGCCCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCAt                                                                                                         <  1:910977/139‑1 (MQ=255)
             ccAAAAATACACAGTACAGCCAGCCCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATAc                                                                                                       <  2:461305/139‑1 (MQ=255)
                             cagccagcCCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCTTATGTTTCa                                                                                       <  1:610486/139‑1 (MQ=255)
                             cagccagcCCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATAGTAAGGCATACCATTCCTTATGTTTCa                                                                                       <  1:1020192/139‑1 (MQ=255)
                              agccagccCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCTTATGTTTCa                                                                                       <  2:51161/138‑1 (MQ=255)
                              agccagccCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCTTATGTTTCa                                                                                       >  1:51161/1‑138 (MQ=255)
                                     cccGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTa                                                                               >  1:112355/1‑139 (MQ=255)
                                     cccGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTa                                                                               >  2:76758/1‑139 (MQ=255)
                                        gcgcAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGcc                                                                            >  1:976939/1‑139 (MQ=255)
                                              ccGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCa                                                                      >  2:940928/1‑139 (MQ=255)
                                                    aGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCaaa                                                                <  1:940928/139‑1 (MQ=255)
                                                      gCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAat                                                              >  1:427967/1‑139 (MQ=255)
                                                               cGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAAt                                                     <  1:895749/139‑1 (MQ=255)
                                                                     ggCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGtt                                               <  1:76758/139‑1 (MQ=255)
                                                                         tCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACa                                           >  2:514825/1‑139 (MQ=255)
                                                                         tCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACa                                           >  2:91821/1‑139 (MQ=255)
                                                                         tCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACa                                           >  1:1068232/1‑139 (MQ=255)
                                                                            cTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACaata                                        <  2:1068232/139‑1 (MQ=255)
                                                                                 aaaTTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGAc                                   <  2:944727/139‑1 (MQ=255)
                                                                                       gcgcCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCTTATGTTTCAGAggg                                                                                  <  1:51118/86‑1 (MQ=255)
                                                                                       gcgcCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCTTATGTTTCAGAggg                                                                                  >  2:51118/1‑86 (MQ=255)
                                                                                       gcgcCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGa                             >  2:142196/1‑139 (MQ=255)
                                                                                             ggCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCcaca                       <  2:427967/139‑1 (MQ=255)
                                                                                                  tACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCa                  >  1:297500/1‑139 (MQ=255)
                                                                                                     gCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCa                                                                                                          >  1:457281/1‑48 (MQ=255)
                                                                                                     gCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCa                                                                                                          <  2:457281/48‑1 (MQ=255)
                                                                                                        gCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTAAGCCACACACCATCACTc            >  1:1010/1‑139 (MQ=255)
                                                                                                         cGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATACCATTcc                                                                                                 >  2:990801/1‑53 (MQ=255)
                                                                                                         cGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATACCATTcc                                                                                                 <  1:990801/53‑1 (MQ=255)
                                                                                                         cGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCTTATGTTTCa                                                                                       <  1:933450/63‑1 (MQ=255)
                                                                                                         cGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCTTATGTTTCa                                                                                       >  2:933450/1‑63 (MQ=255)
                                                                                                          gTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCcg          <  1:91821/139‑1 (MQ=255)
                                                                                                            gCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCcggc        >  1:58533/1‑139 (MQ=255)
                                                                                                             cccGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCcggcg       >  1:510020/1‑139 (MQ=255)
                                                                                                                  cgggcgCTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAAt  <  1:142196/139‑1 (MQ=255)
                                                                                                                           ttCGCTGTATACGTTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTTGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGaaa   <  2:770061/129‑1 (MQ=255)
                                                                                                                           ttCGCTGTATACGTTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGaaa   >  1:770061/1‑129 (MQ=255)
                                                                                                                                        |                                                                                                                    
CGCAAACCAGCCGCCAAAAATACACAGTACAGCCAGCCCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAAT  >  NC_000913/1354369‑1354621

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