BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_007779	1,555,708	T→G	19.8%	D559A (GAC→GCC)	Y75_RS07805 ←	malate dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007779	1,555,708	0	T	G	19.8%	111.7 / 3.2	70	D559A (GAC→GCC)	Y75_RS07805	malate dehydrogenase
Reads supporting (aligned to +/- strand): ref base T (22/29); new base G (13/2); total (36/34)
Fisher's exact test for biased strand distribution p-value = 3.33e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTTCGCCATGTTTCAGTGAGCGCAGTGTTTTATACGTAACGTCAACATGATGATCCTTCGTAACAACTGCAGCCTTCATAGTTCCTCCTTTTCGGATGATGTTCTGCATAGCAGGTGAGGCAAATGAGATTTATTCG‑CCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAG‑CCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGC‑GGTATTCGGC‑TTG‑CCAGAAATTATCGTCAATGGCCTGTTGCA‑GGGCTTCGGCAGA‑GGTTTTCACCGCCACGCCTTGCTGCTGCGCCATTTTGCCAACCGC‑AAACGCAATTGCGCGGGAGACTTTCTGAATATCTTTCAGTTCCGGCAGTA‑CCATACCTTCGCCGTTCAGCACCAATGGTGAATACTGCGCCAGCGTTTCACTTGCCGACATCAGCAT  >  NC_007779/1555473‑1555928

gCTTCGCCATGTTTCAGTGAGCGCAGTGTTTTATACGTAACGTCAACATGATGATCCTTCGTAACAACTGCAGCCTTCATAGTTCCTCCTTTTCGGATGATGTTCTGCATAGCAGGTGAGGCAAATGAGATTTATTCG‑CCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAG‑CCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTcg                                                                                                                                                                                                                                   <  1:1136959/238‑1 (MQ=255)

