Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,631,323	A→T	T197T (ACA→ACT)	slt →	soluble lytic murein transglycosylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,631,323	0	A	T	90.9%	193.2 / 8.5	66	T197T (ACA→ACT)	slt	soluble lytic murein transglycosylase
Reads supporting (aligned to +/- strand):  ref base A (2/4);  new base T (34/26);  total (36/30)
Fisher's exact test for biased strand distribution p-value = 3.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.75e-01

GCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCGAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGCGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACAGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTA‑CCAGACTATCGCCTCGGCAATCATTTCACTGGCGAACAACCCTAATAC‑GGTACTGACCTTCGCGCGTACAACTGGCGCGACCGATTTTACCCGTCAAATGGCGGCGGTGGCGTTTGCCAGTGTGGCGCGGC‑AGGATGCTGAGAATGCACGGCTGATGATCCCATCGCTTGCCCAGGCGCAGCAGCTTAATGAAGATCAGATTC  >  NC_000913/4631077‑4631567                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                        gCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACGATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGcc                                                                                                                                                                                                                                                     <  2:155982/251‑1 (MQ=255)   gTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTg                                                                                                                                                                                                                                                   <  1:360772/251‑1 (MQ=255)     ggCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAATCGCAATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTGGt                                                                                                                                                                                                                                                 <  1:109584/251‑1 (MQ=255)       cTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGCGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTGGTa                                                                                                                                                                                                                                                >  1:310795/1‑250 (MQ=255)          gTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTGGTAACa                                                                                                                                                                                                                                             <  2:235764/250‑1 (MQ=255)           ttAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTGGTAACAGt                                                                                                                                                                                                                                           >  1:361080/1‑251 (MQ=255)                  ttAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTGGTAACAGTGCTggc                                                                                                                                                                                                                                     <  2:361080/250‑1 (MQ=255)                    aGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTGGTAACAGTGCTggcagg                                                                                                                                                                                                                                  >  1:154005/1‑251 (MQ=255)                        cGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTGGTAACAGTGCTggcagggc                                                                                                                                                                                                                                >  1:404282/1‑249 (MQ=255)                          gAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTGGTAACAGTGCTGGCAGGGCAGAt                                                                                                                                                                                                                            <  1:177446/251‑1 (MQ=255)                                ccctgAACTACCGAAGCGCAATGTAATTACTACTATGCGAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGCGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCGGGCGATGAAAGCGGGTAACACAGGCCTGGTAACAGTGCTGGCAGGGCAGAtgcctgc                                                                                                                                                                                                                      <  2:34337/247‑1 (MQ=255)                                       cTACCGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTA‑c                                                                                                                                                                                                              >  1:319224/1‑251 (MQ=255)                                          ccGAAGCGCAATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTGGTAACAGTGCTGGCAGGGCAGAGGCCTGCCGATTA‑ca                                                                                                                                                                                                             >  2:458308/1‑248 (MQ=255)                                               gcgcAATGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTA‑CCAGACTAt                                                                                                                                                                                                      <  1:40670/251‑1 (MQ=255)                                                    aTGTAATTACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGCTTA‑CCAGACTATCGCCt                                                                                                                                                                                                 >  1:31547/1‑251 (MQ=255)                                                          ttACTACTATGCAAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTA‑CCAGACTATCGCCTCGGCa                                                                                                                                                                                            >  1:50377/1‑250 (MQ=255)                                                            actactAGGCAAAATGGAACCCCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTA‑CCAGACTATCGCCTCGGCAAt                                                                                                                                                                                          <  2:58784/250‑1 (MQ=255)                                                                     cGAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGCGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACAGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTA‑CCAGACTATCGCCTCGGCAATCATTTCACTg                                                                                                                                                                                >  1:62487/1‑251 (MQ=255)                                                                         aTGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTA‑CCCGACTATCTCCTCGGCAATCATTTCACTGGCg                                                                                                                                                                             >  1:144105/1‑250 (MQ=255)                                                                                   gggCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTA‑CCAGACTATCGCCTCGGCAATCATTTCACTGGCGAACAACCCTa                                                                                                                                                                   >  1:388873/1‑250 (MQ=255)                                                                                      cAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTA‑CCAGACTATCGCCTCGGCAATCATTTCACTGGCGAACAACCCTAATAc                                                                                                                                                               >  2:444798/1‑251 (MQ=255)                                                                                         agTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTA‑CCAGACTATCGCCTCGGCAATCATTTCACTGGCGAACCACCCTAATAC‑ggg                                                                                                                                                           >  1:349708/1‑250 (MQ=255)                                                                                                         aaGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGGGCCAAGTTATTTAGTGTCTGGCGTGCGTCCGGTAACCAAGATCAGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCATGGGAACAGTGCTGGCAGGGCAGACGCCTTCCGATTAACCA‑ACTATCGCATCGGCAAACATTCCCCTGGCGAACAACCCTAATaaggacctgccctcagccggt                                                                                                                                            >  2:503309/1‑231 (MQ=255)                                                                                                          aGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTA‑CCAGACTATCGCCTCGGCAATCATTTCACTGGCGAACAACCCTAATAC‑GGTACTGACCTTCGCGCGt                                                                                                                                           >  1:321845/1‑250 (MQ=255)                                                                                                             ggCGAACGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTA‑CCAGACTATCGCCTCGGCAATCATTTCACTGGACAACAACCCTAATAC‑GGTACTGACCTTCGCCCGTACa                                                                                                                                        >  2:174226/1‑250 (MQ=255)                                                                                                               cGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTA‑CCAGACTATCGCCTCGGCAATCATTTCACTGGCGAACAACCCTAATAC‑GGTACTGACCTTCGCGCGTACAACt                                                                                                                                     <  1:104878/251‑1 (MQ=255)                                                                                                                gAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGCGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACAGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTA‑CCAGACTATCGCCTCGGCAATCATTTCACTGGCGAACAACCCTAATAC‑GGTACTGACCTTCGCGCGTACAACTg                                                                                                                                    <  2:62487/251‑1 (MQ=255)                                                                                                                 aaaGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTA‑CCAGACTATCGCCTCGGCAATCATTTCACTGGCGAACAACCCTAATAC‑GGTACTGACCTTCGCGCGTACAACTgg                                                                                                                                   >  1:267145/1‑251 (MQ=255)                                                                                                                          tgGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACCAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTGGTAACAGTGCTGGCAGGGCAGAGGCCTGCCGATTA‑CCAGACTATCGCCTCGGCAATCATTTCACTGGCGAACAACCCTAATAC‑GGTACTGACCTTCCCGCGTACAACAGGCGCGACCgg                                                                                                                          >  2:25789/1‑250 (MQ=255)                                                                                                                            gCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTA‑CCAGACTATCGCCTCGGCAATCATTTCACTGGCGAACAACCCTAATAC‑GGTACTGACCTTCGCGCGTACAACTGGCGCGACCGAtt                                                                                                                        >  1:15237/1‑251 (MQ=255)                                                                                                                               aaCCGGCAAGAGCCCTCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGTTAAACAAGATCCGCTGGCTTATTTAGATCGTATCCGTCTGTCGATGAAAGCGGGTAACACTGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTA‑CCAGACTATCGCCTCGGCAATCATTTCACTGTCGAACAACCCTAATAC‑GGTACTGACCTTCGCGCGTACAACTGGCTCGACCGATTTTa                                                                                                                     <  1:91667/251‑1 (MQ=255)                                                                                                                                  cGGCAAGGGACAGCCTAACGCCTGTGACAAGTTATTTAGTGTCGGGCGTGCGGCAGGGAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACCCTGGCCTTGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTA‑CCAGACTATCGCCTCGGCAATCATTTCACTGGCGAACAACCCTAATAC‑GGTACTGACCTTCGCGCGTACAACTGGCGCGACCGATTTTAccc                                                                                                                  <  2:322880/251‑1 (MQ=255)                                                                                                                                     cAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTA‑CCAGACTATCGCCTCGGCAATCATTTCACTGGCGAACAACCCTAATAC‑GGTACTGACCTTCGCGCGTACAACTGGCGCGACCGATTTTACCCGTc                                                                                                               <  1:458308/251‑1 (MQ=255)                                                                                                                                               ccTAACGCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTA‑CCAGACTATCGCCTCGGCAATCATTTCACTGGCGAACAACCCTAATAC‑GGTACTGACCTTCGCGCGTACAACTGGCGCGACCGATTTTACCCGTCAAATggcggg                                                                                                     >  1:22234/1‑250 (MQ=255)                                                                                                                                                  aaCGCCGGTGACAAGTTATTTAGTGTCTGGCGTGCGTCATGTAAACAAGATCCGCTGGCGTATGTAGAGCGTATCCGTCTGGCGAGGAAAGCGGGTAACACTGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTA‑CCAGACTATCGCCTCGGCAAGCATTTCACTGGCGAACAACCCTAATAC‑GGTACTGACCTTCGCGCGTACAACTGGCGCGACCGATTTTACCCGTCAAATGgcggcg                                                                                                    <  2:267145/249‑1 (MQ=255)                                                                                                                                                     gCCTGTGACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACACGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTA‑CCAGACTATCGCCTCGGCAATCATTTCACTGGCGAACAACCCTAATAC‑GGTACTGACCTTCGCGCGTACAACTGGCGCGACCGATTTTACCCGTCAAATGGcggcgggggg                                                                                               >  2:495382/1‑247 (MQ=255)                                                                                                                                                        tgtgACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTA‑CCAGACTATCGCCTCGGCAATCATTTCACTGGCGAACAACCCTAATAC‑GGTACTGACCTTCGCGCGTACAACTGGCGCGACCGATTTTACCCGTCAAATGGCGGCGGCGGCGtt                                                                                            >  1:334620/1‑251 (MQ=255)                                                                                                                                                           gACAAGTTATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTGGTAACAGTGCTGGCAGGGCCGATGCCTGCCGATTA‑ACAGACTATCGCCTCGGCAATCATTTCACTGGCGAACAACCCTAATAA‑GGGACTGACCTTCGCGCGTACAACTGGCGAGACCGATTTTACCCGTCAAAAGGCGGGGGTGGCGttttc                                                                                         >  2:386628/1‑249 (MQ=255)                                                                                                                                                                 ttATTTAGTGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTA‑CCAGACTATCGCCTCGGCAATCATTTCACTGGCGAACAACCCTAATAC‑GGTACTGACCTTCGCGCGTACAACTGGCGCGACCGATTTTACCCGTCAAATGGCGGCGGTgggggtttgcccggg                                                                                   >  2:500479/1‑238 (MQ=255)                                                                                                                                                                    tttAGCGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACAGGCCTGGTAACAGTGCTGGCAGGGCAGAtgcc                                                                                                                                                                                                                         <  1:169280/116‑1 (MQ=255)                                                                                                                                                                    tttAGCGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACAGGCCTGGTAACAGTGCTGGCAGGGCAGAtgcc                                                                                                                                                                                                                         >  2:169280/1‑116 (MQ=255)                                                                                                                                                                         tgTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTA‑CCAGACTATCGCCTCGGCAATCATTTCACTGGCGAACAACCCTAATAC‑GGTACTGACCTTCGCGCGTACAACTGGCGCGACCGATTTTACCCGTCAAATGGCGGCGGTGGCGTTTGCCAGTGTGGCGCggg                                                                           >  1:220396/2‑250 (MQ=255)                                                                                                                                                                          gTCTGGCGTGCGTCAGGTAAACAAGTTCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTA‑CCAGACTATCGCCTCGGCAATCATTTCACTGGCGAACAACCCTAATAC‑GGTACTGACCTGCGCGCGTACAACTGGCGCGACCGATTTTACCCTTCAAATGGCGGCGGTGGCGTTTGCCAGTGTGGCGCGGC‑a                                                                         <  1:200902/251‑1 (MQ=255)                                                                                                                                                                          gTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTA‑CCAGACTATCGCCTCGGCAATCATTTCACTGGCGAACAACCCTAATAC‑GGTACTGACCTTCGCGCGTACAACTGGCGCGACCGATTTTACCCGTCAAATGGCGGCGGTGGCGTTTGCCAGTGTGGCGCGGc                                                                           <  1:362586/250‑1 (MQ=255)                                                                                                                                                                          gTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTA‑CCAGACTATCGCCTCGGCAATCATTTCACTGGCGAACAACCCTAATAC‑GGTACTGACCTTCGCGCGTACAACTGGCGCGACCGATTTTACCCGTCAAATGGCGGCGGTGGCGTTTGCCAGTGTGGCGCGGC‑a                                                                         <  1:495382/251‑1 (MQ=255)                                                                                                                                                                                       cAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGCTTA‑CCAGACTATCGCCTCGGCAATCATTTCACTGGCGAACAACCCTAATAC‑GGTACTGACCTTCGCGCGTACAACTGGCGGGACCGATTTTACCCGTCCAATGGCGGCGGTGGCGTTTGACAGGGTGGCCCGGC‑AGGGTGCTgagac                                                             >  1:305712/1‑249 (MQ=255)                                                                                                                                                                                       cAGGTAAACAAGATCCGCTGGCGGATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTA‑CCAGACTATCGCCTCGGCAATCATTTCACTGGCGAACAACCCTAATAC‑GGTACTGACCTTCGCGCGTACAACTGGCGCGACCGATTTTACCCGTCAAATGGCGGCGGTGGCGTTTGCCAGTGTGGCGCGGC‑AGGATGCTGAGAAt                                                            <  2:63547/251‑1 (MQ=255)                                                                                                                                                                                       cAGGAAAACAAGATCAGCTGGAGAATTTAGAGCGTATCCGTCTGGCGATGACAGCGGGTAACACTGGCCTGGTAACAGTGCTGGCAGGGCCGATGCCTGCCGATTA‑CCAGACTATCGACTCGGCAATCATTTCACTGGCGAACAACCCTCATAC‑GGTAATGTCCTTCGAGCGGACAAATGGCACGACCCATTTTACACGTCACATGGCGGCGGGGCCGGTTGACAGTGCGGcgcgcacgagtcatggggat                                                             >  2:476828/1‑235 (MQ=255)                                                                                                                                                                                                   aTCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTA‑CCAGACTATCGCCTCGGCAATCATTTCACTGGCGAACAACCCTAATAC‑GGTACTGACCTTCGCGCGTACAACTGGCGCGACCGATTTTACCCGTCAAATGGCGGCGGTGGCGTTTGCCAGGGTGGCGCGGCCAGG‑TGCTGAGAAAGCACGGCtggtg                                                >  1:219963/1‑248 (MQ=255)                                                                                                                                                                                                      cGCTGGCGTATTGAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTA‑CCAGACTATCGCCTCGGCAATCATTTCACTGGCGAACAACCCTAATAC‑GGTACTGACCTTCGCGCGTACAACTGGCGCGACCGATTTTACCCGTCAAATGGCGGCGGTGGCGTTTGCCAGTGTGGCGCGGC‑AGGATGCTGAGAATGCACGGCTGATGATc                                             <  1:455979/251‑1 (MQ=255)                                                                                                                                                                                                        cTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTA‑CCAGACTATCGCCTCGGCAATCATTTCACTGGCGAACAACCCTAATAC‑GGTACTGACCTTCGCGCGTACAACTGGCGCGACCGATTTTACCCGTCAAATGGCGGCGGTGGCGTTTGCCAGTGTGGCGCGGC‑CGGGTGCTGAGAATGCACGGCTGGTGGTccc                                           >  1:146938/1‑251 (MQ=255)                                                                                                                                                                                                        cTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTA‑CCAGACTATCGCCTCGGCAATCATTTCACTGGCGAACAACCCTAATAC‑GGTACTGACCTTCGCGCGTACAACTGGCGCGACCGATTTTACCCGTCAAATGGCGGCGGTGGCGTTTGCCAGTGTGGCGCGGC‑AGGATGCTGAGAATGCCCGGCTGATGATccc                                           >  1:370878/1‑251 (MQ=255)                                                                                                                                                                                                           gCGTATTTAGAACGTATCCCTCTGGTGAGTAAAGCGGGTAACCCTGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTA‑CCAGACTATCGCCTCGGCAAGCATTTCACTGGCGAACAACCCTAATAC‑GGTACTGACCTTCGCGCGTACAACTGGCGCGACCGATTTTACCCGTCAAATGGCGGCGGTGGCGTTTGCCAGTGTGGCGCGGC‑AGGATGCTGATAATGCACGGCTGATGACCCCATc                                        <  2:326873/251‑1 (MQ=255)                                                                                                                                                                                                            cgtatttggcttgtgaccttCTGGCGATGAAAGCGGGTAACACTGTCCTGGTAACAGTGCGGGCAGGGCAGACGCCTGCCGATTA‑CCAGACTGTCGCCGCTGCAATCATTTCACTGGCGAACAACCCTAATACAGGT‑CTGACCTTCGCGCGTTCACCTGGCGCGACCGATTTTACCCGTCAAATGGCGGCGGTGGCGTTTGCCAGTGTGGCGCGGA‑AGGATGCTGAGCACGCGCGGCTGATGATCCCATCg                                       <  2:73368/232‑1 (MQ=255)                                                                                                                                                                                                               aTTTAGAGCGTATCCGTCGGGCGAGGAAAGCGGGTAACACAGGCCTGGTAACAGGGCTGGCAGGGCAGATGCCTGCCGATTA‑CCAGACTATCGCCTCGGCAATCATTTCACTGGCGAACAACCCTAATAC‑GGTACTGACCTTCGCGCGTACAACTGGCGCGACCGAGTTTACCCGTCAAATGGCGGCGGTGGCGTTTGCCAGTGTGGCGCGGC‑AGGATGCTGAGAATGCACGGCTGATGATCCCATCGCtt                                    <  2:154005/251‑1 (MQ=255)                                                                                                                                                                                                                      gCGTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTA‑CCAGACTATCGCCTCGGCAATCATTTCACTGGCGAACAACCCTAATAC‑GGTACTGACCTTCGCGCGTACAACTGGCGCGACCGATTTTACCCGTCAAATGGCGGCGGTGGCGTTTGCCAGTGGGGCGCGGC‑AGGGTGCTGGGAATGCCCGGCTGATGATCCCCTCGCCTGCCCacg                             >  1:298463/1‑249 (MQ=255)                                                                                                                                                                                                                       cGTATCCGTCTGGCGATGAAAGCGGGTATCACTGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTA‑CCAGACTATCGCCTCGGCAATCATTTCACTGGCGAACAACCCTAATAC‑GGTACTGACCTTCGCGCGTACAACTGGCGCGACCGATTTTACCCGTCAAATGGCGGCGGTGGCGTTTGCCAGTGTGGCGCGGC‑AGGATGCTGAGAATGCACGGCTGATGATCCCATCGCTTGCCCAGgc                            <  1:155165/251‑1 (MQ=255)                                                                                                                                                                                                                        gTATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTA‑CCAGACTATCGCCTCGGCAATCATTTCACTGGCGAACAACCCTAATAC‑GGTACTGACCTTCGCGCGTACAACTGGCGCGACCGATTTTACCCGTCAAATGGCGGCGGTGGCGTTTGCCAGTGTGGCGCGGC‑AGGATGCTGAGAATGCACGGCTGATGATCCCATCGCTTGCCCAGgcg                           <  2:50377/251‑1 (MQ=255)                                                                                                                                                                                                                         tATCCGTCTGGCGATGAAAGCGGGTAACACTGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTA‑CCAGACTATCGCCTCGGCAATCATTTCACTGGCGAACAACCCTAATAC‑GGTACTGACCTTCGCGCGTACAACTGGCGCGACCGATTTTACCCGTCAAATGGCGGCGGTGGCGTTTGCCAGTGTGGCGCGGC‑AGGATGCTGAGAATGCACGGCCGATGATCCCATCGCTTGCCCAGgcgc                          >  1:125255/1‑251 (MQ=255)                                                                                                                                                                                                                              gTCTGGCGATGAAGGCGGGTAACACTGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTA‑CCAGACTATCGCCTCGGCAATCATTTCACTGGCGAACAACCCTAATAC‑GGTACTGACCTTCGCGCGTACAACTGGCGCGACCGATTTTACCCGTCAAATGGCGGCGGTGGCGTTTGCCAGTGTGGCGCGGC‑AGGATGCTGAGAATGCACGGCTGATGATCCCATCGCTTGCCCAGGCGcagcag                     <  2:388873/251‑1 (MQ=255)                                                                                                                                                                                                                                 tGGCGATGAAAGCGGGTAACACTGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTA‑CCAGACTATCGCCTCGGCAATCATTTCACTGGCGAACAACCCTAATAC‑GGTACTGACCTTCGCGCGTACAACTGGCGCGACCGATTTTACCCGTCAAATGGCGGCGGTGGCGTTTGCCAGTGTGGCGCGGC‑AGGATGCTGAGCATGCACGGGTGATGATCCCCTCGCTTGCCCCGGCGCAGCAGCtt                  >  1:361361/1‑251 (MQ=255)                                                                                                                                                                                                                                  ggCGATGAAAGCGGGTAACACTGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTA‑CCAGACTATCGCCTCGGCAATCATTTCACTGGCGAACAACCCTAATAC‑GGTACTGACCTTCGCGCGTACAACTGGCGCGACCGATTTTACCCGTCAAATGGCGGCGGTGGCGTTTGCCAGTGTGGCGCGGC‑AGGATGCTGAGAATGCACGGCTGATGATCCCATCGCTTGCCCAGGCGCAGCAGCt                   <  1:457390/249‑1 (MQ=255)                                                                                                                                                                                                                                        gAAAGCGGGTAACACTGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTA‑CCAGACTATCGCCTCGGCAATCATTTCACTGGCGAACAACCCTAATAC‑GGTACTGACCTTCGCGCGTACAACTGGCGCGACCGATTTTACCCGTCAAATGGCGGCGGTGGCGTTTGCCAGTGGGGCGCGGC‑AGGATGCTGGGAATGGCCGGGTGGTGATCCCCTCGCTTTCCCCAGCGCAGCAGCTTCATGa             >  2:23027/1‑249 (MQ=255)                                                                                                                                                                                                                                          aaGCGGGTAACACTGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTA‑CCAGACTATCGCCTCGGCAATCATTTCACTGGCGAACAACCCTAATAC‑GGTACTGACCTTCGCGCGTACAACTGGCGCGACTGATTTTACCCGTCAAATGGCGGCGGTGGCGTTTGCCAGTGTGGCGCGGC‑AGGATGCTGAGAATGCACGGCTGATGATCCCATCGCTTGCCCAGGCGCAGCAGCTTAATGAAGAt         <  1:346872/251‑1 (MQ=255)                                                                                                                                                                                                                                                  aaCACTGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTA‑CCAGACTATCGCCTCGGCAATCATTTCACTGGCGAACAACCCTAATAC‑GGTACTGACCTTCGCGCGTACAACTGGCGCGACCGATTTTACCCGTCAAATGGCGGCGGTGGCGTTTGCCAGTGTGGCGCGGC‑AGGATGCTGAGAATGCACGGCTGATGATCCCATCGCTTGCCCAGGCGCAGCAGCTTAATGAAGATCAGATTc  <  1:85779/250‑1 (MQ=255)                                                                                                                                                                                                                                                   acacTGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTA‑CCAGACTATCGCCTCGGCAATCATTTCACTGGCGAACAACCCTAATAA‑GGTACTGACCTTCGCGCGTACAACTGGCGCGACCGATTTTACCCGTCAAATGGCGGCGGAGGAGTTTGCCAGAGTGGCGCGGC‑GGGATGGTGAGGATGCAAGGGTCATGATCCCCTCGCGTGCCCACGCGCAGCAGCTTAATGAAGATcctcttc  >  1:434932/1‑243 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                        GCGTGGCTTGTTAGCCTTTAGCCCGGAAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCGAAATGGAACACCGGGCAGAGTGAAGAAGCCTGGCAAGGGGCGAAAGAGCTGTGGCTAACCGGCAAGAGCCAGCCTAACGCCTGTGACAAGTTATTTAGCGTCTGGCGTGCGTCAGGTAAACAAGATCCGCTGGCGTATTTAGAGCGTATCCGTCTGGCGATGAAAGCGGGTAACACAGGCCTGGTAACAGTGCTGGCAGGGCAGATGCCTGCCGATTA‑CCAGACTATCGCCTCGGCAATCATTTCACTGGCGAACAACCCTAATAC‑GGTACTGACCTTCGCGCGTACAACTGGCGCGACCGATTTTACCCGTCAAATGGCGGCGGTGGCGTTTGCCAGTGTGGCGCGGC‑AGGATGCTGAGAATGCACGGCTGATGATCCCATCGCTTGCCCAGGCGCAGCAGCTTAATGAAGATCAGATTC  >  NC_000913/4631077‑4631567

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