Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_007779	3,777,754	A→C	32.5%	T315P (ACC→CCC)	Y75_RS18735 →	PTS transporter subunit EIIC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007779	3,777,754	0	A	C	32.5%	29.1 / 12.1	43	T315P (ACC→CCC)	Y75_RS18735	PTS transporter subunit EIIC
Reads supporting (aligned to +/- strand):  ref base A (20/8);  new base C (2/12);  total (22/21)
Fisher's exact test for biased strand distribution p-value = 7.78e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.25e-02

GGACAGTTTATCTTTGGTCCGGCAGCTGTTGAAGGCGGCATTCAGATGTACTGGGCGCAGCATCTGCAAGAGTTCAGTTTGAGCGCCGAGCCGCTGAAATCGTTGTTCCCGGAAGGCGGTTT‑TGCCCTGCACGGTAACTCAAAAATCTTTGGTGCC‑GTGGGCATTTCTTT‑AGCGATGTAC‑TTCACTG‑CCGCACCGGAAAATCGGGTAAAA‑GTGGCGGG‑C‑TTG‑CT‑GA‑TTCCCGCAACC‑TTAACCGCCATGCT‑GGTGGGAATT‑ACCGAACCGCTGGAATTTACCTTCCTGTTCATTTCACCGTTGCTGTTTGCGGTACACGCCGTGCTGGCGGCCTCAATGTCGACCGTAATGTATCTCTTTGGTGTGGTGGGCAACATGGGCGGAGGTCTGATTGACTAGGTTTTACCGCAAAACTGGATCCCGATGTTCAGCAACCACGCGGATATGATGCTGACCCAAATCGCCATTGGGTT  >  NC_007779/3777520‑3777992                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                  ggACAGTTTATCTTTGGTCCGGCAGCTGTTGAAGGCGGCATTCAGATGTACTGGGCGCAGCATCTGCAAGAGTTCAGTTTGAGCGCCGAGCCGCTGAAATCGTTGTTCCCGGAAGGCGGTTT‑TGCCCTGCACGGTAACTCAAAAATCTTTGGTGCC‑GTGGGCATTTCTTT‑AGCGATGTAC‑TTCACTG‑CCGCACCGGAAGATCGGGTAAAA‑GTGGCGGG‑G‑TTT‑CT‑GG‑TTCTCCCAACC‑tt                                                                                                                                                                                                                                              >  1:1421049/1‑239 (MQ=255)                 gTCCGGCAGCTGTTGAAGGCGGCATTCAGATGTACTGGGCGCAGCATCTGCAAGAGTTCAGTTTGAGCGCCGAGCCGCTGAAATCGTTGTTCCCGGAAGGCGGTTT‑TGCCCTGCACGGTAACTCAAAAATCTTTGGTGCC‑GTGGGCATTTCTTT‑AGCGATGTAC‑TTCACTG‑CCGCACCGGAAAATCGGGTAAAA‑GTGGCGGG‑G‑TTT‑GT‑GA‑TTTCCGCCACC‑CTAAACGccccgccgggg                                                                                                                                                                                                                              >  1:1413727/1‑230 (MQ=255)                 gTCCGGCAGCTGTTGAAGGCGGCATTCAGATGTACTGGGCGCAGCATCTGCAAGAGTTCAGTTTGAGCGCCGAGCCGCTGAAATCGTTGTTCCCGGAAGGCGGTTT‑TGCCCTGCACGGTAACTCAAAAATCTTTGGTGCC‑GTGGGCATTTCTTT‑AGCGATGTAC‑TTCACTG‑CCGCACCGGAAAATCGGGTAAAA‑GTGGCGGG‑G‑TTG‑CT‑GA‑TTTCCGCCACC‑CTAACCGCCATGCt‑ggtg                                                                                                                                                                                                                             >  1:860311/1‑239 (MQ=255)                 gTCCGGCAGCTGTTGAAGGCGGCATTCAGATGTACTGGGCGCAGCATCTGCAAGAGTTCAGTTTGAGCGCCGAGCCGCTGAAATCGTTGTTCCCGGAAGGCGGTTT‑TGCCCTGCACGGTAACTCAAAAATCTTTGGTGCC‑GTGGGCATTTCTTT‑AGCGATGTAC‑TTCACTG‑CCGCACCGGAAAATCGGGTAAAA‑GTGGCGGG‑C‑TTT‑CT‑GG‑TTTCCCCCACC‑TTAAACGCCATGCt‑gggg                                                                                                                                                                                                                             >  1:1288643/1‑237 (MQ=255)                   cgggCAGTGGTGGAGGGGGGCATTAAGATGTCCTGGGCGCAGAATTTGAAAGATTTCAGTTGGGGCGCGGAGCCGCTGAAATCGTTGTTCCCGGAAGGCGGTTT‑TGCCCTGCACGGTAACTAAAAAATCTTTGGTGCC‑GTGGGCTTTTTTTT‑AGCGATGTAC‑TTCACTG‑CCGCACCGGAAAATCGGGTAAAA‑GTGGCGGG‑C‑TTG‑CT‑GA‑TTCCCGCAACC‑TTAACCGCCATGCT‑GGTggg                                                                                                                                                                                                                           <  1:1420307/237‑1 (MQ=255)                                    cggcATTCAGATGTACTGGGCGCAGCATCTGCAAGAGTTCAGTTTGAGCGCCGAGCCGCTGAAATCGTTGTTCCCGGAAGGCGGTTT‑TGCCCTGCACGGTAACTCAAAAATCTTTGGTGCC‑GTGGGCATTTCTTT‑AGCGATGTAC‑TTCACTG‑CCGCACCGGAAAATCGGGTAAAA‑GTGGCGGG‑C‑GTG‑CTGGT‑TTCCCCCACCC‑TAAACCCCCCT‑CT‑GGTGGGGAAT‑ACCGAACCGcccg                                                                                                                                                                                                         >  1:1733469/1‑236 (MQ=255)                                       cATTCAGATGTACTGGGCGCAGCATCTGCAAGAGTTCAGTTTGAGCGCCGAGCCGCTGAAATCGTTGTTCCCGGAAGGCGGTTT‑TGCCCTGCACGGTAACTCAAAAATCTTTGGTGCC‑GTGGGCATTTCTTT‑AGCGATGTAC‑TTCACTG‑CCGCACCGGAAAATCGGGTAAAA‑GTGGCGGG‑C‑TTG‑CT‑GG‑TTTCCCCCACC‑CTAACCCCCCTGCT‑TGTGGGGAAT‑ACCCAAcccccgggat                                                                                                                                                                                                      >  1:1293316/1‑231 (MQ=255)                                                 tACTGGGCGCAGCATCTGCAAGAGTTCAGTTTGAGCGCCGAGCCGCTGAAATCGTTGTTCCCGGAAGGCGGTTT‑TGCCCTGCACGGTAACTCAAAAATCTTTGGTGCC‑GTGGGCATTTCTTT‑AGCGATGTAC‑TTCACTG‑CCGCACCGGAAAATCGGGTAAAA‑GTGGCGGG‑C‑TTGTCG‑TT‑TTCCCCCAACC‑TAAACC‑CCCTGCT‑GGGGGGGATT‑ACCTAAACGCTGGAATTTattctc                                                                                                                                                                                              >  1:1049148/1‑232 (MQ=255)                                                   cTGGGCGCAGCATCTGCAAGAGTTCAGTTTGAGCGCCGAGCCGCTGAAATCGTTGTTCCCGGAAGGCGGTTT‑TGCCCTGCACGGTAACTCAAAAATCTTTGGTGCC‑GTGGGCATTTCTTT‑AGCGATGTAC‑TTCACTG‑CCGCACCGGAAAATCGGGTAAAA‑GTGGCGGG‑G‑CTT‑CT‑GA‑TTTCCGCCACC‑CTAACCGCCATGCT‑GGGGGGGAGT‑ACCGAACCGCTGGGAGTTACCTTTCTGt                                                                                                                                                                                          >  1:700674/1‑239 (MQ=255)                                                   cTGGGCGCAGCATCTGCAAGAGTTCAGTTTGAGCGCCGAGCCGCTGAAATCGTTGTTCCCGGAAGGCGGTTT‑TGCCCTGCACGGTAACTCAAAAATCTTTGGTGCC‑GTGGGCATTTCTTT‑AGCGATGTAC‑TTCACTG‑CCGCACCGGAAAATCGGGTAAAA‑GTGGCGGG‑C‑TTG‑CT‑GA‑TTCCCCCCACC‑CTAAACCCCCTCCC‑GGGGGGGATT‑ACCCAACCCCTTGAATTttccttcccct                                                                                                                                                                                          >  1:1853613/1‑236 (MQ=255)                                                   cTGGGCGCAGCATCTGCAAGAGTTCAGTTTGAGCGCCGAGCCGCTGAAATCGTTGTTCCCGGAAGGCGGTTT‑TGCCCTGCACGGTAACTCAAAAATCTTTGGTGCC‑GTGGGCATTTCTTT‑AGCGATGTAC‑TTCACTG‑CCGCACCGGAAAATCGGGTAAAA‑GTGGCGGG‑C‑GTG‑CT‑GG‑TTCCCCCAACC‑ATAAACGCCCCGCT‑GGGGGGGAGT‑TCACAAACGCCGGGAATTATCTTCCTGt                                                                                                                                                                                          >  1:1084808/1‑239 (MQ=255)                                                      ggCGCAGCATCTGCAAGAGTTCAGTTTGAGCGCCGAGCCGCTGAAATCGTTGTTCCCGGAAGGCGGTTT‑TGCCCTGCACGGTAACTCAAAAATCTTTGGTGCC‑GTGGGCATTTCTTT‑AGCGATGTAC‑TTCACTG‑CCGCACCGGAAAATCGGGTAAAA‑GTGGCGGG‑C‑GTG‑CTGGT‑TTCCCCCAACC‑TAAACC‑CCCTGCC‑GGGGGGGATTAACCAAACCCCGGGATTTTTcccgtctcttt                                                                                                                                                                                        >  1:1847960/1‑230 (MQ=255)                                                         gcagaaTCTGCAAGAGTTAGGTTTGAGCGCCGGGCCGCTAAAATCGTTGTTCCCGGAAGGCGGTTT‑TCCCCTGCAGGGTAACTAAAAAATCTTTGGTGCG‑GTGGGCTTTTCTTT‑AGCGATGTAT‑TTCACTG‑CCGCACGGAAAAATGGGGTAAAA‑GGGGGGGG‑T‑TTG‑CT‑GT‑TCCCCGCACCC‑TTAACCGCCATGCT‑GGTGGGAATT‑ACCGAACCGCTGGAATTTACCTTCCTGTTCAttt                                                                                                                                                                                    <  1:1722305/239‑1 (MQ=255)                                                           agcaTCTGCAAGAGTTCAGTTTGAGCGCCGAGCCGCTGAAATCGTTGTTCCCGGAAGGCGGTTT‑TGCCCTGCACGGTAACTCAAAAATCTTTGGTGCC‑GTGGGCATTTCTTT‑AGCGATGTAC‑TTCACTG‑CCGCACCGGAAAATCGGGTAAAA‑GTGGCGGG‑G‑TTG‑CT‑GA‑TGTCCCCCACC‑CTAACCCCCCTCGC‑GGGGGGGGTT‑TTCGCACCCCTGGAAATTTCCCTCCCCTTCAtttt                                                                                                                                                                                   >  2:769514/1‑237 (MQ=255)                                                           agcaTCTGCAAGAGTTCAGTTTGAGCGCCGAGCCGCTGAAATCGTTGTTCCCGGAAGGCGGTTT‑TGCCCTGCACGGTAACTCAAAAATCTTTGGTGCC‑GTGGGCATTTCTTT‑AGCGATGTAC‑TTCACTG‑CCGCACCGGAAAATCGGGTAAAA‑GTGGCGGG‑CTTTG‑CT‑‑A‑TTCCCGCCACC‑CTTACCGCCCCGCC‑GGTGGGGATT‑ACCCAACCGCCGGGAATTTCCTTCcctttcctttt                                                                                                                                                                                   >  1:77431/1‑237 (MQ=255)                                                                       agTTCAGTTTGAGCGCCGAGCCGCTGAAATCGTTGTTCCCGGAAGGCGGTTT‑TGCCCTGCACGGTAACTCAAAAATCTTTGGTGCC‑GTGGGCATTTCTTT‑AGCGATGTAC‑TTCACTG‑CCGCACCGGAAAATCGGGTAAAA‑GTGGCGGG‑C‑TTG‑CT‑GA‑TTTCCGCCACC‑TTAACCCCCATGCT‑GGGGGGGATT‑ACCGAACCGCTGGGATTTACCTTCCTGTTCATTTTACCGTTGCTGttt                                                                                                                                                                      >  1:1368718/1‑239 (MQ=255)                                                                         ttCAGTTTGAGCGCCGAGCCGCTGAAATCGTTGTTCCCGGAAGGCGGTTT‑TGCCCTGCACGGTAACTCAAAAATCTTTGGTGCC‑GTGGGCATTTCTTT‑AGCGATGTAC‑TTCACTG‑CCGCACCGGAAAATCGGGTAAAA‑GTGGCGGG‑C‑TTG‑CT‑GA‑TTCCCCCCACC‑TTAACCCCCCCGCT‑GGGGGGGATT‑ATCGAACCGCTGGAATTTACCTTCCTTTTCCTTTCACCGTTGCTGTTTGc                                                                                                                                                                    >  1:1066947/1‑239 (MQ=255)                                                                           cAGTTTGAGCGCCGAGCCGCTGAAATCGTTGTTCCCGGAAGGCGGTTT‑TGCCCTGCACGGTAACTCAAAAATCTTTGGTGCC‑GTGGGCATTTCTTT‑AGCGATGTAC‑TTCACTG‑CCGCACCGGAAAATCGGGTAAAA‑GTGGCGGG‑G‑CTG‑GT‑TG‑TTCCCCCCACC‑TTAACCCCCCAGCTGGGGGGGAGTT‑ACCAAACC‑CCGGGATTTTTCTTCCTGTTCATTTTTCCCTTGTTGTTTGTgg                                                                                                                                                                  >  1:1019850/1‑239 (MQ=255)                                                                              tttGAGCGCCGAGCCGCTGAAATCGTTGTTCCCGGAAGGCGGTTT‑TGCCCTGCACGGTAACTCAAAAATCTTTGGTGCC‑GTGGGCATTTCTTT‑AGCGATGTAC‑TTCACTG‑CCGCACCGGAAAATCGGGTAAAA‑GTGGCGGG‑C‑GTG‑GT‑GG‑TTCCCGCCACC‑CTAAACGCCATGTT‑GGTGGGGGTT‑ATCCCTCCGCCGGGATTTTACTTCCTGTTCATTTTACCGCTGCTGTTTGTGGtgc                                                                                                                                                               >  2:1594096/1‑237 (MQ=255)                                                                                           cggctaaattggttgtgCCGGGAGGGCGGTTT‑CGCCCGGCAGGGTACCTAAAAATTTTTTGGGGCG‑GGGGGCTTTTTTTA‑GGCGTGTTAC‑TTCCCTG‑CCGCCCCGGAAAATGGGGTAAAG‑GGGGGGGG‑C‑TTG‑CT‑GT‑TCCCCGAACCC‑TTACCCGCCATGCT‑GGTGGGAATT‑ACCGAACCGCTGGAATTTACCTTCCTGTTCATTTCACCGTTGCTGTTTGCGGTACACGCCGTGCTgg                                                                                                                                                   <  2:77431/227‑1 (MQ=255)                                                                                            ggataaatgtgttgtcCCCGGAGGGGGGTTT‑T‑CCCTGGCCGGTAAACTAAAAATTGTTGGGGCC‑GGGGGGATTTTTTTGGGTGTTTTCC‑TTCCCGC‑CCGCCCCGAAAAATGGGGAAAAG‑GGGGGGGG‑T‑TTG‑TT‑GT‑TCCCCGCACCT‑TTAACCGCCATGCT‑GGTGGGAATT‑ACCGAACCGCTGGAATTTACCTTCCTGTTCATTTCACCGTTGCTGTTTGCGGTACACGCCGTGCTggc                                                                                                                                                  <  2:1505299/228‑1 (MQ=255)                                                                                                       ttgttCCGGGAGGGCGTTTT‑TCCCCGGCACGGTAACTAAAAATTCTTTGGGGCC‑GGGGGTTTTTTTTT‑AGCGATTTCC‑TTCACTG‑CCGCCCGGAAAAATGGGGAAAAG‑GGGGGGGG‑T‑TTG‑CT‑GA‑TCCCCGCACCT‑TTAACCGCCATGCT‑GGTGGGAATT‑ACCGAACCGCTGGAATTTACCTTCCTGTTCATTTCACCGTTGCTGTTTGCGGTACACGCCGTGCTGGCGGCCTCAATGTc                                                                                                                                      <  1:1655449/239‑1 (MQ=255)                                                                                                       gttttCCCGGGAGGGGTTTTGTGCCTTGC‑CGGTAAACTAAAAATTTTTGGGGCC‑GTGGGGTTTTTTTT‑GGGGGTTTTCTTTCTCTG‑CCGCCCCGAAAACTGGGATAAAG‑GGGGGGGT‑T‑TTG‑CT‑TT‑TCCCCGCAACT‑TTACCCGCCATGCT‑GGTGGGAATT‑ACCGAACCGCTGGAATTTACCTTCCTGTTCATTTCACCGTTGCTGTTTGCGGTACACGCCGTGCTGGCGGCCTCAATGt                                                                                                                                       <  2:1757653/236‑1 (MQ=255)                                                                                                          tccccGGAGGGGGGTTT‑TGCCCTGCAGGGTAACTAAAAATTTTTTGGTGCC‑GTGGGTTTTTTTTT‑AGCGATGTAC‑TTCCCTG‑CCCCCCCGAAAAATGGGGAAAAA‑GTGGCGGG‑T‑TTG‑CT‑GA‑TCCCCGCACCC‑TTACCCGCCATGCT‑GGTGGGAATT‑ACCGAACCGCTGGAATTTACCTTCCTGTTCATTTCACCGTTGCTGTTTGCGGTACACGCCGTGCTGGCGGCCTCAATGTCGAc                                                                                                                                   <  1:313552/237‑1 (MQ=255)                                                                                                            ccgggAGGGGGGTTT‑TCCCCTGCAGGTAAACTAAAAAATCTTTGGTGCC‑GGGGGCTTTTTTTT‑AGCGATGTCC‑TTCCCTG‑CCCCCCCGAAAAATGGGGAAAAA‑GGGGGGGG‑T‑TTG‑CT‑GA‑TCCCCGCACCT‑TAAACCGCCATGCT‑GGTGGGAATT‑ACCGAACCGCTGGAATTTACCTTCCTGTTCATTTCACCGTTGCTGTTTGCGGTACACGCCGTGCTGGCGGCCTCAATGTCGACCg                                                                                                                                 <  1:464106/236‑1 (MQ=255)                                                                                                                    ggttttt‑tCCCCTGCACGTAAACTAAAAATTTTTTGGGGCT‑GT‑GGCTTTTTTTT‑GGGGGTTTTC‑TTTCCCG‑GCGGCCCGGAAAATGGGGGAAAA‑GGGGGGGGTC‑TTG‑CT‑GT‑TCCCCACACCC‑TTACCCGCCATGCT‑GGTGGGAATT‑ACCGAACCGCTGGAATTTACCTTCCTGTTCATTTCACCGTTGCTGTTTGCGGTACACGCCGTGCTGGCGGCCTCAATGTCGACCGTAATGTa                                                                                                                          <  2:1756710/234‑1 (MQ=255)                                                                                                                        ttt‑tGGCCCTCACCGAAACCCAAAATAGTTTGGGGCCAGTGGGTTTTTCTTG‑AGGGATGTCC‑TCCCCCG‑CCGCACGGAAAAACGGGGTAAAG‑GGGGGAGG‑C‑TTG‑AT‑GT‑TCCCCGAATCT‑TCAACCGCCATGCT‑GGTGGGAATT‑ACCGACCCGCTGGAATTTCCCTTCCTGTTCATTTCACCGTTGCTGTTTGAGGTACACGCCGTGCTGGCGGCCTCAATGTCGACCGTAATGTAtct                                                                                                                       <  2:162115/238‑1 (MQ=255)                                                                                                                        gtttgccccGGCGCGTTAACTAAAAAATTTTGGGTGCG‑GTGGGTTTTTTTTT‑GGCGTTGTAC‑TCCCCGC‑CCGCCCCGAAAAATGGGGAAAAG‑GTGGGGGG‑T‑TTG‑TT‑GT‑TCCCCGCACCT‑TTACCCGCCATGCT‑GGTGGGAATT‑ACCGAACCGCTGGAATTTACCTTCCTGTTCATTTCACCGTTGCTGTTTGCGGTACACGCCGTGCTGGCGGCCTCAATGTCGACCGTAATGTAtct                                                                                                                       <  1:777373/232‑1 (MQ=255)                                                                                                                                          aCTCAAAAATTTTGGGTGCG‑GTGGGTTTTTCTTG‑AGCGATGTAC‑TTCCCTG‑CCGCACGGAAAAATGGGGTAAAG‑GGGGCGGG‑C‑TTG‑CT‑GA‑TTCCCGCAACC‑TTAACCGCCATGCT‑GGTGGGAATT‑ACCGAACCGCTGGAATTTACCTTCCTGTTCATTTCACCGTTGCTGTTTGCGGTACACGCCGTGCTGGCGGCCTCAATGTCGACCGTAATGTATCTCTTTGGTGTGGTGGGCAAc                                                                                                    <  1:299324/239‑1 (MQ=255)                                                                                                                                                 aaTTTTTGGGGCG‑GTGGGTTTTTTTTT‑AGCGTTGTCC‑TCCCCTC‑CCCCCCCGAAAAATCGGGAAAAG‑GGGGCGGG‑T‑TTG‑CT‑GA‑TCCCCCCACCT‑TTACCCGCCATGCT‑GGTGGGAATT‑ACCGAACCGCTGGAATTTACCTTCCTGTTCATTTCACCGTTGCTGTTTGCGGTACACGCCGTGCTGGCGGCCTCAATGTCGACCGTAATGTATCTCTTTGGTGTGGTGGGCAACATGGGCg                                                                                             <  1:17282/239‑1 (MQ=255)                                                                                                                                                      ttGGTGCC‑GTGGGCATTTCTTT‑AGCGATGTAC‑TTCACTG‑CCGCACCGGAAAATCGGGTAAAA‑GTGGCGGG‑C‑TTT‑CT‑GA‑TTCCCGCAACC‑TTTAACGCCATGCT‑GGTGGGGATT‑ACCCAACCCCTGGAATTTAACTTCCTTTTCATTTTACCGTTGCTGTTTGCGGGACACGCCGTGCTGGCGGCCTCAATGTCCACCGTAATGTATCTCTTTTGTGTGGTGGGGAACCTGGGCGggggg                                                                                        >  1:1157193/1‑238 (MQ=255)                                                                                                                                                            gcggggggtttttCTTG‑AGCGTTGT‑C‑TTCCCCGCCCGCCCCGGAAAATGGGGTAAAG‑GGGGCGGG‑T‑TTG‑TT‑TA‑TTCCCGCACCT‑TTAACCGCCATGCT‑GGTGGGAATT‑ACCGAACCGCTGGAATTTACCTTCCTGTTCATTTCACCGTTGCTGTTTGCGGTACACGCCGTGCTGGCGGCCTCAATGTCGACCGTAATGTATCTCTTTGGTGTGGTGGGCAACATGGGCGGAGGTCTGa                                                                                    <  1:122977/228‑1 (MQ=255)                                                                                                                                                            cc‑GTGGGATTTTTTTT‑GGGGATGTCC‑TCCCCTC‑CCCCCCGGAAAAATCGGGAAAAG‑GTGGCGGG‑T‑TTG‑CT‑GA‑TCCCCGCAACC‑TTAACCGCCATGCT‑GGTGGGAATT‑ACCGAACCGCTGGAATTTACCTTCCTGTTCATTTCACCGTTGCTGTTTGCGGTACACGCCGTGCTGGCGGCCTCAATGTCGACCGTAATGTATCTCTTtg                                                                                                                  <  1:355937/207‑1 (MQ=255)                                                                                                                                                                                                 ccGCACCGGAAAATCGGGTAAAA‑GTGGCGGG‑G‑TTG‑CT‑TATTTCCCCCACCC‑TCAACC‑CCCTGCT‑TGGGGGGAGT‑ACCGCACCCCCGGGAATTATCTTCCCGTTCTTTTTACCCTTTCTGTTTTCGGGACACGCCGGGCTGGCGGGCTCACTGTCGACCGGTATGTGTCTCTTTGTGGTGGGGGGGTGCATGGGGGGGGGGCCGGGTGAGTAGGTTTTTTCGCGCAACTgtggtccg                                                    >  2:815978/1‑229 (MQ=255)                                                                                                                                                                                                              cgggatAAAG‑GTGGGGGG‑C‑TTG‑CT‑TC‑TTCCCGCACCC‑TAAACCGCCATGCT‑GGTGGGAATT‑ACCGAACCGCTGGAATTTACCTTCCTGTTCATTTCACCGTTGCTGTTTGCGGTACACGCCGTGCTGGCGGCCTCAATGTCGACCGTAATGTATCTCTTTGGTGTGGTGGGCAACATGGGCGGAGGTCTGATTGACTAGGTTTTACCGCAAAACTGGATCCCGATGTTCAGCAACCAcg                                    <  2:860311/234‑1 (MQ=255)                                                                                                                                                                                                               ggggTAAAG‑GGGGCGGG‑T‑TGG‑CT‑AA‑TTCCCGAAACT‑TTAACCGCCATGCT‑GGTGGGAATT‑ACCGAACCGCTGGAATTTACCTTCCTGTTCATTTCACCGTTGCTGTTTGCGGTACACGCCGTGCTGGCGGCCTCAATGTCGACCGTAATGTATCTCTTTGGTGTGGTGGGCAACATGGGCGGAGGTCTGATTGACTAGGTTTTACCGCAAAACTGGATCCCGATGTTCAGCAACCAcgc                                   <  1:1183180/238‑1 (MQ=255)                                                                                                                                                                                                               cGGGTAAAATGTGG‑GGG‑G‑TTG‑CT‑GG‑TTCCCCGAACCTTTAACCGCCATGCT‑GGTGGGAATT‑ACCGAACCGCTGGAATTTACCTTCCTGTTCATTTCACCGTTGCTGTTTGCGGTACACGCCGTGCTGGCGGCCTCAATGTCGACCGTAATGTATCTCTTTGGTGTGGTGGGCAACATGGGCGGAGGTCTGATTGACTAGGTTTTACCGCAAAACTGGATCCCGATGTTCAGCAACCAcg                                    <  2:426232/239‑1 (MQ=255)                                                                                                                                                                                                                gggTAAAA‑GGGGGGGG‑G‑TTT‑GC‑GT‑TTCCCCGAACCTTTAACCGCCATGCT‑GGTGGGAATT‑ACCGAACCGCTGGAATTTACCTTCCTGTTCATTTCACCGTTGCTGTTTGCGGTACACGCCGTGCTGGCGGCCTCAATGTCGACCGTAATGTATCTCTTTGGTGTGGTGGGCAACATGGGCGGAGGTCTGATTGACTAGGTTTTACCGCAAAACTGGATCCCGATGTTCAGCAACCAcg                                    <  2:1288643/238‑1 (MQ=255)                                                                                                                                                                                                                                      g‑CT‑GA‑TTCCCGCAACC‑TTAACCGCCATGCT‑GGTGGGAATT‑ACCGAACCGCTGGAATTTACCTTCCTGTTCATTTCACCGTTGCTGTTTGCGGTACACGCCGTGCTGGCGGCCTCAATGTCGACCGTAATGTATCTCTTTGGTGTGGTGGGCAACATGGGCGGAGGTCTGATTGACTAGGTTTTACCGCAAAACTGGATCCCGATGTTCAGCAACCACGCGGATATGATGCTGACCCaaa               >  2:1203566/1‑239 (MQ=255)                                                                                                                                                                                                                                         t‑GA‑TTCCCCCACCC‑TTAACCGCCATGCT‑GGTGGGAATT‑ACCGAACCGCTGGAATTTACCTTCCTGTTCATTTCACCGTTGCTGTTTGCGGTACACGCCGTGCTGGCGGCCTCAATGTCGACCGTAATGTATCTCTTTGGTGTGGTGGGCAACATGGGCGGAGGTCTGATTGACTAGGTTTTACCGCAAAACTGGATCCCGATGTTCAGCAACCACGCGGATATGATGCTGACCCAAAt              <  1:1692672/238‑1 (MQ=255)                                                                                                                                                                                                                                                   gCAACC‑TTAACCGCCATGCT‑GGTGGGAATT‑ACCGAACCGCTGGAATTTACCTTCCTGTTCATTTCACCGTTGCTGTTTGCGGTACACGCCGTGCTGGCGGCCTCAATGTCGACCGTAATGTATCTCTTTGGTGTGGTGGGCAACATGGGCGGAGGTCTGATTGACTAGGTTTGACCGCAAAACTGGATCCCGATGTTCAGCAACCACGCGGATATGATGCTGACCCAAATCGCCATTgg     >  1:1633706/1‑239 (MQ=255)                                                                                                                                                                                                                                                    cAACC‑TTAACCGCCATGCT‑GGTGGGAATT‑ACCGAACCGCTGGAATTTACCTTCCTGTTCATTTCACCGTTGCTGTTTGCGGTACACGCCGTGCTGGCGGCCTCAATGTCGACCGTAATGTATCTCTTTGGTGTGGTGGGCAACATGGGCGGAGGTCTGATTGACTAGGTTTTACCGCAAAACTGGATCCCGAt                                                  >  1:833698/1‑193 (MQ=255)                                                                                                                                                                                                                                                    cAACC‑TTAACCGCCATGCT‑GGTGGGAATT‑ACCGAACCGCTGGAATTTACCTTCCTGTTCATTTCACCGTTGCTGTTTGCGGTACACGCCGTGCTGGCGGCCTCAATGTCGACCGTAATGTATCTCTTTGGTGTGGTGGGCAACATGGGCGGAGGTCTGATTGACTAGGTTTTACCGCAAAACTGGATCCCGATGTTCAGCAACCACGCGGATATGATGCTGACCCAAATCGCCATTggg    <  1:100227/239‑1 (MQ=255)                                                                                                                                                                                                                                                      aCC‑TTACCCGCCATGCT‑GGTGGGAATT‑ACCGAACCGCTGGAATTTACCTTCCTGTTCATTTCACCGTTGCTGTTTGCGGTACACGCCGTGCTGGCGGCCTCAATGTCGACCGTAATGTATCTCTTTGGTGTGGTGGGCAACATGGGCGGAGGTCTGATTGACTAGGTTTTACCGCAAAACTGGATCCCGATGTTCAGCAACCACGCGGATATGATGCTGACCCAAATCGCCATTGGGtt  <  1:428306/239‑1 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                  GGACAGTTTATCTTTGGTCCGGCAGCTGTTGAAGGCGGCATTCAGATGTACTGGGCGCAGCATCTGCAAGAGTTCAGTTTGAGCGCCGAGCCGCTGAAATCGTTGTTCCCGGAAGGCGGTTT‑TGCCCTGCACGGTAACTCAAAAATCTTTGGTGCC‑GTGGGCATTTCTTT‑AGCGATGTAC‑TTCACTG‑CCGCACCGGAAAATCGGGTAAAA‑GTGGCGGG‑C‑TTG‑CT‑GA‑TTCCCGCAACC‑TTAACCGCCATGCT‑GGTGGGAATT‑ACCGAACCGCTGGAATTTACCTTCCTGTTCATTTCACCGTTGCTGTTTGCGGTACACGCCGTGCTGGCGGCCTCAATGTCGACCGTAATGTATCTCTTTGGTGTGGTGGGCAACATGGGCGGAGGTCTGATTGACTAGGTTTTACCGCAAAACTGGATCCCGATGTTCAGCAACCACGCGGATATGATGCTGACCCAAATCGCCATTGGGTT  >  NC_007779/3777520‑3777992

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