Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	AE015451	1,932,640	+C	intergenic (‑107/+43)	rluA ← / ← minE	23S rRNA (746) and tRNA (32) pseudouridine synthase/Cell division topological specificity factor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	AE015451	1,932,638	1	.	C	100.0%	138.8 / NA	39	intergenic (‑105/+45)	rluA/minE	23S rRNA (746) and tRNA (32) pseudouridine synthase/Cell division topological specificity factor
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base C (21/18);  total (21/18)

GGATGATGCCGGAGGTTTCCCAGTCCAGACGGTGCACGATCAGCGCATCGGGGTAGCCGTTTTCCTGCAGACGGGTGATCAGGCAGTCCTTGTTGTCCTCGGCACGGCCAGGTACCGACAGCAGCAGGGTCGGTTTGTTGATCACCAGGATGGCGGCGTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTG‑CCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACATATTTGCGGATCACTTCCAGCAGTTCTTTCTGCAGCGCGGGCAGGTAGTCCGGTTCGCTGCGTTGGCCGCGCTCGTGCGCCACGATGATCTGTAGACGCTCTTTCGCTACCGACGCGCTGGTCTGTTTCTGTCTGCCACGAAAGAAGTCAAAAAG  >  AE015451/1932341‑1932931                                                                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                                                                      gaagacggcggggggTTTCCCAGTCCGGACGGGGCACGATCGTCGCAGCGGGGTAGCCGTGTTCCTGCAGACGGGTGATCAGGCAGTCCTTGTTGGCCTCGGCACGGCCAGGTACCGACAGCAGCAGGGTCGGTTTGTTGATCACCAGGATGGCGGCGTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTG‑ccc                                                                                                                                                                                                                                                                                                    <  2:370346/288‑2 (MQ=255)          gggAGGTTTCCCAGTCCAGACGGTGCACGATCAGCGCATCGGGGTAGCCGTTTTCCTGCAGACGGGTGATCAGGCAGTCCTTGTTGTCCTCGGCACGGCCAGGTACCGACAGCAGCAGGGTCGGTTTGTTGATCACCAGGATGGCGGCGTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGc                                                                                                                                                                                                                                                                                            <  1:223596/300‑1 (MQ=255)             gggTTTACCAGACCCGACAGTGCACGGTCCGGGCATCGGGGTAGCCGTTGTCCTGCAGACGGGTGAGCAGGCAGGCCTTGTTGTCCTCGGCACGGCCAGGTACCGACAGCAGCAGGGTCGGTTTGTTGATCACCAGGATGGCGGCGTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGc                                                                                                                                                                                                                                                                                         <  2:372261/300‑1 (MQ=255)                  tCCCAGTCCAGACGGTGCACGATCAGCGCATCGGGGTAGCCGTTTTCCTGCAGACGGGTGATCAGGCAGTCCTTGTTGTCCTCGGCACGGCCAGGTACCGACAGCAGCAGGGTCGGTTTGTTGATCACCAGGATGGCGGCGTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCCAggg                                                                                                                                                                                                                                                                                    >  1:208097/1‑300 (MQ=255)                                                           gTTTTCCTGCAGACGGGTGATCAGGCAGTCCTTGTTGTCCTCGGCACGGCCAGGTACCGACAGCAGCAGGGTCGGTTTGTTGATCACCAGGATGGCGGCGTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAgg                                                                                                                                                                                                                                           >  3:171902/1‑301 (MQ=255)                                                             tttCCTGCAGACGGGTGATCAGGCAGTCCTTGTTGTCCTCGGCACGGCCAGGTACCGACAGCAGCAGGGTCGGTTTGTTGATCACCAGGATGGCGGCGTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGt                                                                                                                                                                                                                                         <  4:12925/301‑1 (MQ=255)                                                                                 aGGCAGTCCTTGTTGTCCTCGGCACGGCCAGGTACCGACAGCAGCAGGGTCGGTTTGTTGATCACCAGGATGGCGGCGTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCCGCGATCCGGCAGGGTGATGTTCAGGTCCAggaccc                                                                                                                                                                                                                     >  2:204705/1‑296 (MQ=255)                                                                                     aGTCCTTGTTGTCCTCGGCACGGCCAGGTACCGACAGCAGCAGGGTCGGTTTGTTGATCACCAGGATGGCGGCGTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGagca                                                                                                                                                                                                                 >  3:513470/1‑301 (MQ=255)                                                                                                               ggTACCGACAGCAGCAGGGTCGGTTTGTTGATCACCAGGATGGCGGCGTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGatat                                                                                                                                                                                       <  3:208099/301‑1 (MQ=255)                                                                                                                     gcagggagaagGGTCGGTTTGTTGATCACCAGGATGGCGGCGTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGACCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAAGTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCAcat                                                                                                                                                                                 <  2:626330/292‑1 (MQ=255)                                                                                                                                  gcGGTTTGTTGATCACCAGGATGGCGGCGTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGAt                                                                                                                                                                    <  1:171896/300‑1 (MQ=255)                                                                                                                                               cACCAGGATGGCGGCGTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACATATTTg                                                                                                                                                       <  1:513456/301‑1 (MQ=255)                                                                                                                                                       tGGCGGCGTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTACAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACATATTTGCGGATCAc                                                                                                                                               >  1:567965/1‑301 (MQ=255)                                                                                                                                                                  tCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACATATTTGCGGATCACTTCCAGCAGtt                                                                                                                                    <  4:204743/301‑1 (MQ=255)                                                                                                                                                                      agagGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACATATTTGCGGATCACTTCCAGCAGttcttt                                                                                                                                >  3:317463/1‑301 (MQ=255)                                                                                                                                                                        gggATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACATATTTGCGGATCACTTCCAGCAGttctttct                                                                                                                              <  1:317455/300‑1 (MQ=255)                                                                                                                                                                                 aCATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACATATTTGCGGATCACTTCCAGCAGTTCTTTCTGCATCGCgg                                                                                                                     >  3:520837/1‑301 (MQ=255)                                                                                                                                                                                          cttccgatctcAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACATATTTGCGGATCACTTCCAGCAGTTCTTTCTGCAGCGCGGGCAGGTAGt                                                                                                            <  4:703360/291‑1 (MQ=255)                                                                                                                                                                                            cGGCATCGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACATATTTGCGGATCACTTCCAGCAGTTCTTTCTGCAGCGCGGGCAGGTAGTcc                                                                                                          <  4:716282/301‑1 (MQ=255)                                                                                                                                                                                             ggCATTGCAGGGCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACATATTTGCGGATCACTTCCAGCAGTTCTTTCTGCAGCGCGGGCAGGTAGTCCg                                                                                                         <  1:157030/301‑1 (MQ=255)                                                                                                                                                                                               cATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACATATTTGCGGATCACTTCCAGCAGTTCTTTCTGCAGCGCGGGCAGGTAGTCCGGt                                                                                                       >  4:159104/1‑301 (MQ=255)                                                                                                                                                                                                    cAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACATTTTTGCGGATCACTTCCAGCAGTTCTTTCTGCAGCGCGGGCAGgtcgtcgaccggaag                                                                                                  >  2:703178/1‑292 (MQ=255)                                                                                                                                                                                                      ggTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACATATTTGCGGATCACTTCCAGCAGTTCTTTCTGCAGCGCGGGCAGGTAGTCCGGTTCGCTGc                                                                                                >  3:280024/1‑301 (MQ=255)                                                                                                                                                                                                               aaaaTAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACATATTTGCGGATCACTTCCAGCAGTTCTTTCTGCAGCGCGGGCAGGTAGTCCGGTTCGCTGCGTTGGCcgc                                                                                       >  4:569800/1‑301 (MQ=255)                                                                                                                                                                                                                               cTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACATATTTGCGGATCACTTCCAGCAGTTCTTTCTGCAGCGCGGGCAGGTAGTCCGGTTCGCTGCGTTGGCCGCGCTCGTGCGCCACgat                                                                       <  1:280019/301‑1 (MQ=255)                                                                                                                                                                                                                                              ggCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACATATTTGCGGATCACTTCCAGCAGTTCTTTCTGCAGCGCGGGCAGGTAGTCCGGTTCGCTGCGTTGGCCGCGCTCGTGCGCCACGATTATCTGTAGACGctc                                                        >  1:531312/1‑301 (MQ=255)                                                                                                                                                                                                                                                cccTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACATATTTGCGGATCACTTCCAGCAGTTCTTTCTGCAGCGCGGGCAGGTAGTCCGGTTCGCTGCGTTGGCCGCGCTCGTGCGCCACGCGGATCTTTAGaagccttt                                                      >  2:764239/1‑294 (MQ=255)                                                                                                                                                                                                                                                 cctgcgcgGGAAAGCCCGCTCCTACAAAGGATATTCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACATATTTGCGGATCACTTCCAGCAGTTCTTTCTGCAGCGCGGGCAGGTAGTCCGGTTCGCTGCGTTGGCCGCGCTCGTGCGCCACGATGATCTGTAGACGCTCttt                                                     <  2:159074/297‑1 (MQ=255)                                                                                                                                                                                                                                                              gCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACATATTTGCGGATCACTTCCAGCAGTTCTTTCTGCAGCGCGGGCAGGTAGTCCGGTTCGCTGCGTTGGCCGCGCTCGTGCGCCACGATGATCTGTAGACGCTCTTTCGCTACCGAcgcg                                        >  3:446935/1‑301 (MQ=255)                                                                                                                                                                                                                                                                           aaGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACATATTTGCGGATCACTTCCAGCAGTTCTTTCTGCAGCGCGGGCAGGTAGTCCGGTTCGCTGCGTTGGCCGCGCTCGTGCGCCACGATGAACTGTAGACGCTCTTTCGCTACCGACGCGCTGGTCTTTTtct                           >  2:354508/1‑301 (MQ=255)                                                                                                                                                                                                                                                                            aGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACATATTTGCGGATCACTTCCAGCAGTTCTTTCTGCAGCGCGGGCAGGTAGTCCGGTTCGCTGCGTTGGCCGCGCTCGTGCGCCACGATGATCTGTAGACGCTCTTTCGCTACCGACGCGCTGGTCTGTTtctg                          <  1:520823/301‑1 (MQ=255)                                                                                                                                                                                                                                                                            aGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACATATTTGCGGATCACTTCCAGCAGTTCTTTCTGCAGCGCGGGCAGGTAGTCCGGTTCGCTGCGTTGGCCGCGCTCGTGCGCCACGATGATCTGTAGACGCTCTTTCGCTACCGACGCGCTGGTCTGTTtctg                          <  4:354576/301‑1 (MQ=255)                                                                                                                                                                                                                                                                                 atGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACATATTTGCGGATCACTTCCAGCAGTTCTTTCTGCAGCGCGGGCAGGTAGTCCGGTTCGCTGCGTTGGCCGCGCTCGTGCGCCACGATGATCTGTAGACGCTCTTTCGCTACCGACGCGCTGGTCTGTTTCTGTCTGc                     >  3:372492/1‑301 (MQ=255)                                                                                                                                                                                                                                                                                    cGTCCCCGGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACATATTTGCGGATCACTTCCAGCAGTTCTTTCTGCAGCGCGGGCAGGTAGTCCGGTTCGCTGCGTTGGCCGCGCTCGTGCGCCACGATGATCTGTAGACGCTCTTTCGCTACCGACGCGCTGGTCTGTTTCTGTCTGCCat                  <  1:372482/301‑2 (MQ=255)                                                                                                                                                                                                                                                                                      tCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACATATTTGCGGATCACTTCCAGCAGTTCTTTCTGCAGCGCGGGCAGGTAGTCCGGTTCGCTGCGTTGGCCGCGCTCGTGCGCCACGATGATCTGTAGACGCTCTTTCGCTACCGACGCGCTGGTCTGTTTCTGTCTGCCACgg                >  1:67762/1‑300 (MQ=255)                                                                                                                                                                                                                                                                                        cccTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACATATTTGCGGATCACTTCCAGCAGTTCTTTCTGCAGCGCGGGCAGGTAGTCCGGTTCGCTGCGTTGGCCGCGCTCGTGCGCCACGATGATCTGTAGACGCTCTTTCGCTACCGACGCGCTGGTCTGTTTCTGTCTGCCACGaaa              <  2:569657/301‑1 (MQ=255)                                                                                                                                                                                                                                                                                           tGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACATATTTGCGGATCACTTCCAGCAGTTCTTTCTGCAGCGCGGGCAGGTAGTCCGGTTCGCTGCGTTGGCCGCGCTCGTGCGCCACGATGATCTGTAGACGCCCTTTCGCTACCGACGCCATGGTCTGTTTCGTTTTGCCACGaaatga           >  2:618691/1‑298 (MQ=255)                                                                                                                                                                                                                                                                                               ggAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACATATTTGCGGATCACTTCCAGCAGTTCTTTCTGCAGCGCGGGCAGGTAGTCCGGTTCGCTGCGTTGGCCGCGCTCGTGCGCCACGATGATCTGTAGACGCTCTTTCGCTACCGACGCGCTGGTCTGTTTCTGTCTGCCACGAAAGAAGTCa       >  1:84112/1‑301 (MQ=255)                                                                                                                                                                                                                                                                                                    gggCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACATATTTGCGGATCACTTCCAGCAGTTCTTTCTGCAGCGCGGGCAGGTAGTCCGGTTCGCTGCGTTGGCCGCGCTCGTGCGCCACGATGATCTGTAGACGCTCTTTCGCTACCGACGCGCTGGTCTGTTTCTGTCTGCCACGAAAGAAGTCAAAAAg  >  3:453734/1‑301 (MQ=255)                                                                                                                                                                                                                                                                                                    gggCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACATATTTGCGGATCACTTCCAGCAGTTCTTTCTGCAGCGCGGGCAGGTAGTCCGGTTCGCTGCGTTGGCCGCGCTCGTGCGCCACGATGATCTGTAGACGCTCTTTCGCTACCGACGCGCTGGCCTGTTTCTTTCTGCCACGCAAGAAGTCaaaata  >  2:197763/1‑299 (MQ=255)                                                                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                                                                      GGATGATGCCGGAGGTTTCCCAGTCCAGACGGTGCACGATCAGCGCATCGGGGTAGCCGTTTTCCTGCAGACGGGTGATCAGGCAGTCCTTGTTGTCCTCGGCACGGCCAGGTACCGACAGCAGCAGGGTCGGTTTGTTGATCACCAGGATGGCGGCGTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTG‑CCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACATATTTGCGGATCACTTCCAGCAGTTCTTTCTGCAGCGCGGGCAGGTAGTCCGGTTCGCTGCGTTGGCCGCGCTCGTGCGCCACGATGATCTGTAGACGCTCTTTCGCTACCGACGCGCTGGTCTGTTTCTGTCTGCCACGAAAGAAGTCAAAAAG  >  AE015451/1932341‑1932931

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