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.8%	111.8 / ‑5.7	31	L180P (CTG→CCG)	ldtC	L,D‑transpeptidase LdtC
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (16/14);  minor base T (0/1);  total (16/15)
Fisher's exact test for biased strand distribution p-value = 4.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

TGTACTTCAACCAGACGCGCACCGTTTGGTTCGGCAGAGACTTTTATCGGAGTG‑TTGATGATATTCACTTTGGTGCCTGGGGTGACCTGGCTAAAGAGTGTTTTGATATCGTCATCCCGCAGACGAATACAGCCAGAACTTACCCGCATGCCAATGCCGAAATCGGCGTTCGTACCATGAAGCAAATAAACGCCGCCATAGGCCGCCAGACGAATCGCATGATGGCCCATTGGGTTATCCAGTCCAGCCGGCACTACCGCAGGCAACTCAATTCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCCGGCGGGTAGTAATAGAGACGCAGCTCCGCAATGTTGATCACAATGCCTTCGCGCGGCGCA  >  NC_000913/1169596‑1170081                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                      tGTACTTCAACCAGACGCGCACCGTTTGGTTCGGCAGAGACTTTTATCGGAGTG‑TTGATGATATTCACTTTGGTGCCTGGGGTGACCTGGCTAAAGAGTGTTTTGATATCGTCATCCCGCAGACGAATACAGCCAGAACTTACCCGCATGCCAATGCCGAAATCGGCGTTCGTACCATGAAGCAAATAAACGCCGCCATAGGCCGCCAGACGAATCGCATGATGGCCCATTGGGTTATCCGGTCCAGCCgg                                                                                                                                                                                                                                             >  2:207937/1‑251 (MQ=255)  gTACTTCAACCAGACGCGCACCGTTTGGTTCGGCAGAGACTTTTATCGG‑GTGGTTGATGATATTCGCGTGGGTGCTTGGGGTGACCTGGCTAAAGAGTGTTTTGAGATCTTCATCCCTCAGACGAATACAGCCAGAACTTACCCGCATGCCAATGCCGAAATCGGCGTTCGTACTATGAAGCAAATAAACGCCGCCATAGGCCGCCAGACGAATCGCATGATGGCCCATGGGGTTATCCTGTCGAGCCgg                                                                                                                                                                                                                                             <  2:41734/250‑1 (MQ=255)          aCCAGACGCGCACCGTTTGGTTCGGCAGAGACTTTTATCGGAGTG‑TTGATGATATTCACTTTGGTGCCTGGGGTGACCTGGCTAAAGAGTGTTTTGATATCGTCATCCCGCAGACGAATACAGCCAGAACTTACCCGCATGCCAATGCCGAAATCGGCGTTCGTACCATGAAGCAAATAAACGCCGCCATAGGCCGCCAGACGAATCGCATGATGGCCCATTGGGTTATCCGGTCCAGCCGGCACTACCGc                                                                                                                                                                                                                                    >  1:381876/1‑251 (MQ=255)                                    agagACTTTTATCGGAGTG‑TTGATGACATCCACTCTGGTGCCTGGGGTGACCTGGCTAAAGAGTGTTTTGATATCGTCATCACGCAGACGAATACAGCCAGACCTTACCCGCATGCCAATGCCGAACTCGGCGTTCGTACCATGAAGCAAATAAACGCCGCCATAGGCCGCCAGACGAATCGCATGATGGACAATTGGGTTATCCGGTCCAGCCGGAACTACCGCAGGCAACtaactgacctttgctttat                                                                                                                                                                                                          >  2:211833/1‑233 (MQ=255)                                     gagaCTTTTATCGGAGTG‑TTGATGATATTCCCTTTGGTGCCTGGGGTGACCTGGCTAAAGAGTGTTTTGATATCGTCATCCCGCAGACGAATACAGCCAGAACTTACCCGCATGCCAATGCCGAAATCGGCGTTCGTACCATGAAGCAAATAAACGCCGCCATAGGCCGCCAGACGAATCGCATGATGGCCCATTGGGTTATCCGGTCCAGCCGGCACTACCGCAGGCAACTCAATTCCCTGTGCTTtata                                                                                                                                                                                                         <  1:262161/251‑1 (MQ=255)                                                   gaggggtgatgaTATTCACATTGGTGCCTGGGGTGACCTGGCTAAAGAGTGTTTTGATATCGTCATCCCGCAGACGAATACAGCCAGAACTTAACCGCATGCCAATGCCGAAATCGGCGTTCGAACCATGAAGCAACTAAACGCCCCCATATGCCGCCAGACGAATCACATGCTGGCCCATTGGGTTATACGGTCCAGCCGGCACTACCGCTGGCAACTCACTTCCCTGTGCTTTCTAACGGGCGCGGATg                                                                                                                                                                                            >  2:423325/7‑251 (MQ=255)                                                                               tGGGGTGACCTGGCTAAAGAGTGTTTTGATATCGTCATCCCGCAGACGAATACAGCCAGAACTTACCCGCATGCCAATGCCGAAATCGGCGTTCGTACCATGAAGCAAATAAACGCCGCCATAGGCCGCCAGACGAATCGCATGATGGCCCATTGGGTTATCCGGTCCAGCCGGCACTACCGCAGGCAACTCAATTCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTg                                                                                                                                                                >  1:263371/1‑251 (MQ=255)                                                                                          ggCTAAAGAGTGTTTGTATATCGTCATCCCGCAGACGAATACAGCCAGAACTTACCCGCATGCCAATGCCGAAATCGGCGTTCGTACCATGAAGCAAATAAACGCCGCCATAGGCCGCCAGACGAATCGCATGATGGCCCATTGGGTTATCCGGTCCAGCCGGCACTACCGCAGGCAACTCAATTCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTc                                                                                                                                                      <  1:94735/250‑1 (MQ=255)                                                                                          ggCTAAAGAGTGATTAGATATCGCCATCCCGCAGAGGAATAGAGCCAGAACTTACCCGCCTGCCAATGCCGAAATCGGCGTTCGTACAATGAAGCAAATAAACGCCGCAATAGGCCGACAGACGAATCGAACGATGGCCCCTTGGGGTATCCGGTCAAGCGGGAACTACCGCACGCCACTCAATTCCCTGTGCTTTATAACGGCCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGCTTGTCt                                                                                                                                                     >  2:387780/1‑251 (MQ=255)                                                                                           gCTAAAGAGTGTTTTGATATCGTCATCCCGCAGACGAATACAGCCAGAACTTACCCGCATGCCAATGCCGAAATCGGCGTTCGTACCATGAAGCAAATAAACGCCGCCATAGGCCGCCAGACGAATCGCATGATGGCCCATTGGGTTATCCGGTCCAGCCGGCACTACCGCAGGCAACTCAATTCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTg                                                                                                                                                    <  2:172762/251‑1 (MQ=255)                                                                                                                  cATCCCGCAGACGAATACAGCCAGAACTTACCCGCATGCCAATGCCGAAATCGGCGTTCGTACCATGAAGCAAATAAACGCCGCCATAGGCCGCCAGACAAATCGCATGATGGCCCATTGGGTTATCCGGTCCAGCCGGCACTACCGCAGGCAACTCAATTCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGgt                                                                                                                             <  2:334654/251‑1 (MQ=255)                                                                                                                   aTCCCGCAGACGAATACAGCCAGAACTTACCCGCATGCCAATGCCGAAATCGGCGTTCGTACCATGAAGCAAATAAACGCCGCCATAGGCCGCCAGACGAATCGCATGATGGCCCATTGGGTTATCCGGTCCAGCCGGCACTACCGCAGGCAACTCAATTCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGgtg                                                                                                                            >  1:133038/1‑251 (MQ=255)                                                                                                                                     gCCAGAACTTACCCGCATGCCAATGCCGAAATCGGCGTTCGTACCATGAAGCAAATAAACGCCGCCATAGGCCGCCAGACGAATCGCATGATGGCCCATTGGGTTATCCGGTCCAGCCGGCACTACCGCAGGCAACTCAATTCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTa                                                                                                          <  2:383948/251‑1 (MQ=255)                                                                                                                                                      tGCCAATGCCGAAATCGGCGTTCGTACCATGAAGCAAATAAACGCCGCCATAGGCCGCCAGACGAATCGCATGATGGCCCATTGGGTTATCCGGTCCAGCCGGCACTACCGCAGGCAACTCAATTCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTAtt                                                                                         >  2:143787/1‑251 (MQ=255)                                                                                                                                                                     cGGCGTTCGTACCATGAAGCAAATAAACGCCGCCATAGGCCGCCAGACGAATCGCATGATGGCCCATTGGGTTATCCGGTCCAGCCGGCACTACCGCAGGCAACTCAATTCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTAcc                                                                          <  2:435069/251‑1 (MQ=255)                                                                                                                                                                      ggCGTTCGTACCATGAAGCAAATAAACGCCGCCATAGGCCGCCAGACGAATCGCATGATGGCCCATTGGGTTATCCGGTCCAGCCGGCACTACCGCAGGCAACTCAATTCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCg                                                                         <  1:389350/251‑1 (MQ=255)                                                                                                                                                                             gTACCATGAAGCAAATAAACGCCGCCATAGGCCGCCAGACGAATCGCATGATGGCCCATTGGGTTATCCGGTCCAGCCGGCACTACCGCAGGCAACTCAATTCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTACCTGACCAATACCTATTGGATACCCGGTTACCGAAttttt                                                                  >  2:130581/1‑251 (MQ=255)                                                                                                                                                                              tACCATGAAGCAAATAAACGCCGCCATAGGCCGCCAGACGCATCGCATGATGGCCCATTGGGTTATCCGGTCCAGCCGGCACTACCGCAGGCAACTCAATTCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAAttttt                                                                  >  2:262800/1‑250 (MQ=255)                                                                                                                                                                                    gAAGCAAATAAACGCCGCCATAGGCCGCCAGACGAATCGCATGATGGCCCATTGGGTTATCCGGTCCAGCCGGCACTACCGCAGGCAACTCAATTCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCcggc                                                           >  1:467560/1‑251 (MQ=255)                                                                                                                                                                                          aaTAAACGCCGCCATAGGCCGCCAGACGAATCGCATGATGGCCCATTGGGTTATCCGGTCCAGCCGGCACTACCGCAGGCAACTCAATTCCCTGTGCTGTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATAGCTATTGGATACACGGTTACCGAATTTTGACCcggcggataa                                                     <  2:22280/251‑5 (MQ=255)                                                                                                                                                                                                   cgccATAGGCCGCCAGACGAATCGCATGATGGCCCATTGGGTTATCCGGTCCAGCCGGCACTACCGCAGGCAACTCAATTCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCCGGCGGGTAGTAATAgaga                                            <  2:288284/251‑1 (MQ=255)                                                                                                                                                                                                      cATAGGCCGCCAGACGAATCGCATGATGGCCCATTGGGTTATCCGGTCCAGCCGGCACTACCGCAGGCAACTCAATTCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACTGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTAtt                                                                                         <  2:415890/203‑1 (MQ=255)                                                                                                                                                                                                      cATAGGCCGCCAGACGAATCGCATGATGGCCCATTGGGTTATCCGGTCCAGCCGGCACTACCGCAGGCAACTCAATTCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTAtt                                                                                         >  1:415890/1‑203 (MQ=255)                                                                                                                                                                                                      cATAGGCCGCCAGACGAATCGCATGATGGCCCATTGGGTTATCCGGTCCAGCCGGCACTACCGCAGGCAACTCAATTCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCCGGCGGGTAGTAATAGAGACGc                                         >  1:173965/1‑251 (MQ=255)                                                                                                                                                                                                       aTAGGCCGCCAGACGAATCGCATGATGGCCCATTGGGTTATCCGGTACAGCCGGCACTACCGCAGGCAACTCAATTCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTAGGTGTCCGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGGTACCGAATTTTTACCCGGCGGGTAGTAATAGAGACGcc                                        >  1:201940/1‑250 (MQ=255)                                                                                                                                                                                                                       aTCGCATGATGGCCAATTGGTTTATCCGGTCCATCCGGCACTACCGCAGGCAACTCAATTCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATTCCTATTGGATACACGGTTACCGAATTTTTACCCGGCGGGTAGTAATAGAGACTCAGCTCCGCAATGTTGAt                        <  1:254839/251‑1 (MQ=255)                                                                                                                                                                                                                         cGCATGATGGCCCATTGGGTTATCCGGTCCAGCCGGCACTACCGCAGGCAACTCAATTCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGATACCGAATTTTTACCCGGCGGGTAGTAATAGAGACGCAGCTCCGCAATGTTGATc                       >  2:247227/1‑250 (MQ=255)                                                                                                                                                                                                                              gatgGCCCGTTGGGTTATCCGGTCCAGCCGGCACTACCGCAGGCAACTCAATTCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCCGGCGGGTAGTAATAGAGACGCAGCTCCGCAATGTTGATCACAATg                 <  2:401417/251‑1 (MQ=255)                                                                                                                                                                                                                              gatgGCCCATTGGGTTATCCGGTCCAGCCGGCACTACCGCAGGCAACTCAATTCCCTGTGCTTTATAACGGGCGCGGATGTTTTCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCAtt                                                                                                                                   <  2:262280/137‑1 (MQ=255)                                                                                                                                                                                                                              gatgGCCCATTGGGTTATCCGGTCCAGCCGGCACTACCGCAGGCAACTCAATTCCCTGTGCTTTATAACGGGCGCGGATGTTTTCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCAtt                                                                                                                                   >  1:262280/1‑137 (MQ=255)                                                                                                                                                                                                                                              aTCCGGTCCAGCCGGCACTACCGCAGGCAACTCAATTCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCCGGCGGGTAGTAATAGAGACGCAGCTCCGCAATGTTGATCACAATGCCTTCGCGCGGCGCa  <  2:294526/250‑1 (MQ=255)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                      TGTACTTCAACCAGACGCGCACCGTTTGGTTCGGCAGAGACTTTTATCGGAGTG‑TTGATGATATTCACTTTGGTGCCTGGGGTGACCTGGCTAAAGAGTGTTTTGATATCGTCATCCCGCAGACGAATACAGCCAGAACTTACCCGCATGCCAATGCCGAAATCGGCGTTCGTACCATGAAGCAAATAAACGCCGCCATAGGCCGCCAGACGAATCGCATGATGGCCCATTGGGTTATCCAGTCCAGCCGGCACTACCGCAGGCAACTCAATTCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCCGGCGGGTAGTAATAGAGACGCAGCTCCGCAATGTTGATCACAATGCCTTCGCGCGGCGCA  >  NC_000913/1169596‑1170081

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