Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,353,062	C→A	M189I (ATG→ATT)	sapD ←	antimicrobial peptide ABC transporter ATPase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,353,062	0	C	A	100.0%	210.3 / NA	52	M189I (ATG→ATT)	sapD	antimicrobial peptide ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (25/27);  total (25/27)

CTGTTTGACCGCAGTAAAGCACGTTAATTTTATCCGCCCATTGGCTAAGCATTTGTAAGTCATGGCTGATAAGCAAAATAGTGGTATTGCTGTTTTGGTTGAGACGCGTCAGCAGGCGAAAGATTTGCGCCTGGGTTGTTGGCTCCATTGAGTTGGTCGGTTCGTCAGCAATCAGCAGACGCGGTTGATTCGCCAGTGCAATGGCTATCATCACTTTCTGACATTCACCTTCGGTCAACTCATAGGGAAAACTGCGCATCGCATCTTTGTGATCTTTAATCCCCA  >  NC_000913/1352917‑1353201                                                                                                                                                  |                                                                                                                                            cTGTTTGACCGCAGTAAAGCACGTTAATTTTATCCGCCCATTGGCTAAGCATTTGTAAGTCATGGCTGATAAGCAAAATAGTGGTATTGCTGTTTTGGTTGAGACGCGTCAGCAGGCGAAAGATTTGCGCCTGGGTTGTTGGCTCAATTg                                                                                                                                         >  1:877845/1‑150 (MQ=255)     ttGACCGCAGTAAAGCACGTTAATTTTATCCGCCCATTGGCTAAGCATTTGTAAGTCATGGCTGATAAGCAAAATAGTGGTATTGCTGTTTTGGTTGAGACGCGTCAGCAGGCGAAAGATTTGCGCCTGGGTTGTTGGCTCAATTGAGtt                                                                                                                                     <  1:525676/150‑1 (MQ=255)       gACCGCAGTAAAGCACGTTAATTTTATCCGCCCATTGGCTAAGCATTTGTAAGTCATGGCTGATAAGCAAAATAGTGGTATTGCTGTTTTGGTTGAGACGCGTCAGCAGGCGAAAGATTTGCGCCTGGGTTGTTGGCTCAATTGAGTTGGt                                                                                                                                  >  2:472672/1‑151 (MQ=255)         ccGCAGTAAAGCACGTTAATTTTATCCGCCCATTGGCTAAGCATTTGTAAGTCATGGCTGATAAGCAAAATAGTGGTATTGCTGTTTTGGTTGAGACGCGTCAGCAGGCGAAAGATTTGCGCCTGGGTTGTTGGCTCAATTGAGTTGGTc                                                                                                                                 <  1:128564/150‑1 (MQ=255)         ccGCAGTAAAGCACGTTAATTTTATCCGCCCATTGGCTAAGCATTTGTAAGTCATGGCTGATAAGCAAAATAGTGGTATTGCTGTTTTGGTTGAGACGCGTCAGCAGGCGAAAGATTTGCGCCTGGGTTGTTGGCTCAATTGAGTTGGTCg                                                                                                                                >  1:523583/1‑151 (MQ=255)            cAGTAAAGCACGTTAATTTTATCCGCCCATTGGCTAAGCATTTGTAAGTCATGGCTGATAAGCAAAATAGTGGTATTGCTGTTTTGGTTGAGACGCGTCAGCAGGCGAAAGATTTGCGCCTGGGTTGTTGGCTCAATTGAGTTGGTCGGtt                                                                                                                             >  2:1060921/1‑151 (MQ=255)                 aaGCACGTTAATTTTATCCGCCCATTGGCTAAGCATTTGTAAGTCATGGCTGATAAGCAAAATAGTGGTATTGCTGTTTTGGTTGAGACGCGTCAGCAGGCGAAAGATTTGCGCCTGGGTTGTTGGCTCAATTGAGTTGGTCGGTTCGTCa                                                                                                                        <  2:849742/151‑1 (MQ=255)                     aCGTTAATTTTATCCGCCCATTGGCTAAGCATTTGTAAGTCATGGCTGATAAGCAAAATAGTGGTATTGCTGTTTTGGTTGAGACGCGTCAGCAGGCGAAAGATTTGCGCCTGGGTTGTTGGCTCAATTGAGTTGGTCGGTTCGTCAGCa                                                                                                                     >  2:504072/1‑150 (MQ=255)                       gTTAATTTTATCCGCCCATTGGCTAAGCATTTGTAAGTCATGGCTGATAAGCAAAATAGTGGTATTGCTGTTTTGGTTGAGACGCGTCAGCAGGCGAAAGATTTGCGCCTGGGTTGTTGGCTCAATTGAGTTGGTCGGTTCGTCAGCAATc                                                                                                                  >  1:37089/1‑151 (MQ=255)                             tttATCCGCCCATTGGCTAAGCATTTGTAAGTCATGGCTGATAAGCAAAATAGTGGTATTGCTGTTTTGGTTGAGACGCGTCAGCAGGCGAAAGATTTGCGCCTGGGTTGTTGGCTCAATTGAGTTGGTCGGTTCGTCAGCAATcagc                                                                                                               >  1:257574/1‑148 (MQ=255)                              ttATCCGCCCATTGGCTAAGCATTTGTAAGTCATGGCTGATAAGCAAAATAGTGGTATTGCTGTTTTGGTTGAGACGCGTCAGCAGGCGAAAGATTTGCGCCTGGGTTGTTGGCTCAATTGAGTTGGTCGGTTCGTCAGCAATCAGCAGAc                                                                                                           >  2:831546/1‑151 (MQ=255)                                aTCCGCCCATTGGCTAAGCATTTGTAAGTCATGGCTGATAAGCAAAATAGTGGTATTGCTGTTTTGGTTGAGACGCGTCAGCAGGCGAAAGATTTGCGCCTGGGTTGTTGGCTCAATTGAGTTGGTCGGTTCGTCAGCAATCAGCAGAcg                                                                                                          <  1:605239/150‑1 (MQ=255)                                  ccGCCCATTGGCTAAGCATTTGTAAGTCATGGCTGATAAGCAAAATAGTGGTATTGCTGTTTTGGTTGAGACGCGTCAGCAGTCGAAAGATTTGCGCCTGGGTTGTTGGCTCAATTTAGTTGGTCGGTTCGTCAGCAATCAGCAGACGCgg                                                                                                       <  1:927332/151‑1 (MQ=255)                                    gCCCATTGGCTAAGCATTTGTAAGTCATGGCTGATAAGCAAAATAGTGGTATTGCTGTTTTGGTTGAGACGCGTCAGCAGGCGAAAGATTTGCGCCTGGGTTGTTGGCTCAATTGAGTTGGTCGGTTCGTCAGCAATCAGCAGACGCGGtt                                                                                                     <  1:43645/151‑1 (MQ=255)                                                aGCATTTGTAAGTCATGGCTGATATGCAAAATAGTGGTATTGCTGTTTTGGTTGAGACGCGTCAGCAGGCGAAAGATTTGCGCCTGGGTTGTTGGCTCAATTGAGTTGGTCGGTTCGTCAGCAATCAGCAGACGCGGTTGATTCGCCAGTg                                                                                         <  1:226339/151‑1 (MQ=255)                                                         aaGTCATGGCTGATAAGCAAAATAGTGGTATTGCTGTTTTGGTTGAGACGCGTCAGCAGGCGAAAGATTTGCGCCTGGGTTGTTGGCTCAATTGAGTTGGTCGGTTCGTCAGCAATCAGCAGACGCGGTTGATTCGCCAGTGCAATGGCTa                                                                                <  1:908750/151‑1 (MQ=255)                                                            tCATGGCTGATAAGCAAAATAGTGGTATTGCTGTTTTGGTTGAGACGCGTCAGCAGGCGAAAGATTTGCGCCTGGGTTGTTGGCTCAATTGAGTTGGTCGGTTCGTCAGCAATCAGCAGACGCGGTTGATTCGCCAGTGCAATGGCTatc                                                                              >  1:493338/1‑150 (MQ=255)                                                            tCATGGCTGATAAGCAAAATAGTGGTATTGCTGTTTTGGTTGAGACGCGTCAGCAGGCGAAAGATTTGCGCCTGGGTTGTTGGCTCAATTGAGTTGGTCGGTTCGTCAGCAATCAGCAGACGCGGTTGATTCGCCAGTGCAATGGCTa                                                                                >  2:440693/1‑148 (MQ=255)                                                             cATGGCTGATAAGCAAAATAGTGGTATTGCTGTTTTGGTTGAGACGCGTCAGCAGGCGAAAGATTTGCGCCTGGGTTGTTGGCTCAATTGAGTTGGTCGGTTCGTCAGCAATCAGCAGACGCGGTTGATTCGCCAGTGCAATGGCTatcat                                                                            <  1:891031/151‑1 (MQ=255)                                                                          cAAAATAGTGGTATTGCTGTTTTGGTTGAGACGCGTCAGCAGGCGAAAGATTTGCGCCTGGGTTGTTGGCTCAATTGAGTTGGTCGGTTCGTCAGCAATCAGCAGACGCGGTTGATTCGCCAGTGCAATGGCTATCATCACTTTCTGACAt                                                               >  2:1071350/1‑151 (MQ=255)                                                                              aTAGTGGTATTGCTGTTTTGGTTGAGACGCGTCAGCAGGCGAAAGATTTGCGCCTGGGTTGTTGGCTCAATTGAGTTGGTCGGTTCGTCAGCAATCAGCAGACGCGGTTGATTCGCCAGTGCAATGGCTATCATCACTTTCTGACATTCAc                                                           >  2:1030319/1‑151 (MQ=255)                                                                               tAGTGGTATTGCTGTTTTGGTTGAGACGCGTCAGCAGGCGAAAGATTTGCGCCTGGGTTGTTGGCTCAATTGAGTTGGTCGGTTCGTCAGCAATCAGCAGACGCGGTTGATTCGCCAGTGCAATGGCTATCATCACTTTCTGACATTCAcc                                                          >  2:670448/1‑151 (MQ=255)                                                                               tAGTGGTATTGCTGTTTTGGTTGAGACGCGTCAGCAGGCGAAAGATTTGCGCCTGGGTTGTTGGCTCAATTGAGTTGGTCGGTTCGTCAGCAATCAGCAGACGCGGTTGATTCGCCAGTGCAATGGCTATCATCACTTTCTGACATTCAcc                                                          >  2:354764/1‑151 (MQ=255)                                                                                aGTGGTATTGCTGTTTTGGTTGAGACGCGTCAGCAGGCGAAAGATTTGCGCCTGGGTTGTTGGCTCAATTGAGTTGGTCGGTTCGTCAGCAATCAGCAGACGCGGTTGATTCGCCAGTGCAATGGCTATCATCACTTTCTGACATTCACCt                                                         >  2:904605/1‑151 (MQ=255)                                                                                   ggTATTGCTGTTTTGGTTGAGACGCGTCAGCAGGCGAAAGATTTGCGCCTGGGTTGTTGGCTCAATTGAGTTGGTCGGTTCGTCAGCAATCAGCAGACGCGGTTGATTCGCCAGTGCAATGGCTATCATCACTTTCTGACATTCACCTTc                                                       <  2:776081/150‑1 (MQ=255)                                                                                   ggTATTGCTGTTTTGGTTGAGACGCGTCAGCAGGCGAAAGATTTGCGCCTGGGTTGTTGGCTCAATTGAGTTGGTCGGTTCGTCAGCAATCAGCAGACGCGGTTGATTCGCCAGTGCAATGGCTATCATCACTTTCTGACATTCACCTTCg                                                      <  2:866875/151‑1 (MQ=255)                                                                                        tGCTGTTTTGGTTGAGACGCGTCAGCAGGCGAAAGATTTGCGCCTGGGTTGTTGGCTCAATTGAGTTGGTCGGTTCGTCAGCAATCAGCAGACGCGGTTGATTCGCCAGTGCAATGGCTATCATCACTTTCTGACATTCACCTTCGGTCa                                                  <  2:1043970/150‑1 (MQ=255)                                                                                        tGCTGTTTTGGTTGAGACGCGTCAGCAGGCGAAAGATTTGCGCCTGGGTTGTTGGCTCAATTGAGTTGGTCGGTTCGTCAGCAATCAGCAGACGCGGTTGATTCGCCAGTGCAATGGCTATCATCACTTTCTGACATTCACCTTCGGTCa                                                  <  1:635590/150‑1 (MQ=255)                                                                                              tttGGTTGAGACGCGTCAGCAGGCGAAAGATTTGCGCCTGGGTTGTTGGCTCAATTGAGTTGGTCGGTTCGTCAGCAATCAGCAGACGCGGTTGATTCGCCAGTGCAATGGCTATCATCACTTTCTGACATTCACCTTCGGTCAACTCAt                                            >  2:81606/1‑150 (MQ=255)                                                                                              tttGGTTGAGACGCGTCAGCAGGCGAAAGATTTGCGCCTGGGTTGTTGGCTCAATTGAGTTGGTCGGTTCGTCAGCAATCAGCAGACGCGGTTGATTCGCCAGTGCAATGGCTATCATCACTTTCTGACATTCACCTTCGGTCAACTCAt                                            <  1:168811/150‑1 (MQ=255)                                                                                                tGGTTGAGACGCGTCAGCAGGCGAAAGATTTGCGCCTGGGTTGTTGGCTCAATTGAGTTGGTCGGTTCGTCAGCAATCAGCAGACGCGGTTGATTCGCCAGTGCAATGGCTATCATCACTTTCTGACATTCACCTTCGGTCAACTCATAg                                          <  2:190804/150‑1 (MQ=255)                                                                                                 ggTTGAGACGCGTCAGCAGGCGAAAGATTTGCGCCTGGGTTGTTGGCTCAATTGAGTTGGTCGGTTCGTCAGCAATCAGCAGACGCGGTTGATTCGCCAGTGCAATGGCTATCATCACTTTCTGACATTCACCTTCGGTCAACTCATAggg                                        <  2:521338/151‑1 (MQ=255)                                                                                                    tGAGACGCGTCAGCAGGCGAAAGATTTGCGCCTGGGTTGTTGGCTCAATTGAGTTGGTCGGTTCGTCAGCAATCAGCAGACGCGGTTGAGTCGCCAGTGCAATGGCTATCATCACTGTCTGACATTCACCTTCGGTCAACTCATAGGGaa                                      <  2:681362/150‑1 (MQ=255)                                                                                                       gaCGCGTCAGCAGGCGAAAGATTTGCGCCTGGGTTGTTGGCTCAATTGAGTTGGTCGGTTCGTCAGCAATCAGCAGACGCGGTTGATTCGCCAGTGCAATGGCTATCATCACTTTCTGACATTCACCTTCGGTCAACTCATAGGGAAAACt                                  <  2:590283/151‑1 (MQ=255)                                                                                                          gcgTCAGCAGGCGAAAGATTTGCGCCTGGGTTGTTGGCTCAATTGAGTTGGTCGGTTCGTCAGCAATCAGCAGACGCGGTTGATTCGCCAGTGCAATGGCTATCATCACTTTCTGACATTCACCTTCGGTCAACTCATAGGGAAAACTgcg                               >  2:1045941/1‑151 (MQ=255)                                                                                                           cgTCAGCAGGCGAAAGATTTGCGCCTGGGTTGTTGGCTCAATTGAGTTGGTCGGTTCGTCAGCAATCAGCAGACGCGGTTGATTCGCCAGTGCAATGGCTATCATCACTTTCTGACATTCACCTTCGGTCAACTCATAGGGAAAACTgcgc                              <  2:37088/151‑1 (MQ=255)                                                                                                              cagcagGCGAAAGATTTGCGCCTGGGTTGTTGGCTCAATTGAGTTGGTCGGTTCGTCAGCAATCAGCAGACGCGGTTGATTCGCCAGTGCAATGGCTATCATCACTTTCTGACATTCACCTTCGGTCAACTCATAGGGAAAACTGcgcatc                           >  1:872082/1‑151 (MQ=255)                                                                                                              cagcagGCGAAAGATTTGCGCCTGGGTTGTTGGCTCAATTGAGTTGGTCGGTTCGTCAGCAATCAGCAGACGCGGTTGATTCGCCAGTGCAATGGCTATCATCACTTTCTGACATTCACCTTCGGTCAACTCATAGGGAAAACTGcgcatc                           <  1:831523/151‑1 (MQ=255)                                                                                                               agcagGCGAAAGATTTGCGCCTGGGTTGTTGGCTCAATTGAGTTGGTCGGTTCGTCAGCAATCAGCAGACGCGGTTGATTCGCCAGTGCAATGGCTATCATCACTTTCTGACATTCACCTTCGGTCAACTCATAGGGAAAACTGcgcatcg                          <  1:472674/151‑1 (MQ=255)                                                                                                                  agaataaaGATTTGCCCCTGGGTTGGTGGCTCAATTGAGTTGGTCGGTTCGTGAACAATCAGCAGACGCGGCTGATTCGCCAGTGCAATGGCCATCATCACTTTCTGACATTCACCTTCGGTCAACTCCTAGGGCAACCTGcgcatc                           >  2:376187/6‑147 (MQ=255)                                                                                                                    gCGAAAGATTTGCGCCTGGGTTGTTGGCTCAATTGAGTTGGTCGGTTCGTCAGCAATCAGCAGACGCGGTTGATTCGCCAGTGCAATGGCTATCATCACTTTCTGACATTCACCTTCGGTCAACTCATAGGGAAAACTGCGCATCGCATCt                     <  1:1108944/151‑1 (MQ=255)                                                                                                                      gAAAGATTTGCGCCTGGGTTGTTGGCTCAATTGAGTTGGTCGGTTCGTCAGCAATCAGCAGACGCGGTTGATTCGCCAGTGCAATGGCTATCATCACTTTCTGACATTCACCTTCGGTCAACTTATAGGGAAAACTGCGCATCGCATCttt                   <  1:880378/151‑1 (MQ=255)                                                                                                                      gAAAGATTTGCGCCTGGGTTGTTGGCTCAATTGAGTTGGTCGGTTCGTCAGCAATCAGCAGACGCGGTTGATTCGCCAGTGCAATGGCTATCATCACTTTCTGACATTCACCTTCGGTCAACTCATAGGGAAAACTGCGCATCGCATCttt                   >  1:510777/1‑151 (MQ=255)                                                                                                                            tttGCGCCTGGGTTGTTGGCTCAATTGAGTTGGTCGGTTCGTCAGCAATCAGCAGACGCGGTTGATTCGCCAGTGCAATGGCTATCATCACTTTCTGACATTCACCTTCGGTCAACTCATAGGGTAAACTGCGCATCGCATCTTTGTGa               >  1:852755/1‑149 (MQ=255)                                                                                                                             ttGCGCCTGGGTTGTTGGCTCAATTGAGTTGGTCGGTTCGTCAGCAATCAGCAGACGCGGTTGATTCGCCAGTGCAATGGCTATCATCACTTTCTGACATTCACCTTCGGTCAACTCATAGGGAAAACTGCGCATCGCATCTTTGTGATc             <  2:877874/150‑1 (MQ=255)                                                                                                                              tGCGCCTGGGTTGTTGGCTCAATTGAGTTGGTCGGTTCGTCAGCAATCAGCAGACGCGGTTGATTCGCCAGTGCAATGGCTATCATCACTTTCTGACATTCACCTTCGGTCAACTCATAGGGAAAACTGCGCATCGCATCTTTGTGATCt            <  2:257586/150‑1 (MQ=255)                                                                                                                                cgcCTGGGTTGTTGGCTCAATTGAGTTGGTCGGTTCGTCAGCAATCAGCAGACGCGGTTGATTCGCCAGTGCAATGGCTATCATCACTTTCTGACATTCACCTTCGGTCAACTCATAGGGAAAACTGCGCATCGCATCTTTGTGATCTTTa         <  2:523597/151‑1 (MQ=255)                                                                                                                                  ccTGGGTTGTTGGCTCAATTGAGTTGGTCGGTTCGTCAGCAATCAGCAGACGCGGTTGATTCGCCAGTGCAATGGCTATCATCACTTTCTGACATTCACCTTCGGTCAACTCATAGGGAAAACTGCGCATCGCATCTTTGTGATCTTTAAt       >  1:65926/1‑151 (MQ=255)                                                                                                                                   cTGGGTTGTTGGCTCAATTGAGTTGGTCGGTTCGTCAGCAATCAGCAGACGCGGTTGATTCGCCAGTGCAATGGCTATCATCACTTTCTGACATTCACCTTCGGTCAACTCATAGGGAAAACTGCGCATCGCATCTTTGTGATCTTTAATc      >  2:29923/1‑151 (MQ=255)                                                                                                                                   cTGGGTTGTTGGCTCAATTGAGTTGGTCGGTTCGTCAGCAATCAGCAGACGCGGTTGATTCGCCAGTGCAATGGCTATCATCACTTTCTGACATTCACCTTCGGTCAACTCATAGGGAAAACTGCGCATCGCATCTTTGTGATCTTTAATc      >  2:1064150/1‑151 (MQ=255)                                                                                                                                     gggTTGTTGGCTCAATTGAGTTGGTCGGTTCGTCAGCAATCAGCAGACGCGGTTGATTCGCCAGTGCAATGGCTATCATCACTTTCTGACATTCACCTTCGGTCAACTCATAGGGAAAACTGCGCATCGCATCTTTGTGATCTTTAATccc    >  2:1074025/1‑151 (MQ=255)                                                                                                                                       gttgttGGCTCAATTGAGTTGGTCGGTTCGTCAGCAATCAGCAGACGCGGTTGATTCGCCAGTGCAATGGCTATCATCACTTTCTGACATTCACCTTCGGTCAACTCATAGGGAAAACTGCGCATCGCATCTTTGTGATCTTTAATCCCCa  <  2:102435/151‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                            CTGTTTGACCGCAGTAAAGCACGTTAATTTTATCCGCCCATTGGCTAAGCATTTGTAAGTCATGGCTGATAAGCAAAATAGTGGTATTGCTGTTTTGGTTGAGACGCGTCAGCAGGCGAAAGATTTGCGCCTGGGTTGTTGGCTCCATTGAGTTGGTCGGTTCGTCAGCAATCAGCAGACGCGGTTGATTCGCCAGTGCAATGGCTATCATCACTTTCTGACATTCACCTTCGGTCAACTCATAGGGAAAACTGCGCATCGCATCTTTGTGATCTTTAATCCCCA  >  NC_000913/1352917‑1353201

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