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.4%	85.6 / ‑5.9	28	L180P (CTG→CCG)	ldtC	L,D‑transpeptidase LdtC
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (14/13);  total (14/14)
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.89e-01

GCTGATGTACTTCAACCAGACGCGCACCGTTTGGTTCGGCAGAGACTTTTATCGGAGTGTTGATGATATTCACTTTGGTGCCTGGGGTGACCTGGCTAAAGAGTGTTTTGATATCGTCATCCCGCAGACGAATACAGCCAGAACTTACCC‑GCATGCCAATGCCGAAATCGGCGTTCGTACCATGAAGCAAATAAACGCCGCCATAGGCCGCCAGACGAATCGCATGATGGCCCATTGGGTTATCCAGTCCAGCCGGCACTACCGCAGGCAACTCAATTCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCCGGCGGGTAGTAATAGAGACGCAGCTCCGCAATGTTGATCACAATGCCT  >  NC_000913/1169591‑1170069                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                          gCTGATGTACTTCAACCAGACGCGCACCGTTTGGTTCGGCAGAGACTTTTATCGGAGTGTTGATGATATTCACTTTGGTGCCTGGGGTGACCTGGCTAAAGAGTGTTTTGATATCATCATCCCGCAGACGAATACAGCCAGAACTTACCC‑GCATGCCAATGCCGAAATCGGCGTTCGTACCATGAAGCAAATAAACGCCGCCATAGGCCGCCAGACGAATCGCATGATGGCCCATTGGGTTAtccggtcc                                                                                                                                                                                                                                       >  1:300063/1‑250 (MQ=255)                 cAGAAGCGCACCGTGTGTGTCGGCAGAGACTTTTATCGGAGTGTTGATGATATTCACTTGGGTGCCTGGGGGGACCTGGCTAAAGAGTGTTTTGATACGGTCAACCCTCAGACGAATACAGCCAGAACTTACCC‑GCATGCCAATGCCGAAAGCGGCGTTCGTACCAGGAAGCAAATAAACGCCGCCATAGGCCGCCAGCCGAATCGCATTATGGCCCATTGGGTTATCCGGTCCAGATGGCACTACCGCAg                                                                                                                                                                                                                      <  2:161979/251‑1 (MQ=255)                                            gatttttATCGTACTGTTGATTATAATCACTTTGGTGACTGGGGTGTCCTGGCTTAAGAGGGTTTTGATATCGGCACCCCGCAGACGAACACAACCAGAACTTATCC‑GCATTCCAATGCCGAAATCGGCGTTCGTACCATGAAGCATATAAACGCCGCCATAGGCCGCCAGACTAATCGCATGATTGCCCATTGGTTTACCCGTTTCAGCGGGCACTACCGCAGGCAGCTCAAGTTCCtgttgttgataac                                                                                                                                                                                           <  2:267011/248‑11 (MQ=255)                                            gaCTTTTATCGGAGTGTTGATGATATTCACTTTGGTGCCTGGGGTGACCTGGCTAAAGAGTGTTTTGATATCGTCATCCCGCAGACGAATACAGCCAGAACTTACCC‑GCATGCCAATGCCGAAATCGGCGTTCGTACCATGAAGCAAATAAACGCCGCCATAGGCCGCCAGACGAATCGCATGATGGCCCATTGGGTTATCCGGTCCAGCCGGCACTACCGCAGGCAACTCAATTCCCTGTGCTTTATAAc                                                                                                                                                                                           >  1:151576/1‑251 (MQ=255)                                             aCTTTTATCGGAGTGTTGATGATATTCACTTTGGTGCCTGGGGTGACCTGGCTAAAGAGTGTTTTGATATCGTCATCCCGCAGACGAATACAGCCAGAACTTACCC‑GCATGCCAATGCCGAAATCGGCGTTCGTACCATGAAGCAAATAAACGCCGCCATAGGCCGCCAGACGAATCGCATGATGGCCCATTGGGTTATCCGGTCCAGCCGGCACTACCGCAGGCAACTCAATTCCCTGTGCTTTATAACg                                                                                                                                                                                          >  2:83124/1‑251 (MQ=255)                                                                            tGGTGCCTGGGGTGACCTGGCTAAAGAGTGTTTTGATATCGTCATCCCGCAGACGAATACAGCCAGAACTTACCC‑GCATGCCAATGCCGAAATCGGCGTTCGTACCATGAAGCAAATAAACGCCGCCATAGGCCGCCAGACGAATCGCATGATGGCCCATTGGGTTATCCGGTCCAGCCGGCACTACCGCAGGCAACTCAATTCCCTGTGCttt                                                                                                                                                                                                >  2:84568/1‑214 (MQ=255)                                                                            tGGTGCCGGGGGTGACCTGGCTAAAGAGTGTTTTGATATCGTCATCCCGCAGACGAATACAGCCAGAACTTACCC‑GCATGCCAATGCCGAAATCGGCGTTCGTACCATGAAGCAAATAAACGCCGCCATAGGCCGCCAGACGAATCGCATGATGGCCCATTGGGTTATCCGGTCCAGCCGGCACTACCGCAGGCAACTCAATTCCCTGTGCttt                                                                                                                                                                                                <  1:84568/214‑1 (MQ=255)                                                                             ggTGCGTGGGGGGACCTGGCTAAAGAGTGTTTTGATATCGTACTCCGGCAGACTAATAGAACTAGAACT‑ACCCTGCAAGCCAATGCGGAGATGGGCGTTCGAACCATGAAGCATATAATCGCCGCCATAGGGCGCCAGACGAATCGCGTGATGGCCCATGGGGTTATCCGGTCCAGCCGGCACTACCGCAGGCAACTCAATTCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGGCCAGGTTg                                                                                                                                                          <  2:455391/251‑1 (MQ=255)                                                                             ggTGCCTGGGGTGCCCGGGCTAAAGAGTGTTTGGATATCGGCATCCCGCAGACGAATACCGCCAGAACTTACCC‑GCATGCCAATGCCGAAATCGGCGTTCGTACCATGAAGCAAATAAACGCCGCCATAGGCCGCCAGACGAATCGCATGATGGCCCATTGGGTTATCCGGTCCAGCCGGCACTTCCGCAGGCAACTCAATTCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGtt                                                                                                                                                           <  2:268169/250‑1 (MQ=255)                                                                                      ggTGACCTGGCTAAAGAGTGTTTTGATATCGTCATCCCGCAGACGAATACAGCCAGAACTTACCC‑GCATGCCAATGCCGAAATCGGCGTTCGTACCATGAAGCAAATAAACGCCGCCATAGGCCGCCAGACGAATCGCATGCTGGCCCATTGGGTTATCCGGTCCAGCCGGCACTACCGCAGGCAACTCAATTCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCAc                                                                                                                                                 >  1:95718/1‑251 (MQ=255)                                                                                              ggCTAAAGAGTGTTTTGATATCGTCATCCCGCAGACGAATACAGCCAGAACTTACCC‑GCATGCCAATGCCGAAATCGGCGTTCGTACCATGAAGCAAATAAACGCCGCCATAGGCCGCCAGACGAATCGCATGATGGCCCATTGGGTTATCCGGTCCAGCCGGCACTACCGCAGGCAACTCAATTCCCt                                                                                                                                                                                                       >  1:114224/1‑189 (MQ=255)                                                                                              ggCTAAAGAGTGTTTTGATATCGTCATCCCGCAGACGAATACAGCCAGAACTTACCC‑GCATGCCAATGCCGAAATCGGCGTTCGTACCATGAAGCAAATAAACGCCGCCATAGGCCGCCAGACGAATCGCATGATGGCCCATTGGGTTATCCGGTCCAGCCGGCACTACCGCAGGCAACTCAATTCCCt                                                                                                                                                                                                       <  2:114224/189‑1 (MQ=255)                                                                                                                tatCGTCATCCCGCAGACGAATACAGCCAGAACTTACCC‑GCATGCCACTGCCGAAATCGGCGTTCGTACCATGAAGCAAATAAACGCCGCCATAGGCCGCCAGCCGACTCGCATGATGGCCCATTGGGTTATCCGGTCCAGCCGGCACTACCGCAGGCAACTCAATTCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAACCGGTGGTCACCAtt                                                                                                                       >  1:306282/1‑251 (MQ=255)                                                                                                                     tCATCCCGCAGACGAATACAGCCAGAACTTACCC‑GCATGCCAATGCAGAAATCGGCGTTCGTACCATGAAGCAAATAAACGCCGCCATCGGCCGCCAGACGAATCGCATGATGGCCCATTGGGTTATCCGGTCCAGCCGGCACTACCGCAGGCAACTCAATTCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCGTAAACGGTGGTCCCCCTTGTCgg                                                                                                                  >  2:149227/1‑251 (MQ=255)                                                                                                                           cccAGACGAATACAGCCAGAACTTACCC‑GCATGCCAATGCCGAAAGTGGCGTTCGTACCATGAAGCAAATAAACGCCGCCATAGGCCGCCAGACGAATCGCATGATGGCCCATTGTGTTATCCGGTCCAGCCGGCACTACCGCAGTCAACTCAATTCCCTGTGCTTTATAGCGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGGCTGAAACGGTGGTGACCATTGGCGGTGGCAg                                                                                                            <  2:267703/249‑1 (MQ=255)                                                                                                                                                          aTGCCAATGCCGAAATCGGCGTTCGTACCATGAAGCAAATAAACGCCGCCATAGGCCGCCAGACGAATCGCATGATGGCCCATTGGGTTATCCGGTCCAGCCGGCACTACCGCAGGCAACTCAATTCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTca                                                                                                     <  1:191711/228‑1 (MQ=255)                                                                                                                                                          aTGCCAATGCCGAAATCGGCGTTCGTACCATGAAGCAAATAAACGCCGCCATAGGCCGCCAGACGAATCGCATGATGGCCCATTGAGTTATCCGGTCCAGCCGGCACTACCGCAGGCAACTCAATTCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTca                                                                                                     >  2:191711/1‑228 (MQ=255)                                                                                                                                                              cAATGCCGAAATCGGCGTTCGTACCATGAAGCAAATAAACGCCGCCATAGGCCGCCAGACGAATCGCATGATGGCCCATTGGGTTATCCGGTCCAGCCGGCACTACCGCAGGCAACTCAATTCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTgg                                                                           >  2:17523/1‑250 (MQ=255)                                                                                                                                                                aTGCCGAAATCGGCGTTCGTACCATGAAGCAAATAAACGCCGCCATAGGCCGCCAGACGAATCGCATGATGGCCCATTGGGTTATCCGGTCCAGCCGGCACTACCGCAGGCAACTCAATTCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGAt                                                                         >  1:188818/1‑250 (MQ=255)                                                                                                                                                                          cGGCGTTCGTACCATGAAGCAAATAAACGCCGCCATAGGCCGCCAGACGAATCGCATGATGGCCCATTGGGATATCCGGTCCAGCCGGCACTACCGCAGGCCACTCAATTCCCTGTGCTTTATCACGGGCGCGGATGTTTGCCGTTGGCGTCCCGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTAcc                                                              >  1:414427/1‑251 (MQ=255)                                                                                                                                                                                tCGTACCATGAAGCAAATAAACGCCGCCATAGGCCGCCAGACGAATCGCATGATGGCCCATTGGGTTATCCGGTCCAGCCGGCACTACCGCAGGCAACTCAATTCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAAttt                                                        >  1:116407/1‑251 (MQ=255)                                                                                                                                                                                    aCCATGAAGCAAATAAACGCCGCCATAGGCCGCCAGACGAATCGCATGATGGCCCATTGGGTTATCCGGTCCAGCCGGCACTACCGCAGGCAACTCAATTCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTAc                                                    <  1:453954/251‑1 (MQ=255)                                                                                                                                                                                      cATGAATCAAATAAACGCCGCCATAGGCCGCCAGACGAATCGCATGATGGCCCATTGGGTTATCCGGTCCAGCCGGCTCTACCGCAGGCAACTCAATTCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTAccc                                                  <  2:368499/251‑1 (MQ=255)                                                                                                                                                                                                 tAAACGCCGCCATAGGCCGCCAGACGAATCGCATGATGGCCCATTGGGTTATCCGGTCCAGCCGGCACTACCGCAGGCAACTCAATTCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCCGGCGGgtagt                                        >  1:224933/1‑250 (MQ=255)                                                                                                                                                                                                    aCGCCGCCATAGGCCGCCATACGAAGCGCATGATGGCCCATTGGGTTTTCCGGTCCAGCCGGCACTACCGCAGGCAACTCAATTCCCTGTGCGTTATAACGGGCGCGGATGTGTGCCGTTGGCTTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGTTGTTCACCATTGTCGGTGTCAGCGTGTCACCACCGAACTGACCAATACCTATGGGATACACGTTTACCGAATTTTTTCCCGGCTGGTAGTAATa                                    <  2:306282/251‑1 (MQ=255)                                                                                                                                                                                                      gccgcaATAGGCCGCCAGACGAATCGCCTGATGGCCCATTGGGTTATCCGGTCCAGCCGGCACTACCGCAGGCAACTCAATTCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACATATTGGATACACGGTTACCGAATTTTTACCCGGCGGGTAGTAAt                                     >  2:276720/1‑248 (MQ=255)                                                                                                                                                                                                                          gAATCGCATGAGTGCCCATTGGGTTATCCGGGCCAGCCGGCACTACCGCAGGCAACTCAATTCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATGTTTACCCGGCGGGTAGTAATAGAGACGCAGCTCCGCAATGTTg              <  2:95718/251‑1 (MQ=255)                                                                                                                                                                                                                              cGCATGATGGCCCATTGGGTTATCCGGTCCAGCCGGCACTACCGCAGGCAACTCAATTCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCCGGCGGGTAGTAATAGAGACGCAGCTCCGCAATGTTGATca          <  2:18677/251‑1 (MQ=255)                                                                                                                                                                                                                                      ggCCCATCGGGTTATACAGTCCAGCCGGCACTACCGCAGTCAACTCAATTCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCATGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGTTGTCATCGTGTCCCCTCCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCCGGCGTGTAGTAATAGAGACGCAGCTCCGCAATTTTGATCCCAATGCCt  <  1:181310/251‑1 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                          GCTGATGTACTTCAACCAGACGCGCACCGTTTGGTTCGGCAGAGACTTTTATCGGAGTGTTGATGATATTCACTTTGGTGCCTGGGGTGACCTGGCTAAAGAGTGTTTTGATATCGTCATCCCGCAGACGAATACAGCCAGAACTTACCC‑GCATGCCAATGCCGAAATCGGCGTTCGTACCATGAAGCAAATAAACGCCGCCATAGGCCGCCAGACGAATCGCATGATGGCCCATTGGGTTATCCAGTCCAGCCGGCACTACCGCAGGCAACTCAATTCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCCGGCGGGTAGTAATAGAGACGCAGCTCCGCAATGTTGATCACAATGCCT  >  NC_000913/1169591‑1170069

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