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%	107.1 / NA	30	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 (18/12);  total (18/12)

GGATGATGCCGGAGGTTTCCCAGTCCAGACGGTGCACGATCAGCGCATCGGGGTAGCCGTTTTCCTGCAGACGGGTGATCAGGCAGTCCTTGTTGTCCTCGGCACGGCCAGGTACCGACAGCAGCAGGGTCGGTTTGTTGATCACCAGGATGGCGGCGTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTG‑CCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACATATTTGCGGATCACTTCCAGCAGTTCTTTCTGCAGCG‑CGGGCAGGTAGTCCGGTTCGCTGCGTTGGCCGCGCTCGTGCGCCACGATGATCTGTAGACGCTCTTTCGCTACCGACGCGCTGGTCTGTTTCTGTCTGCCACGAAAG  >  AE015451/1932341‑1932920                                                                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                                                            ggATGATGCCGGAGGTTTCCCAGTCCAGACGGTGCACGATCAGCGCATCGGGGTAGCCGTTTTCCTGCAGACGGGTGATCAGGCAGTCCTTGTTGTCCTCGGCACGGCCAGGTACCGACAGCAGCAGGGTCGGTTTGTTGATCACCAGGATGGCGGCGTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTG‑ccc                                                                                                                                                                                                                                                                                          <  3:135315/301‑2 (MQ=255)                                    aCGATCAGCGCATCGGGGTAGCCGTTTTCCTGCAGACGGGTGATCAGGCAGTCCTTGTTGTCCTCGGCACGGCCAGGTACCGACAGCAGCAGGGTCGGTTTGTTGATCACCAGGATGGCGGCGTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCAGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACt                                                                                                                                                                                                                                                        >  3:147145/1‑301 (MQ=255)                                     ggacaggggcATCGGGGTAGCCGGTTTCCTGCAGAAGGGTGATCAGGCAGTCCTTGTTGTCCTCGGCACGGCCAGGTACCGACAGCAGCAGGGTCGGTTTGTTGATCACCAGGATGGCGGCGTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTa                                                                                                                                                                                                                                                       <  2:464272/293‑1 (MQ=255)                                                     tcagacgtgtgctcttccgatcttGATCAGGCAGTCCTTGTTGTCCTCGGCACGGCCAGGTACCGACAGCAGCAGGGTCGGTTTGTTGATCACCAGGATGGCGGCGTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATcc                                                                                                                                                                                                                                       <  4:533174/278‑1 (MQ=255)                                                     gTAGCCGTTTTCCTGCAGACGGGTGATCAGGCAGTCCTTGTTGTCCTCGGCACGGCCAGGTACCGACAGCAGCAGGGTCGGTTTGTTGATCACCAGGATGGCGGCGTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGAttc                                                                                                                                                                                                                                       >  1:309030/1‑299 (MQ=255)                                                           gTTTTCCTGCAGACGGGTGATCAGGCAGTCCTTGTTGTCCTCGGCACGGCCAGGTACCGACAGCAGCAGGGTCGGTTTGTTGATCACCAGGATGGCGGCGTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGCCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAgg                                                                                                                                                                                                                                 >  4:472203/1‑301 (MQ=255)                                                                                 aGGCAGTCCTTGTTGTCCTCGGCACGGCCAGGTACCGACAGCAGCAGGGTCGGTTTGTTGATCACCAGGATGGCGGCGTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCg                                                                                                                                                                                                           >  4:758508/1‑301 (MQ=255)                                                                                           tgttgtCCTCGGCACGGCCAGGTACCGACAGCAGCAGGGTCGGTTTGTTGATCACCAGGATGGCGGCGTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGcc                                                                                                                                                                                                 <  3:309029/301‑1 (MQ=255)                                                                                                          ggCCAGGTACGGACAGCAGCAGGTTCGGTTTGTTGATCACCAGGATGGCGGCGTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTc                                                                                                                                                                                  <  2:758229/301‑1 (MQ=255)                                                                                                                      aCAGCAGCAGGGTCGGTTTGTTGATCACCAGGATGGCGGCGTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCAcatc                                                                                                                                                                      >  3:116022/1‑301 (MQ=255)                                                                                                                        ggcagcagGGTCGGTTTGTTGATCACCAGGATGGCGGCGTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCAcatcat                                                                                                                                                                    <  1:116020/300‑1 (MQ=255)                                                                                                                        agcagcagGGTCGGTTTGTTGATCACCAGGATGGCGGCGTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGACCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCCAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGATCGATATGCAcataat                                                                                                                                                                    >  1:409465/1‑298 (MQ=255)                                                                                                                                  gcGGTTTGTTGATCACCAGGATGGCGGCGTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGACCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGAt                                                                                                                                                          <  3:409467/300‑1 (MQ=255)                                                                                                                                      ctcttccgatctcAGGATGGCGGCGTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTc                                                                                                                                                      <  1:395855/289‑1 (MQ=255)                                                                                                                                                  cAGGATGGCGGCGTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTcagatcggaagag                                                                                                                                          >  3:395855/1‑290 (MQ=255)                                                                                                                                                             cttCTTAGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCAGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACATATTTGCGGATCACTTCcag                                                                                                                               <  1:147143/299‑1 (MQ=255)                                                                                                                                                                                     tCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGCCTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGCATCGAGCAGCTGCACTGGTTTTACAGTTCGATCTGCCCATCATCGTTGCCGATGTTCACATATTTGCGGATCACTTCCAGCAGGTCTTTTCGCAGCG‑CGGGccg                                                                                                      >  1:318483/1‑299 (MQ=255)                                                                                                                                                                                                       gTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACATATTTGCGGATCACTTCCAGCAGTTCTTTCTGCAGCG‑CGGGCAGGTAGTCCGGTTCGCTGCg                                                                                    >  3:542888/1‑301 (MQ=255)                                                                                                                                                                                                               aaaaTAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTTCCGATGTTCCCATCTTTCCGGATCACTTCCATCCTTTTTTTCTGCA‑CGCCCGGCCTGTAGTCCGGTTCGCTTCGTTTGCcgc                                                                            >  2:257126/1‑301 (MQ=255)                                                                                                                                                                                                                              gCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCACGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCCGGGAGCTGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACATATTTGCGGATCAATTCCAGCAGTTCTTTCTGCAGCG‑CGGGCAGGTAGTACGGTTCGCTGCGTCGGCCGCGCTCGTCCGCCACgg                                                             >  3:331593/1‑300 (MQ=255)                                                                                                                                                                                                                                      ctacacgacgctcttccgatctaaGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACATATTTGCGGATCACTTCCAGCAGTTCTTTCTGCAGCG‑CGGGCAGGTAGTCCGGTTCGCTGCGTTGGCCGCGCTCGTGCGCCACGATGATCTGt                                                     <  1:206460/279‑1 (MQ=255)                                                                                                                                                                                                                                                        ctttccctacacgacgctcttccgatctcGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACATATTTGCGGATCACTTCCAGCAGTTCTTTCTGCAGCG‑CGGGCAGGTAGTCCGGTTCGCTGCGTTGGCCGCGCTCGTGCGCCACGATGATCTGTAGACGCTCTTTCGCTAcc                                   <  1:543143/273‑1 (MQ=255)                                                                                                                                                                                                                                                            aaGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACATATTTGCGGATCACTTCCAGCAGTTCTTTCTGCAGCG‑CGGGCAGGTAGTCCGGTTCGCTGCGTTGGCCGCGCTCGTGCGCCACGATGATCTGTAGatcggaagagcacacgtctg                               >  3:206458/1‑282 (MQ=255)                                                                                                                                                                                                                                                             ttccgatctcTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACATATTTGCGGATCACTTCCAGCAGTTCTTTCTGCAGCG‑CGGGCAGGTAGTCCGGTTCGCTGCGTTGGCCGCGCTCGTGCGCCACGATGATCTGTAGACGCTCTTTCGCTACCGAcgc                              <  3:266098/292‑1 (MQ=255)                                                                                                                                                                                                                                                                 cGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACATATTTGCGGATCACTTCCAGCAGTTCTTTCTGCAGCG‑CGGGCAGGTAGTCCGGTTCGCTGCGTTGGCCGCGCTCGTTCGCCACGATGATCTGTAGACGCTCTTTCGCTACCGACGCGCTg                          >  1:403817/1‑301 (MQ=255)                                                                                                                                                                                                                                                                     ccTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACATATTTGCGGATCACTTCCAGCAGTTCTTTCTGCAGCG‑CGGGCAGGTAGTCCGGTTCGCTGCGTTGGCCGCGCTCGTGCGCCACGATGATCTGTAGACGCTCTTTCGCTACCGACGCGCTGGTCt                      >  3:163875/1‑301 (MQ=255)                                                                                                                                                                                                                                                                      cTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACATATTTGCGGATCACTTCCAGCAGTTCTTTCTGCAGCG‑CGGGCAGGTAGTCCGGTTCGCTGCGTTGGCCGCGCTCGTGCGCCACGATGATCTGTAGACGCTCTTTCGCTACCGAcgcggatcggaa                     >  1:266099/1‑293 (MQ=255)                                                                                                                                                                                                                                                                             ggATATGCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACATATTTGCGGATCACTTCCAGCAGTTCTTTCTGCAGCG‑CGGGCAGGTAGTCCGGTTCGCTGCGTTGGCCGCGCTCGTGCGCCACGATGATCTGTAGACGCTCTTTCGCTACCGACGCGCTGGTCTGTTtctgt              >  1:441070/1‑301 (MQ=255)                                                                                                                                                                                                                                                                                 ggcgtccccggagggacctgagtggg‑CCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACATATTTGCGGATCACTTCCAGCAGTTCTTTCTGCAGCG‑CGGGCAGGTAGTCCGGTTCGCTGCGTTGGCCGCGCTCGTGCGCCACGATGATCTGTAGACGCTCTTTCGCTACCGACGCGCTGGTCTGTTTCTGTCTGcc         <  2:472069/276‑1 (MQ=255)                                                                                                                                                                                                                                                                                  ggCGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACATATTTGCGGATCACTTCCAGCAGTTCTTTCTGCAGCG‑CGGGCAGGTAGTCCGGTTCGCTGCGTTGGCCGCGCTCGTGCGCCACGATGATCTGTAGACGCTCTTTCGCTACCGACGCGCTGGTCTGTTTCTGTCTGcc         <  3:441072/300‑1 (MQ=255)                                                                                                                                                                                                                                                                                    cGTCCCCTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACATATTTGCGGATCACTTCCAGCAGTTCTTTCTGCAGCG‑CGGGCAGGTAGTCCGGTTCGCTGCGTTGGCCGCGCTCGTGCGCCACGATGATCTGTAGACGCTCTTTCGCTAccagatcggaagagcacacttctgaactcc       >  3:543145/1‑273 (MQ=255)                                                                                                                                                                                                                                                                                         ccTGTAGGAGCGGGCTTGCCCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTACAGCTCGATATGCACATCATCGTTGCCGATGTTCACATATTTGCGGATCACTTCCAGCAGTTCTTTCTGCAGCG‑CGGGCAGGTAGTCCGGTTCGCTGCGTTGGCCGCGCTCGTGCGCCACGATGATCTGTAGACGCTCTTTCGCTACCGACGCGCTGGTCTGTTTCTGTCTGCCACGAaag  >  4:245321/1‑301 (MQ=255)                                                                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                                                            GGATGATGCCGGAGGTTTCCCAGTCCAGACGGTGCACGATCAGCGCATCGGGGTAGCCGTTTTCCTGCAGACGGGTGATCAGGCAGTCCTTGTTGTCCTCGGCACGGCCAGGTACCGACAGCAGCAGGGTCGGTTTGTTGATCACCAGGATGGCGGCGTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTTG‑CCGCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACATATTTGCGGATCACTTCCAGCAGTTCTTTCTGCAGCG‑CGGGCAGGTAGTCCGGTTCGCTGCGTTGGCCGCGCTCGTGCGCCACGATGATCTGTAGACGCTCTTTCGCTACCGACGCGCTGGTCTGTTTCTGTCTGCCACGAAAG  >  AE015451/1932341‑1932920

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