Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	981,850	(G)6→5	coding (804/1848 nt)	ldtD →	L,D‑transpeptidase LdtD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	981,845	0	G	.	96.2%	98.6 / ‑2.3	26	coding (799/1848 nt)	ldtD	L,D‑transpeptidase LdtD
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base . (14/11);  total (15/11)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATGAATCGTTACTGGCCTTACTCTGTGACACCAAACCGTGGCCCCAACTGACCGGCAAAGCAACGTTGCGCCCAGGGCAGTGGAGTAACGACGTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGACGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGACGCGGTAGTCAGCCCATCCGCTGTTACTGTTGAAACAGCAGAAACTAAGCCGATGGATAAGCAAACGACGTCTCGTAGTAAACC  >  NC_000913/981708‑981972                                                                                                                                          |                                                                                                                                aTGAATCGTTACTGGCCTTACTCTGTGACACCAAACCGTGGCCCCAACTGACCGGCAAAGCAACGTTGCGCCCAGGGCAGTGGAGTAACGACGTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGACggg                                                                                                                               <  1:660799/140‑1 (MQ=255)          tACTGGCCTTACTCTGTGACACCAAACCGTGGCCCCAACTGACCGGCAAAGCAACGTTGCGCCCAGGGCAGTGGAGTAACGACGTACCGGCTTTGCGCGATATATTGCTACGCACAGTCATGTTGGAC‑GGGGGCCGaaaa                                                                                                                     >  2:219942/1‑140 (MQ=255)          tACTGGCCTTACTCTGTGACACCAAACCGTGGCCCCAACTGACCGGCAAAGCAACGTTGCGCCCAGGGCAGTGGAGTAACGACGTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGAC‑GGGGGCCGaaaa                                                                                                                     >  2:664981/1‑140 (MQ=255)          tACTGGCCTTACTCTGTGACACCAAACCGTGGCCCCAACTGACCGGCAAAGCAACGTTGCGCCCAGGGCAGTGGAGTAACGACGTACCGGCCTTGCGCGAAATAGTGCGACGCACAGACATGCTGGAC‑GGGGGCCGaaaa                                                                                                                     >  2:195266/1‑140 (MQ=255)               gCCTTACTCTGTGACACCAAACCGTGGCCCCAACTGACCGGCAAAGCAACGTTGCGCCCAGGGCAGTGGAGTAACGACGTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGAC‑GGGGGCCGAAAATTAct                                                                                                                >  2:140196/1‑140 (MQ=255)                    aCTCTGTGACACCAAACCGTGGCCCCAACTGACCGGCAAAGCAACGTTGCGCCCAGGGCAGTGGAGTAACGACGTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGAC‑GGGGGCCGAAAATTACTCTAcc                                                                                                           >  2:145516/1‑140 (MQ=255)                       ctGTGACACCAAACCGTGGCCCCAACTGACCGGCAAAGCAACGTTGCGCCCAGGGCAGTGGAGTAACGACGTACCGGCGTTGCGCGAAATATTGCAACGCACGGGCATGTTGGA‑GGGGGGCCGAGAATTACTCTACCTgg                                                                                                        >  1:170538/1‑140 (MQ=255)                                cAAACCGTGGCCCCAACTGACCGGCAAAGCAACGTTGCGCCCAGGGCAGTGGAGTAACGACGTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGAC‑GGGGGCCGAAAATTACTCTACCTGGCGATGacac                                                                                               >  2:416337/1‑140 (MQ=255)                                  aaCCGTGGCCCCAACTGACCGGCAAAGCAACGTTGCGCCCAGGGCAGTGGAGTAACGACGTACCGGCGTTGCGCGAAATATTGCGACGCACAGGCATGTTGGAC‑GGGGGCCGAAAATTACTCTACCTGGCGATGACACGc                                                                                             >  2:368655/1‑140 (MQ=255)                                  aaCCGTGGCCCCAACTGACCGGCAAAGCAACGTTGCGCCCAGGGCAGTGGAGTAACGACGTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGAC‑GGGGGCCGAAAATTACTCTACCTGGCGATGACACGc                                                                                             >  2:662045/1‑140 (MQ=255)                                          ccccAACTGACCGGCAAAGCAACGTTGCGCCCAGGGCAGTGGAGTAACGACGTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGAC‑GGGGCCCGAAAAGTACTCTACATGGCCATGACACGCCAACTGAc                                                                                     >  2:284756/1‑140 (MQ=255)                                            ccAACTGACCGGCAAAGCAACGTTGCGCCCAGGGCAGTGGAGTAACGACGTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGAC‑GGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGAcgc                                                                                   <  1:145516/140‑1 (MQ=255)                                            ccAACTGACCGGCAAAGCAACGTTGCGCCCAGGGCAGTGGAGTAACGACGTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGAC‑GGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGAcgc                                                                                   <  2:556187/140‑1 (MQ=255)                                                         aaaGCAACGTTGCGCCCAGGGCAGTGGAGTAACGACGTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGAC‑GGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGACGCGGTAGTCAGCCCa                                                                      <  1:599331/140‑1 (MQ=255)                                                                  ttGCGCCCAGGGCAGTGAAGTAACGACGTACTGGCCTTGCGCGAAATATTGCAACGCACAGGCATGTTGGAC‑GGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGACGCGGTAGTCAGCCCATCCGCTGtt                                                             <  2:250220/140‑1 (MQ=255)                                                                     cgcCCAGGGCAGTGGAGTAACGACGTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGAC‑GGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGACGCGGTAGTCAGCCCATCCGCTGTTACt                                                          >  2:444507/1‑140 (MQ=255)                                                                            ggCAGTGGAGTAACGACGTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGAC‑GGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGACGCGGTAGTCAGCCCATCCGCTGTTACTGTTGaaa                                                   <  1:444507/140‑1 (MQ=255)                                                                                gTGGAGTAACGACGTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGAC‑GGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGACGCGGTAGTCAGCCCATCCGCTGTTACTGTTGAAAcagc                                               <  1:416337/140‑1 (MQ=255)                                                                                 tGGAGTAACGACGTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGAC‑GGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGACGCGGTAGTCAGCCCATCCGCTGTTACTGTTGAAAcagca                                              <  1:140196/140‑1 (MQ=255)                                                                                      tAACGACGTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGAC‑GGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGACGCGGTAGTCAGCCCATCCGCTGTTACTGTTGAAACAGCAGAAAc                                         >  1:159102/1‑140 (MQ=255)                                                                                           acgTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGAC‑GGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGACGCGGTAGTCAGCCCATCCGCTGTTACTGTTGAAACAGCAGAAACTAAGc                                    >  2:3841/1‑140 (MQ=255)                                                                                              tACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGAC‑GGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGACGCGGTAGTCAGCCCATCCGCTGTTACTGTTGAAACAGCAGAAACTAAGCCGa                                 >  1:357393/1‑140 (MQ=255)                                                                                              tACCGGCGTTGCGCGAAATATTCCAACGCACTGGCTTGTGGGAC‑GGGGTCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGACGCGGTAGTCAGCCCATCCGCTGTTACTGTTGAAACAGCAGAAACTAAGCCGa                                 <  2:170538/140‑1 (MQ=255)                                                                                                      ttGCGCGAAATATTGCAACGCAGAGACATCTTGGAC‑GTGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGACGCGGTAGTCAGCCCATCCGCTGTTACTGTTGAAACAGCAGAAACTAAGCCGATGGATAAg                         <  1:186066/140‑1 (MQ=255)                                                                                                                  ttGCAACGCACAGGCCTGTTGGTC‑GTGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGACGCGGTAGTCAGCCCATCCGCTGTTACTGTTGAAACAGCAGAAACTAAGCCGATGGATAAGCAAACGACGtct             <  1:219942/140‑1 (MQ=255)                                                                                                                  ttGCAACGCACAGGCATGTTGGAC‑GGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGACGCGGTAGTCAGCCCATCCGCTGTTACTGTTGAAACAGCAGAAACTAAGCCGATGGATAAGCAAACGACGtct             <  1:664981/140‑1 (MQ=255)                                                                                                                        cGCACAGGCATGTTGGAC‑GGGGGCCGAAAATTACTCTACCTGGCGACGACACGCCAACTGACGCGGTAGTCAGCCCATCCGCTGTTACTGTTGAAACAGCAGAAACTAAGCCGATGGATAAGCAAACGACGTCTCgtagt       >  1:514879/1‑140 (MQ=255)                                                                                                                             aGGCATGTTGGAC‑GGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGACGCGGTAGTCAGCCCATCCGCTGTTACTGTTGAAACAGCAGAAACTAAGCCGATGGATAAGCAAACGACGTCTCGTAGTAAAcc  >  2:367651/1‑140 (MQ=255)                                                                                                                                          |                                                                                                                                ATGAATCGTTACTGGCCTTACTCTGTGACACCAAACCGTGGCCCCAACTGACCGGCAAAGCAACGTTGCGCCCAGGGCAGTGGAGTAACGACGTACCGGCGTTGCGCGAAATATTGCAACGCACAGGCATGTTGGACGGGGGGCCGAAAATTACTCTACCTGGCGATGACACGCCAACTGACGCGGTAGTCAGCCCATCCGCTGTTACTGTTGAAACAGCAGAAACTAAGCCGATGGATAAGCAAACGACGTCTCGTAGTAAACC  >  NC_000913/981708‑981972

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