Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,354,596 C→T 6.6% G292D (GGC→GAC)  sapB ← 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,5960CT 6.6% 185.2 / 2.6 60G292D (GGC→GAC) sapBantimicrobial peptide transport ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base C (26/30);  new base T (2/2);  total (28/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

AAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCCTGCTGGCGGATGGCGTTAATTAACCAGCGCCCTAAACCCGGCCAGCTAAAAACCATTTCGGTGATCATCGCCAGCGTCAGCATGGTGGAAAACTGTAAGCCCAGGCG  >  NC_000913/1354449‑1354741
                                                                                                                                                   |                                                                                                                                                 
aaaaTTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGc                                                                                                                                                   >  5:43572/1‑148 (MQ=255)
      tGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGCGCcacaac                                                                                                                                             >  6:22373/1‑146 (MQ=255)
         gcCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCa                                                                                                                                          <  6:43572/148‑1 (MQ=255)
            aGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCCTGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCa                                                                                                                                       >  3:138137/1‑148 (MQ=255)
             ggCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATc                                                                                                                                        >  4:272778/1‑146 (MQ=255)
             ggCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATc                                                                                                                                        <  3:272778/146‑1 (MQ=255)
                       aGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGaa                                                                                                                            <  6:327210/148‑1 (MQ=255)
                         cGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAAt                                                                                                                          <  7:435564/148‑1 (MQ=255)
                                         ttttCGCTGTATACGCTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGTCACACACCATCACTCCTGCGGaa                                                                                                                            >  1:199835/1‑130 (MQ=255)
                                         ttttCGCTGTATACGCTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGTCACACACCATCACTCCTGCGGaa                                                                                                                            <  2:199835/130‑1 (MQ=255)
                                         ttttCGCTGTATACGCTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGTCACACACCATCACTCCGGCGGaa                                                                                                                            <  2:205927/130‑1 (MQ=255)
                                         ttttCGCTGTATACGCTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGTCACACACCATCACTCCAGCGGaa                                                                                                                            >  1:205927/1‑130 (MQ=255)
                                          tttCGCTGTATACGCTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCctgct                                                                                                         >  3:2180/1‑148 (MQ=255)
                                           ttCGCTGTATACGCTATCGTAAGGCATACCATTCCTTATGTTTCATAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCctgctg                                                                                                        >  6:233648/1‑148 (MQ=255)
                                               cTGTATACGCTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCCTGCTGGCg                                                                                                     >  7:173089/1‑147 (MQ=255)
                                                  tataCGCTTTAGTTAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCAGCCAAAATATCATAAATAAAGTTAACAATAATGACCAGTTAGCCACTCACAATCACTCCGGCGGAAATGGCTGCGTAATCCTGCTGGCGGa                                                                                                   <  6:426936/146‑1 (MQ=255)
                                                  tataCGCTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATTACGTTAACAATAATGACCAGTGAGCCACACACCATCACt                                                                                                                                     >  8:445823/1‑112 (MQ=255)
                                                  tataCGCTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATTACGTTAACAATAATGACCAGTGAGCCACACACCATCACt                                                                                                                                     <  7:445823/112‑1 (MQ=255)
                                                  tataCGCTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCCTGCTGGCGGATg                                                                                                 <  8:192488/148‑1 (MQ=255)
                                                  tataCGCTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGGGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCCTGCTGGCGGa                                                                                                   <  3:140152/146‑1 (MQ=255)
                                                     aCGCTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCTTGCTGGCGGATGGCg                                                                                              >  8:424304/1‑148 (MQ=255)
                                                      cGCTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCCTGCTGGCGGATGGCg                                                                                              <  6:91159/147‑1 (MQ=255)
                                                          aTCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCAc                                                                                                                                      <  2:169064/103‑1 (MQ=255)
                                                          aTCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCAc                                                                                                                                      >  1:169064/1‑103 (MQ=255)
                                                             gTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGCTCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCATGCTGGCGGATGGCGttaattaa                                                                                      >  3:243004/1‑148 (MQ=255)
                                                             gTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCCTGCTGGCGGATGGCGttaatta                                                                                       <  4:138137/147‑1 (MQ=255)
                                                                  gCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCCTGCTGGCGGATGGCGTTAATTAACCAgc                                                                                 <  5:22373/148‑1 (MQ=255)
                                                                      aCCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGaa                                                                                                                            <  3:208342/101‑1 (MQ=255)
                                                                      aCCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGaa                                                                                                                            >  4:208342/1‑101 (MQ=255)
                                                                      aCCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCCTGCTGGCGGATGGCGttaattaa                                                                                      <  4:169108/139‑1 (MQ=255)
                                                                      aCCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCCTGCTGGCGGATGGCGttaattaa                                                                                      >  3:169108/1‑139 (MQ=255)
                                                                      aCCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCCTGCTGGCGGATGGCGtta                                                                                           <  7:299591/134‑1 (MQ=255)
                                                                      aCCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCCTGCTGGCGGATGGCGtta                                                                                           >  8:299591/1‑134 (MQ=255)
                                                                       ataTTCCTTATGTTTCAGAGGTTTATCCATGGCACACAAAATTTCAGAAATCACGTTAACAATAATGACCAGTGAGCCAGACACCATAACTCCGGCGGATATGGCTGCGTTTTCCTGCTGGCGGATGGCGTTATTTAAGCAGCGccg                                                                             <  7:79636/145‑2 (MQ=255)
                                                                        cATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATAAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTTCTTAATCCTGCTGGCGGATGGCTTTAATTAACCAGAGccat                                                                            >  4:253354/1‑145 (MQ=255)
                                                                          ttCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCCTGCTGGCGGATGGCGTTAATTAACCAGCGCCCTaaa                                                                         >  1:269285/1‑148 (MQ=255)
                                                                             cTTATGTTTCAGAGGGTTAGCAATGGCACCCAAAATATAAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCCTGCTGGAGGATGGCGTTAATTAACCAGCGCCCTAAAcc                                                                       <  3:253354/147‑1 (MQ=255)
                                                                              ttATGTTTCAGAGGGTAAGCCATGCCACCCAAAATATCAGAAATCACGTAAACAATAATGACCAGAGAGCCAAAGACCATCCCTCCGTCGTCAAAGGCTACCTAATCCTCCTGGCCCCACGCGGAAATTAACAAGCGCCCTaacaccg                                                                     >  4:285394/1‑143 (MQ=255)
                                                                               tATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCCTGCTGGCGGATGGCGTTAATTAACCAGCGCCCTAAACCCgg                                                                    >  2:269312/1‑148 (MQ=255)
                                                                                   tttCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCAGCACTCCGGCGGAAATGGCTGCGTAATCCTGCTGGCGGATGGCGTTAATTAACCAGCGCCCTAAACCCGGCCAg                                                                <  3:72023/148‑1 (MQ=255)
                                                                                                gCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCCTGCTGGCGGATGGCGTTAATTAACCAGCGCCCTAAACCCGGCCAGCTAAAAACCAtt                                                    <  5:233648/147‑1 (MQ=255)
                                                                                                    tGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCCTGCTGGCGGATGGCGTTAATTAACCAGCGCCCTAAACCCGGCCAGCTAAAAACCATTTCGGt                                               <  8:173089/148‑1 (MQ=255)
                                                                                                    tGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCCTGCTGGCGGATGGCGTTAATTAACCAGCGCCCTAAACCCGGCCAGCTAAAAACCATTTCGGt                                               >  4:9651/1‑148 (MQ=255)
                                                                                                         cccAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACCCCATCACTCCGGCGGAAATGGCTGCGTAATCCTGCTGGCGGATGGCGTTAATTAACCAGCGCCCTACACCCGGCCAGCTAAAACCCATTTCGGTGatca                                          >  1:396015/1‑148 (MQ=255)
                                                                                                              aaTATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCCTGCTGGCGGATGGCGTTAATTAACCAGCGCCCTAAACCCGGCCAGCTAAAAACCAtt                                                    >  4:186329/1‑133 (MQ=255)
                                                                                                              aaTATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCCTGCTGGCGGATGGCGTTAATTAACCAGCGCCCTAAACCCGGCCAGCTAAAAACCAtt                                                    <  3:186329/133‑1 (MQ=255)
                                                                                                                 atCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGc                                                                                                                       >  7:391833/1‑63 (MQ=255)
                                                                                                                 atCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGc                                                                                                                       <  8:391833/63‑1 (MQ=255)
                                                                                                                  tCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCCTGCTGGCGGATGGCGTTAATTAACCAGCGCCCTAAACCCGGCCAGCTAAAAACCATTTCGGTGATCATCGCCAGCg                                 >  8:447056/1‑148 (MQ=255)
                                                                                                                         tCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCTTGCTGGCGGATGGCGTTAATTAACCAGCGCCCTAAACCCGGCCAGCTAAAAACCATTTCGGTGATCATCGCCAGCGTCAGCAt                          <  7:424304/148‑1 (MQ=255)
                                                                                                                          cACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCCTGCTGGCGGATGGCGTTAATTAACCAGCGCCCTAAACCCGGCCAGCTAAAAACCATTTCGGTGATCATCGCCAGCGTCAGc                            >  5:413099/1‑145 (MQ=255)
                                                                                                                          cACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCCTGCTGGCGGATGGCGTTAATTAACCAGCGCCCTAAACCCGGCCAGCTAAAAACCATTTCGGTGATCATCGCCAGCGTCAGc                            <  6:413099/145‑1 (MQ=255)
                                                                                                                           aCGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCCTGCTGGCGGATGGCGTTAATTAACCAGCGCCCTAAACCCGGCCAGCTAAAAACCATTTCGGTGATCATCGCCAGCGTCAGCAtg                         <  2:269285/147‑1 (MQ=255)
                                                                                                                            cGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCCTGCTGGCGGATGGCGTTAATTAACCAGCGCCCTAAACCCGGCCAGCTAAAAACCATTTCGGTGATCATCGCCAGCGTCAGCAtggt                       >  7:52497/1‑148 (MQ=255)
                                                                                                                                aaCAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCCTGCTGGCGGATGGCGATAATTAACCAGCGCCCTAAACCCGGCCAGCTAAAAACCATTTCGGTGATCATCGCCAGCGTCAGCATGGTGGa                    >  7:228/1‑147 (MQ=255)
                                                                                                                                 aCAATAATGACCTTTTAGCCACACTGCATCACTCCGGCGGAAATGGCTGCGTAATCCTGCTGGCGGATTGCTTTAATTAACCAGCTCACTAAACCCTGCCAGCTAAAAACGATTTCGGTGATCATCGCCAGCGTCAGCATGGTGGaaa                  <  2:403348/148‑1 (MQ=255)
                                                                                                                                    ataatGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCCTGCTGGCGGATGGCGttaatta                                                                                       <  5:167219/76‑1 (MQ=255)
                                                                                                                                    ataatGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCCTGCTGGCGGATGGCGttaatta                                                                                       >  6:167219/1‑76 (MQ=255)
                                                                                                                                          aCCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCCTGCTGGCGGATGGCGTTAATTAACCAGCGCCCTAAACCCGGCCAGCTAAAAACCATTTCGGTGATCATCGCCAGCGTCAGCATGGTGGAAAACTGTAAGc         <  1:269312/148‑1 (MQ=255)
                                                                                                                                            cAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCCTGCTGGCGGATGGCGTTAATTAACCAGCGCCCTAAACCCGGCCAGCTAAAAACCATTTCGGTGATCATCGCCAGCGTCAGCATGGTGGAAAACTGTAAGccc       <  8:52497/148‑1 (MQ=255)
                                                                                                                                                  gCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCCTGCTGGCGGATGGCGTTAATTAACCAGCGCCCTAAACCCGGCCAGCTAAAAACCATTTCGGTGATCATCGCCAGCGTCAGCATGGTGGAAAACTGTAAGCCCAGGCg  <  3:9651/147‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
AAAATTTGCGCCAGGCGGTACGCAGCGTGCCCGGCGGGCGCTTTTCGCTGTATACGCTATCGTAAGGCATACCATTCCTTATGTTTCAGAGGGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCCTGCTGGCGGATGGCGTTAATTAACCAGCGCCCTAAACCCGGCCAGCTAAAAACCATTTCGGTGATCATCGCCAGCGTCAGCATGGTGGAAAACTGTAAGCCCAGGCG  >  NC_000913/1354449‑1354741

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