Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	CP006584	4,178,979	T→C	noncoding (1009/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,979	0	T	C	97.6%	62.4 / NA	42	noncoding (1009/1541 nt)	LY180_19945	16S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (16/25);  total (17/25)
Fisher's exact test for biased strand distribution p-value = 4.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

CGGTGGCGAAGGCGGCCCCCTGGACGAAGACTGACGCTCAGGTGCGAAAGCGTGGGGAGCAAACAGGATTAGATACCCTGGTAGTCCACGCCGTAAACGATGTCGACTTGGAGGTTGTGCCCTTGAGGCGTGGCTTCCGGAGCTAACGCGTTAAGTCGACCGCCTGGGGAGTACGGCCGCAAGGTTAAAACTCAAATGAATTGACGGGGGCCCGCACAAGCGGTGGAGCATGTGGTTTAATTCGATGCAACGCGAAGAACCTTACCTGGTCTTGACATCCACAGAACTTTTCAGAGATGAGAAT‑GTG‑CCTTCGGGAACCGT‑GAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGTGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCCTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACGACCAGGGCTACACACGTGCTACAATGGCGCATACAAAGAGAAGCGACCTCGCGAGAGCAAGCGGACCTCATAA  >  CP006584/4178689‑4179257                                                                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                                                          cGGTGGCGAAGGCGGCCCCCTGGACGAAGACTGACGCTCAGGTGCGAAAGCGTGGGGAGCAAACAGGATTAGATACCCTGGTAGTCCACGCCGTAAACGATGTCGACTTGGAGGTTGTGCCCTTGAGGCGTGGCTTCCGGAGCTAACGCGTTAAGTCGACCGCCTGGGGAGTACGGCCGCAAGGTTAAAACTCAAATGAATTGACGGGGGCCCGCACAAGCGGTGGAGCATGTGGTTTAATTCGATGCAACGCGAAGAACCTTACCTGGTCTTGACATCCACAGAACTTTCCAGAGAt                                                                                                                                                                                                                                                                                    <  1:20856/298‑1 (MQ=17)        gagggcggcCCCCTGAAGGAGGAGGGGCGCTCAGGGGGGAAAGCGTGGGGAGCAAACAGGTTTGGATACCCTGGTAGTCCACGCCGTAAACGATGTCGACTTGGAGGTTGGGCCCTTGAGGCGTGGCTTCCGGAGCTAACGCGTTAAGTCGACCGCCTGGGGAGTACGGCCGCAAGGTTAAAACTCAAATGAATTGACGGGGGCCCGCACAAGCGGTGGAGCATGTGGTTTAATTCGATGCAACGCGAAGAACCTTACCTGGTCTTGACATCCACAGAACTTTCCAGAGATggattgg                                                                                                                                                                                                                                                                             <  2:386048/295‑7 (MQ=11)               gcCCCCTGGACGAAGACTGACGCTCAGGTGCGAAAGCGTGGGGAGCAAACAGGATTAGATACCCTGGTAGTCCACGCCGTAAACGATGTCGACTTGGAGGTTGTGCCCTTGAGGCGTGGCTTCCGGAGCTAACGCGTTAAGTCGACCGCCTGGGGAGTACGGCCGCAAGGTTAAAACTCAAATGAATTGACGGGGGCCCGCACAAGCGGTGGAGCATGTGGTTTAATTCGATGCAACGCGAAGAACCTTACCTGGTCTTGACATCCACAGAACTTTCCAGAGATggattggtgccttc                                                                                                                                                                                                                                                                      <  1:136983/298‑14 (MQ=12)                                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                                                                                                                                                                                                                                                      <  1:97117/230‑1 (MQ=11)                                                                                                aaCGATGTCGACTTGGAGGTTGTGCCCTTGAGGCGTGGCTTCCGGAGCTAACGCGTTAAGTCGACCGCCTGGGGAGTACGGCCGCAAGGTTAAAACTCAAATGAATTGACGGGGGCCCGCACAAGCGGTGGAGCATGTGGTTTAATTCGATGCAACGCGAAGAAACTTACCTGGTCTTGACATCCACAGAACTTTCCAGAGATG‑GATTGGTG‑CCTTCGGGAACTGT‑gaga                                                                                                                                                                                                                                                      >  2:97117/1‑230 (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)                                                                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                                                          CGGTGGCGAAGGCGGCCCCCTGGACGAAGACTGACGCTCAGGTGCGAAAGCGTGGGGAGCAAACAGGATTAGATACCCTGGTAGTCCACGCCGTAAACGATGTCGACTTGGAGGTTGTGCCCTTGAGGCGTGGCTTCCGGAGCTAACGCGTTAAGTCGACCGCCTGGGGAGTACGGCCGCAAGGTTAAAACTCAAATGAATTGACGGGGGCCCGCACAAGCGGTGGAGCATGTGGTTTAATTCGATGCAACGCGAAGAACCTTACCTGGTCTTGACATCCACAGAACTTTTCAGAGATGAGAAT‑GTG‑CCTTCGGGAACCGT‑GAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGTGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCCTTTGTTGCCAGCGGTCCGGCCGGGAACTCAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACGACCAGGGCTACACACGTGCTACAATGGCGCATACAAAGAGAAGCGACCTCGCGAGAGCAAGCGGACCTCATAA  >  CP006584/4178689‑4179257

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