Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	AE015451	4,741,258	A→G	intergenic (+495/+258)	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,258	0	A	G	96.9%	110.2 / ‑6.0	32	intergenic (+495/+258)	gltA/yeiW	citrate synthase/putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (16/15);  total (16/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

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: ‑