Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,792,267	C→T	*155* (TAG→TAA)	nlpC ←	putative C40 clan peptidase lipoprotein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,792,267	0	C	T	100.0%	131.6 / NA	41	*155* (TAG→TAA)	nlpC	putative C40 clan peptidase lipoprotein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (20/21);  total (20/21)

TTTTTCAGGGGTCCTTTTTTTAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCACTCCCTTGCTGGTAGAGGCGTGGATAAATTGGTTGTTGGTATCGTAAATACCTACATGCAAACCATTTTG  >  NC_000913/1792131‑1792401                                                                                                                                         |                                                                                                                                       tttttCAGGGGTCCTTTTTTTAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATATTa                                                                                                                                      >  1:222202/1‑139 (MQ=255)             ccTTTTTTTAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATATTAGATTCGTCTCGc                                                                                                                          >  2:823891/1‑139 (MQ=255)                                           aataatCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATATTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACAtt                                                                                            >  1:1321886/1‑139 (MQ=255)                                             taatCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATATTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTAt                                                                                          <  2:759930/139‑1 (MQ=255)                                             taatCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATATTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTAt                                                                                          <  1:1607685/139‑1 (MQ=255)                                              aatCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATATTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATc                                                                                         <  2:563833/139‑1 (MQ=255)                                                               tttGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATATTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCa                                                                        <  2:222202/139‑1 (MQ=255)                                                                ttGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATATTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCAc                                                                       >  2:160363/1‑139 (MQ=255)                                                                    cTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATATTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCACTccc                                                                   <  1:160363/139‑1 (MQ=255)                                                                                  tGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATATTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCACTCCCTTGCTGGTAGAGGc                                                     <  2:839691/139‑1 (MQ=255)                                                                                      ttCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATATTAGATTCGTCTCGCCTGCCAGAAAtttt                                                                                                            >  2:673534/1‑80 (MQ=255)                                                                                      ttCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATATTAGATTCGTCTCGCCTGCCAGAAAtttt                                                                                                            <  1:673534/80‑1 (MQ=255)                                                                                         aGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATATTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCa                                                                        >  2:1539640/1‑113 (MQ=255)                                                                                         aGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATATTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCa                                                                        <  1:1539640/113‑1 (MQ=255)                                                                                           gAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATATTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCACTCCCTTGCTGGTAGAGGCGTGGATaaa                                            >  1:17622/1‑139 (MQ=255)                                                                                              tAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATATTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCACTCCCTTGCTGGTAGAGGCGTGGATAAATTg                                         >  2:1579577/1‑139 (MQ=255)                                                                                                       ggCGATGTATCATCGCCGTCCATTTTTCGACATATTAGATTCGTCTCGCCTGcc                                                                                                                     <  1:280415/54‑1 (MQ=255)                                                                                                       ggCGATGTATCATCGCCGTCCATTTTTCGACATATTAGATTCGTCTCGCCTGcc                                                                                                                     >  2:280415/1‑54 (MQ=255)                                                                                                        gCGATGTATCATCGCCGTCCATTTTTCGACATATTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCACTCCCTTGCTGGTAGAGGCGTGGATAAATTGGTTGTTGGTa                               <  2:331185/139‑1 (MQ=255)                                                                                                        gCGATGTATCATCGCCGTCCATTTTTCGACATATTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCACTCCCTTGCTGGTAGAGGCGTGGATAAATTGGTTGTTGGTa                               <  2:37134/139‑1 (MQ=255)                                                                                                        gCGATGTATCATCGCCGTCCATTTTTCGACATATTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCACTCCCTTGCTGGTAGAGGCGTGGATAAATTGGTTGTTGGTa                               <  2:438561/139‑1 (MQ=255)                                                                                                               atcatcGCCGTCCATTTTTCGACATATTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCACTCCCTTGCTGGTAGAGGCGTGGATAAATTGGTTGTTGGTATCGTaaa                        <  1:104217/139‑1 (MQ=255)                                                                                                                   tcGCCGTCCATTTTTCGACATATTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCACTcc                                                                    <  2:1575272/91‑1 (MQ=255)                                                                                                                   tcGCCGTCCATTTTTCGACATATTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCACTcc                                                                    >  1:1575272/1‑91 (MQ=255)                                                                                                                      ccGTCCATTTTTCGACATATTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGAc                                                                                               >  2:1347015/1‑61 (MQ=255)                                                                                                                      ccGTCCATTTTTCGACATATTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGAc                                                                                               <  1:1347015/61‑1 (MQ=255)                                                                                                                      ccGTCCATTTTTCGACATATTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCa                                                                        <  2:1427220/84‑1 (MQ=255)                                                                                                                      ccGTCCATTTTTCGACATATTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCa                                                                        >  1:1427220/1‑84 (MQ=255)                                                                                                                      ccGTCCATTTTTCGACATATTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCACTCCCTTGCTGGTAGAGGCGTGGATAAATTGGTTGTTGGTATCGTAAATACCTAc                 >  2:797819/1‑139 (MQ=255)                                                                                                                      ccGTCCATTTTTCGACATATTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCACTCCCTTGCTGGTAGAGGCGTGGATAAATTGGTTGTTGGTATCGTAAATACCTAc                 >  2:744693/1‑139 (MQ=255)                                                                                                                       cGTCCATTTTTCGACATATTAGATTCGTCTCGCCTGCCAGAAATTTTTCt                                                                                                         <  1:1336761/50‑1 (MQ=255)                                                                                                                       cGTCCATTTTTCGACATATTAGATTCGTCTCGCCTGCCAGAAATTTTTCt                                                                                                         >  2:1336761/1‑50 (MQ=255)                                                                                                                        gTCCATTTTTCGACATATTAGATTCGTCTCGCCTGCCAGAAATTTTTCTg                                                                                                        >  1:16122/1‑50 (MQ=255)                                                                                                                        gTCCATTTTTCGACATATTAGATTCGTCTCGCCTGCCAGAAATTTTTCTg                                                                                                        <  2:16122/50‑1 (MQ=255)                                                                                                                        gTCCATTTTTCGACATATTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCACTCCCTTGCTGGTAGAGGCGTGGATAAATTGGTTGTTGGTATCGTAAATACCTACAt               >  1:218090/1‑139 (MQ=255)                                                                                                                        gTCCATTTTTCGACATATTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCACTCCCTTGCTGGTAGAGGCGTGGATAAATTGGTTGTTGGTATCGTAAATACCTACAt               >  1:1625645/1‑139 (MQ=255)                                                                                                                        gTCCATTTTTCGACATATTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCACTCCCTTGCTGGTAGAGGCGTGGATAAATTGGTTGTTGGTATCGTAAATACCTACAt               <  1:744693/139‑1 (MQ=255)                                                                                                                        gTCCATTTTTCGACATATTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCACTCCCTTGCTGGTAGAGGCGTGGATAAATTGGTTGTTGGTATCGTAAATACCTACAt               <  2:1625645/139‑1 (MQ=255)                                                                                                                          ccATTTTTCGACATATTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCACTccc                                                                   <  2:1008454/85‑1 (MQ=255)                                                                                                                          ccATTTTTCGACATATTAGATTCGTCTCGCCTGCCAGAAAATTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCACTccc                                                                   >  1:1008454/1‑85 (MQ=255)                                                                                                                                     cATATTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCACTCCCTTGCTGGTAGAGGCGTGGATAAATTGGTTGTTGGTATCGTAAATACCTACATGCAAACCATTTTg  >  2:1349780/1‑139 (MQ=255)                                                                                                                                         |                                                                                                                                       TTTTTCAGGGGTCCTTTTTTTAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCACTCCCTTGCTGGTAGAGGCGTGGATAAATTGGTTGTTGGTATCGTAAATACCTACATGCAAACCATTTTG  >  NC_000913/1792131‑1792401

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