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,5050CT98.4% 192.3 / NA 61S5N (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);  major base T (28/32);  minor base A (1/0);  total (29/32)
Fisher's exact test for biased strand distribution p-value = 4.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCAAACCAGCCGCCAAAAATACACAGTACAGCCAGCCCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCCTGCT  >  NC_000913/1354369‑1354638
                                                                                                                                        |                                                                                                                                     
cgcAAACCAGCCGCCAAAAATACACAGTACAGCCAGCCCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTa                                                                                                                                     >  1:1274149/1‑139 (MQ=255)
           cgccAAAAATACACAGTACAGCCAGCCCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCAt                                                                                                                          <  1:567673/139‑1 (MQ=255)
              cAAAAATACACAGTACAGCCAGCCCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATAcc                                                                                                                       <  1:1151520/139‑1 (MQ=255)
                           tACAGCCAGCCCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAgg                                                                                                                             <  1:1762530/120‑1 (MQ=255)
                           tACAGCCAGCCCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAgg                                                                                                                             >  2:1762530/1‑120 (MQ=255)
                             cagccagcCCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCTTATGTTTCa                                                                                                        <  1:1014555/139‑1 (MQ=255)
                             cagccagcCCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCTTATGTTTCa                                                                                                        >  2:149037/1‑139 (MQ=255)
                             cagccagcCCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCTTATGTTTCa                                                                                                        <  2:1274149/139‑1 (MQ=255)
                                ccagccCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCTTATGTTTCagag                                                                                                     <  1:149037/139‑1 (MQ=255)
                                 cagccCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCTTATGTTTCAGAgg                                                                                                    <  1:1739675/139‑1 (MQ=255)
                                   gccCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCTTATGTTTCAGATGGt                                                                                                  >  2:292179/1‑139 (MQ=255)
                                     cccGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGATATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTa                                                                                                >  2:1377616/1‑139 (MQ=255)
                                          gcAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCAt                                                                                           <  2:1533984/139‑1 (MQ=255)
                                                      gCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAat                                                                               >  1:33771/1‑139 (MQ=255)
                                                      gCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAat                                                                               >  1:1455141/1‑139 (MQ=255)
                                                           ccATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCg                                                                                                                                  >  2:1040071/1‑83 (MQ=255)
                                                           ccATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCg                                                                                                                                  <  1:1040071/83‑1 (MQ=255)
                                                           ccATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATa                                                                                                                         <  1:1559851/92‑1 (MQ=255)
                                                           ccATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATa                                                                                                                         >  2:1559851/1‑92 (MQ=255)
                                                                  agagGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCTTATGTTTCAGATGGTTAGCCATGGCACCCAAAATATCAGAAATCAc                                                                   <  1:292179/139‑1 (MQ=255)
                                                                  agagGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCAc                                                                   <  2:860761/139‑1 (MQ=255)
                                                                       cGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTaa                                                              <  1:1379131/139‑1 (MQ=255)
                                                                         tCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAgg                                                                                                                             <  2:1017015/74‑1 (MQ=255)
                                                                         tCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAgg                                                                                                                             >  1:1017015/1‑74 (MQ=255)
                                                                         tCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAg                                                                                                                              <  2:1648806/73‑1 (MQ=255)
                                                                         tCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAg                                                                                                                              >  1:1648806/1‑73 (MQ=255)
                                                                         tCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCTTATGTTTCa                                                                                                        >  2:1384249/1‑95 (MQ=255)
                                                                         tCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCTTATGTTTCa                                                                                                        <  1:1384249/95‑1 (MQ=255)
                                                                         tCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACa                                                            >  2:97997/1‑139 (MQ=255)
                                                                                aaaaTTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCTTAt                                                                                                              >  1:1482786/1‑82 (MQ=255)
                                                                                aaaaTTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCTTAt                                                                                                              <  2:1482786/82‑1 (MQ=255)
                                                                                     ttGCGCCAGGCGGTACGCAGCGTGCCCGGCTGGCGTTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGt                                                <  1:1377616/139‑1 (MQ=255)
                                                                                       gcgcCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGc                                                                                                                            >  2:367023/1‑61 (MQ=255)
                                                                                       gcgcCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGc                                                                                                                            <  1:367023/61‑1 (MQ=255)
                                                                                       gcgcCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCa                                                                                                                           >  2:1766291/1‑62 (MQ=255)
                                                                                       gcgcCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCa                                                                                                                           <  1:1766291/62‑1 (MQ=255)
                                                                                       gcgcCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCAcc                                                                                     <  1:1623439/100‑1 (MQ=255)
                                                                                       gcgcCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCAcc                                                                                     >  2:1623439/1‑100 (MQ=255)
                                                                                        cgcCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAg                                             <  2:33771/139‑1 (MQ=255)
                                                                                         gcCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTa                                                               <  2:1504453/120‑1 (MQ=255)
                                                                                         gcCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTa                                                               >  1:1504453/1‑120 (MQ=255)
                                                                                           cAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCca                                          <  2:1501935/139‑1 (MQ=255)
                                                                                             ggCGGTACGCAGCGTGTCCGGCGGGCGCTTTTCGCTGTATACGTTATCg                                                                                                                                  <  2:688404/49‑1 (MQ=255)
                                                                                             ggCGGTACGCAGCGTGTCCGGCGGGCGCTTTTCGCTGTATACGTTATCg                                                                                                                                  >  1:688404/1‑49 (MQ=255)
                                                                                                  tACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTa                                                                                                                                <  1:549762/46‑1 (MQ=255)
                                                                                                  tACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTa                                                                                                                                >  2:549762/1‑46 (MQ=255)
                                                                                                            gCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCAcc                                                                                     >  1:894465/1‑79 (MQ=255)
                                                                                                            gCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCAcc                                                                                     <  2:894465/79‑1 (MQ=255)
                                                                                                            gCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCAc                                                                   <  1:393869/97‑1 (MQ=255)
                                                                                                            gCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCAc                                                                   >  2:393869/1‑97 (MQ=255)
                                                                                                            gCCCGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCcggc                         >  2:1808564/1‑139 (MQ=255)
                                                                                                              ccGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAAt                                                                      >  1:17626/1‑92 (MQ=255)
                                                                                                              ccGGCGGGCGCTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAAt                                                                      <  2:17626/92‑1 (MQ=255)
                                                                                                                 gcgggcgCTTTTCGCTGTATACGTTATCGTAAGGc                                                                                                                            >  2:1723078/1‑35 (MQ=255)
                                                                                                                 gcgggcgCTTTTCGCTGTATACGTTATCGTAAGGc                                                                                                                            <  1:1723078/35‑1 (MQ=255)
                                                                                                                   gggcgCTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATg                  >  1:36094/1‑139 (MQ=255)
                                                                                                                    ggcgCTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATgg                 <  2:36094/139‑1 (MQ=255)
                                                                                                                      cgcTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCTTATGTTTCa                                                                                                        >  2:308078/1‑50 (MQ=255)
                                                                                                                      cgcTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCTTATGTTTCa                                                                                                        <  1:308078/50‑1 (MQ=255)
                                                                                                                        cTTTTCGCTGTATACGTTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGc             >  2:606272/1‑139 (MQ=255)
                                                                                                                                   ataCGTTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCctgct  >  1:1714371/1‑139 (MQ=255)
                                                                                                                                        |                                                                                                                                     
CGCAAACCAGCCGCCAAAAATACACAGTACAGCCAGCCCCGCGCAGCCGTACAGGCCGACCATCGCAGAGGCGTCACTATAAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCCTGCT  >  NC_000913/1354369‑1354638

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: