Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	AE015451	4,741,231	C→G	80.8%	intergenic (+468/+285)	gltA → / ← yeiW	citrate synthase/putative oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	AE015451	4,741,231	0	C	G	80.8%	189.5 / 46.6	73	intergenic (+468/+285)	gltA/yeiW	citrate synthase/putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base C (12/2);  new base G (26/33);  total (38/35)
Fisher's exact test for biased strand distribution p-value = 6.64e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.36e-04

ACTTTTGCCACCCGAACAAAACCCTTCTCTTCCGGCAGGTAAGTCGGCTCTTCTTTCAGGTTACGGGTCAGGCTGCCATCCGCCTCCTTGA‑TGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑AACCTCC‑GCCAAC‑TGCAATTCGTGAA‑CCAAGCAAAAGGCCCCGAAAAGGCCCTC‑C‑GCAGCCCCCCCAACCTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAGGGACACCTTGATCCCGTCCTCGTCCTTGGCCACTTGCCAGTTTTCAGCCCAGGCCACAGGCGACAGCAGAACACCCACACC  >  AE015451/4741044‑4741449                                                                                                                                                                                              |                                                                                                                                                                                                                                aCTTTTGCCACCCGAACAAAACCCTTCTCTTCCGGCAGGTAAGTCGGCTCTTCTTTCAGGTTACGGGTCAGGCTGCCATCCGCCTCCTTGA‑TGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CAGCTTC‑TGCGAC‑TTGCATTCGTGAA‑TCCAGGCACAGGCCCCggacacgtc                                                                                                                                                                              >  1:1042414/1‑228 (MQ=255)                             ctTCCGGCAGGTAAGTCGGCTCTTCTTTCAGGTTACGGGTCAGGCTGCCATCCGCCTCCTTGA‑TGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CAGCTTC‑TGCGAC‑TTGCATTCGTGAA‑TCCAGGCACAGGCCCCggaca                                                                                                                                                                                  >  1:355626/1‑200 (MQ=255)                             ctTCCGGCAGGTAAGTCGGCTCTTCTTTCAGGTTACGGGTCAGGCTGCCATCCGCCTCCTTGA‑TGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CAGCTTC‑TGCGAC‑TTGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTC‑CTGCAGCGCCTCGACCTTGGCCAg                                                                                                                                                  >  2:1011584/1‑235 (MQ=255)                              ttCCGGCAGCTCAGTCGGCTCTTCTTTCAGGTTACGGGTCAGGCTGCCATCCGCCTCCTTGA‑TGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CAGCTTC‑TGCGAC‑TTGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTC‑CTGCAGCGCCTCGACCTTGGCCa                                                                                                                                                   <  1:834158/233‑1 (MQ=255)                                      ggTAAGTCGGCTCTTCTTTCAGGTTACGGGTCAGGCTGCCATCCGCCTCCTTGA‑TGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CAGCTTC‑TGCGAC‑TTGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTC‑CTGCAGCGCCTCGACCTTGGCCAGCGgcgc                                                                                                                                            >  2:673713/1‑232 (MQ=255)                                          attCGGCTCTTCTTTCAGGTTACGGGTCAGGCTGCCATCCGCCTCCTTGA‑TGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CAGCTTC‑TGCGAC‑TTGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTC‑CTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATc                                                                                                                                     <  1:491475/233‑1 (MQ=255)                                            tCGGCTCTTCTTTCAGGTTACGGGTCAGGCTGCCATCCGCCTCCTTGA‑TGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CAGCTTC‑TGCGAC‑TTGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTC‑CTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCAc                                                                                                                                   <  1:182985/235‑1 (MQ=255)                                              ggCTCGTCTTTCAGGTTACGGGTCAGGCTGCCATCCGCCTCCTTGA‑TGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CAGCTTC‑TGCGAC‑TTGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGtcctc‑ct                                                                                                                                                                        <  1:65049/196‑2 (MQ=255)                                               gCTCTTCTTTCAGGTTACGGGTCAGGCTGCCATCCGCCTCCTTGA‑TGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CAGCTTC‑TGCGACTTGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTC‑CTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTc                                                                                                                               <  1:673716/236‑1 (MQ=255)                                                 tcttctTCCAGGTTACGGGTCAGGCTGCCATCCGCCTCCTTGA‑TGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CAGCTTC‑TGCGAC‑TTGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTC‑CTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTca                                                                                                                              <  1:254810/235‑1 (MQ=255)                                                         cAGTTTACGGGTCAGGCTGCCATCCGCCTCCTTGA‑TGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CAGCTTC‑TGCGAC‑TTGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTC‑CTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGa                                                                                                                      >  2:633773/1‑235 (MQ=255)                                                         cAGGTTACGGGTCAGGCTGCCATCCGCCTCCTTGA‑TGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATCGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CAGCTTC‑TGCGAC‑TTGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTC‑CTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTca                                                                                                                              <  2:536016/227‑1 (MQ=255)                                                             ttACGGGTCAGGCTGCCATCCGCCTCCTTGATTGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CAGCTTC‑TGCGAC‑TTGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTC‑CTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAgg                                                                                                                  >  1:420541/1‑236 (MQ=255)                                                                cGGGTCAGGCTGCCATCCGCCTCCTTGA‑TGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CAGCTTC‑TGCGAC‑TCGCATTCGTGAA‑TCCAGGCACAGGCCCCgg                                                                                                                                                                                     <  1:1184436/166‑2 (MQ=255)                                                                         cTGCCATCCGCCTCCTTGA‑TGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCT‑GCGACT‑TGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTC‑CTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTg                                                                                                       >  2:991576/1‑234 (MQ=255)                                                                             cctCCGCCTCCTTGA‑TGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CAGCTTCTGCGACT‑TGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTC‑CTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGc                                                                                                                 <  2:941281/218‑1 (MQ=255)                                                                                      cctTGA‑TGGTCGTAATGTGCATCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTTGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CAGCTTC‑TGCGAC‑TTGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTC‑CTGCAGCGCCTCGACCTTGGCCAGCGgcgc                                                                                                                                            <  1:457545/184‑1 (MQ=255)                                                                                          gA‑TGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCT‑GCGACT‑TGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTC‑CTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGa                                                                                    >  1:895503/1‑236 (MQ=255)                                                                                           a‑TGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CAGCTTC‑TGCGAC‑TTGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTC‑CTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAg                                                                                   >  1:65812/1‑236 (MQ=255)                                                                                           a‑TGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCATGCTCGAG‑CAGCTTC‑TGCGAC‑TTGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTC‑CTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAg                                                                                   <  1:652571/236‑1 (MQ=255)                                                                                                tCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGATCACCTTCGTGCTTGAG‑CAGCTTC‑TGCGACTTGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCAT‑GCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTg                                                                                     <  2:650861/230‑1 (MQ=255)                                                                                                     aTGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCT‑GCGACT‑TGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTC‑CTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGAccc                                                                                                  >  1:992474/1‑212 (MQ=255)                                                                                                     aTGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCT‑GCGACT‑TGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTC‑CTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGAccc                                                                                                  <  2:992471/212‑1 (MQ=255)                                                                                                            gcaCGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCT‑GCGACT‑TGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTC‑CTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTg                                                                                                       <  2:929013/200‑1 (MQ=255)                                                                                                             caCGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CAGCTTC‑TGCGAC‑TTGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTC‑CTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCa                                                                                                                                    >  2:69224/1‑170 (MQ=255)                                                                                                              aCGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CAGCTTCTGCGACT‑TGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTC‑CTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAGGGACACCTTGATCCcgt                                                                  <  1:1011587/235‑1 (MQ=255)                                                                                                                ggAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCT‑GCGACT‑TGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTC‑CTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAGGGACACCTTGATCCcgtc                                                                 >  2:1106805/1‑234 (MQ=255)                                                                                                                ggAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CAGCTTC‑TGCGAC‑TTGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTC‑CTGCAGCGCCTCGACCTTGGCTAGCGGCGCCTTGATCACGGTCACACcgc                                                                                                                        >  1:356505/1‑179 (MQ=255)                                                                                                                 gAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCT‑GCGACT‑TGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTC‑CTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAGGGACACCTTGATCCcgtcct                                                               >  1:712245/1‑235 (MQ=255)                                                                                                                     cACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCT‑GCGACT‑TGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTC‑CTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAGGGACACCTTGATCCcgtcctcgtcc                                                          >  1:425249/1‑236 (MQ=255)                                                                                                                       cgcgGCGTAACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CAGCTTC‑TGCGAC‑TTGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTC‑CTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACcgc                                                                                                                        <  2:594219/172‑1 (MQ=255)                                                                                                                         cggcgTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CAGCTTC‑TGCGAC‑TTGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTC‑CTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCgg                                                                                                <  2:1205137/194‑1 (MQ=255)                                                                                                                          ggcgTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCCCCTTCGTGCTTGAG‑CAGCTTC‑TGCGAC‑TTGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTC‑CTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCg                                                                                                 <  2:186512/192‑1 (MQ=255)                                                                                                                          ggcgTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCT‑GCGACT‑TGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTC‑CTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAGGGACACCTTGATCCCGTCCTCGTCCTTGGc                                                     >  1:614084/1‑236 (MQ=255)                                                                                                                          ggcgTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCT‑GCGACT‑TGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTC‑CTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAGGGACACCTTGATCCCGTCCTCGTCCTTGGc                                                     >  1:473379/1‑236 (MQ=255)                                                                                                                          ggcgTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CAGCTTC‑TGCGAC‑TTGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTC‑CTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGAc                                                                                                    >  1:6783/1‑189 (MQ=255)                                                                                                                          ggcgTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CAGCTTC‑TGCGAC‑TTGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTC‑CTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAGGGa                                                                                >  1:483893/1‑209 (MQ=255)                                                                                                                           gcgTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CAGCTTC‑TGCGAC‑TTGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTC‑CTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCa                                                                                             <  2:483890/195‑1 (MQ=255)                                                                                                                               cACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCT‑GCGACT‑TGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTC‑CTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAGGGACACCTTGATCCCGTCCTCGTCCTTGGCCACt                                                 >  1:484650/1‑235 (MQ=255)                                                                                                                               cACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CAGCTTC‑TGCGAC‑TTGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTC‑CTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTg                                                                                     <  1:634214/199‑1 (MQ=255)                                                                                                                                aCCGGACATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CAGCTTC‑TGCGAC‑TTGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTC‑CTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCAc                                                                                       <  2:1073924/196‑1 (MQ=255)                                                                                                                                 ccGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CAGCTTC‑TGCGAC‑TTGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTC‑CTGCAGCGCCTCGACCTTGGCCAGCGGCGCCtt                                                                                                                                         <  2:141701/145‑1 (MQ=255)                                                                                                                                   ggCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CAGCTTC‑TGCGAC‑TTGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTC‑CTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCt                                                                                          <  1:25080/190‑1 (MQ=255)                                                                                                                                        tGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTTCTTGAG‑CAGCTTC‑TGCGAC‑TTGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTC‑CTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCa                                                                                             <  2:519745/182‑1 (MQ=255)                                                                                                                                         ggCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CAGCTTC‑TGCGAC‑TTGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTC‑CTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAGGGACACCTTGATCCCGTCCTCGTCCTTGGCCACTTGCCAGttt                                        >  2:548128/1‑234 (MQ=255)                                                                                                                                            gTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCT‑GCGACT‑TGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTC‑CTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAGGGACACCTTGATCCCGTCCTCGTCCTTGGCCCCTTGCCAGt                                          >  2:936479/1‑229 (MQ=255)                                                                                                                                            gTATTGAACT‑GGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCT‑GCGACT‑TGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTC‑CTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAGGGACACCTTGATCCCGTCCTCGTCCTTGGCCACTTGCCAGTTTTCAGcc                                  >  1:823086/1‑236 (MQ=255)                                                                                                                                               ttGACCTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CAGCTTC‑TGCGAC‑TTGCATTCGTGAA‑TCCAGCCACAGGCCCCGGACACGTCCTC‑CTGCAGCGCCTCGACCTTGGCCAgcggcg                                                                                                                                             <  2:81494/127‑1 (MQ=255)                                                                                                                                                  aaCTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CAGCTTC‑TGCGAC‑TTGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTC‑CTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACCGTCACACCGCGAt                                                                                                                     >  2:108258/1‑148 (MQ=255)                                                                                                                                                         tataGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CAGCTTC‑TGCGAC‑TTGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTC‑CTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAGGGACACCTTGATCCCGTCCTCGTCCTTGGCCACTTGCCAGTTTTCAGCCCAGGCcaca                         <  2:425246/233‑1 (MQ=255)                                                                                                                                                          ataGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCT‑GCGACT‑TGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTC‑CTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAGGGACACCTTGATCCCGTCCTCGTCCTTGGCCACTTGCCAGTTTTCAGCCCAGGCCACAg                        >  2:930868/1‑233 (MQ=255)                                                                                                                                                                 cAGCTCTGGTCACCTTCGTGCTTGAG‑CAGCTTC‑TGCGACTTGC‑ATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTC‑CTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAGGGACACCTTGATCCCGTCCTCGTCCTTGGCCACTTGCCAGTTTTCAGCCCAGGCCACAGGCGAc                   <  1:1220645/231‑1 (MQ=255)                                                                                                                                                                 cAGCTCTGGTCACCTTCGTGCTTGAG‑CAGCTTC‑TGCGAC‑TTGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTC‑CTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAGGGACACCTTGATCCCGTCCTCGTCCTTGGCCACTTGCCAGTTTTCAGCCCAGGCCACAGGCGAc                   >  2:1220642/1‑231 (MQ=255)                                                                                                                                                                  aGCTCTGGTCACCTTCGTGCTTGAG‑CAGCTTC‑TGCGAC‑TTGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTC‑CTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAGGGACACCTTGATCCcgt                                                                  <  1:1136058/183‑1 (MQ=255)                                                                                                                                                                    ctctGGTCACCTTCGTGCTTGAG‑CAGCTTC‑TGCGACTTGC‑ATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTC‑CTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAGGGACACCTTGATCCCGTCCTCGTCCTTGGCCACTTGCCAGTTTTCAGCCCAGGCCACAGGCGAcag                 >  2:451747/1‑230 (MQ=255)                                                                                                                                                                    ctctGGTCACCTTCGTGCTTGAG‑CAGCTTC‑TGCGAC‑TTGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTC‑CTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAGGGACACCTTGATCCCGTCCTCGTCCTTGGcc                                                    >  2:764612/1‑195 (MQ=255)                                                                                                                                                                    ctctGGTCACCTTCGTGCTTGAG‑CAGCTTC‑TGCGAC‑TTGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTC‑CTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAGGGACACCTTGATCCCGTCCTCGTCCTTGGCCACTTGCCAGTTTTCAGCCCa                                >  1:16568/1‑215 (MQ=255)                                                                                                                                                                    ctctGGTCACCTTCGTGCTTGAG‑CAGCTTC‑TGCGAC‑TTGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTC‑CTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAGGGACACCTTGATCCCGTCCTCGTCCTTGGCCACTTGCCAGTTTTCAGCCCAGGCCACAGGCGACAGCAGaa            >  1:70371/1‑235 (MQ=255)                                                                                                                                                                    ctctGGTCACCTTCGTGCTTGAG‑CAGCTTC‑TGCGAC‑TTGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTC‑CTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAGGGACACCTAGATCCCGTCCTCGTCCTTGGCCACTTGCCAGTTTTCAGCCCAGGCCACAGGCGAc                   >  1:978627/1‑228 (MQ=255)                                                                                                                                                                    ctctGGTCACCTTCGTGCTTGAG‑CAGCTTC‑TGCGAC‑TTGCATCCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTC‑CTGCAGCGCCTCGACCTTGGCCAGCGGCGCCtag                                                                                                                                        >  2:357670/1‑109 (MQ=255)                                                                                                                                                                     tctGGTCACCTTCGTGCTTGAG‑CAGCTTC‑TGCGAC‑TTGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTC‑CTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAGGGACCCCTTGATCCCGTCCTCGTCCTTGGCCACTTGCCAGTTTTCAGCCCAGGCCACAgg                       >  2:777113/1‑223 (MQ=255)                                                                                                                                                                     tctGGTCACCTTCGTGCTTGAG‑CAGCTTC‑TGCGAC‑TGTCACTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTC‑CTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAGGGACACCTTGATCCCGTCCTCGTCCTTGGCCACTTGCCAGTTTTCAGCCCAGGCCACAGGCGAcag                 <  1:369432/229‑1 (MQ=255)                                                                                                                                                                      ctGGTCACCTTCGTGCTTGAG‑CAGCTTC‑TGCGAC‑TTGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTC‑CTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAGGGACACCTTGATCCCGTCCTCGTCCTTGGCCACTTGCCAGTTTTCAGCCCAGGCCACAGGCGAcagcag              >  1:1147505/1‑231 (MQ=255)                                                                                                                                                                      ctGGTCACCTTCGTGCTTGAG‑CAGCTTC‑TGCGAC‑TTGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTC‑CTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAGGGACACCTTGATCCCGTCCTCGTCCTTGGCCACTTGCCAGTTTTCAGCCCAGGCCACAGGCGAcagc                >  1:284103/1‑229 (MQ=255)                                                                                                                                                                      ctGGTCACCTTCGTGCTTGAG‑CAGCTTC‑TGCGAC‑TTGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTC‑CTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAGGGACACCTTGATCCCGTCCTCGTCCTTGGCCACTTGCCAGTTTTCAGCCCAGGCCACAGGCGACAGCAGAacac         >  1:1190014/1‑236 (MQ=255)                                                                                                                                                                       tGGTCACCTTCGTGCTTGAG‑CAGCTTC‑TGCGAC‑TTGCATTCGTGCA‑TCCAGGCACAGGCCCCGGACACGTCCTC‑CTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAGGGACACCTTGATCCCGTCCTCGTCCTTGGCCACTTGCCAGTTTTCAGCCCAGGCCACAGGCGACAGCAGAacac         <  1:137528/235‑1 (MQ=255)                                                                                                                                                                       tGGTCACCTTCGTGCTTGAG‑CAGCTTC‑TGCGAC‑TTGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTC‑CTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAGGGACACCTTGATCCCGTCCTCGTCCTTGGCCACTTGCCAGTTTTCAGCCCAGGCCACAGGCGACAGCAGAacac         <  1:737887/235‑1 (MQ=255)                                                                                                                                                                       tGGTCACCTTCGTGCTTGAG‑CAGCTTC‑TGCGAC‑TTGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTC‑CTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAGGGACACCTTGATCCCGTCCTCGTCCTTGGCCACTTGCCAGTTTTCAGCCCAGGCCACAGGCGACAGCAGAacac         <  1:474295/235‑1 (MQ=255)                                                                                                                                                                        ggTCACCTTCGTGCTTGAG‑CAGCTTC‑TGCGAC‑TTGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTC‑CTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAGGGACACCTTGATCCCGTCCTCGTCCTTGGCCACTTGCCAGTTTTCAGCCCAGGCCACAGGCGAcag                 <  2:887522/226‑1 (MQ=255)                                                                                                                                                                            aCCTTCGTGCTTGAG‑CAGCTTC‑TGCGAC‑TTGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTC‑CTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAGGGGCACCTTGATCCCGTCCTCGTCCTTGGCCACTTGCCAGTTTTCAGCCCAGGCCACAGGCGAcagc                >  1:332747/1‑223 (MQ=255)                                                                                                                                                                            aCCTTCGTGCTTGAG‑CAGCTTC‑TGCGAC‑TTGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTC‑CTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAGGGACACCTTGATCCCGTCCTCGTCCTTGGCCACTTGCCAGTTTTCAGCCCAGGCCACAgg                       <  2:528270/216‑1 (MQ=255)                                                                                                                                                                            aCCTTCGTGCTTGAG‑CAGCTTC‑TGCGAC‑TTGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTC‑CTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGAGATAGGCCTTGTACTTGGAc                                                                                                    >  2:1270573/1‑139 (MQ=255)                                                                                                                                                                                tCGTGCTTGAG‑CAGCTTC‑TGCGAC‑TTGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTC‑CTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAGGGACACCTTGATCCCGTCCTCGTCCTTGGCCACTTGCCAGTTTTCAGCCCAGGCCACAGGCGACAGCAGAACACccacacc  <  2:823083/233‑1 (MQ=255)                                                                                                                                                                                              |                                                                                                                                                                                                                                ACTTTTGCCACCCGAACAAAACCCTTCTCTTCCGGCAGGTAAGTCGGCTCTTCTTTCAGGTTACGGGTCAGGCTGCCATCCGCCTCCTTGA‑TGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑AACCTCC‑GCCAAC‑TGCAATTCGTGAA‑CCAAGCAAAAGGCCCCGAAAAGGCCCTC‑C‑GCAGCCCCCCCAACCTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAGGGACACCTTGATCCCGTCCTCGTCCTTGGCCACTTGCCAGTTTTCAGCCCAGGCCACAGGCGACAGCAGAACACCCACACC  >  AE015451/4741044‑4741449

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