Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	CP006584	4,178,988	Δ1 bp	noncoding (1018/1541 nt)	LY180_19945 →	16S ribosomal RNA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP006584	4,178,988	0	A	.	97.4%	57.1 / ‑2.1	39	noncoding (1018/1541 nt)	LY180_19945	16S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base . (16/21);  total (17/22)
Fisher's exact test for biased strand distribution p-value = 4.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

TGACGCTCAGGTGCGAAAGCGTGGGGAGCAAACAGGATTAGATACCCTGGTAGTCCACGCCGTAAACGATGTCGACTTGGAGGTTGTGCCCTTGAGGCGTGGCTTCCGGAGCTAACGCGTTAAGTCGACCGCCTGGGGAGTACGGCCGCAAGGTTAAAACTCAAATGAATTGACGGGGGCCCGCACAAGCGGTGGAGCATGTGGTTTAATTCGATGCAACGCGAAGAACCTTACCTGGTCTTGACATCCACAGAACTTTTCAGAGATGAGAAT‑GTG‑CCTTCGGGAACCGT‑GAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGTGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCCTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACGACCAGGGCTACACACGTGCTACAATGGCGCATACAAAGAGAAGCGACCTCGCGAGAGCAAGCGGACCTCATAA  >  CP006584/4178720‑4179257                                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                                                 tGACGCTCAGGTGCGAAAGCGTGGGGAGACAACAGGATTAGATACCCTGGTAGTCCACGCCGTAAACGATGTCGACTTGGAGGTTGTGCCCTTGAGGCGTGGCTTCCGGAGCTAACGCGTTAAGTCGACCGCCTGGGGAGTACGGCCGCAAGGTTAAAACTCAAATGAATTGACGGGGGCCCGCACAAGCGGTGGAGCATGTGGTTTAATTCGATGCAACGCGAAGAACCTTACCTGGTCTTGACATCCACAGAACTTTCCAGAGATG‑GATTGGTG‑CCTTCGGGAACTGT‑GAGACAgg                                                                                                                                                                                                                                                  <  1:382299/298‑1 (MQ=255)         aGGTCCGAAAGCGTGGGGAGCAAACAGGATTAGATACCCTGGTAGTCCACGCCGTAAACGATGTCGACTTGGAGGTTGTGCCCTTGAGGCGCGGCTTCCGGAGCTAACGCGTTAAGTCGACCGCCTGGGGAGCACGGCCGCAAGGTTAAAACTCAACTGACTTGACGGGGGCCCGACCAAGCGGTGGCGCAAGTCGTTTAATCAGATGCAACCACAAGAACATTACCTGGTCTGGACATCCACCGTACTTTTCAGAGATGAAGAT‑GTG‑CCTTCGGGAACCGT‑GGAACCGGTGCCGCac                                                                                                                                                                                                                                          >  2:119326/1‑297 (MQ=11)                     gTGGGGAGCAAACAGGATTAGATACCCTGGTAGTCCACGCCGTAAACGATGTCGACTTGGAGGTTGTGCCCTTGAGGCGTGGCTTCCGGAGCTAACGCGTTAAGTCGACCGCCTGGGGAGTACGGCCGCAAGGTTAAACATCAAATGAATTGACGGGGGCCCGCACAAGCGGTGGAGCATGTGGTTTAATTCGATGCAACGCGAACAACCTTACCTGGTCTTGACATCCACAGAACTTTCCCCAGATG‑GATTGGTG‑CCTTCGGGCCCTGT‑GAGACAGGTCCGGAGTGGCGgtcgtccg                                                                                                                                                                                                                              >  2:261439/1‑296 (MQ=11)                     gTGGGGAGCAAACAGGATTAGATACCCTGGTAGTCCACGCCGTAAACGATGTCGACTTGGAGGTTGTGCCCTTGAGGCGTGGCTTCCGGAGCTAACGCGTTAAGTCGACCGCCTGGGGAGTACGGCCGCAAGGTTAAAACTCAAATGAATTGACGGGGGCCCGCACAAGCGGTGGAGCATGTGGTTTAATTCGATGCAACGCGAAGAACCTTACCTGGTCTTGACATCCACAGAACTTTCCAGAGATG‑GATTGGTG‑CCTTCGGGAACTGT‑GAGACAGGTGCTGCATGGCTGTCGTCAg                                                                                                                                                                                                                              >  1:357158/1‑298 (MQ=9)                      tGGGGAGCAAACAGGATTAGATACCCTGGTAGTCCACGCCGTAAACGATGTCGACTTGGAGGTTGTGCCCTTGAGGCGTGGCTTCCGGAGCTAACGCGTTAAGTCGACCGCCTGGGGAGTACGGCCGCAAGGTTAAAACTCAAATGAATTGACGGGGGCCCGCACAAGCGGTGGAGCATGTGGTTTAATTCGATGCAACGCGAAGAACCTTACCTGGTCTTGACATCCACAGAACTTTCCAGGCATG‑GAAT‑GGGTCCTTCGGGAACTGT‑GAGACAGGGTATGCATGGCttgtggtagc                                                                                                                                                                                                                             >  2:19429/1‑289 (MQ=11)                       ggggAGCAAACAGGATTAGATACCCTGGTAGTCCACGCCGTAAACGATGTCGACTTGGAGGTTGTGCCCTTGAGGCGTGGCTTCCGGAGCTAACGCGTTAAGTCGACCGCCTGGGGAGTACGGCCGCAAGGTTAAAACTCAAATGAATTGACGGGGGCCCGCACAAGCGGTGGAGCATGTGGTTTAATTCGATGCAACGCGAAGAACCTTACCTGGTCTTGACATCCACAGAACTTTCCAGAGATG‑GATTGGTG‑CCTTCGGGAACTGT‑GAGACAGGTGCTGCATGGCTGTCGTCAGCt                                                                                                                                                                                                                            <  2:74231/298‑1 (MQ=9)                        gggAGCAAACAGGATTAGATACCCTGGTAGTCCACGCCGTAAACGATGTCGACTTGGAGGTTGTGCCCTTGAGGCGTGGCTTCCGGAGCTAACGCGTTAAGTCGACCGCCTGGGGAGTACGGCCGCAAGGTTAAAACTCAAATGAATTGACGGGGGCCCGCACAAGCGGTGGAGCATGTGGTTTAATTCGATGCAACGCGAAGAACCTTACCTGGTCTTGACATCCACAGAACTTTCCAGAGATG‑GATTGGTG‑CCTTCGGGAACTGT‑GAGACAGGTGCTGCATGGCTGTCGTCAGCTc                                                                                                                                                                                                                           <  1:425162/298‑1 (MQ=9)                          gAGCAAACAGGATTAGATACCCTGGTAGTCCACGCCGTAAACGATGTCGACTTGGAGGTTGTGCCCTTGAGGCGTGGCTTCCGGAGCTAACGCGTTAAGTCGACCGCCTGGGGAGTACGGCCGCAAGGTTAAAACTCAAATGAATTGACGGGGGCCCGCACAAGCGGTGGAGCATGTGGTTTAATTCGATGCAACGCGAAGAACCTTACCTGGTCTTGACATCCACAGAACTTTCCAGAGATG‑GATTGGTG‑CCTTCGGGAACTGT‑GAGACAGGTGCTGCATGGCTGTCGTCAGCTCgt                                                                                                                                                                                                                         >  1:419438/1‑298 (MQ=9)                            gCAAACAGGATTAGATACCCTGGTAGTCCACGCCGTAAACGATGTCGACTTGGAGGTTGTGCCCTTGAGGCGTGGCTTCCGGAGCTAACGCGTTAAGTCGACCGCCTGGGGAGTACGGCCGCAAGGTTAAAACTCAAATGAATTGACGGGGGCCCGCACAAGCGGTGGAGCATGTGGTTTAATTCGATGCAACGCGAAGAACCTTACCTGGTCTTGACATCCACAGAACTTTCCAGAGATG‑GATTGGTG‑CCTTCGGGAACTGT‑GAGACAGGTGCTGCATGGCTGTCGTCAGCTCgtgt                                                                                                                                                                                                                       <  1:329784/298‑1 (MQ=9)                                  aGGATTAGATACCCTGGTAGTCCACGCCGTAAACGATGTCGACTTGGAGGTTGTGCCCTTGAGGCGTGGCTTCCGGAGCTAACGCGTTAAGTCGACCGCCTGGGGAGTACGGCCGCAAGGTTAAAACTCAAATGAATTGACGGGGGCCCGCACAAGCGGTGGAGCATGTGGTTTAATTCGATGCAACGCGAAGAACCTTACCTGGTCTTGACATCCACAGAACTTTCCAGAGATG‑GATTGGTT‑CCTTCGGGAACCGT‑GAGACAGGTGCTGCCTGGCTGTCGCCAGCCCgtggtggtga                                                                                                                                                                                                                 >  2:146584/1‑294 (MQ=11)                                      ttAGATACCCTGGTAGTCCACGCCGTAAACGATGTCGACTTGGAGGTTGTGCCCTTGAGGCGTGGCTTCCGGAGCTAACGCGTTAAGTCGACCGCCTGGGGAGTACGGCCGCAAGGTTAAAACTCAAATGAATTGACGGGGGCCCGCACAAGCGGTGGAGCATGTGGTTTAATTCGATGCAACGCGAAGAACCTTACCTGGTCTTGACATCCACAGAACTTTCCAGAGATG‑GATTGGTG‑CCTTCGGGAACTGT‑GAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGTGAAATGt                                                                                                                                                                                                             >  1:240375/1‑298 (MQ=9)                                        aGATACCCTGGTAGTCCACGCCGTAAACGATGTCGACTTGGAGGTTGTGCCCTTGAGGCGTGGCTTCCGGAGCTAACGCGTTAAGTCGACCGCCTGGGGAGTACGGCCGCAAGGTTAAAACTCAAATGAATTGACGGGGGCCCGCACAAGCGGTGGAGCATGTGGTTTAATTCGATGCAACGCGAAGAACCTTACCTGGTCTTGACATCCACAGAACTTTCCAGAGATG‑GATTGGTG‑CCTTCGGGAACTGT‑GAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGTGAAATGTTg                                                                                                                                                                                                           <  1:181319/298‑1 (MQ=9)                                                                aaaccgtggcGACTTTGAGTCTTTACCATTGGGGCGTGGCTTCCGGAGCTAACGCTTTAAGTCTGCCGCCTGGGGAGTACGGCCGCAAGGTTAAAACTCAAATGAATTGGCGGGGCCCCGCACAAGCGGTGGAGCATGTGGTTTCATTCGATGCAACGCGAAGAACCTTACCTGGTCTTGACATCCACAGAACTTTCCAGAGATG‑GATTGGTG‑CCTTCGGGAACTGT‑GAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGTGAAATGTTGGGTTAAGTCCCGCAACGAGCGCaa                                                                                                                                                                                   <  2:89846/289‑1 (MQ=11)                                                                 aaCGATGTCGACTTGGAGGTTGTGCCCTTGAGGCGTGGCTTCCGGAGCTAACGCGTTAAGTCGACCGCCTGGGGAGTACGGCCGCAAGGTTAAAACTCAAATGAATTGACGGGGGCCCGCACAAGCGGTGGAGCATGTGGTTTAATTCGATGCAACGCGAAGAAACTTACCTGGTCTTGACATCCACAGAACTTTCCAGAGATG‑GATTGGTG‑CCTTCGGGAACTGT‑gaga                                                                                                                                                                                                                                                      >  2:97117/1‑230 (MQ=11)                                                                 aaCGATGTCGACTTGGAGGTTGTGCCCTTGAGGCGTGGCTTCCGGAGCTAACGCGTTAAGTCGACCGCCTGGGGAGTACGGCCGCAAGGTTAAAACTCAAATGAATTGACGGGGGCCCGCACAAGCGGTGGAGCATGTGGTTTAATTCGATGCAACGCGAAGAAACTTACCTGGTCTTGACATCCACAGAACTTTCCAGAGATG‑GATTGGTG‑CCTTCGGGAACTGT‑gaga                                                                                                                                                                                                                                                      <  1:97117/230‑1 (MQ=11)                                                                  aCGATGTCGACTTGGAGGTTGTGCCCTTGAGGCGTGGCTTCCGGAGCTAACGCGTTAAGTCGACCGCCTGGGGAGTACGGCCGCAAGGTTAAAACTCAAATGAATTGACGGGGGCCCGCACAAGCGGTGGAGCATGTGGTTTAATTCGATGCAACGCGAAGAACCTTACCTGGTCTTGACATCCACAGAACTTTCCAGAGATG‑GATTGGTG‑CCTTCGGGAACTGT‑GAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGGGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAAcc                                                                                                                                                                                 >  1:263327/1‑298 (MQ=9)                                                                      tGTCGACTTGGAGGTTGTGCCCTTGAGGCGTGGCTTCCGGAGCTAACGCGTGAAGTCGACCGCCTGGGGAGTACGGCCGCAAGGTTAAAACTCAAATGAATTGACGGGGGCCCGCACAAGCGGTGGAGCATGTGGTTTAATTCGATGCAACGCGAAGAACCTTACCTGGTCTTGACATCCACAGAACTTTCCAGAGATG‑GATTGGTG‑CCTTCGGGAACTGT‑GAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGTGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTa                                                                                                                                                                             <  1:199863/298‑1 (MQ=9)                                                                          gACTTGGAGGTTGTGCCCTTGAGGCGTGGCTTCCGGAGCTAACGCGTTAAGTCGACCGCCTGGGGAGTACGGCCGCAAGGTTAAAACTCAAATGAATTGACGGGGGCCCGCACAAGCGGTGGAGCATGTGGTTTAATTCGATGCAACGCGAAGAACCTTACCTGGTCTTGACATCCACAGAACTTTCCAGAGATG‑GATTGGTG‑CCTTCGGGAACTGT‑GAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGTGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATcat                                                                                                                                                                         >  1:156240/1‑296 (MQ=9)                                                                              tGGAGGTTGTGCCCTTGAGGCGTGGCTTCCGGAGCTAACGCGTTAAGTCGACCGCCTGGGGAGTACGGCCGCAAGGTTAAAACTCAAATGAATTGACGGGGGCCCGCACAAGCGGTGGAGCATGTGGTTTAATTCGATGCAACGCGAAGAACCTTACCTGGTCTTGACATCCACAGAACTTTCCAGAGATG‑GATTGGTG‑CCTTCGGGAACTGT‑GAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGTGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCCTttgt                                                                                                                                                                     <  1:421265/298‑1 (MQ=9)                                                                                  ggTTGTGCCCTTGAGGCGTGGCTTCCGGAGCTAACGCGTTAAGTCGACCGCCTGGGGAGTACGGCCGCAAGGTTAAAACTCAAATGAATTGACGGGGGCCCGCACAAGCGGTGGAGCATGTGGTTTAATTCGATGCAACGCGAAGAACCTTACCTGGTCTTGACATCCACAGAACTTTCCAGAGATG‑GATTGGTG‑CCTTCGGGAACTGT‑GAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGTGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCCTTTGTTGcc                                                                                                                                                                 <  1:421483/298‑1 (MQ=9)                                                                                  ggTTGTGCCCTTGAGGCGTGGCTTCCGGAGCTAACGCGTTAAGTCGACCGCCTGGGGAGTACGGCCGCAAGGTTAAAACTCAAATGAATTGACGGGGGCCCGCACAAGCGGTGGAGCATGTGGTTTAATTCGATGCAACGCGAAGAACCTTACCTGGTCTTGACATCCACAGAACTTTCCAGAGATG‑GATTGGTG‑CCTTCGGGAACTGT‑GAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGTGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCCTTTGTTGcc                                                                                                                                                                 <  1:216156/298‑1 (MQ=9)                                                                                   ggtgggCCCTTGGGGCGTGGCTTCCGGAGCTAACGCGTTAAGTCGACCGCCTGGGGAGTACGGCCGCAAGGTTAAAACTCAAATGAATTGACGGGGGCCCGCACAAGCGGTGGAGCATGTGGTTTAATTCGATGCAACGCGAAGAACCTTACCTGGGCTTGACATCCACAGAACTTTCCAGAGATG‑GATTGGTG‑CCTTCGGGAACTGT‑GAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGTGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCCTTTGTTGCCa                                                                                                                                                                <  2:200829/293‑1 (MQ=11)                                                                                      gtgCCCTTGGGGCGTGGCTTCGGGGGCTGACGGGTTAAGTCGACCGCCTGGGGAGTACGGGCGCAAGGTTAAAACTCAAATGAATTGACGGGGGCCCGAACAAGCGGTGGAGCATGTGGTTTAAGTCGATGCAACGCGAAGAACCTTACCTGGTCTTGACATCCACAGAACTTTCCAGAGATG‑GATTGGTG‑CCTTCGGGAACTGT‑GAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGTGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCCTTTGTTGCCAGCg                                                                                                                                                             <  2:142989/298‑1 (MQ=11)                                                                                         cgcggggggtttGGCGTCCGGATTTAAAGAGTTAAGAGGACTTGCTGGGGAGTACGGACGCAAGTTTAACACTCATATTAAGTGACGGGGGCCCGCACAAGCGGGGGAGCATGTGTTTTTATTCGATGCACCGGGAAGAGCCTTACCTGGTATTGACATCCACAGAACTTTCCAGAGATG‑GATTGGTG‑CCTTCGGGTACTGT‑GAGACAGGTGCTGCATGGCTGTCGTCGGCTCGTGTTGTGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCCTTTGTTGCCAGCGGTc                                                                                                                                                          <  2:247354/287‑1 (MQ=11)                                                                                          ccTTGAGGCGTGGCTTCCGGAGCTAACGCGTTAAGTCGACCGCCTGGGGAGTACGGCCGCAAGGTTAAAACTCAAATGAATTGACGGGGGCCCGCACAAGCGGTGGAGCATGTGGTTTAATTCGATGCAACGCGAAGAACCTTACCTGGTCTTGACATCCACAGAACTTTCCAGAGATG‑GATTGGTG‑CCTTCGGGAACTGT‑GAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGTGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCCTTTGTTGCCAGCGGTcc                                                                                                                                                         <  1:146584/298‑1 (MQ=9)                                                                                                                  aaCGCGTTAAGTCGACCGCCTGGGGAGTACGACCGCAAGGTTAAAACTCAAATGAATTGACGGGGGCCCGCACAAGCGGTGGAGCATGTGGTTTAATTCGATGCAACGCGAAGAACCTTACCTGGTCTTGACATCCACAGAACTTTCCAGAGATG‑GATTGGTG‑CCTTCGGGAACTGT‑GAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGTGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCCTTTGTTGCCAGCGGTCCGGCCGGTAACTCAAAGGAGACTGc                                                                                                                                 >  1:360924/1‑298 (MQ=255)                                                                                                                          aaGTGGCCCGCCGGGGGAGTGCGGCCTCAAGGTTAAACATAAAATGAATTGCCGGGGGACCGCCAAAGCTGGGGAGCATGTGGTTTAATTCGATGCAACGCGACGAACCTTACCTGGTCTTGACATCCACAGAACTTTCCAGGGATG‑GATTGGTG‑CCTTCGGGAACTGT‑GAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGTGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCCTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATa                                                                                                                         <  2:240375/298‑1 (MQ=11)                                                                                                                                                  ccGCAAGGTTAAAACTCAAATGAATTGACGGGGGCCCGCACAAGCGGTGGAGCTTGTGGTTTAATTCGATGCAACGCGAAGAACCTTACCTGGTCTTGACATCCACAGAACTTTCCAGAGATG‑GATTGGTG‑CCTTCGGGAACTGT‑GAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGTGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCCTTTGTTGCCAGCGGTCCGGCCGGGGACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGGTGGCg                                                                                                 >  2:23514/1‑298 (MQ=11)                                                                                                                                                         gTTAAAACTCAAATGAATTGACGGGGGCCCGCACAAGCGGTGGAGCATGTGGTTTAATTCGATGCAACGCGAAGAACCTTACCTGGTCTTGACATCCACAGAACTTTCCAGAGATG‑GATTGGTG‑CCTTCGGGAACTGT‑GAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGTGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCCTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGtc                                                                                          >  1:321285/1‑298 (MQ=9)                                                                                                                                                                           tGACGGGGGCCCGCACAAGCGGTGGAGCATGTGGTTTAATTCGATGCAACGCGAAGAACCTTAGCTGGTCTTGACATCCACAGAACTTTCCAGAGATG‑GATTGGTG‑CCTTCGGGAACTGT‑GAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGTGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCCTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACGAcc                                                                        <  1:287206/298‑1 (MQ=9)                                                                                                                                                                                 gggCCCGCACAAGCGGTGGAGCATGTGGTTTAATTCGATGCAACGCGAAGAACCTTACCTGGTCTTGACATCCACAGAACTTTCCAGAGATG‑GATTGGTG‑CCTTCGGGAACTGT‑GAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGTGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCCTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACGACCAGGGCt                                                                  <  1:30972/298‑1 (MQ=9)                                                                                                                                                                                  ggCCCGCACAAGCGGTGGAGCATGTGGTTTAATTCGATGCAACGCGAAGAACCTTACCTGGTCTTGACATCCACAGAACTTTCCAGAGATG‑GATTGGTG‑CCTTCGGGAACTGT‑GAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGTGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCCTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACGACCAGGGCTa                                                                 >  1:13738/1‑298 (MQ=9)                                                                                                                                                                                       gCACAAGCGGTGGAGCATGTGGTTTAATTCGATGCAACGCGAAGAACCTTACCTGGTCTTGACATCCACAGAACTTTCCAGAGATG‑GATTGGTG‑CCTTCGGGAACTGT‑GAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGTGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCCTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACGACCAGGGCTacacac                                                            <  1:261439/298‑1 (MQ=9)                                                                                                                                                                                                gTGGAGCATGTGGTTTAATTCGATGCAACGCGAAGAACCTTACCTGGTCTTGACATCCACAGAACTTTCCAGAGATG‑GATTGGTG‑CCTTCGGGAACTGT‑GAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGTGAAATGTTGGGTTAAGTCCCACAACGAGCGCAACCCTTATCCTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACGACCAGGGCTACACACGTGCTACaa                                                   <  1:325735/298‑1 (MQ=9)                                                                                                                                                                                                  ggAGCATGTGGTTTAATTCGATGCAACGCGAAGAACCTTACCTGGTCTTGACATCCACAGAACTTTCCAGAGATG‑GATTGGTG‑CCTTCGGGAACTGT‑GAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGTGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCCTTTGTTGCCAGCGGTCCAGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACGACCAGGGCTACACACGTGCTACAATg                                                 <  1:88530/298‑1 (MQ=9)                                                                                                                                                                                                                aaTTCGATGCAACGCGAAGAACCTTACCTGGTCTTGACATCCACAGAACTTTCCAGAGATG‑GATTGGTG‑CCTTCGGGAACTGT‑GAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGTGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCCTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACGACCAGGGCTACACACGTGCTACAATGGCGCATACAAagag                                   >  1:133946/1‑298 (MQ=9)                                                                                                                                                                                                                               gaagaaCCTTACCTGGTCTTGACATCCACAGAACTTTCCAGAGATG‑GATTGGTG‑CCTTCGGGAACTGT‑GAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGTGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCCTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACGACCAGGGCTACACACGTGCTACAATGGCGCATACAAAGAGAAGCGACCTcgcgcg                    >  1:65864/1‑296 (MQ=9)                                                                                                                                                                                                                                      cTTACCTGGTCTTGACATCCCAAGAACTGTCCAGAGAGGGGTTG‑GGG‑CGTTCGGG‑ACCGTGGAGACAGGTGCTGCATGGCTGGCGTCAGCTTGTGTTGTGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCCTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACGACCAGGGCTACACACGTGCTACAATGGCGCATACAAAGAGAAGCGACCTCGCGAGAGCAAGc             <  2:419438/298‑1 (MQ=11)                                                                                                                                                                                                                                                 ttGACATCCACAGAACTTTCCAGAGATG‑GATTGGTG‑CCTTCGGGAACTGT‑GAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGTGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCCTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACGACCAGGGCTACACACGTGCTACAATGGCGCATACAAAGAGAAGCGACCTCGCGAGAGCAAGCGGACCTCATaa  >  1:51883/1‑298 (MQ=9)                                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                                                 TGACGCTCAGGTGCGAAAGCGTGGGGAGCAAACAGGATTAGATACCCTGGTAGTCCACGCCGTAAACGATGTCGACTTGGAGGTTGTGCCCTTGAGGCGTGGCTTCCGGAGCTAACGCGTTAAGTCGACCGCCTGGGGAGTACGGCCGCAAGGTTAAAACTCAAATGAATTGACGGGGGCCCGCACAAGCGGTGGAGCATGTGGTTTAATTCGATGCAACGCGAAGAACCTTACCTGGTCTTGACATCCACAGAACTTTTCAGAGATGAGAAT‑GTG‑CCTTCGGGAACCGT‑GAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGTGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCCTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACGACCAGGGCTACACACGTGCTACAATGGCGCATACAAAGAGAAGCGACCTCGCGAGAGCAAGCGGACCTCATAA  >  CP006584/4178720‑4179257

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