Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,169,836	A→G	L180P (CTG→CCG)	ldtC ←	L,D‑transpeptidase LdtC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,169,836	0	A	G	96.9%	107.7 / NA	32	L180P (CTG→CCG)	ldtC	L,D‑transpeptidase LdtC
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (14/17);  minor base T (1/0);  total (15/17)
Fisher's exact test for biased strand distribution p-value = 4.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

CGTTTGGTTCGGCAGAGACTTTTATCGGAGTGTTGATGATATTCACTTTGGTGCCTGGGGTGACCTGGCTAAAGAGTGTTTTGATATCGTCATCCCGCAGACGAATACAGCCAGAACTTACCCGCATGCCAATGCCGAAATCGGCGTTCGTACCATGAAGCAAATAAACGCCGCCATAGGCCGCCAGACGAATCGCATGATGGCCCATTGGGTTATCCAGTCCAGCCGGCACTACCGCAGGCAACTCAATTCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCCGGCGGGTAGTAATAGAGACGCAGCTCCGCAATGTTGATCACAATGCCTTCGCGCG  >  NC_000913/1169618‑1170076                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                 cGTTTGGTTCGGCAGAGACTTTTATCGGAGTGTTGATGATATTCACTTTGGTGCCTGGGGTGACCTGGCTAAAGAGTGTTTTGATATCGTCATCCCGCAGACGAATACAGCCAGAACTTACCCGCATGCCAATGCCGAAATCGGCGTTCGTACCATGAAGCAAATAAACGCCGCCATCGGCCGCCAGACGACTCGCATGATGGCCCATTGGGTTATCCGGTCCAGCCGGCACTACCGCAGGCCACTCAAtt                                                                                                                                                                                                                  >  2:280096/1‑251 (MQ=255)               gagaCTTTTATCGGAGTGTTGATGATATTCACTTTGGTGCCTGGGGTGACCTGGCTAAAGAGTGTTTTGATATCGTCATCCCGCAGACGAATACAGCCAGAACTTACCCGCATGCCAATGCCGAAATCGGCGTTCGTACCATGAAGCAAATAAACGCCGCCATAGGCCGCCAGACGAATCGCATGATGGCCCATTGGGTTATCCGGTCCAGCCGGCACTACCGCAGGCAACTCAATTCCCTGTGCTTtat                                                                                                                                                                                                     >  1:156557/1‑250 (MQ=255)                agaCTTTTATCGGAGTGTTGATGATATTCACTTTGGTGCCTGGGGTGACCTGGCTAAAGAGTGTTTTGATATCGTCATCCCGCAGACGAATACAGCCAGAACTTACCCGCATGCCAATGCCGAAATCGGCGTTCGTACCATGAAGCAAATAAACGCCGCCATAGGCCGCCAGACGAATCGCATGATGGCCCATTGGGTTATCCGGTCCAGCCGGCACTACCGCAGGCAACTCAATTCCCTGTGCTTTATaa                                                                                                                                                                                                   <  1:328603/251‑1 (MQ=255)                       tATCGGAGTGTTGATGATATTCACTTTGGTGCCTGGGGTGACCTGGCTCAAGAGTGCTTTGATATCGTCATCCCGCAGACGAATACAGCCAGAACTTACCCGCATGCCAATGCCGAAATCGGCGTTCGTACCATGAAGCAAATAAACGCCGCCATAGGCCGCCAGACGAATCGCATGATGGCCCATTGGGTTATCCGGTCCCGCCGGCACTACCGCAGGCAACTCAATTCCCTGTGCTTTCTAACGGgcgc                                                                                                                                                                                            >  1:40936/1‑251 (MQ=255)                              gtgtTGATGATATTCACTTTGGTGCCTGGGGTGACCTGGCTAAAGAGTGTTTTGATATCGTCATCCCGCAGACGAATACAGCCAGAACTTACCCGCATGCCAATGCCGAAATCGGCGTTCGTACCATGAAGCAAATAAACGCCGCCATAGGCCGCCAGACGAATCGCATGATGGCCCATTGGGTTATCCGGTCCAGCCGGCACTACCGCAGGCAACTCAATTCCCTGTGCTTTATAACGGGCGCGGATGtt                                                                                                                                                                                     <  1:569558/251‑1 (MQ=255)                                      gatATTCACTTTGGTGCCTGGGGTGACCTGGCTAAAGAGTGTTTTGATATCGTCATCCCGCAGACGAATACAGCCAGAACTTACCCGCATGCCAATGCCGAAATCGGCGTTCGTACCATGAAGCAAATAAACGCCGCCATAGGCCGCCAGACGAATCGCATGATGGCCCATTGGGTTATCCGGTCCAGCCGGCACTACCGCAGGCAACTCAATTCCCTGTGCTTTATaa                                                                                                                                                                                                   <  1:136570/229‑1 (MQ=255)                                      gatATTCACTTTGGTGCCTGGGGTGACCTGGCTAAAGAGTGTTTTGATATCGTCATCCCGCAGACGAATACAGCCAGAACTTACCCGCATGCCAATGCCGAAATCGGCGTTCGTACCATGAAGCAAATAAACGCCGCCATAGGCCGCCAGACGAATCGCATGATGGCCCATTGGGATATCCGGTCCAGCCGGCACTACCGCAGGCAACTCAATTCCCTGTGCTTTATaaagaacggaagagcgtcggga                                                                                                                                                                               >  2:136570/1‑229 (MQ=255)                                         atTCACTTTGGTGCCTGGGGTGACCTGGCTAAAGAGTGTTTTGATATCGTCATCCCGCAGACGAATACAGCCAGAACTTACCCGCATGCCAATGCCGAAATCGGCGTTCGTACCATGAAGCAAATAAACGCCGCCATAGGCCGCCAGACGAATCGCATGATGGCCCATTGGGTTATCCGGTCCAGCCGGCACTACCGCAGGCAACTCAATTCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGc                                                                                                                                                                           <  1:114677/250‑1 (MQ=255)                                              cTTTGGTGCCTGGGGTGACCTGGCTAAAGAGTGTTTTGATATCGTCATCCCGCAGACGAATACAGCCAGAACTTACCCGCATGCCAATGCCGAAATCGGCGTTCGTACCATGAAGCAAATAAACGCCGCCATAGGCCGCCAGACGAATCGCATGATGGCCCATTGGGTTATCCGGTCCAGCCGGCAc                                                                                                                                                                                                                                     >  1:499048/1‑187 (MQ=255)                                                  ggTGCCTGGGGTGACCTGGCTAAAGAGTGTTTTGATATCGTCATCCCGCAGACGACTACAGCCAGAACTTACCCGCATGCCAATGCCGAAATCGGCGTTCGTACCATGAAGCAAATAAACGCCGCCATAGGCCGCCAGACGAATCGCATGATGGCCCATTGGGTTATCCGGTCCAGCCGGCACTACCGCAGGCAACTCAATTCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTg                                                                                                                                                                 >  1:231799/1‑251 (MQ=255)                                                    ggCCTGGGGTTCCCTGGCTAAAGAGTTGTTTGATAACGTCTGCCCGCAGACGAATACAGCCAGAACCTACCCGCATGCAAATGCCGAAATCGGCGTTCGTCCCATGAAGCAAATAAACGCCGCCATAGGCCGCCAGACGAATCGCATGATGGCCCATTGGGTTATCCGGTCCAGCCGGCACTACCGCAGGCAACTCAATTCCCTGTGCTTTATAACGGGCGCGGAGGTTGGCCGTTGGCGTCCAGGTTggg                                                                                                                                                               <  2:284632/250‑1 (MQ=255)                                                                       aaaGAGTGTTTTGATATCGTCATCCCGCAGACGAATACAGCCAGAACTTACCCGCATGCCAATGCCGAAATCGGCGTTCGTACCATGAAGCAAATAAACGCCGCCATAGGCCGCCAGACGAATCGCATGATGGCCCATTGGGTTATCCGGTCCAGCCGGCACTACCGCAGGCAACTCAATTCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGa                                                                                                                                              <  1:54430/249‑1 (MQ=255)                                                                        aaGAGTGTTTTGATATCGTCATCCCGCAGACGAATACAGCCAGAACTTACCCGCATGCCAATGCCGAAATCGGCGTTCGTACCATGAAGCAAATAAACGCCGCCCTAGGCCGCCAGACGAATCGCATGCAGGCCCATTGGGTTATCCTGTCCAGCCGGCACTACCGCAGGCAACTCAATTCCCTGTGCTTTATAACAGGCTCGGATGTTTGCCGTTGGCCTCCATGTTGGGTTTGCACGTTTTTCTGacac                                                                                                                                           >  2:468/1‑248 (MQ=255)                                                                         agagTGTTTTGATATCGTCATCCCGCAGACGAATACAGCCAGAACTTACCCGCATGCCAATGCCGAAATCGGCGTTCGTACCATGAAGCAAATAAACGCCGCCATAGGCCGCCAGACGAATCGCATGATGGCCCATTGGGTTATCCGGTCCAGCCGGCACTACCGCAGGCAACTCAATTCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACg                                                                                                                                          <  1:458993/251‑1 (MQ=255)                                                                             gtttttGATATCGTCACCCCGCAGACGAATACAGCCAGAACGTACCCGCATGCCAATGCCGAAATCGTCGTGCGTACCATGAAGCAAATAAACGCCGCCATAGGCCGCCAGACGAATCGCATGATGGCCCATTGGGTTATCCGGTCCAGCCGGCACTACCGCAGGCAACTCAATTCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACggtgg                                                                                                                                      <  2:233968/249‑1 (MQ=255)                                                                                                      cGAATACAGCCAGAAATTACCCGCATGCCAATGCCGAAATCGGCGTTCGAACAATGAAGCAAATAAACGCCGCCATAGGCCGCCAGACGAATCGCATGATGGCCCATTGGGTTATCCGGTCCAGCCGGCACTACCGCAGGCAACTCAATTCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTc                                                                                                             >  1:447562/1‑251 (MQ=255)                                                                                                             aGCCAGAACTTACCCGCATGCCAATGCCGAAATCGGCGTTCGTACCATGAAGCAAATAAACGCCGCCATAGGCCGCCAGACGAATCGCATGATGGCCCATTGGGTTATCCGGTCCAGCCGGCACTACCGCAGGCAACTCAATTCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCt                                                                                                      <  1:280096/251‑1 (MQ=255)                                                                                                                                 ccAATGCCGAAATCGGCGTTCGTACCATGAAGCAAATAAACGCCACCATAGGCCGCCAGACGAATCGCATGATGGCCCATTGGGTTATCCGGTCCAGCCGGCACTACCGCAGGCAACTCAATTCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTgg                                                                                  <  1:549862/251‑1 (MQ=255)                                                                                                                                     tGCCGAAATCGGCGTTCGTACCATGAAGCAAATAAACGCCGCCATAGGCCGCCAGACGAATCGCATGATGGCCCATTGGGTTATCCGGTCCAGCCGGCACTACCGCAGGCAACTCAATTCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCCGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATac                                                                              >  1:6090/1‑251 (MQ=255)                                                                                                                                          aaaTCGGCTGTCGTCCCAGGAAGAAATTAAACGCCGCCATAGGCCGCCAGACGAATCGCATGAGGGCCCATTGGGTTATCCGGTCCAGCCGGCACTACCGCAGGCAACTCAATTCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACTTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGt                                                                         <  2:380760/251‑1 (MQ=255)                                                                                                                                          aaaTCGGCGTTCGTACCATGAAGCAAATAAACGCCGCCATAGGCCGCCAGACGAATCGCATGATGGCCCATTGGGTTATCCGGTCCAGCCGGCACTACCGCAGGCAACTCAATTCCCTGTGc                                                                                                                                                                                                          >  2:150261/1‑122 (MQ=255)                                                                                                                                          aaaTCGGCGTTCGTACCATGAAGCAAATAAACGCCGCCATAGGCCGCCAGACGAATCGCATGATGGCCCATTGGGTTATCCGGTCCAGCCGGCACTACCGCAGGCAACTCAATTCCCTGTGc                                                                                                                                                                                                          <  1:150261/122‑1 (MQ=255)                                                                                                                                                                gCAAATAACCGCCGCCATAGGCCTCCAGAAGAATCGCTTGAGGGCCAATTGGGTTATCCGGTCCAGCCGGCACTACCGCAGGCAACTCAATTCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCCGGCggg                                                   <  2:182851/251‑1 (MQ=255)                                                                                                                                                                    aTAAACGCCGCCATAGGCCGCCAGACGAATCGCATGATGGCCCATTGGGTTATCCGGTCCAGCCGGCACTACCGCCGGCAACTCAATTCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAAATGACCAATACCTATTGGATCCACGGTTACCGAATTTTTACCCGGCGGgtagt                                               >  2:98602/1‑251 (MQ=255)                                                                                                                                                                    aTAAACGCCGCCATAGGCCGCCAGACGAATCGCATGATGGCCCATTGGGTTATCCGGTCCAGCCGGCACTACCGCAGGCAACTCAATTCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGGGTCACCCCCTAACTTACCAATCACTATTGGATACCCGGTTACCGAATTTTTACCCGGCGGgtggg                                               >  2:438489/1‑248 (MQ=255)                                                                                                                                                                          gccgccATAGGCCGCCAGACGAATCGCATGATGGCCCATTGGGTTATCCGGTCCAGCCGGCACTACCGCAGGCAACTCAATTCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCCGGCGGGTAGTAATaga                                         <  1:280547/251‑1 (MQ=255)                                                                                                                                                                               cATAGGCCGCCAGACGAATCGCATGATGGCCCATTGGGTTATCCGGTCCAGCCGGCACTACCGCAGGCAACTCAATTCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCCGGCGGGTAGTAATAGAGACGc                                    >  1:104734/1‑251 (MQ=255)                                                                                                                                                                                                  cGCATGATGGCCCATTGGGTTATCCGGTCCAGCCGGCACTACCGCAGGCAACTCAATTCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCCGGCGGGTAGTAATAGAGACGCAGCTCCGCAATGTTGATc                  <  1:359912/250‑1 (MQ=255)                                                                                                                                                                                                    cATGATGGCCCATTGGGTTATCCGGTCCAGCCGGCACTACCGCAGGCAACTCAATTCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCCGGCGGGTAGTAATAGAGACGCAGCTCCGCAATGTTGATcac                >  1:490756/1‑250 (MQ=255)                                                                                                                                                                                                        atgGCCCATTGGGTTATCCGGTCCAGCCGGCACTACCGCAGGCAACTCAATTCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCCCCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCCGGCGGGTAGTAATAAAGACGCAGCTCCGCAATGTTGAACAAAATGc           >  2:103491/1‑251 (MQ=255)                                                                                                                                                                                                             ccATGGGGTTTTCCGGTCCAGCCGGCATTACCGCAGGCAACGCAATTCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCCGGCGGGTAGTAATAGAGACGCAGCTCCGCAATGTTGATCACAATGCCTTcg      <  2:305738/251‑1 (MQ=255)                                                                                                                                                                                                                  gggTTATCCGGTCCAGCCGGCACTACCGCAGGCAACTCAATTCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCCGGCGGGTAGTAATAGAGACGCAGCTCCGCAATGTTGATCACAATGCCTTcgcgcg  <  1:125580/250‑1 (MQ=255)                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                 CGTTTGGTTCGGCAGAGACTTTTATCGGAGTGTTGATGATATTCACTTTGGTGCCTGGGGTGACCTGGCTAAAGAGTGTTTTGATATCGTCATCCCGCAGACGAATACAGCCAGAACTTACCCGCATGCCAATGCCGAAATCGGCGTTCGTACCATGAAGCAAATAAACGCCGCCATAGGCCGCCAGACGAATCGCATGATGGCCCATTGGGTTATCCAGTCCAGCCGGCACTACCGCAGGCAACTCAATTCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCCGGCGGGTAGTAATAGAGACGCAGCTCCGCAATGTTGATCACAATGCCTTCGCGCG  >  NC_000913/1169618‑1170076

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