Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,630,912	G→A	K60K (AAG→AAA)	slt →	soluble lytic murein transglycosylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,630,912	0	G	A	94.7%	108.0 / ‑3.7	38	K60K (AAG→AAA)	slt	soluble lytic murein transglycosylase
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base A (18/18);  total (19/19)
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

ATAAATCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGTCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAGGATTATCCGCTTTATCCCTACCTGGAATACCGCCAGATCACCGATGATCTGATGAATCAACCGGCGGTGACGGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGG‑AAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCGAAATGGAACAC  >  NC_000913/4630674‑4631157                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                       aaaaatttaatacgcggatagACGCTGCATGGATGTATTTCAACTTAGAGGATGCGCTTGGGGAAAAAGCCAAAAAAGTTCCCTGTCGGGTGGTGGCTGCCGGTGTCTTTCTGCTGTCGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGACGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGGGGGCGAACAAATGATGCCTGGACTGAAAGAGTATCCACtt                                                                                                                                                                                                                                            <  2:156863/231‑1 (MQ=255)    aaTCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTAt                                                                                                                                                                                                                                         >  1:18737/1‑251 (MQ=255)      tCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGATGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGGCAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATcc                                                                                                                                                                                                                                       >  1:362976/1‑251 (MQ=255)             cgcgGCATGACGCTGCATTGATGTATTTACACTTAGAGGATTCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTg                                                                                                                                                                                                                                >  1:152035/1‑251 (MQ=255)               cgGCATGACGCTTCATTGATGTATTTACACGTAGAGGAGGCGCTTGTGGAAAAAGCCAAACAAGTTACCGGGCGGCTGTTGGCTGCCGGTGTCTGTCGTCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGa                                                                                                                                                                                                                              <  1:383609/251‑1 (MQ=255)                    tGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAATCCAAACAAGTTACCTGACGGCTGTTGGCTGCCGGTGTCAGTCTGCTTAAGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAACAGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGGTGTGGTCGAACAACTGATGCCTGGAATGAAAGATTATCCACTTTATCCCTAACTGGAATAcc                                                                                                                                                                                                                         >  2:287179/1‑251 (MQ=255)                                gATGTATTTACACTTAGAGGATGCGCTTGTTGAAAAAGCCAAACAAGTTACCAGGCGGCTGTTGGCTGCCGGTGTCTGGCTGCTGACGGTCAGCAGCTTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCg                                                                                                                                                                                                             <  2:61042/251‑1 (MQ=255)                                   gTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAATTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGAGCGAGCCGACTCACTGGATGAGCAGCGTAGTCGTTACGCGCAAAACAAGCAGGCCTGGGATAATCGACAAATGGCTGCGGTCGCACAAATGATGCCTGGACTGAAGGATTATCCGCTTTACCCCTACCTGGAATACCGCCAGATCACCgatg                                                                                                                                                                                                          >  2:166629/1‑251 (MQ=255)                                   gTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCgat                                                                                                                                                                                                           <  1:303670/250‑1 (MQ=255)                                        taaacTTTGAGTATGCGCTGGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGTTGTCTGTTTGCTGACGGTCAGCAGCGGGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCgatgacct                                                                                                                                                                                                      <  2:338386/247‑4 (MQ=255)                                          aaattAGATGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCtgat                                                                                                                                                                                                   <  1:257291/248‑1 (MQ=255)                                               agagGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATc                                                                                                                                                                                              >  1:76005/1‑251 (MQ=255)                                                       cgcTTGTGGAAAACGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACcggcg                                                                                                                                                                                      <  1:232798/251‑1 (MQ=255)                                                           tggggAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCTACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTCCCTGGAATACCGTCAGAGCACCGATGACTTGATGAATCAACCGGCGGt                                                                                                                                                                                    <  1:41860/246‑1 (MQ=255)                                                              ggAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGActg                                                                                                                                                                               >  1:93071/1‑249 (MQ=255)                                                                       cAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGAAGCGTAGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGATACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGCGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACATGATGAATCAACCGGCGGTGACTGTCTCTAACt                                                                                                                                                                      >  1:182575/1‑251 (MQ=255)                                                                            aaGTTTCCTGGCGGCTGTTGTCGGCCGGTGTCTGGCTGCTGACGGTCAGCAGCGTGGCGCGAGCCTACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACGTTATCCCTACCTGGAATACCGTCAGATCACCGTTGGCCTGATGAATCAACCGGCGGTGATTGTCACGAACTTTGtt                                                                                                                                                                 <  2:192554/251‑1 (MQ=255)                                                                                           tGTTGGCTGCCGGTGACTGTCTGCTGACGGTCAGAAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGGGACTGTCACTAACTTTGTTCGCGCTAAACCCACg                                                                                                                                                  >  1:109147/1‑251 (MQ=255)                                                                                                  tGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTtcctc                                                                                                                                           >  1:291028/1‑251 (MQ=255)                                                                                                          tctgtctgCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGc                                                                                                                                   >  1:33422/1‑251 (MQ=255)                                                                                                                    gACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCa                                                                                                                           <  1:344879/249‑1 (MQ=255)                                                                                                                         tCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCACTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGGCTGTCACTAACTTTGTTCGCGCAAACCCCACGCTTCCTCCCACTCGCACGCTGCAATCTCGt                                                                                                                    >  1:114909/1‑251 (MQ=255)                                                                                                                                        gAGCCGACGCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACtt                                                                                                     <  1:235210/251‑2 (MQ=255)                                                                                                                                                     ggATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAAtctc                                                                                                                      >  1:148927/1‑221 (MQ=255)                                                                                                                                                       aTGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGAc                                                                                      <  1:146205/251‑1 (MQ=255)                                                                                                                                                        tGAGCAGCGTAGCCGTTCCGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCa                                                                                                                                                                                             >  2:171816/1‑147 (MQ=255)                                                                                                                                                        tGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCa                                                                                                                                                                                             <  1:171816/147‑1 (MQ=255)                                                                                                                                                        tGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAAtctc                                                                                                                      >  1:347843/1‑218 (MQ=255)                                                                                                                                                        ggAGCAGCGTAGCCGTGACGCGCAAATCAAGCAGGCCTGGTATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCTCTTTATCCCTACCTGGAATACCGTCAGATCTCCGATGACCTGATGAAGCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAAtctc                                                                                                                      <  2:347843/217‑1 (MQ=255)                                                                                                                                                                gTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGTTAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCAACTTTATCCCTACCTGGAATACCATCAGATCACAGATGACCTGATGAATCAACCGGCGGTGACTGTCCCTAACTTTTTTCGCGCTAAACCCACGCTTCCCCACACCCGCCCTCCGAAATCCCCTTTCGTCACTGAACTGGCGCGGCGTGAAGaacagcgtggg                                                                             >  2:160982/1‑241 (MQ=255)                                                                                                                                                                                 aTCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGccc                                                            <  1:47889/251‑1 (MQ=255)                                                                                                                                                                                               ggATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGgc                                                                                                   <  1:52515/198‑1 (MQ=255)                                                                                                                                                                                                        aCAAATGGATGTGGTCGAACAAATGATGTCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCCTGTTAGCCTTTAGCCCGGGCAAAGCC‑GGAACTACCGa                                      >  1:379634/1‑249 (MQ=255)                                                                                                                                                                                                           aaTGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTTGCCTTTAGCCCGG‑GAAAGCCCGGAAATACCCAAgcgc                                 >  1:238817/1‑251 (MQ=255)                                                                                                                                                                                                                         aacaaatatgtcttGGACTGAAGGATGTTTCGCTCTACCCTTACCTGGAATACCGCCAGTTCCCCGATGTTCTGATGAATCAACCGGTGGCGACGTTCACTAACTTTGTGCGCGCTAACCCCACGCTTCCTCCCGCTCGCAGGCTGCAATCTCGTTTCGTCAATGATCTGGCGCGGCGTGAAGACTGTCGTGTCTTGTTAGCCTTTATCCCGG‑AAAAGCCCGGAACTACCGAAGCGCAATGTAATtac                      <  2:131987/235‑1 (MQ=255)                                                                                                                                                                                                                           aaaaTGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATTCCGTCAGATCACCGAGGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTGTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACTCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGG‑AAAAGCCCGGAACTACCGAAGCGCAATGTAATTActacta                 <  1:160982/250‑1 (MQ=255)                                                                                                                                                                                                                                       ggACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACAGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGG‑AAAAGCCCGGAACTACCGAAGCGCTATGTAATTACTACTATGCAAAATGGa      <  1:339407/250‑1 (MQ=255)                                                                                                                                                                                                                                          cTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGG‑AAAAACCCGGAACTACCGAAGCGCCATGTAATTAATACTATGCAAAATGGaaaac  >  1:237233/1‑248 (MQ=255)                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                       ATAAATCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGTCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAGGATTATCCGCTTTATCCCTACCTGGAATACCGCCAGATCACCGATGATCTGATGAATCAACCGGCGGTGACGGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGCTTGTTAGCCTTTAGCCCGG‑AAAAGCCCGGAACTACCGAAGCGCAATGTAATTACTACTATGCGAAATGGAACAC  >  NC_000913/4630674‑4631157

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: ‑