Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	CP006584	3,545,697	A→G	noncoding (2803/2931 nt)	LY180_16870 ←	23S ribosomal RNA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP006584	3,545,697	0	A	G	88.6%	49.7 / 2.5	35	noncoding (2803/2931 nt)	LY180_16870	23S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (14/17);  minor base . (4/0);  total (18/17)
Fisher's exact test for biased strand distribution p-value = 1.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

ACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTAAGGCCGCCAGGCAAATTCTGTTTTATCAGACCGCTTCTGCGTTCTGATTTAATCTGTATCAGGCTGAAAATCTTCTCTCATCCGCCAAAACACCTTCGGCGTTGTAAGGTTAAGCCTCACGGTTCATTAGTACCGGTTAGCTCAACGCATCGCTGCGCTTACACACCCGGCCTATCAACGTCGTCGTCTTCAACGTTCCTTCAGGA‑CTC‑TCAAGGAGTCA‑GGGAGAACTCATCTCGGGGCAAGTTTCGTGCTTAGATGCTTTCAGCACTTATCTCTTCCGCATTTAGCTACCGGGCA  >  CP006584/3545413‑3545778                                                                                                                                                                                                                                                                                               |                                                                                   acacTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTAAGGCCGCCAGGCAAATTCTGTTTTATCAGACCGCTTCTGCGTTCTGATTTAATCTGTATCAGGCTGAAAATCTTCTCTCATCCGCCAAAACACCTTCGGCGTTGTAAGGTTAAGCCTCACGGTTCATTAGTACCGGTTAGCTCAACGCATCGCTGCGCTTACACACCCGGCCTATCAACGTCGTCGTCTTCAACGTTCCTTCAGGA‑CCC‑TTAAAGGGTCA‑GGGAGAACt                                                                      >  1:742140/1‑298 (MQ=14) acacTACCATCGGCGCTACGGCGTTTCACTGCTGAGTTCGGCATGGGGTCAGGCGGGACCACCGCGCTAAGGCCGCCAGGCCAATTCTGTTTTATCAGACCGCCTCTGCGTTCTGATTTAATCTGTATCAGGCTGAAAATCTCCTCTCATCCGCCCAAACACCTTCGGCGTTGTCAGGTTAAGCCTCACGGTCCATTAGTACCGGTTAGCTCAACGCATCGCTGCGCTTACCACCCCGGCCTCCCAACCTCTTCCTACTCCACGTTCCTTCACGACCCC‑TCAAGG‑GTCA‑CGGAGAAcg                                                                      >  2:130502/1‑297 (MQ=11)       cattgggggcTCCGCCGTTTGACTTCTGGGTTTGGAATGGGGTCTGTGGTGACCACCGCGCTAAGCCGGCCAGGCAAATTCTGTGTTATCAGACCGCTTCGCCGTTCTGATTTAATCTGTATCAGGCTGAAAATCTTCTATCATCAGCCAAAACACCTTCGGCGTTGTAAGGTTAAGCCTCACGGTTCATTAGTACCGGTTAGCTCAACGCATCGCTGCGCTTACACACCCGGCCTATCAACGTCGTCGTCTTCAACGTTCCTTCAGGA‑CCC‑TTAAAGGGTCA‑GGGAGAACTCAtctc                                                                <  2:405609/290‑1 (MQ=11)         aTCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTAAGGCCGCCAGGCAAATTCTGTTTTATCAGACCGCTTCTGCGTTCTGATTTAATCTGTATCAGGCTGAAAATCTGCTCTCATCCGCCAAAACAACTTCGGCGTTGTAAGGTTAAGCCTCACGGTTCATTAGTACCGGTTAGCTCAACGCATCGCTGCGCTTACACAACCTGCCTATCAACGTCGTCGTCTTCAACGTTCCTTCAGGA‑CCC‑TTAAAGGGTCA‑GGGAGAACTCAtctccg                                                              >  1:614328/1‑296 (MQ=9)                 tACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTAAGGCCGCCAGGCAAATTCTGTTTTATCAGACCGCTTCTGCGTTCTGATTTAATCTGTATCAGGCTGAAAATCTTCTCTCATCCGCCAAAACACCTTCGGCGTTGTAAGGTTAAGCCTCACGGTTCATTAGTACCGGTTAGCTCAACGCATCGCTGCGCTTACACACCCGGCCTATCAACGTCGTCGTCTTCAACGTTCCTTCAGGA‑CCC‑TTAAAGGGTCA‑GGGAGAACTCATCTCGGGGCAAGtt                                                      >  1:25476/1‑298 (MQ=14)                 gaCGGCGTTCCCCTGCTGAGTTCGGCATGGGGGCAAGGGGGACCACCGCGCTAAGGCCGCCAGGCAAATTCTGTTGTATGCGACCGCTTTTGCGTTCTGATTTAATCTGTATCAGGCTGAAAATCTTCCCGCCTCCGCCAAAACACCTTCGGCGTTGTAAGGTTAAGCCTCACGGTTCATTAGTACCGGTTAGCTCAACGCATCGCTGCGCTTACACACCCGGCCTATCAACGTCGTCGTCTTCAACGTTCCTTCAGGA‑CCC‑TTAAAGGGTCA‑GGGAGAACTCATCTCGGGGCAAGtt                                                      <  2:311633/297‑1 (MQ=11)                  aCGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTAAGGCCGCCAGGCAAATTCTGTTTTATCAGACCGCTTCTGCGTTCTGATTTAATCTGTATCAGGCTGAAAATCTTCTCTCATCCGCCAAAACACCTTCGGCGTTGTAAGGTTAAGCCTCACGGTTCATTAGTACCGGTTAGCTCAACGCATCGCTGCGCTTACACACCCGGCCTATCAACGTCGTCGTCTTCAACGTTCCTTCAGGA‑CCC‑TCAACGGCTCA‑GGGAGAACTCATCCACGGCCAAGttt                                                     >  2:1000851/1‑298 (MQ=9)                    ggCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTAAGGCCGCCAGGCAAATTCTGTTTTATCAGACCGCTTCTGCGTTCTGATTTAATCTGTATCAGGCTGAAAATCTTCTCTCATCCGCCAAAACACCTTCGGCGTTGTAAGGTTAAGCCTCACGGTTCATTAGTACCGGTTAGCTCAACGCATCGCTGCGCTTACACACCCGGCCTATCAACGTCGTCGTCTTCAACGTTCCTTCAGGA‑CCC‑TTAAAGGGTCA‑GGGAGAACTCATCTCGGGGCAAGTTTCg                                                   <  1:918694/298‑1 (MQ=14)                     tgttttCACTTCTGGGTTGGGCAGGGGGTCCGGTGGGACCAACGCGCTAAGGCCGCCAGGCAAATTCTGTTTTATCAGACCGCTTCTGCGTTCGGATTTAATCTGTATCAGGCTGAAAAGCTTCTCTCATCCGCCAAAACACCTTCGGCGTTGTAAGGTTAAGCCTCACGGTTCATTAGTACCGGTTAGCTCAACGCATCGCTGCGCTTACACACCCGGCCTATCAACGTCGTCGTCTTCAACGTTCCTTCAGGA‑CCC‑TTAAAGGGTCA‑GGGAGAACTCATCTCGGGGCAAGTTTCGt                                                  <  2:681468/295‑1 (MQ=11)                     gCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTAAGGCCGCCAGGCAAATTCTGTTTTATCAGACCGCTTCTGCGTTCTGATTTAATCTGTATCAGGCTGAAAATCTTCTCTCATCCGCCAAAACACCTTCGGCGTTGTAAGGTTAAGCCTCACGGTTCATTAGTACCGGTTAGCTCAACGCATCGCTGCGCTTACACACCCGGCCTATCAACGTCGTCGTCTTCAACGTTCCTTCAGGA‑CCC‑TTAAAGGGTCA‑GGGAGAACTCATCTCGGGGCAAGTTTCGt                                                  <  1:520185/298‑1 (MQ=14)                        tttCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTAAGGCCGCCAGGCAAATTCTGTTTTATCAGACCGCTTCTGCGTTCTGATTTAATCTGTATCAGGCTGAAAATCTTCTCTCATCCGCCAAAACACCTTCGGCGTTGTAAGGTTAAGCCTCACGGTTCATTAGTACCGGTTAGCTCAACGCATCGCTGCGCTTACACACCCGGCCTATCAACGTCGTCGTCTTCAACGTTCCTTCAGGA‑CCC‑TTAAAGGGTCA‑GGGAGAACTCATCTCGGGGCAAGTTTCGTggt                                               >  1:521226/1‑296 (MQ=14)                          tCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTAAGGCCGCCAGGCAAATTCTGTTTTATCAGACCGCTTCTGCGTTCTGATTTAATCTGTATCAGGCTGAAAATCTTCTCTCATCCGCCAAAACACCTTCGGCGTTGTAAGGTTAAGCCTCACGGTTCATTAGTACCGGTTAGCTCAACGCATCGCTGCGCTTACACACCCGGCCTATCAACGTCGTCGTCTTCAACGTTCCTTCCGGACCTG‑TAAAGG‑CTTA‑GGGATAAACCCTTCCGGGGCAGGTTTCGTCCTTa                                             >  2:944196/1‑298 (MQ=11)                          ggccttgtggatttggCATGGGGTCAGGTGGGCCCACCGCGCTAAGGCCGCAAGGCAAATTCTGTTTTATCAGACCGCTGCTGCGTTCTGATTTAATCTGTATCATGCTGAAAATCTTCTCTCCTCCGCCAAAACACCTTCGGCGTTGTAAGGTTAAGCCTCACGGTTCATTAGTACCGGTTAGCTCAACGCATCGCTGCGCTTACACACCCGGCCTATCAACGTCGTCGTCTTCAACGTTCCTTCAGGA‑CCC‑TTAAAGGGTCA‑GGGAGAACTCATCTCGGGGCAAGTTTCGTGCTTa                                             <  2:902104/284‑1 (MQ=11)                             cTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTAAGGCCGCCAGGCAAATTCTGTTTTATCAGACCGCTTCTGCGTTCTGATTTAATCTGTATCAGGCTGAAAATCTTCTCTCATCCGCCAAAACACCTTCGGCGTTGTAAGGTTAAGCCTCACGGTTCATTAGTACCGGTTAGCTCAACGCATCGCTGCGCTTACACACCCGCCCTATCAACGCCCTCGTCTTCACCGTTCCTTCAGGC‑CCC‑TTAACGGGCCC‑GGCAGCCCCCACACCGGGACACGTTTCGTGCTTAGAt                                          >  2:118763/1‑298 (MQ=11)                               tgtgtggtgggCGTGGGGTCTGGGGGGACCCACGGGTTAAGGCCGCCAGGCAAATTCTGTTGTATCAGACCGCTTCTGCGTTCTGCTTTAATCTGGTTCAGGCTGAAAATCTTCTCTCATCCGCCAAAACACCTTCGGCGTTGTAAGGTTAAGCCTCACGGTTCATTAGTACCGGTTAGCTCAACGCATCGCTGCGCTTACACACCCGGCCTATCAACGTCGTCGTCTTCAACGTTCCTTCAGGA‑CCC‑TTAAAGGGTCA‑GGGAGAACTCATCTCGGGGCAAGTTTCGTGCTTAGATGc                                        <  2:977286/289‑1 (MQ=11)                               tCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTAAGGCCGCCAGGCAAATTCTGTTTTATCAGACCGCTTCTGCGTTCTGATTTAATCTGTATCAGGCTGAAAATCTTCTCTCATCCGCCAAAACACCTTCGGCGTTGTAAGGTTAAGCCTCACGGTTCATTAGTACCGGTTAGCTCAACGCATCGCTGCGCTTACACACCCGGCCTATCAACGTCGTCGTCTTCAACGTTCCTTCAGGA‑CCC‑TTAAAGGGTCA‑GGGAGAACTCATCTCGGGGCAAGTTTCGTGCTTAGATGc                                        >  1:1444/1‑298 (MQ=14)                                  gAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTAAGGCCGCCAGGCAAATTCTGTTTTATCAGACCGCTTCTGCGTTCTGATTTAATCTGTATCAGGCTGAAAATCTTCTCTCATCCGCCAAAACACCTTCGGCGTTGTAAGGTTAAGCCTCACGGTTCATTAGTACCGGTTAGCTCAACGCATCGCTGCGCTTACACACCCGGCCTATCAACGTCGTCGTCTTCAACGTTCCTTCAGGA‑CCC‑TTAAAGGGTCA‑GGGAGAACTCATCTCGGGGCAAGTTTCGTGCTTAGATGCttt                                     >  1:659915/1‑298 (MQ=14)                                   aGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTAAGGCCGCCAGGCAAATTCTGTTTTATCAGACCGCTTCTGCGTTCTGATTTAATCTGTATCAGGCTGAAAATCTTCTCTCATCCGCCAAAACACCTTCGGCGTTGTAAGGTTAAGCCTCACGGTTCATTAGTACCGGTTAGCTCAACGCATCGCTGCGCTTACACACCCGGCCTATCAACGTCGTCGTCTTCAACGTTCCTTCAGGA‑CCC‑TTAAAGGGTCA‑GGGAGAACTCATCTCGGGGCAAGTTTCGTGCTTAGATGCTTTc                                    <  1:609758/298‑1 (MQ=14)                                           aTGGGGTCAGGTGGGACCACCGCGCTAAGGCCGCCAGGCAAATTCTGTTTTATCAGACCGCTTCTGCGTTCTGATTTAATCTGTATCAGGCTGAAAATCTTCTCTCATCCGCCAAAACACCTTCGGCGTTGTAAGGTTAAGCCTCACGGTTCATTAGTACCGGTTAGCGCAACGCATCGCTGCGCTTACACACCCGGCCTATCAACGTCGTCGTCTTCAACGTTCCTTCAGGA‑CCC‑TTAAAGGGTCA‑GGGAGAACTCATCTCGGGGCAAGTTTCGTGCTTAGATGCTTTCAGCACTTa                            >  1:993909/1‑298 (MQ=9)                                            tGGGGTCAGGTGGGACCACCGCGCTAAGGCCGCCAGGCAAATTCTGTTTTATCAGACCGCTTCTGCGTTCTGATTTAATCTGTCTCAGGCTGAAAATCTTCTCTCATCCGCCAAAACACCTTCGGCGTTGTAAGGTTAAGCCTCACGGTTCATTAGTCCCGGTTAGCCCCACGCATCGCTGCGCTTACACACCCGGCCTACCAACGTCGTCGTCTTCCACGTTCCTTACGGA‑CCC‑TTACAGGGCTA‑GGGAGAACTCAACCCGGGGCAAGGTTCCTGCCTAAATGCCTTCACCACttca                           >  2:17807/1‑296 (MQ=11)                                            tGGGGTCAGGTGGGACCACCGCGCTAAGGCCGCCAGGCAAATTCTGTTTTATCAGACCGCTTCTGCGTTCTGATTTAATCTGTATCAGGCTGAAAATCTTCTCTCATCCGCCAAAACACCTTCGGCGTTGTAAGGTTAAGCCTCACGGTTCATTAGTACCGGTTAGCTCAACGCATCGCTGCGCTTACACACCCGGCCTATCAACGTCGTCGTCTTCAACGTTCCTTCAGGA‑CCC‑TTAAAGGGTCA‑GGGAGAACTCATCTCGGGGCAAGTTTCGTGCTTAGATGCTTTCAGCACTTAt                           >  1:676720/1‑298 (MQ=14)                                              gggtagcgTGGGACCACCGGGCTGAGGCCGCCAGGCAAATTCTGTTTTTTCAGACCGCTTCGGCGTTCGGATTTAATTTGTATCAGGCTGAAGATCTTCTCTCATCCGCCAAAACACCGGCGGCGTTGTAATGTTAAGCCTCACGGTTCATTAGTACCGGTTAGCGCAACGCATCGCTGCGCTTACACACCCGGCCTATCAACGGCGTCGTCTTCAACGTTGCTTCAGGA‑CTC‑TTAAAGGGTCA‑GGGAGAACTCATCTCGGGGCAAGTTTCGTGCTTAGATGCTTTCAGCACTTAtct                         <  2:183754/291‑1 (MQ=11)                                                gTCAGGTGGGACCACCGCGCTAAGGCCGCCAGGCAAATTCTGTTTTATCAGACCGCTTCTGCGTTCTGATTTAATCTGTATCAGGCTGAAAATCTTCTCTCATCCGCCAAAACACCTTCGGCGTTGTAAGGTTAAGCCTCACGGTTCATTAGTACCGGTTAGCTCAACGCATCGCTGCGCTTACACACCCGGCCTATCAACGTCGTCGTCTTCAACGTTCCTTCAGGA‑CCC‑TTAAAGGGTCA‑GGGAGAACTCATCTCGGGGCAAGTTTCGTGCTTAGATGCTTTCAGCACTTATctct                       >  1:374075/1‑298 (MQ=14)                                                gTCAGGGGGGACCACCGCGCTAAGGCCGCCAGGCAAATTCTGTTTTATCAGACCGCTTCTGCGTTCTGATTTAATCTGTATCAGGCTGAAAATCTTCTCTCATCCGCCAAAACACCTTCGGCGTTGTAAGGTTAAGCCTCACGGTTCATTAGTACCGGTTAGCTCAACGCATCGCTGCGCTTACACACCCGGCCTATCAACGTCGTCGTCTTCAACGTTCCTTCAGGA‑CCC‑TTAAAGGGTCA‑GGGAGAACTCATCTCGGGGCAAGTTTCGTGCTTAGATGCTTTCAGCACTTATctct                       <  1:264386/298‑1 (MQ=9)                                                   aGGTGGGACCACCGCGCTAAGGCCGCCAGGCAAATTCTGTTTTATCAGACCGCTTCTGCGTTCTGATTTAATCTGTATCAGGCTGAAAATCTTCTCTCATCCGCCAAAACACCTTCGGCGTTGTAAGGTTAAGCCTCACGGTTCATTAGTACCGGTTAGCTCAACGCATCGCTGCCCTTACACACCCGGCCGATCCACGTCGTCGTCTTCAACGTTCCTTCCGGA‑CCC‑TTAAAGGGTCA‑GGGAGAACTCATCTCGGGGCAAGTTTCCTGCCTAAATGCTTTCACCACTTAGATCTTcc                    >  2:906559/1‑298 (MQ=11)                                                    ggTGGGACCACCGCGCTAAGGCCGCCAGGCAAATTCTGTTTTATCAGACCGCTTCTGCGTTCTGATTTAATCTGTATCAGGCTGAAAATCTTCTCTCATCCGCCAAAACACCTTCGGCGTTGTAAGGTTAAGCCTCACGGTTCATTAGTACCGGTTAGCTCAACGCATCGCTGCGCTTACACACCCGGCCTATCAACGTCGTCGTCTTCAACGTTCCTTCAGGA‑CCC‑TTAAAGGGTCA‑GGGAGAACTCATCTCGGGGCAAGTTTc                                                    <  1:224580/265‑1 (MQ=14)                                                    ggTGGGACCACCGCGCTAAGGCCGCCAGGCAAATTCTGTTTTATCAGACCGCTTCTGCGTTCTGATTTAATCTGTATCAGGCTGAAAATCTTCTCTCATCCGCCAAAACACCTTCGGCGTTGTAAGGTTAAGCCTCACGGTTCATTAGTACCGGTTAGCTCAACGCATCGCTGCGCTTACACACCCGGCCTATCAACGTCGTCGTCTTCAACGTTCCTTCAGGA‑CCC‑TTAAAGGGTCA‑GGGAGAACTCATCCCGGGGCAAGTTTc                                                    >  2:224580/1‑265 (MQ=9)                                                     gTGGGACCACCGCGCTAAGGCCGCCAGGCAAATTCTGTTTTATCAGACCGCTTCTGCGTTCTGATTTAATCTGTATCAGGCTGAAAATCTTCTCTCATCCGCCAAAACACCTTCGGCGTTGTAAGGTTAAGCCTCACGGTTCATTAGTACCGGTTAGCTCAACGCATAGACTCGCTTACACACCAGGCCTCTCACCGGCTTAGTCTTCACAGTTCCTTCCGGA‑CCC‑TTAAAG‑GACACGGGAGACCCCATCTCGGGTCCAGTTACGGGCCTTAAGTCTTTCAGCGTTTATATCTCCCGc                  >  2:1014295/1‑298 (MQ=11)                                                       gggACCACCGCGCTAAGGCCGCCAGGCAAATTCTGTTTTATCAGACCGCTTCTGCGTTCTGATTTAATCTGTATCAGGCTGAAAATCTTCTCTCATCCGCCAAAACACCTTCGGCGTTGTAAGGTTAAGCCTCACGGTTCATTAGTACCGGTTAGCTCAACGCATCGCTGCGCTTACACACCCGGCCTATCAACGTCGTCGTCTTCAACGTTCCTTCAGGA‑CCC‑TTAAAGGGTCA‑GGGAGAACTCATCTCGGGGCAAGTTTCGTGCTTAGATGCTTTCAGCACTTATCTCTTCCGCAt                <  1:626756/298‑1 (MQ=14)                                                          accaccGCGCTAAGGCCGCCAGGCAAATTCTGTTTTATCAGACCGCTTCTGCGTTCTGATTTAATCTGTATCAGGCTGAAAATCTTCTCTCATCCGCCAAAACACCTTCGGCGTTGTAAGGTTAAGCCTCACGGTTCATTAGTACCGGTTAGCTCAACGCATCGCTGCGCTTACACACCCGGCCTATCAACGTCGTCGTCTTCAACGTTCCTTCAGGA‑CCC‑TTAAAGGGTCA‑GGGAGAACTCATCTCGGGGCAAGTTTCGTGCTTAGATGCTTTCAGCACTTATCTCTTCCGCATTTa             <  2:116526/298‑1 (MQ=14)                                                           ccaccGCGCTAAGGCCGCCAGGCAAATTCTGTTTTATCAGACCGCTTCTGCGTTCTGATTTAATCTGTATCAGGCTGAAAATCTTCTCTCATCCGCCAAAACACCTTCGGCGTTGTAAGGTTAAGCCTCACGGTTCATTAGTACCGGTTAGCTCAACGCATCGCTGCGCTTACACACCCGGCCTATCAACGTCGTCGTCTTCAACGTTCCTTCAGGA‑CCC‑TTAAAGGGTCA‑GGGAGAACTCATCTCGGGGCAAGTTTCGTGCTTAGATGCTTTCAGCACTTATCTCTTCCGCATTTAg            <  1:128254/298‑1 (MQ=14)                                                            caccGCGCTAAGGCCGCCAGGCAAATTCTGTTTTATCAGACCGCTTCTGCGTTCTGATTTAATCTGTATCAGGCTGAAAATCTTCTCTCATCCGCCAAAACACCTTCGGCGTTGTAAGGTTAAGCCTCACGGTTCATTAGTACCGGTTAGCTCAACGCATCGCTGCGCTTACACACCCGGCCTATCAACGTCGTCGTCTTCAACGTTCCTTCAGGA‑CCC‑TTAAAGGGTCA‑GGGAGAACTCATCTCGGGGCAAGTTTCGTGCTTAGATGCTTTCAGCACTTATCTCTTCCGCATTTAGc           <  1:286816/298‑1 (MQ=14)                                                              ccGCGCTAAGGCCGCCAGGCAAATTCTGTTTTATCAGACCGCTTCTGCGTTCTGATTTAATCTGTATCAGGCTGAAAATCTTCTCTCATCCGCCAAAACACCTTCGGCGTTGTAAGGTTAAGCCTCACGGTTCATTAGTACCGGTTAGCTCAACGCATCGCTGCGCTTACACACCCGGCCTATCAACGTCGTCGTCTTCAACGTTCCTTCAGGA‑CCC‑TTAAAGGGTCA‑GGGAGAACTCATCTCGGGGCAAGTTTCGTGCTTAGATGCTTTCAGCACTTATCTCTTCCGCATTTAGCTa         <  1:282244/298‑1 (MQ=14)                                                                     aaGGCCGCCAGGCAAATTCTGTTTTATCAGACCGCTTCTGCGTTCTGATTTAATCTGTATCAGGCTGAAAATCTTCTCTCATCCGCCAAAACACCTTCGGCGTTGTAAGGTTAAGCCTCACGGTTCATTAGTACCGGTTAGCTCAACGCATCGCTGCGCTTACACACCCGGCCTATCAACGTCGTCGTCTTCAACGTTCCTTCAGGA‑CCCTTCAAGG‑GTCA‑GGGAGAACTCATCCCGGGGCAAGTTTCGCGCTTAGATGCTTTCAGCACTTATCACTTCCGCATTGAGCTACCgcgac  >  2:365374/1‑294 (MQ=11)                                                                     aaGGCCCCCAGGCAAATTCTGTTTTATCAGGCCCCTTTTGCGTTCTGATTTAATCTGTATCAGGCTGAAAACCTTCTCTCAGCCGCTAAAACACCTTCGGCGTTGTAAGGTTAAGCCTCACGGTTCATTAGTACCGGTTAGCTCAACGCATCGCTGCGCTTACACACCCGGCCTATCAACGTCGTCGTCTTCAACGTTCCTTCAGGA‑CCC‑TTAAAGGGTCA‑GGGAGAACTCATCTCGGGGCAAGTTTCGTGCTTAGATGCTTTCAGCACTTATCTCTTCCGCATTTAGCTACCGGGCa  <  2:452168/298‑1 (MQ=11)                                                                                                                                                                                                                                                                                               |                                                                                   ACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTAAGGCCGCCAGGCAAATTCTGTTTTATCAGACCGCTTCTGCGTTCTGATTTAATCTGTATCAGGCTGAAAATCTTCTCTCATCCGCCAAAACACCTTCGGCGTTGTAAGGTTAAGCCTCACGGTTCATTAGTACCGGTTAGCTCAACGCATCGCTGCGCTTACACACCCGGCCTATCAACGTCGTCGTCTTCAACGTTCCTTCAGGA‑CTC‑TCAAGGAGTCA‑GGGAGAACTCATCTCGGGGCAAGTTTCGTGCTTAGATGCTTTCAGCACTTATCTCTTCCGCATTTAGCTACCGGGCA  >  CP006584/3545413‑3545778

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