Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	AE015451	4,741,283	+T	intergenic (+520/+233)	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,283	1	.	T	96.8%	103.7 / ‑6.9	32	intergenic (+520/+233)	gltA/yeiW	citrate synthase/putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base . (0/1);  major base T (16/14);  minor base G (0/1);  total (16/16)
Fisher's exact test for biased strand distribution p-value = 4.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

CGTCGACTTTTGCCACCCGAACAAAACCCTTCTCTTCCGGCAGGTAAGTCGGCTCTTCTTTCAGGTTACGGGTCAGGCTGCCATCCGCCTCCTTGATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑AACCTCCGCCAAC‑TGCAATTCGTG‑AA‑CCAAGCAAAAGGCCCC‑GAAAAGGCCCTCC‑GCAGCCCCCCCAACCTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAGGGACACCTTGATCCCGTCCTCGTCCTTGGCCACTTGCCAGTTTTCAGCCCAGGCCACAGGCGACAGCAGAACACCCACACCACACAGCAGTGCAATACGCTTGAACGATCTCATCATGTTTATCCTTGTTGTCGA  >  AE015451/4741039‑4741503                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                             cGTCGACTTTTGCCACCCGAACAAAACCCTTCTCTTCCGGCAGGTAAGTCGGCTCTTCTTTCAGGTTACGGGTCAGGCTGCCATCCGCCTCCTTGATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACT‑TGC‑ATTCGTG‑AATCCAGGC‑ACAGGCCCC‑GGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGTTCACACCCCGATAGGc                                                                                                                                                                       >  1:517204/1‑301 (MQ=255)                cccGAACAAAACCCTTCTCTTCCGGCAGGTAAGTCGGCTCTTCTTTCAGGTTACGGGTCAGGCTGCCATCCGCCTCCTTGATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACT‑TGC‑ATTCGTG‑AATCCAGGC‑ACAGGCCCC‑GGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGAccc                                                                                                                                                        <  1:88250/301‑1 (MQ=255)                  cGAACAAAACCCTTCTCTTCCGGCAGGTAAGTCGGCTCTTCTTTCAGGTTACGGGTCAGGCTGCCATCCGCCTCCTTGATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACT‑TGC‑ATTCGTG‑AATCCAGGC‑ACAGGCCCC‑GGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCgg                                                                                                                                                      >  1:712222/1‑301 (MQ=255)                                ctctTCCGGCAGGTAAGTCGGCTCTTCTTTCAGGTTACGGGTCAGGCTGCCATCCGCCTCCTTGATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACT‑TGC‑ATTCGTG‑AATCCAGGC‑ACAGGCCCC‑GGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTCCTTGGACTCGGCCACCTCACTGAgg                                                                                                                                        >  1:10017/1‑301 (MQ=255)                                  ctTCCGGCAGGTAAGTCGGCTCTTCTTTCAGGTTACGGGTCAGGCTGCCATCCGCCTCCTTGATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACT‑TGC‑ATTCGTG‑AATCCAGGC‑ACAGGCCCC‑GGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAGGGa                                                                                                                                      <  1:809742/301‑1 (MQ=255)                                   gtCCGGCAGGTAAGTCGGCGCTTCTTTCAGGTTACGGGTCAGTCTGCCATCCGCCTCCTTGAGGGTCGTAATGTGCAGAACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACT‑TGC‑ATTCGTG‑AATCCAGGC‑ACAGGCCCC‑GGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAGGGac                                                                                                                                     <  2:295781/300‑1 (MQ=255)                                                    gttctacgtgcAGGGTACGGGTCAGGCTGCCATCCGGCTCCTGGATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTAGTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACT‑TGC‑ATTCGTG‑AATCCAGGC‑ACAGGCCCC‑GGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAGGGACACCTTGATCCcgtcctc                                                                                                                    <  2:551576/291‑1 (MQ=255)                                                    cttcttctTTCAGGTTACGGGTCAGGCTGCCATCCGCCTCCTTGATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACT‑TGC‑ATTCGTG‑AATCCAGGC‑ACAGGCCCC‑GGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAGGGACACCTTGATCCcgtcctc                                                                                                                    <  2:117851/299‑1 (MQ=255)                                                       cttcttTCAGGTTACGGGTCAGGCTGCCATCCGCCTCCTTGATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACT‑TGC‑ATTCGTG‑AATCCAGGC‑ACAGGCCCC‑GGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAGGGACACCTTGATCCcgtcctcgtc                                                                                                                 >  1:299071/1‑301 (MQ=255)                                                           tttCAGGTTACGGGTCAGGCTGCCATCCGCCTCCTTGATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACT‑TGC‑ATTCGTG‑AATCCAGGC‑ACAGGCCCC‑GGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAGGGACACCTTGATCCCGTCCTCGTCCttt                                                                                                             >  1:787824/1‑300 (MQ=255)                                                                  ttcggggtcaggatcacattcGCCTCATTGATGGTCGGAATGTGCCGCACGGAAGCACGCGGCGTCACCGGACATGGCGTATTGAACGGGGGATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACT‑TGC‑ATTAGTG‑AATCCAGGC‑ACAGGCCCC‑GGACACGTCCTCCTGCAGCGCCTCGAACTTGGCCAGCGGCGCCTTGATCACGGACACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAGGGACACCTTGAGCCCGTCCTCGTCCTTGGCCACtt                                                                                                      <  2:290075/281‑1 (MQ=255)                                                                                    ccacgacgcgcttacgatttaaggtgCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATGGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACT‑TGC‑ATTCGTG‑AATCCAGGC‑ACAGGCCCC‑GGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAGGGACACCTTGATCCCGTCCTCGTCCTTGGCCACTTGCCAGTTTTCAGCCCAgg                                                                                    <  2:102914/278‑1 (MQ=255)                                                                                         ctcctTGATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACT‑TGC‑ATTCGTG‑AATCCAGGC‑ACAGGCCCC‑GGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAGGGACACCTTGATCCCGTCCTCGTCCTTGGCCACTTGCCAGTTTTCAGCCCAGGCcaca                                                                               <  1:618357/301‑1 (MQ=255)                                                                                         ctcctTGATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACT‑TGC‑ATTCGTG‑AATCCAGGC‑ACAGGCCCC‑GGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAGGGACACCTTGATCCCGTCCTCGTCCTGAGCCACTTCCCAGTTTTCAGccccgcacgac                                                                               >  2:578964/1‑293 (MQ=255)                                                                                          tgcttgagggtcgataggtgaagcaCGGAATCACCAGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTAGTCGCATTCGTGCTTGAGCAGCTTCCGCGCCT‑TGC‑ATTCGTG‑AATCCAGGC‑ACAGGCCCCAGACACGTCCTCCC‑GCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAGGGACACCTTGATCCCGTCCTCGTCCTTGGCCACTTGCCAGTTTTCAGCCCAGGCCACAg                                                                              <  2:10016/280‑1 (MQ=255)                                                                                          tcctTGATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACT‑TGC‑ATTCGTG‑AATCCAGGC‑ACAGGCCCC‑GGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAGGGACACCTTGATCCCGTCCTCGTCCTTGGCCACTTGCCAGTTTTCAGCCCAGGCCACAg                                                                              <  2:873313/301‑1 (MQ=255)                                                                                            ctTGATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACT‑TGC‑ATTCGTG‑AATCCAGGC‑ACAGGCCCC‑GGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAGGGACACCTTGATCCCGTCCTCGTCCTTGGCCACTTGCCAGTTTTCAGCCCAGGCCACAGGc                                                                            >  2:374789/1‑301 (MQ=255)                                                                                                    tCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACT‑TGC‑ATTCGTG‑AATCCAGGC‑ACAGGCCCC‑GGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAGGGACACCTTGATCCCGTCCTCGTCCTTGGCCACTTGCCAGTTTTCAGCCCAGGCCAAGGGCGAcagcaa                                                                    >  2:824229/1‑300 (MQ=255)                                                                                                        aaTGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACT‑TGC‑ATTCGTG‑AATCCAGGC‑ACAGGCCCC‑GGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAGGGACACCTTGATCCCGTCCTCGTCCTTGGCCACTTGCCAGTTTTCAGCCCAggagatcggaagagcacccgtc                                                                >  1:102944/1‑281 (MQ=255)                                                                                                              cagcaCGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACT‑TGC‑ATTCGTG‑AATCCAGGC‑ACAGGCCCC‑GGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAGGGACACCTTGATCCCGTCCTCGTCCTTGGCCACTTGCCAGTTTTCAGCCCAGGCCACAGGCGACAGCAGAACACccaca                                                          >  1:131506/1‑301 (MQ=255)                                                                                                                  aCGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACT‑TGC‑ATTCGTG‑AATCCAGGC‑ACAGGCCCC‑GGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAGGGACACCTTGATCCCGTCCTCGTCCTGTGCCACTTGCCAGTTTTCAGCCCCGGCCCCCGGCGACACCAGAACAcccaccccac                                                      >  1:194127/1‑301 (MQ=255)                                                                                                                   gggaatgaagcggcgTAACCGGCCATGGCGGATTCAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGGG‑CAGCTGCTGCAACTTGC‑ATTCGTGCAT‑CCAGGC‑ACAGGCCCC‑GGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAGGGACACCTTGATCCCGTCCTCGTCCTTGGCCACTTGCCAGTTTTCAGCCCAGGCCACAGGCGACAGCAGAACACccacaccaca                                                     <  2:10277/292‑1 (MQ=255)                                                                                                                                cgTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACT‑TGC‑ATTCGTG‑AATCCAGGC‑ACAGGCCCT‑GGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAGGGACACCTTGATCCCGTCCTCGTCCTTGGCCACTTGCCAGGTTTCAGCCCAGGCCACAGGCGACGGCAGAACACCCACACCACACAGCAGTGCacaa                                        >  2:307977/1‑298 (MQ=255)                                                                                                                                         ccATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACT‑TGC‑ATTCGTG‑AATCCAGGC‑ACAGGCCCC‑GGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAGGGACACCTTGATCCCGTCCTCGTCCTTGGCCACTTGCCAGTTTTCAGCCCTGGTCACCGTCCCCAGCAGAACACCCCCCCAACAcagcacgcccatccccttgaac                               >  2:638228/1‑284 (MQ=255)                                                                                                                                         ccATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACT‑TGC‑ATTCGTG‑AATCCAGGC‑ACAGGCCCC‑GGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAGGGACACCTTGATCCCGTCCTCGTCCTCGGCCACTTGCCAGTTTTCAGCCCAGGCCACAGGCGACAGCAGAACACCCACACCACACAGCAGTGCAATACGCTTGAAc                               <  1:218434/301‑1 (MQ=255)                                                                                                                                              gCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACT‑TGC‑ATTCGTG‑AATCCAGGC‑ACAGGCCCC‑GGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAGGGACACCTTGATCCCGTCCTCGTCCTTGGCCACTTGCCAGTTTTCAGCCCAGGCCACAGGCGACAGCAGAACACCCACACCACACAGCAGTGCAATACGCTTGAACGaacc                          >  2:848508/1‑298 (MQ=255)                                                                                                                                                   ttGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACT‑TGC‑ATTCGTG‑AATCCAGGC‑ACAGGCCCC‑GGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAGGGACACCTTGATCCCGTCCTCGTCCTTGGCCACTTGCCAGTTTTCAGCCCAGGCCACAGGCGACAGCAGAACACCCACACCACACAGCAGTGCAATACGCTTGAACGATCtcatca                     <  2:194058/301‑1 (MQ=255)                                                                                                                                                      aaCTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACT‑TGC‑ATTCGTG‑AATCCAGGC‑ACAGGCCCC‑GGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAGGGACACCTTGATCCCGTCCTCGTCCTTGGCCACTTGCCAGTTTTCAGCCCAGGCCACAGCCGCCAGCAGAACACCCACACACCACATCAGTCCACTACTCTTGAACGATCTCATCATGt                  >  2:607450/1‑301 (MQ=255)                                                                                                                                                       aCTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACT‑TGC‑ATTCGTG‑AATCCAGGC‑ACAGGCCCC‑GGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAGGGACACCTTGATCCCGTCCTCGTCCTTGGCCACTTGCCAGTTTTCAGCCCAGGCCACAGGCGACAGCAGAACACCCACACCACACAGCAGTGCAATACGCTTGAACGATCTCATCATGtt                 <  1:824519/301‑1 (MQ=255)                                                                                                                                                          gggTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACT‑TGC‑ATTCGTG‑AATCCAGGC‑ACAGGCCCC‑GGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAGGGACACCTTGATCCCGTCCTCGTCCTTGGCCACTTGCCAGTTTTCAGCCCAGGCCACAGGCGACAGCAGAACACCCACACCACACAGCAGTGCAATACGCTTGAACGATCTCATCATGTTTAt              <  1:194358/301‑1 (MQ=255)                                                                                                                                                                    ccAGCTCTGGTCACCTTCGTGCTTGAGCAGCTTCTGCGACT‑TGC‑ATTCGTG‑AATCCAGGC‑ACAGGCCCC‑GGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAGGGACACCTTGATCCCGTCCTCGTCCTTGGCCACTTGCCAGTTTTCAGCCCAGGCCACCGGCGACAGCAGAACAACCACACCACACAGCAGTGCAATAAGCTTGAACGATCTCATCATGTTTATCCTtgttgtt    >  1:625315/1‑300 (MQ=255)                                                                                                                                                                      ggcgctGGTCACCTTCGGGCTTGGG‑CAGCTTCGCCGAC‑TGTGATTCGTG‑AA‑TCAAGGCACAGGCCCC‑GGAAACGTCCTCCGGCAGGGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCAACCTAACTGAGGGACACCTTGATCCCGTCCTCGTCCTTGGCCACTGGCCAGTTTTCAGCCCAGGCCACAGGCGACAGCAGAACACCCACACCACACAGCAGTGCAATACGCTTGAACGATCTCATCATGTTTATCCTTGTTGTCGa  <  2:711976/297‑1 (MQ=255)                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                             CGTCGACTTTTGCCACCCGAACAAAACCCTTCTCTTCCGGCAGGTAAGTCGGCTCTTCTTTCAGGTTACGGGTCAGGCTGCCATCCGCCTCCTTGATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑AACCTCCGCCAAC‑TGCAATTCGTG‑AA‑CCAAGCAAAAGGCCCC‑GAAAAGGCCCTCC‑GCAGCCCCCCCAACCTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACTGAGGGACACCTTGATCCCGTCCTCGTCCTTGGCCACTTGCCAGTTTTCAGCCCAGGCCACAGGCGACAGCAGAACACCCACACCACACAGCAGTGCAATACGCTTGAACGATCTCATCATGTTTATCCTTGTTGTCGA  >  AE015451/4741039‑4741503

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