Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,352,926 C→T G235S (GGT→AGT)  sapD ← antimicrobial peptide ABC transporter ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,352,9260CT98.0% 177.6 / ‑4.2 49G235S (GGT→AGT) sapDantimicrobial peptide ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (18/30);  total (18/31)
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.92e-01

GTGGGATAGCGCCGGGCAGCGTATTGAGGCGACTTTTATGCGGCATCGCGCTGCCGAAGTCTGGTATCGCGCGGATCAGCGCCTGGGTATAAGGATGATGTGGCATCGTCACCAACTCCTTACTCGGCGCGGTTTCCACTGTTTGACCGCAGTAAAGCACGTTAATTTTATCCGCCCATTGGCTAAGCATTTGTAAGTCATGGCTGATAAGCAAAATAGTGGTATTGCTGTTTTGGTTGAGACGCGTCAGCAGGCGAAAGATTTGCGCCTGGGTTGTTGGCTCCATTGAGT  >  NC_000913/1352779‑1353069
                                                                                                                                                   |                                                                                                                                               
gTGGGATAGCGCCGGGCAGCGTATTGAGGCGACTTTTATGCGGCATCGCGCTGCCGAAGTCTGGTATCGCGCGGATCAGCGCCTGGGTATAAGGATGATGTGGCATCGTCACCAACTCCTTACTCGGCGCGGTTTCCACTGTTTGACTGCa                                                                                                                                              <  2:799691/151‑1 (MQ=255)
  gggATAGCGCCGGGCAGCGTATTGAGGCGACTTTTATGCGGCATCGCGCTGCCGAAGTCTGGTATCGCGCGGATCAGCGCCTGGGTATAAGGATGATGTGGCATCGTCACCAACTCCTTACTCGGCGCGGTTTCCACTGTTTGACTGCAGt                                                                                                                                            <  2:501425/151‑1 (MQ=255)
   ggATAGCGCCGGGCAGCGTATTGAGGCGACTTTTATGCGGCATCGCGCTGCCGAAGTCTGGTATCGCGCGGATCAGCGCCTGGGTATAAGGATGATGTGGCATCGTCACCAACTCCTTACTCGGCGCGGTTTCCACTGTTTGACTGCAGTa                                                                                                                                           <  1:120202/151‑1 (MQ=255)
   ggATAGCGCCGGGCAGCGTATTGAGGCGACTTTTATGCGGCATCGCGCTGCCGAAGTCTGGTATCGCGCGGATCAGCGCCTGGGTATAAGGATGATGTGGCATCGTCACCAACTCCTTACTCGGCGCGGTTTCCACTGTTTGACTGCAGTa                                                                                                                                           <  1:155986/151‑1 (MQ=255)
     aTAGCGCCGGGCAGCGTATTGAGGCGACTTTTATGCGGCATCGCGCTGCCGAAGTCTGGTATCGCGCGGATCAGCGCCTGGGTATAAGGATGATGTGGCATCGTCACCAACTCCTTACTCGGCGCGGTTTCCACTGTTTGACTGCAGTaaa                                                                                                                                         <  1:238287/151‑1 (MQ=255)
     aTAGCGCCGGGCAGCGTATTGAGGCGACTTTTATGCGGCATCGCGCTGCCGAAGTCTGGTATCGCGCGGATCAGCGCCTGGGCATAAGGATGATGTGGCATCGTCACCAACTCCTTACTCGGCGCGGTTTCCACTGTTTGACTGCAGTa                                                                                                                                           >  1:776768/1‑149 (MQ=255)
            cGGGCAGCGTATTGAGGCGACTTTTATGCGGCATCGCGCTGCCGAAGTCTGGTATCGCGCGGATCAGCGCCTGGGTATAAGGATGATGTGGCATCGTCACCAACTCCTTACTCGGCGCGGTTTCCACTGTTTGACTGCAGTAAAGCACGtt                                                                                                                                  <  1:630652/151‑1 (MQ=255)
            cGGGCAGCGTATTGAGGCGACTTTTATGCGGCATCGCGCTGCCGAAGTCTGGTATCGCGCGGATCAGCGCCTGGGTATAAGGATGATGTGGCATCGTCACCAACTCCTTACTCGGCGCGGTTTCCACTGTTTGACTGCAGTAAAGCACGtt                                                                                                                                  <  1:75146/151‑1 (MQ=255)
               gCAGCGTATTGAGGCGACTTTTATGCGGCATCGCGCTGCCGAAGTCTGGTATCGCGCGGATCAGCGCCTGGGTATAAGGATGATGTGGCATCGTCACCAACTCCTTACTCGGCGCGGTTTCCACTGTTTGACTGCAGTAAAGCACGTTAAt                                                                                                                               <  1:361022/151‑1 (MQ=255)
                    gTATTGAGGCGACTTTTATGCGGCATCGCGCTGCCGAAGTCTGGTATCGCGCGGATCAGCGCCTGGGTATAAGGATGATGTGGCATCGTCACCAACTCCTTACTCGGCGCGGTTTCCACTGTTTGACTGCAGTAAAGCACGTTAATTTTAt                                                                                                                          >  2:325432/1‑151 (MQ=255)
                      aTTGAGGCGACTTTTATGCGGCATCGCGCTGCCGAAGTCTGGTATCGCGCGGATCAGCGCCTGGGTATAAGGATGATGTGGCATCGTCACCAACTCCTTACTCGGCGCGGTTTCCACTGTTTGACTGCAGTAAAGCATGTTAATTTTATcc                                                                                                                        <  2:495682/151‑1 (MQ=255)
                       ttGAGGCGACTTTTATGCGGCATCGCGCTGCCGAAGTCTGGTATCGCGCGGATCAGCGCCTGGGAATAAGGATGATGTGGCATCGTCACCAACTCCTTACTCGGCGCGGTTTCCACTGTTTGACTGCAGTAAAGCACGTTAATTTTATCCg                                                                                                                       >  1:388082/1‑151 (MQ=255)
                         gAGGCGACTTTTATGCGGCATCGCGCTGCCGAAGTCTGGTATCGCGCGGATCAGCGCCTGGGTATAAGGATGATGTGGCATCGTCACCAACTCCTTACTCGGCGCGGTTTCCACTGTTTGACTGCAGTAAAGCACGTTAATTTTATCCGcc                                                                                                                     >  2:950428/1‑151 (MQ=255)
                          aGGCGACTTTTATGCGGCATCGCGCTGCCGAAGTCTGGTATCGCGCGGATCAGCGCCTGGGTATAAGGATGATGTGGCATCGTCACCAACTCCTTACTCGGCGCGGTTTCCACTGTTTGACTGCAGTAAAGCACGTTAATTTTATCCGccc                                                                                                                    <  2:243929/151‑1 (MQ=255)
                               aCTTTTATGCGGCATCGCGCTGCCGAAGTCTGGTATCGCGCGGATCAGCGCCTGGGTATAAGGATGATGTGGCATCGTCACCAACTCCTTACTCGGCGCGGTTTCCACTGTTTGACTGCAGTAAAGCACGTTAATTTTATCCGCCCATTgg                                                                                                               >  1:904289/1‑151 (MQ=255)
                                  tttATGCGGCATCGCGCTGGCGAAGTCTGGTGTCGCGCGGATCATCGCCTGGGTATAAGGATGATGTGGCTTCGTCACCAACTCCTTACTCGGCGCGGTTTCCACTGTTTGACCGCAGTAAATCACGTTAATTTTATCCGCCCATTGGc                                                                                                              <  2:213944/149‑1 (MQ=255)
                                    tATGCGGCATCGCGCTGCCGAAGTCTGGTATCGCGCGGATCAGCGCCTGGGTATAAGGATGATGTGGCATCGTCACCAACTCCTTACTCGGCGCGGTTTCCACTGTTTGACTGCAGTAAAGCACGTTAATTTTATCCGCCCATTGGCTa                                                                                                            >  1:598689/1‑149 (MQ=255)
                                     aTGCGGCATCGCGCTGCCGAAGTCTGGTATCGCGCGGATCAGCGCCTGGGTATAAGGATGATGTGGCATCGTCACCAACTCCTTACTCGGCGCGGTTTCCACTGTTTGACTGCAGTAAAGCACGTTAATTTTATCCGCCCATTGGCTAAg                                                                                                          <  1:280882/150‑1 (MQ=255)
                                                 gcTGCCGAAGTCTGGTATCGCGCGGATCAGCGCCTGGGTATAAGGATGATGTGGCATCGTCACCAACTCCTTACTCGGCGCGGTTTCCACTGTTTGACTGCAGTAAAGCACGTTAATTTTATCCGCCCATTGGCTAAGCATTTGTAAGTc                                                                                              >  1:589611/1‑150 (MQ=255)
                                                    gCCGAAGTCTGGTATCGCGCGGATCAGCGCCTGGGTATAAGGATGATGTGGCATCGTCACCAACTCCTTACTCGGCGCGGTTTCCACTGTTTGACTGCAGTAAAGCACGTTAATTTTATCCGCCCATTGGCTAAGCATTTGTAAGTCATgg                                                                                          <  2:738042/151‑1 (MQ=255)
                                                    gCCGAAGTCTGGTATCGCGCGGATCAGCGCCTGGGTATAAGGATGATGTGGCATCGTCACCAACTCCTTACTCGGCGCGGTTTCCACTGTTTGACTGCAGTAAAGCACGTTAATTTTATCCGCCCATTGGCTAAGCATTTGTAAGTCATg                                                                                           <  1:139169/150‑1 (MQ=255)
                                                      cGAAGTCTGGTATCGCGCGGATCAGCGCCTGGGTATAAGGATGATGTGGCATCGTCACCAACTCCTTACTCGGCGCGGTTTCCACTGTTTGACTGCAGTAAAGCACGTTAATTTTATCCGCCCATTGGCTAAGCATTTGTAAGTCATGGCt                                                                                        <  2:803346/151‑1 (MQ=255)
                                                                 aTCGCGCGGATCAGCGCCTGGGTATAAGGATGATGTGGCATCGTCACCAACTCCTTACTCGGCGCGGTTTCCACTGTTTGACTGCAGTAAAGCACGTTAATTTTATCCGCCCATTGGCTAAGCATTTGTAAGTCATGGCTGATAAGCaa                                                                               >  2:121016/1‑149 (MQ=255)
                                                                     cgcgGATCAGCGCCTGGGTATAAGGATGATGTGGCATCGTCACCAACTCCTTACTCGGCGCGGTTTCCACTGTTTGACTGCAGTAAAGCACGTTAATTTTATCCGCCCATTGGCTAAGCATTTGTAAGTCATGGCTGATAAGCAAAATAg                                                                          <  2:946416/150‑1 (MQ=255)
                                                                       cgGATCAGCGCGTGGGTATAAGGATGATGTGTCATCGTCACCACCTCCTTACTCGGCGCGGTTTCCACTGTTTGCCTGCAGTAAAGCACGTTAATTTTATCCGCCCATAGGCTAAGCATTTGTAAGTCATGGCTAATAAGCAAAATAGTg                                                                        <  2:829943/150‑1 (MQ=255)
                                                                                    gggTATAAGGATGATGTGGCATCGTCACCAACTCCTTACTCGGCGCGGTTTCCACTGTTTGACTGCAGTAAAGCACGTTAATTTTATCCGCCCATTGGCTAAGCATTTGTAAGTCATGGCTGATAAGCAAAATAGTGGTATTGCTGttt                                                            <  2:108716/149‑1 (MQ=255)
                                                                                       tataAGGATGATGTGGCATCGTCACCAACTCCTTACTCGGCGCGGTTTCCACTGTTTGACTGCAGTAAAGCACGTTAATTTTATCCGCCCATTGGCTAAGCATTTGTAAGTCATGGCTGATAAGCAAAATAGTGGTATTGCTGTTTTGGtt                                                       >  1:474705/1‑151 (MQ=255)
                                                                                        ataAGGATGATGTGGCATCGTCACCAACTCCTTACTCGGCGCGGTTTCCACTGTTTGACTGCAGTAAAGCACGTTAATTTTATCCGCCCATTGGCTAAGCATTTGTAAGTCATGGCTGATAAGCAAAATAGTGGTATTGCTGTTTTGGtt                                                       <  2:725019/150‑1 (MQ=255)
                                                                                        ataAGGATGATGTGGCATCGTCACCAACTCCTTACTCGGCGCGGTTTCCACTGTTTGACTGCAGTAAAGCACGTTAATTTTATCCGCCCATTGGCTAAGCATTTGTAAGTCATGGCTGATAAGCAAAATAGTGGTATTGCTGTTTTGGtt                                                       <  2:598089/150‑1 (MQ=255)
                                                                                              atgatgTGGCATCGTCACCAACTCCTTACTCGGCGCGGTTTCCACTGTTTGACTGCAGTAAAGCACGTTAATTTTATCCGCCCATTGGCTAAGCATTTGTAAGTCATGGCTGATAAGCAAAATAGTGGTATTGCTGTTTTGGTTGAGAcgc                                                >  1:499554/1‑151 (MQ=255)
                                                                                                gatgTGGCATCGTCACCAACTCCTTACTCGGCGCGGTTTCCACTGTTTGACTGCAGTAAAGCACGTTAATTTTATCCGCCCATTGGCTAAGCATTTGTAAGTCATGGCTGATAAGCAAAATAGTGGTATTGCTGTTTTGGTTGAGAcgcg                                               <  1:92908/150‑1 (MQ=255)
                                                                                                gatgTGGCATCGTCACCAACTCCTTACTCGGCGCGGTTTCCACTGTTTGACTGCAGTAAAGCACGTTAATTTTATCCGCCCATTGGCTAAGCATTTGTAAGTCATGGCTGATAAGCAAAATAGTGGTATTGCTGTTTTGGTTGAGACGCGt                                              >  1:494635/1‑151 (MQ=255)
                                                                                                 atgTGGCATCGTCACCAACTCCTTACTCGGCGCGGTTTCCACTGTTTGACTGCAGTAAAGCACGTTAATTTTATCCGCCCATTGGCTAAGCATTTGTAAGTCATGGCTGATAAGCAAAATAGTGGTATTGCTGTTTTGGTTGAGACGCGTc                                             <  1:756209/151‑1 (MQ=255)
                                                                                                             cACCAACTCCTTACTCGGCGCGGTTTCCACTGTTTGACTGCAGTAAAGCACGTTAATTTTATCCGCCCATTGGCTAAGCATTTGTAAGTCATGGCTGATAAGCAAAATAGTGGTATTGCTGTTTTGGTTGAGACGCGTCAGCAGGCGaaa                                  <  1:410734/150‑1 (MQ=255)
                                                                                                             cACCAACTCCTTACTCGGCGCGGTTTCCACTGTTTGACTGCAGTAAAGCACGTTAATTTTATCCGCCCATTGGCTAAGCATTTGTAAGTCATGGCTGATAAGCAAAATAGTGGTATTGCTGTTTTGGTTGAGACGCGTCAGCAGGCGaa                                   >  1:483464/1‑149 (MQ=255)
                                                                                                              aCCAACTCCTTACTCGGCGCGGTTTCCACTGTTTGACTGCAGTAAAGCACGTTAATTTTATCCGCCCATTGGCTAAGCATTTGTAAGTCATGGCTGATAAGCAAAATAGTGGTATTGCTGTTTTGGTTGAGACGCGTCAGCAGGCGaa                                   >  2:629590/1‑148 (MQ=255)
                                                                                                                 aaCTCCTTACTCGGCGCGGTTTCCACTGTTTGACTGCAGTAAAGCACGTTAATTTTATCCGCCCATTGGCTAAGCATTTGTAAGTCATGGCTGATAAGCAAAATAGTGGTATTGCTGTTTTGGTTGAGACGCGTCAGCAGGCGAAAGAttt                             <  2:350753/151‑1 (MQ=255)
                                                                                                                     ccTTACTCGGCGCGGTTTCCACTGTTTGACTGCAGTAAAGCACGTTAATTTTATCCGCCCATTGGCTAAGCATTTGTAAGTCATGGCTGATAAGCAAAATAGTGGTATTGCTGTTTTGGTTGAGACGCGTCAGCAGGCGAAAGATTTgc                           <  2:54913/149‑1 (MQ=255)
                                                                                                                         aCTCGGCGCGGTTTCCACTGTTTGACTGCAGTAAAGCACGTTAATTTTATCCGCCCATTGGCTAAGCATTTGTAAGTCATGGCTGATAAGCAAAATAGTGGTATTGCTGTTTTGGTTGAGACGCGTCAGCAGGCGAAAGATTTGCGCCTgg                     >  2:672383/1‑151 (MQ=255)
                                                                                                                            cGGCGCGGTTTCCACTGTTTGACTGCAGTAAAGCACGTTAATTTTATCCGCCCATTGGCTAAGCATTTGTAAGTCATGGCTGATAAGCAAAATAGTGGTATTGCTGTTTTGGTTGAGACGCGTCAGCAGGCGAAAGATTTGCGCCTGGgtt                  <  1:149557/151‑1 (MQ=255)
                                                                                                                            cGGCGCGGTTTCCACTGTTTGACTGCAGTAAAGCACGTTAATTTTATCCGCCCATTGGCTAAGCATTTGTAAGTCATGGCTGATAAGCAAAATAGTGGTATTGCTGTTTTGGTTGAGACGCGTCAGCAGGCGAAAGATTTGCGCCTGGgtt                  <  1:8414/151‑1 (MQ=255)
                                                                                                                            cGGCGCGGTTTCCACTGTTTGACTGCAGTAAAGCACGTTAATTTTATCCGCCCATTGGCTAAGCATTTGTAAGTCATGGCTGATAAGCAAAATAGTGGTATTGCTGTTTTGGTTGAGACGCGTCAGCAGGCGAAAGATTTGCGCCTGGgtt                  >  2:311083/1‑151 (MQ=255)
                                                                                                                                gcgGTTTCCACTGTTTGACTGCAGTAAAGCACGTTAATTTTATCCGCCCATTGGCTAAGCATTTGTAAGTCATGGCTGATAAGCAAAATAGTGGTATTGCTGTTTTGGTTGAGACGCGTCAGCAGGCGAAAGATTTGCGCCTGGGttgttg              >  2:819210/1‑151 (MQ=255)
                                                                                                                                 cgGTTTCCACTGTTTGACTGCAGTAAAGCACGTTAATTTTATCCGCCCATTGGCTAAGCATTTGTAAGTCATGGCTGATAAGCAAAATAGTGGTATTGCTGTTTTGGTTGAGACGCGTCAGCAGGCGAAAGATTTGCGCCTGGGttgttg              <  1:378214/150‑1 (MQ=255)
                                                                                                                                 cgGTTTCCACTGTTTGACTGCAGTAAAGCACGTTAATTTTATCCGCCCATTGGCTAAGCATTTGTAAGTCATGGCTGATAAGCAAAATAGTGGTATTGCTGTTTTGGTTGAGACGCGTCAGCAGGCGAAAGATTTGCGCCTGGGTTGTTgg             >  2:846360/1‑151 (MQ=255)
                                                                                                                                     ttCCACTGTTTGACTGCAGTAAAGCACGTTAATTTTATCCGCCCATTGGCTAAGCATTTGTAAGTCATGGCTGATAAGCAAAATAGTGGTATTGCTGTTTTGGTTGAGACGCGTCAGCAGGCGAAAGATTTGCGCCTGGGTTGTTGGCTcc         >  1:421112/1‑151 (MQ=255)
                                                                                                                                       ccACTGTTTGACTGCAGTAAAGCACGTTAATTTTATCCGCCCATTGGCTAAGCATTTGTAAGTCATGGCTGATAAGCAAAATAGTGGTATTGCTGTTTTGGTTGAGACGCGTCAGCAGGCGAAAGATTTGCGCCTGGGTTGTTGGCTCCa        <  2:196503/150‑1 (MQ=255)
                                                                                                                                       ccACTGTTTGACTGCAGTAAAGCACGTTAATTTTATCCGCCCATTGGCTAAGCATTTGTAAGTCATGGCTGATAAGCAAAATAGTGGTATTGCTGTTTTGGTTGAGACGCGTCAGCAGGCGAAAGATTTGCGCCTGGGTTGTTGGCTCCa        <  1:846735/150‑1 (MQ=255)
                                                                                                                                             tttGACTGCAGTAAAGCACGTTAATTTTATCCGCCCATTGGCTAAGCATTTGTAAGTCATGGCTGATAAGCAAAATAGTGGTATTGCTGTTTTGGTTGAGACGCGTCAGCAGGCGAAAGATTTGCGCCTGGGTTGTTGGCTCCATTGAGt  <  2:869168/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                               
GTGGGATAGCGCCGGGCAGCGTATTGAGGCGACTTTTATGCGGCATCGCGCTGCCGAAGTCTGGTATCGCGCGGATCAGCGCCTGGGTATAAGGATGATGTGGCATCGTCACCAACTCCTTACTCGGCGCGGTTTCCACTGTTTGACCGCAGTAAAGCACGTTAATTTTATCCGCCCATTGGCTAAGCATTTGTAAGTCATGGCTGATAAGCAAAATAGTGGTATTGCTGTTTTGGTTGAGACGCGTCAGCAGGCGAAAGATTTGCGCCTGGGTTGTTGGCTCCATTGAGT  >  NC_000913/1352779‑1353069

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