gCTTCGCCATGTTTCAGTGAGCGCAGTGTTTTATACGTAACGTCAACATGATGATCCTTCGTAACAACTGCAGCCTTCATAGTTCCTCCTTTTCGGATGATGTTCTGCATAGCAGGTGAGGCAAATGAGATTTATTCG‑CCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAG‑CCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTcg                                                                                                                                                                                                                                   <  1:1786680/238‑1 (MQ=255)

           tttCAGTGAGCGCAGTGTTTTATACGTAACGTCAACATGATGATCCTTCGTAACAACTGCAGCCTTCATAGTTCCTCCTTTTCGGATGATGTTCTGCATAGCAGGTGAGGCAAATGAGATTTATTCG‑CCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAG‑CCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGC‑GGGTTTCGGC‑t                                                                                                                                                                                                                     >  1:391325/1‑239 (MQ=255)

           tttCAGTGAGCGCAGTGTTTTATACGTAACGTCAACATGATGATCCTTCGTAACAACTGCAGCCTTCATAGTTCCTCCTTTTCGGATGATGTTCTGCATAGCAGGTGAGGCAAATGAGATTTATTCG‑CCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAG‑CCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGGTGC‑GGTATTCGGC‑t                                                                                                                                                                                                                     >  1:489745/1‑239 (MQ=255)

           tttCAGTGAGCGCAGTGTTTTATACGTAACGTCAACATGATGATCCTTCGTAACAACTGCAGAGTTCATAGTTCCTCCTTTTCGGATGATGTTCTGCATAGCAGGTGAGGCAAATGAGATTTATTCG‑CCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAG‑CGCGGTAGCCTTCACAACCGGGCGCAGGCGTAGATGGAGGTACGGCGGTAAGCGC‑CGGAATCGGcg                                                                                                                                                                                                                      >  2:1508079/1‑238 (MQ=255)

                  tggcgcAGTGTTTTATACGGAACGTCAACATGGTGATCCTTCGTAACAACTGCAGCCTTCATAAGTCCTCCTTTTCGGATGATGTTCGGCATAGCAGGTGAGGCAAATGAGATTTATTCG‑CCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAG‑CCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGC‑GGTATTCGGC‑TTG‑cc                                                                                                                                                                                                                <  2:1908617/234‑1 (MQ=255)

                  gAGCGCAGTGTTTTATACGTAACTTCAACATGATGATCCTTCGTAACAACTGCAGCCTTCATAGTTCCTCCTTTTCGGATGATGTTCTGCATAGCAGGTGAGGCAAATGAGATTTATTCG‑CCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAG‑GCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGC‑GGTATTCGGC‑TTG‑CCAg                                                                                                                                                                                                              <  1:1968438/238‑1 (MQ=255)

                        ggggtTTTATACGTAACGCTAACATGGTGACCCTTCGTAACAACTGCAGCCTTGAGAGTTCCTCCTTTTCGCATAATGTTCTGCATAGCAGTGGAGACAAATGGGATTTATTCG‑CCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAG‑CCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGC‑GGTATTCGGC‑TTG‑CCAGAAATTAt                                                                                                                                                                                                       <  2:1201758/236‑1 (MQ=255)

                           gtTTTTTACGTAACGTCAACATTATGTTCCTTCGTAACAACTGCAGCCTTCATAGTTCCTCCTTTTCGGATGATGTTCTGCATAGCAGGTGAGGCAAATGAGATTTATTCG‑CCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAG‑CCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGC‑GGTATTCGGC‑TTG‑CCAGAAATTATCg                                                                                                                                                                                                     <  2:489745/238‑1 (MQ=255)

                             tttATCCGTATCGTCAACATGTTGCTCGTTCGTAACAACTGCAGCCTTCATAGTTCGTCCTTTTCGGATGATGTTCTGCATAGCAGGGGAGGCAAAGGAGATTTATTCG‑CCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAG‑CCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGC‑GGTATTCGGC‑TTG‑CCAGAAATTATCGTCa                                                                                                                                                                                                  <  2:1724700/239‑1 (MQ=255)

                             tttATACGTAACGTCAACATGATGATCCTTCGGAACAACTGCAGCCTTCATAGTTCCTCCTTTTCGGATGATGTTCTGCATAGCAGGTGAGGCAAATGAGATTTATTCG‑CCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAG‑CCCGGTAGCCGTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGC‑GGTATTCGGC‑TTG‑CCAGAAATTATCGTCa                                                                                                                                                                                                  <  1:1874485/239‑1 (MQ=255)

                               tataCGTAACGTCAACATGATGATCCTTCGTAACAACTGCAGCCTTCATAGTTCCTCCTTTTCGGATGATGTTCTGCATAGCAGGTGAGGCAAATGAGATTTATTCG‑CCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAG‑CCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGCGGGTATCGGGC‑TT‑‑CCAGGAAATATCgctcat                                                                                                                                                                                                >  2:266751/1‑234 (MQ=255)

                                   cttAACGTAAACATGAGGATCCTTCGTAACAACTGCAGCCTTCATGGTTCCTCCTTTTCGGATGATGTTCTGCATGGCAGGTAGGGAAAATGAGATTTATTCG‑CCACTCCCCAGTATGGATGAGATCTGAAAAAGAGAGAGGGAAATAG‑CCCGGTAGCCTCCCCTACCCGGCGCCGGCTTAGATGGATGTACGGCGTGAGTCGC‑GTTATCCGGC‑TTG‑CCAGAAATTATCGTCAATCGca                                                                                                                                                                                            <  2:2094314/237‑2 (MQ=255)

                                      aaCGTCAACATGATGATCCTTCGTAACAACTGCAGCCTTCATAGTTCCTCCTTTTCGGATGATGTTCTGCATAGCAGGTGAGGCAAATGAGATTTATTCG‑CCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAG‑CCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGCTGGAGGTACGGCGGTAGTCCG‑GGGATACGGC‑TTG‑CCACAAATTCTTGTccgtgggccgt                                                                                                                                                                                         >  1:697357/1‑229 (MQ=255)

                                      aaCGTCAACATGATGATCCTTCGTAACAACTGCAGCCTTCATAGTTCCTCCTTTTCGGATGATGTTCTGCATAGCAGGTGAGGCAAATGAGATTTATTCG‑CCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAG‑CCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTGGGCGC‑GGGATTTGGC‑TTG‑CCCGAAATTATTGTCAATGGCCCGt                                                                                                                                                                                         >  1:96172/1‑239 (MQ=255)

                                      aaCGTCAACATGATGATCCTTCGTAACAACTGCAGCCTTCATAGTTCCTCCTTTTCGGATGATGTTCTGCATAGCAGGTGAGGCAAATGAGATTTATTCG‑CCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAG‑CCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGC‑GGGATTCGGC‑TTG‑CCAGAAATTATTGTCCATGGCCTGt                                                                                                                                                                                         >  1:1478000/1‑239 (MQ=255)

                                      aaCGTCAACATGATGATCCTTCGTAACAACTGCAGCCTTCATAGTTCCTCCTTTTCGGATGATGTTCTGCATAGCAGGTGAGGCAAATGAGATTTATTCG‑CCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAG‑CCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGGCGC‑GGGATTCGGC‑TTG‑CCAGAAATTATTGTCACTGGCCTgg                                                                                                                                                                                         >  1:1562826/1‑238 (MQ=255)

                                      aaCGTCAACATGATGATCCTTCGTAACAACTGCAGCCTTCATAGTTCCTCCTTTTCGGATGATGTTCTGCATAGCAGGTGAGGCAAATGAGATTTATTCG‑CCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAG‑CCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGGTGGTGC‑GGGGATCGGC‑TTG‑CCAGAAAATATATGCAATgggccgt                                                                                                                                                                                         >  1:1394154/1‑234 (MQ=255)

                                      aaCGTCAACATGATGATCCTTCGTAACAACTGCAGCCTTCATAGTTCCTCCTTTTCGGATGATGTTCTGCATAGCAGGTGAGGCAAATGAGATTTATTCG‑CCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAG‑CCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGAACGGCGGGAGCTCC‑GGGATTCGGG‑TCG‑TCCGAAAATAACGTGACTgggccgt                                                                                                                                                                                         >  2:356953/1‑234 (MQ=255)

                                            aaCATGATGATCCTTCGTAACAACTGCAGCCTTCATAGTTCCTCCTTTTCGGATGATGTTCTGCATAGCAGGTGAGGCAAATGAGATTTATTCG‑CCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAG‑CCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGC‑GGGATTCGGG‑TTT‑GCAGAAATTATCGTCAAAGGGCTGTTTCA‑gg                                                                                                                                                                                  >  1:1007200/1‑239 (MQ=255)

                                                   tgatCCTTCGTAACAACTGCAGCCTTCATAGTTCCTCCTTTTCGGATGATGTTCTGCATAGCAGGTGAGGCAAATGAGATTTATTCG‑CCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAG‑CCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGGAGTCGC‑GGGTATCGGC‑CTG‑GCCAAAATTATAGTTAAAGGCCTGTTgtcggggatcg                                                                                                                                                                             >  1:1571768/1‑228 (MQ=255)

                                                      tCCTTCGTAACAACTGCAGCCTTCATAGTTCCTCCTTTTCGGATGATGTTCTGCATAGCAGGTGAGGCAAATGAGATTTATTCG‑CCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAG‑CCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGC‑GGTATTCGGC‑TTG‑CCAGAAATTATCGTCAATGGCCTGTTGCA‑GGGCTTCGGCa                                                                                                                                                                         <  1:939060/238‑1 (MQ=255)

                                                      tCCTTCGTAACAACTGCAGCCTTCATAGTTCCTCCTTTTCGGATGATGTTCTGCATAGCAGGTGAGGCAAATGAGATTTATTCG‑CCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAG‑CCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGGAGTCCC‑GGGATTCGGC‑CTG‑GCAGAAATTTTTGTCAATAGCCTGGTGGC‑GGGCTTCGGCag                                                                                                                                                                        >  1:1547703/1‑239 (MQ=255)

                                                       ccTTCGTAACAACTGCAGCCTTCATAGTTCCTCCTTTTCGGATGATGTTCTGCATAGCAGGTGAGGCAAATGAGATTTATTCG‑CCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAG‑CCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTGGC‑GGGATTCGGC‑CTG‑GCCGCAAATATCGTCAATGGGCTGGTGCA‑GGGCGTGGGCaga                                                                                                                                                                       >  1:1055036/1‑239 (MQ=255)

                                                        cTTCGTAACAACTGCAGCCTTCATAGTTCCTCCTTTTCGGATGATGTTCTGCATAGCAGGAGAGGCAAATGAGATTTATTCG‑CCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAG‑CCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGGAGGCGC‑GGGGTTCGGG‑GTT‑GCAAACATTATTGTTAATGGCCTGTTGGC‑GGGCCTCCGCcga‑g                                                                                                                                                                     >  1:96326/1‑239 (MQ=255)

                                                         ttCGTAACAACTGCAGCCTTCATAGTTCCTCCTTTTCGGATGATGTTCTGCATAGCAGGTGAGGCAAATGAGATTTATTCG‑CCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAG‑CCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTATC‑GGGATTTGTC‑TTG‑CCCGAAATTTGCTGTAATGGTCTGTTGGT‑GGGcgcgggcgggag                                                                                                                                                                     >  2:776906/1‑228 (MQ=255)

                                                         ttCGTAACAACTGCAGCCTTCATAGTTCCTCCTTTTCGGATGATGTTCTGCATAGCAGGTGAGGCAAATGAGATTTATTCG‑CCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAG‑CCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGAAGGCGC‑GGTACTCCGG‑TTT‑TCAGGAAATATAGTCCATTGGCCGTTTGA‑GGGCTTCGGGAgg‑gg                                                                                                                                                                    >  1:486100/1‑239 (MQ=255)

                                                           cGTAACAACTGCAGCCTTCATAGTTCCTCCTTTTCGGATGATGTTCTGCATAGCAGGTGAGGCAAATGAGATTTATTCG‑CCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAG‑CCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTcg                                                                                                                                                                                                                                   >  1:536078/1‑179 (MQ=255)

                                                           cGTAACAACTGCAGCCTTCATAGTTCCTCCTTTTCGGATGATGTTCTGCATAGCAGGTGAGGCAAATGAGATTTATTCG‑CCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAG‑CCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTcg                                                                                                                                                                                                                                   <  2:536078/179‑1 (MQ=255)

                                                           cGTAACAACTGCAGCCTTCATAGTTCCTCCTTTTCGGATGATGTTCTGCATAGCAGGTGAGGCAAATGAGATTTATTCG‑CCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAG‑CCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGC‑GGGATTCGGC‑TTG‑CCAGAAAGTATTGGCAATGGGCTGTTGCA‑GGGCTTCGGCAgg‑gggt                                                                                                                                                                  >  1:2085820/1‑237 (MQ=255)

                                                           cGTAACAACTGCAGCCTTCATAGTTCCTCCTTTTCGGATGATGTTCTGCATAGCAGGTGAGGCAAATGAGATTTATTCG‑CCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAG‑CCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGCTGGCGC‑GGGATTTGGG‑TTG‑GCCGAAAATATTCGCCATGGGCTGTTGTC‑GGGCTTCGGGGGA‑gggt                                                                                                                                                                  >  1:2173385/1‑237 (MQ=255)

                                                                      cAGCCTTCATAGTTCCTCCTTTTCGGATGATGTTCTGCATAGCAGGTGAGGCAAATGAGATTTATTCG‑CCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAG‑CCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGGCGC‑GGTATTTGGC‑TTT‑CCCGAAATTATCGTCAAATGTCTGCTGCC‑GGGCGTCGGCGGA‑GGGTTTCatccccc                                                                                                                                                        >  2:2100158/1‑232 (MQ=255)

                                                                        gCCTTCATAGTTCCTCCTTTTCGGATGATGTTCTGCATAGCAGGTGAGGCAAATGAGATTTATTCG‑CCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAG‑CCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTGGTAGC‑GGGGATCGGC‑GTGTCCGGAAATT‑TCGGTAATGGCCTGTTGCA‑GGGCGTCCGGAGG‑GGGTTTCACCGCCAca                                                                                                                                                      >  2:313788/1‑237 (MQ=255)

                                                                        gCCTTCATAGTTCCTCCTTTTCGGATGATGTTCTGCATAGCAGGTGAGGCAAATGAGATTTATTCG‑CCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAG‑CCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAAGCCC‑GGGATTTGGC‑CTG‑CCAGAAATTATTCTTAATGGGCTGTTGGA‑GGGGTTCGGGGGG‑GGGTTTCACCGCCAcc                                                                                                                                                      >  1:435926/1‑237 (MQ=255)

                                                                         ccTTCATAGTTCCTCCTTTTCGGATGATGTTCTGCATAGCAGGTGAGGCAAATGAGATTTATTCG‑CCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAG‑CCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGC‑GGTATTCGGC‑TTG‑CCAGAAATTAt                                                                                                                                                                                                       <  1:1626394/190‑1 (MQ=255)

                                                                             cATAGTTCCTCCTTTTCGGATGATGTTCTGCATAGCAGGTGAGGCAAATGAGATTTATTCG‑CCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAG‑CCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAATCGC‑GGGATTCGGC‑GTG‑CCAGGAATTATCGTCAATGGGCCGGTGCA‑GGGGTTCGGCAGA‑GGGTTTCACCGCCACGCCtt                                                                                                                                                  >  1:1079628/1‑237 (MQ=255)

                                                                                    cctcctTTTCGGATGATGTGCTGTATAGCAGTGGAGGCAAATGATATTTATTCG‑CCCCCACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAG‑CCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGC‑GGTATTCGGC‑TTG‑CCCGAAATTATCGTCAATGGCCTGTTGCA‑GGGCTTCGGCAGA‑GGTTTCCACCGCCACGCCTTGCTGCTgcg                                                                                                                                         <  1:238117/239‑1 (MQ=255)

                                                                                     ctcctTTTCGGATGATGTTCTGCATAGCAGGTGAGGCAAATGAGATTTATTCG‑CCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAG‑CCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGCAGTCGC‑GGTAATTGGCTTTC‑CCAGAAA‑TATCTGCAATGGCCCTTTTCC‑GGGCCTTGGGCCA‑GGTTTTCACCCCCCCCCCCtgcctgcgcg                                                                                                                                         >  2:1936913/1‑231 (MQ=255)

                                                                                          tttCGGATGATGTTCTGCATAGCAGGTGAGGCAAATGAGATTTATTCG‑CCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAG‑CCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGC‑GGGATTCGGC‑TTG‑CCAGAAATTATCGTCAATGGGCTGTTGCA‑GGGGTTCGGCAGA‑GGGTTTCACCGCCACGCCTTGCTGCTGCGCCAttt                                                                                                                                   >  1:663350/1‑239 (MQ=255)

                                                                                                    tgTTCTGAATGGAAGGTGAGGCAAAGGGGATTTATTCG‑CCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGAAAATAG‑CCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGC‑GGTATTCGGC‑TTG‑CCAGAAATTATCGTCAATGGCCTGTTGCA‑GGGCTTCGGCAGA‑GGTTTTCACCGCCACGCCTTGCTGCTGCGCCATTTTGCCAACCGc                                                                                                                         <  1:266751/239‑1 (MQ=255)

                                                                                                              aGCAGGTGAGGCAAATGAGATTTATTCG‑CCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAG‑CCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGC‑GGTATTCGGC‑GTG‑TCAGCAAATATCGTCAATGGCCTGTTGCA‑GGGGTTCGCCAGA‑GGTGTTCACCGCCACGCCTTGCTTCTGCGCCATTTTGCCAACCGC‑AACAGCAAtt                                                                                                              >  1:2232886/1‑239 (MQ=255)

                                                                                                                   gTGAGGCAAATGAGATTTATTCG‑CCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAG‑CCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCAGTAGTCGC‑GGTATTCGGC‑TTG‑CCAGAAATTATCGTCAATGGCCTGTTGCA‑GGGCTTCGGCAGA‑GGTTTTCACCGCCACGCCTTGCTGCTGCGCCATTTTGCCAACCGC‑AAACGCAATTGcgcg                                                                                                         >  1:1250382/1‑239 (MQ=255)

                                                                                                                     gAGGCAAATGAGATTTATTCG‑CCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAG‑CCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGGCGC‑GGGGTTCGGC‑TTT‑ACACAAATTTTCACCATTAGCCCCTTGCT‑GTGTCTCTGGAGG‑AGCTTTTACCAGCCCGCCCTTCCGCGGTGCCCCTTTGCCCACCGC‑ACACACAACTGCGCggg                                                                                                       >  2:1687959/1‑239 (MQ=255)

                                                                                                                     gAGGCAAATGAGATTTATTCG‑CCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAG‑CCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGGAGTCGC‑GGGATTCGGC‑TTG‑CCAGAAATTATTGTCAATGGGCTGTTGCA‑GGGCTTCGGCGGA‑GGTTTTCACCGCCACGCCTTGCTGCTGCGCCCTTTTTCCCACCGC‑AAACGCCCATGCGCggg                                                                                                       >  1:278175/1‑239 (MQ=255)

                                                                                                                        aaaaaaGAGATTTATTCG‑CCACTCCCCATGAGGGATGGAATCAGAAAAAGGGCGGGGGAAATGG‑CGCGTGACCCTTCACTCCGGGGGGCAGGCTAGGGGGGGGGGGCGGGGGTGGCCGG‑GTTTTTCGGT‑TTG‑CAAAAAATTATCGTCAGTGGGCTGTTGCA‑GGGTTTCGGCAGG‑TTTTTTCACCGCCACGCCTTGCTGCTGCGCCATTTTGCCAACCGC‑AAACGCAATTGCGCGGgag                                                                                                     <  2:476617/236‑1 (MQ=255)

                                                                                                                        aaaaaGGAGTTTTTTTCG‑CCCCTACCCAGTAGGGATGGGTTGTGAAAAGGGGGGGGGGAAATAG‑CCGGGTGGCCTTCACTACCGGGCGCAGGTTGAGATGGGGGTGCGGCGGTGGTCGG‑GGTATTCGGC‑TTG‑CCAGAAATTATCGCAAATGGCCTGTGGCA‑GGGTTTCGGCAGA‑GTTTTCCACCGCCACGCCTTGCTGCTGCGCCATTTTGCCAACCGC‑AAACGCAATTGCGCGGgag                                                                                                     <  2:391325/236‑1 (MQ=255)

                                                                                                                            aggagaTTTATCCC‑CCACTACCAGGTATGGATGAGATATGAAAAGGGGGAGGGGAAATAG‑CCCGTTGGCCTTCACTACCGGGCGCGGGTTTAGATGGGGGTAGGGCGGTGGTGGC‑GGTTTGCGGC‑TTC‑CAAGAAATTATCGTCAATGGCTTGTTGCG‑GGGTTTCGGCAGA‑GTTTTTCACCGCCACGCCTTGCTGCTGCGCCATTTTGCCAACCGC‑AAACGCAATTGCGCGGGAGACttt                                                                                                <  2:1007200/237‑1 (MQ=255)

                                                                                                                                   ttATTCGCCCCCTACCTAGGATGAGTGGGAGGTGAAAATGGGGGAGGAAAAGGG‑CCTGGGAGCCGTCACTCCCGGGCGGGGGTTTAGGGGGGGGTAGGGGGGTGGTGGT‑GTTTTGCGGT‑TGG‑CCAGAAATTTGCGTCAATGGCCGGTTGGG‑GGCGTTCGGCGGA‑GGTTTTCGCCGCCACGCCTTGCTGCTGCGCCATTTTGCCAACCGC‑AAACGCAATTGCGCGGGAGACTTTCTGaa                                                                                           <  2:96326/238‑1 (MQ=255)

                                                                                                                                               tACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAG‑CCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGC‑GGGATTCGGC‑TTG‑CCAGAAATTATCGTCAATGGCCTGTTTCA‑GGGCTTCGGCAGA‑GGGTTTTACCGCCACGCCTTGCTGCTGCGGCACTTTTCCAACCGC‑CAACGCAATTGCCCGGGGGGCTTTCTGCATATTTTTCAGtt                                                                               >  1:1459972/1‑238 (MQ=255)

                                                                                                                                                  ccagagtGGATGATATCAGAAAGGGGGGGAGAGAAAGAC‑CCCTGTACCCTTCCTTACCGGACGCGTGCTAAGATGGGGGTACGGCTGTAGCCGC‑GGTATTCGGC‑TTG‑CCAAAATTTCTCTTCATTGGCCTGTTGCG‑GGGCTTCGGCAGA‑TGTTTACACCGCCACGCCTTGCTGCTGCGCCATTTTGCCAACCGC‑AAACGCAATTGCGCGGGAGACTTTCTGAATATCTTTCAGTTCCgg                                                                           <  2:2173385/233‑1 (MQ=255)

                                                                                                                                                    aGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAG‑CCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGC‑GGGATTCGGG‑CTT‑CCAGAAAATATCGTCAATGGGCCGTTGCT‑GGGGTTCGGGAGA‑GGGTTTTCCCCCCACGCCCTGCTGCCGGGCCCTTTTTTCCACCGC‑CAAAGCAACTTCTCGGGAGAGATTTTGGATAACTTTTAGTTCCGGCa                                                                         >  2:470943/1‑239 (MQ=255)

                                                                                                                                                         ggtggagaTCGAAAAAAGGGAGAGGAAAGTAG‑CGCGGTAGCTTTCACTACCGGGCGCAGGCTTAGATGGAGGGGCGGCGGTAGTCGG‑GGTTTTCGGG‑TGC‑CCAGAAATTGGCGGCAAGGGCCTGTTGCG‑GGGTTTCGGCAGA‑GTTTTTCACCGCCACGCCTTGCTGCTGCGCCATTTTGCCAACCGC‑AAACGCAATTGCGCGGGAGACTTTCTGAATATCTTTCAGTTCCGGCAGTA‑c                                                                    <  2:96172/234‑1 (MQ=255)

                                                                                                                                                               tgtCTGGAAAAGGGGACGGGAAAGCGCCCCGGTAGCCGTCAACAGCGGGGGGTGGCTGAGATGGAGGTCCGGCGGGAGTCGC‑GTTATTCGGA‑TTG‑CCACAAATTAGCGTCACTGTCCTGTTGCG‑GGGGTTCGGAATA‑GTTTTCCTCCGCCACGCCTTGCTGCTGCGCCATTTTGCCAACCGC‑AAACGCAATTGCGCGGGAGCCTTTCTGAATATCTTTCAGTTCCGGCAGTA‑CCATAc                                                               <  2:1558325/236‑1 (MQ=255)

                                                                                                                                                                aTCTGAAAAAGGGAGAGGGAAATAG‑CCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGC‑GGTATTCGGG‑TTT‑GCAGAAAATATCGGCAAAGGGCCGTTGGC‑GGGGTTTGGCCGAGGGTTTTCCCCCCCCCCCCTT‑CCCCTGCGCCACTTTTCCAACCCCAAAACCCAATTGCCCGGG‑GGCCTTCTGGAATTCTTTTATGTCCGGCAGAA‑CCACAccc                                                             >  2:1888213/1‑237 (MQ=255)

                                                                                                                                                                   tGAAAAAGGGAGAGGGAAATAG‑CCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGC‑GGTATTCGGT‑TTG‑CCAGAAATTATCGTCAATGGCCTGTTGCA‑GGGCTTCGGCAGA‑GGTTTTCACCGCCACGCCTTGCTGCTGCGCCATTTTGCCAACCGC‑AAACGCAATTGCGCGGGAGACTTTCTGAATATCTTTCAGTTCCGGCAGTA‑CCATACCTTCGc                                                         <  1:1625841/239‑1 (MQ=255)

                                                                                                                                                                           gtgtggggAAATAG‑GCAGGGAGCACTACCTACCGGGCGAAGGCTTAGATGGGGGTCGGGCTGTGGTGGG‑GGTATGCGGC‑TGG‑CCAGAATTTATCGGCATTGGCGTGTTGCG‑GGGTTTCGGAAGA‑GGTTTCCACCGCCACGCCTTGCTGCTGCGCCATTTTGCCAACCGC‑AAACGCAATTGCGCGGGAGACTTTCTGAATATCTTTCAGTTCCGGCAGTA‑CCATACCTTCGCCGTTCAg                                                  <  2:1562826/233‑1 (MQ=255)

                                                                                                                                                                              gagGGAAATAG‑CCCGGTAGCCTTCACTTCCGGGCGCAGGCTTAGATGGAGGTCCGGCGGTAGTCGC‑GGTTTTCGGC‑TTG‑CCAAAAATTATCGTCAATGGCCTGTTGCA‑GGGCTTCGGCAGA‑GTTTTCCACCGCCACGCCTTGCTGCTGCGCCATTTTGCCAACCGC‑AAACGCAATTGCGCGGGAGACTTTCTGAATATCTTTCAGTTCCGGCAGTA‑CCATACCTTCGCCGTTCAGCAc                                               <  1:2100158/238‑1 (MQ=255)

                                                                                                                                                                              gagGGAAAACG‑CCCGTTGGCTTTCCCTCCGGGGGGCGGGTTTAAAGGGAGGTGGGGCGGTGTGCGC‑GGTATTGGGC‑TTG‑CCAAAATTTATCGTCAA‑GGTGTGTTTCAGGGGTTTGGGCGGG‑GGTTTTCACCGCCACGCCTTGCTGCTGCGCCATTTTGCCAACCGC‑AAACGCAATTGCGCGGGAGACTTTCTGAATATCTTTCAGTTCCGGCAGTA‑CCATACCTTCGCCGTTCAGCAc                                               <  2:278175/238‑1 (MQ=255)

                                                                                                                                                                                gggAAATAG‑CCCGGTGGCCTTCCCTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGC‑GGTTTTCGGC‑TTG‑CCAGAAATTATCGTAAATGGCCTTTTGCA‑GGGCTTCGGCAGA‑GGTTTTCACCGCCACGCCTTGCTGCTGCGCCATTTTGCCAACCGC‑AAACGCAATTGCGCGGGAGACTTTCTGAATATCTTTCAGTTCCGGCAGTA‑CCATACCTTCGCCGTTCAGCACCa                                             <  1:1936913/238‑1 (MQ=255)

                                                                                                                                                                                              gTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGGAGGTCC‑GGGAATCGGC‑TTG‑CCAGAAACTATTTCTAATTGCCCGGTGCA‑GGGCTTCGGGACA‑GGGTTTCCCCCCCCCGCCTTTCCGCTGCGCCATTTTTCCCAACGC‑AAAACCAATTCCGCGGGGGGCCTTTTTAATATATTTCTGTTTCGGCAAGT‑CCACACCTTCCCCGTTCCGCTCCAATGTTGacaatcac                               >  2:2147896/1‑232 (MQ=255)

                                                                                                                                                                                                    ttCACTCCCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGC‑GGTATTCGGC‑TTG‑CCAAAAATTATCGTCAATGGCCTGTTGCA‑GGGCTTCGGCAGA‑GGTTTTCACCGCCACGCCTTGCTGCTGCGCCATTTTGCCAACCGC‑AAACGCAATTGCGCGGGAGACTTTCTGAATATCTTTCAGTTCCGGCAGTA‑CCATACCTTCGCCGTTCAGCACCa                                             <  1:438380/219‑1 (MQ=255)

                                                                                                                                                                                                    ttCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGGAGGCGC‑GGTATTTGGC‑TTG‑CCCGAAAATATCGTTAAAGGGCTGTTGCC‑GGGGTTTGGCCGA‑GGGTTTTACCCCCCCGCCCTTCTGCCGCGCCATTTTTCCCACCGC‑AAAACCAAATTCGCGGGGGGCCTTCTTAATTTCTTTTCTGTCCCGCAGTACCCATACCCTCCCCGTTAGCCCCCAATGGGGaaacctcccccc                           >  2:2115969/1‑228 (MQ=255)

                                                                                                                                                                                                     tCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGC‑GGTATTCGGC‑TTG‑CCCGAAAATATCGTCAATGGCCTGTTGCA‑GGGCTTCGGCCGA‑GGGTTTTACCCCCACGCCTTGCTGCTGCGCCATTTTTCCCAACGC‑CAACGCCATTTCGCGGGGGGCTTTCTGAATATCTTTCAGTTCCCGCCGTA‑CCATACCTTCGCCGTTCAGCAACAATGGGGAATACTGCGCCAGCg                        >  1:722146/1‑239 (MQ=255)

                                                                                                                                                                                                               gcgcAGGCTTAGAGGGGGGTACGGCGGTGGTCGC‑GGTATTCGGC‑TGG‑CCAGAAATTATCGTCAATGGCCTGTTGCA‑GGGCTTCGGCAGA‑GTTTTTCACCGCCACGCCTTGCTGCTGCGCCATTTTGCCAACCGC‑AAACGCAATTGCGCGGGAGACTTTCTGAATATCTTTCAGTTCCGGCAGTA‑CCATACCTTCGCCGTTCAGCACCAATGGTGAATACTGCGCCAGCGTTTCACTTg               <  1:548040/238‑1 (MQ=255)

                                                                                                                                                                                                                 gcagctttAGATGGAGGTCCGGCGGTGGTCGC‑GGTGTTCGGC‑TTG‑CCAGAAATTATCGTCAATGCCCTGTTGCA‑GGGCTTCGGCAGA‑GGTTTTCACCGCCACGCCTTGCTGCTGCGCCATTTTGCCAACCGC‑AAACGCAATTGCGCGGGAGACTTTCTGAATATCTTTCAGTTCCGGCAGTA‑CCATACCTTCGCCGTTCAGCACCAATGGTGAATACTGCGCCAGCGTTTCACTTGcc             <  1:356953/232‑1 (MQ=255)

                                                                                                                                                                                                                 gcAGGATTAGGTGGGGTTACGGCGGTAGTCGC‑GGTAGTCGGC‑TTG‑CCAGAATTTATCGTCAATGGCCTGTTGCA‑GGGTTTCGCCAGA‑GTTTTTCACCGCCACGCCTTGCTGCTGCGCCATTTTGCCAACCGC‑AAACGCAATTGCGCGGGAGACTTTCTGAATATCTTTCAGTTCCGGCAGTA‑CCATACCTTCGCCGTTCAGCACCAATGGTGAATACTGCGCCAGCGTTTCACTTGcc             <  1:2171355/238‑1 (MQ=255)

                                                                                                                                                                                                                      tttggATGGAGGTCCGGCGGTAGTCGC‑GGTATTCGGC‑TTC‑CCAGAAATTATCGCCAATGGCCTGTTGCG‑GGGCTCCGGCAGA‑GGTTTCCACCGCCACGCCTTGCTGCTGCGCCATTTTGCCAACCGC‑AAACGCAATTGCGCGGGAGACTTTCTGAATATCTTTCAGTTCCGGCAGTA‑CCATACCTTCGCCGTTCAGCACCAATGGTGAATACTGCGCCAGCGTTTCACTTGCCGACAt        <  1:1849369/234‑1 (MQ=255)

                                                                                                                                                                                                                      cTTAGATGGAGGTACGGCGGTGGCCGC‑GGTTTTCGGT‑TTG‑CCGGAAATTATCGTCATTGGCCTGTTGCA‑GGCCTTCGGCAGA‑GGTTTTCACCGCCACGCCTTGCTGCTGCGCCATTTTGCCAACCGC‑AAACGCAATTGCGCGGGAGACTTTCTGAATATCTTTCAGTTCCGGCAGTA‑CCATACCTTCGCCGTTCAGCACCAATGGTGAATACTGCGCCAGCGTTTCACTTGCCGACAt        <  1:1634814/238‑1 (MQ=255)

                                                                                                                                                                                                                        tAGAGGGAGGTAGGGCGGTAGTCGC‑GGTATTCGGC‑TTG‑CCAGAAATTACCGTCAATGGCCTGTTGCA‑GGGTTTCGGCAGA‑GGTTTTCACCGCCACGCCTTGCTGCTGCGCCATTTTGCCAACCGC‑AAACGCAATTGCGCGGGAGACTTTCTGAATATCTTTCAGTTCCGGCAGTA‑CCATACCTTCGCCGTTCAGCACCAATGGTGAATACTGCGCCAGCGTTTCACTTGCCGACATCa      <  1:313788/238‑1 (MQ=255)

                                                                                                                                                                                                                            ggggggTCCGGCGGTAGGCGG‑GGTATTCGGC‑TTG‑CCAGAAATTATCGTCAATGGCCTGTTGCG‑GGGCTTCGGCAGA‑GTTTTTCACCGCCACGCCTTGCTGCTGCGCCATTTTGCCAACCGC‑AAACGCAATTGCGCGGGAGACTTTCTGAATATCTTTCAGTTCCGGCAGTA‑CCATACCTTCGCCGTTCAGCACCAATGGTGAATACTGCGCCAGCGTTTCACTTGCCGACATCAGCAt  <  1:150310/234‑1 (MQ=255)

GCTTCGCCATGTTTCAGTGAGCGCAGTGTTTTATACGTAACGTCAACATGATGATCCTTCGTAACAACTGCAGCCTTCATAGTTCCTCCTTTTCGGATGATGTTCTGCATAGCAGGTGAGGCAAATGAGATTTATTCG‑CCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAG‑CCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGC‑GGTATTCGGC‑TTG‑CCAGAAATTATCGTCAATGGCCTGTTGCA‑GGGCTTCGGCAGA‑GGTTTTCACCGCCACGCCTTGCTGCTGCGCCATTTTGCCAACCGC‑AAACGCAATTGCGCGGGAGACTTTCTGAATATCTTTCAGTTCCGGCAGTA‑CCATACCTTCGCCGTTCAGCACCAATGGTGAATACTGCGCCAGCGTTTCACTTGCCGACATCAGCAT  >  NC_007779/1555473‑1555928

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