Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_007779	1,555,715	A→T	14.9%	Y557N (TAC→AAC)	Y75_RS07805 ←	malate dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007779	1,555,715	0	A	T	14.9%	116.8 / 2.0	66	Y557N (TAC→AAC)	Y75_RS07805	malate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (29/23);  new base T (2/7);  total (35/31)
Fisher's exact test for biased strand distribution p-value = 8.07e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTTCAGTGAGCGCAGTGTTTTATACGTAACGTCAACATGATGATCCTTCGTAACAACTGCAGCCTTCATAGTTCCTCCTTTTCGGATGATGTTCTGCATAGCAGGTGAGGCAAATGAGATTTATTCG‑CCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAG‑CCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGC‑GGTATTCGGC‑TTG‑CCAGAAATTATCGTCAATGGCCTGTTGCA‑GGGCTTCGGCAGA‑GGTTTTCACCGCCACGCCTTGCTGCTGCGCCATTTTGCCAACCGC‑AAACGCAATTGCGCGGGAGACTTTCTGAATATCTTTCAGTTCCGGCAGTA‑CCATACCTTCGCCGTTCAGCACCAATGGTGAATACTGCGCCAGCGTTTCACTTGCCGACATCAGCAT  >  NC_007779/1555484‑1555928                                                                                                                                                                                                                                           |                                                                                                                                                                                                                            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‑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)                                                                                                                                                                                                                                           |                                                                                                                                                                                                                            TTTCAGTGAGCGCAGTGTTTTATACGTAACGTCAACATGATGATCCTTCGTAACAACTGCAGCCTTCATAGTTCCTCCTTTTCGGATGATGTTCTGCATAGCAGGTGAGGCAAATGAGATTTATTCG‑CCACTACCCAGTATGGATGAGATCTGAAAAAGGGAGAGGGAAATAG‑CCCGGTAGCCTTCACTACCGGGCGCAGGCTTAGATGGAGGTACGGCGGTAGTCGC‑GGTATTCGGC‑TTG‑CCAGAAATTATCGTCAATGGCCTGTTGCA‑GGGCTTCGGCAGA‑GGTTTTCACCGCCACGCCTTGCTGCTGCGCCATTTTGCCAACCGC‑AAACGCAATTGCGCGGGAGACTTTCTGAATATCTTTCAGTTCCGGCAGTA‑CCATACCTTCGCCGTTCAGCACCAATGGTGAATACTGCGCCAGCGTTTCACTTGCCGACATCAGCAT  >  NC_007779/1555484‑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: ‑