Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,640,045	+C	coding (13/534 nt)	rrrQ ←	putative lysozyme

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,640,045	1	.	C	96.5%	99.5 / ‑5.1	30	I5V (ATC→GTC)	rrrQ	putative lysozyme
Reads supporting (aligned to +/- strand):  ref base . (0/1);  major base C (14/14);  minor base A (1/0);  total (15/15)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.19e-01

GCCAGCGCCTTATCACGCTCAATGGCGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑T‑TTTTGTGTTCATAGATATTTATGATGAGGACGCTCGTGCTTATTGGCAGGATTTTCAATCTTAAAGGAGTACTGATGCTGCAGATAAGACTCAACTTTTTCTGACAATTTTTCTGCTACTTCCAGGAAGACTTGCCGGACGCTCCTTCTGGCTGCTGCCTCATAAAACTCCAGCGCAGCTCCTTCAACACGGTCCATGGCGACATCCAGGTCAAAAATTTCACCGTCA  >  NC_000913/1639798‑1640273                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                     gCCAGCGCCTTATCACGCTCAATGGCGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAGCGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTTCATATGCCAGAACCATCGCGGTATGCCATTGTGGGGTTACCTTCTTGTTCGTCCAGAAACTTGTCAAGTATCTTAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGTCCGTATCTGA‑tctt                                                                                                                                                                                                                                    <  1:320214/251‑4 (MQ=255)                 gCTCAATGGCGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgccc                                                                                                                                                                                                                    <  1:83313‑M1/251‑8 (MQ=255)                 gCTCAATGGCGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTACATATGCCAGAACCATCGCTGTATGCCATTGTGTGGTTACCTTCTTTTGCGTCCAGAAAATGGTCCAGTATCTGAGGAGTAGATGCGCCAGCACCACTCAGCGCCAGAACGGCAGCCGACAGGCCGTCACTGA‑TCTTTTGTGTTCAcccgccc                                                                                                                                                                                                                    >  2:442023‑M1/1‑244 (MQ=255)                        ggCGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTAATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgt                                                                                                                                                                                                             >  1:161104‑M1/1‑237 (MQ=255)                          cGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtct                                                                                                                                                                                                           >  1:123205‑M1/1‑235 (MQ=255)                              gACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttac                                                                                                                                                                                                        <  1:144433‑M1/250‑20 (MQ=255)                                  tGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGCCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatat                                                                                                                                                                                                   >  1:186920‑M1/1‑227 (MQ=255)                                   ggTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatatt                                                                                                                                                                                                  <  1:272683‑M1/251‑26 (MQ=255)                                              tcctccGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTTGGTGGGTCCCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCGcccgccctcttcgtcttatgatattcatcataaac                                                                                                                                                                                        <  2:312272‑M1/250‑36 (MQ=255)                                               ccTTCGACAGTTTCATATTGGGAAAAACGGTTTTCCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCAGCGCGGTATGCCATTGTGTGGTGACCTTCTTTTTCGTCCAGAAACTGGTCAAGGATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaact                                                                                                                                                                                       <  2:113492‑M1/250‑37 (MQ=255)                                                  tCGACAGTTGCATATTGGGAAAAACGGTTTTTACATCCACCACTGTGGCACCCCGACATATTGTCCATATGCCAGAACCATCGCTGTATGCCATTGGGTGGTTACCTTCTTTTCCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatagtcatcataaacttt                                                                                                                                                                                     <  2:186920‑M1/249‑39 (MQ=255)                                                   cGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttggg                                                                                                                                                                                  >  2:314786‑M1/1‑210 (MQ=255)                                                            cATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGTTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatact                                                                                                                                                                          <  2:372395‑M1/250‑50 (MQ=255)                                                                                     tccaCCACTGTGCCACCCCGACAGATGGTGCCTATGCCAGACCCCTCGCGGTATGCCATTGCGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGAGCCAGAACGGCCGCCGACAGGCCGCATCTGA‑TCTTTTGTGTTCAcccgcccgcctcgtcttacgatattcatcatcaactttgggatcatactgaagctccccgccagatgcctcctgt                                                                                                                                                >  2:166788‑M1/1‑176 (MQ=255)                                                                                                   cccccgACAGATGGCCCATATGCCAGAACCATCGCGGTCTGCCATTGGGTGGTTACCTTCTATTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCACAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcg                                                                                                                                  <  2:42889‑M1/250‑90 (MQ=255)                                                                                                       cGACAGATGTTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgacct                                                                                                                              <  2:161104‑M1/251‑94 (MQ=255)                                                                                                         aCAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttc                                                                                                                            >  1:9960‑M1/1‑156 (MQ=255)                                                                                                                      atGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGCCAGGCCGTATCCGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcc                                                                                                               >  2:347739‑M1/1‑143 (MQ=255)                                                                                                                          cAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGCCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCATTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacccatcgcgcgaccttcgggaactaaatcccctt                                                                                                           >  2:455207‑M1/1‑139 (MQ=255)                                                                                                                           aGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatccccttt                                                                                                          >  1:375895‑M1/1‑138 (MQ=255)                                                                                                                            gAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttc                                                                                                         >  1:214191‑M1/1‑137 (MQ=255)                                                                                                                             aaCCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttcc                                                                                                        >  2:216493‑M1/1‑136 (MQ=255)                                                                                                                                            ccATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaag                                                                                         <  1:301051‑M1/251‑131 (MQ=255)                                                                                                                                                                                            gaggagCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgcaagttttgtttttgaaccgaaat                                         >  1:432369‑M1/1‑73 (MQ=255)                                                                                                                                                                                                  cagagagcccagcaccaaTCAGCGCATGAAAGGCAGCCTACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatagtcatcataaactttgggatcatactgaagctccccgccaggtgccgcctggagacgcatcgcgcgaccgtcgggaactaaatcccctttccagcgataaagcgaagccaaactaaggcctgctgcttgtgcaagttttgtttttgaaccgaaatacaaaa                                   <  2:287180‑M1/243‑185 (MQ=255)                                                                                                                                                                                                              cACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgcaagttttgtttttgaaccgaaatacaaaagagcgtcagttt                       <  1:216493‑M1/251‑197 (MQ=255)                                                                                                                                                                                                               accaaccaGCGCCAGATCGGCAGCCGACATGCCGTATCTGAGC‑TGTTGGGTTCAcccgcccgcttcgtcttacgatattcatcataaacgttgggatcatactgaagctccccgccagatgcctcctgtagacgcatc                                                                                                                                       <  1:461737‑M1/138‑85 (MQ=255)                                                                                                                                                                                                               aCCAATCAGCGCCAGAACGGCAGCAGACGGGCCGTATCTGA‑TATTTTGTGCTCAcccgcccgcttagtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatc                                                                                                                                       >  2:461737‑M1/1‑54 (MQ=255)                                                                                                                                                                                                                 cAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgcaagttttgtttttgaaccgaaatacaaaagagcgtcagtttta                     >  2:216471‑M1/1‑52 (MQ=255)                                                                                                                                                                                                                      aGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑TCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctagaaagcgaagccaaacgaatacctgctgcttgtgcaagttgtgtttttgaaccgaaatacaaaagagcgtcagttttaagcatt               <  1:166788‑M1/251‑205 (MQ=255)                                                                                                                                                                                                                                   cgtccgacAGGCGGTATCTGA‑TCTTTGGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgcaagttttgtttttgaaccgaaatacaaaagagcgtcagttttaagcatttaaaacaccttta  <  2:9960‑M1/248‑218 (MQ=255)                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                     GCCAGCGCCTTATCACGCTCAATGGCGTTGACCTGGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGA‑T‑TTTTGTGTTCATAGATATTTATGATGAGGACGCTCGTGCTTATTGGCAGGATTTTCAATCTTAAAGGAGTACTGATGCTGCAGATAAGACTCAACTTTTTCTGACAATTTTTCTGCTACTTCCAGGAAGACTTGCCGGACGCTCCTTCTGGCTGCTGCCTCATAAAACTCCAGCGCAGCTCCTTCAACACGGTCCATGGCGACATCCAGGTCAAAAATTTCACCGTCA  >  NC_000913/1639798‑1640273

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