Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	1,790,705	C→T	intergenic (+41/‑57)	unknown → / → unknown	tRNA 2‑thiouridine synthesis protein TusE/Manganese ABC transporter, periplasmic‑binding protein SitA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	1,790,705	0	C	T	100.0%	197.8 / NA	70	intergenic (+41/‑57)	unknown/unknown	tRNA 2‑thiouridine synthesis protein TusE/Manganese ABC transporter, periplasmic‑binding protein SitA
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (33/37);  total (33/37)

CCCATTGGGAAATTATCCATTTTGTGCGTGATTTTTATCAAGAATATAAAACCTCTCCTGCCATTCGTATGTTGGTGAAAGCCACCGCGCAAAAGCTGGGAGAAGATAAAGGCAATAGCCGTTATTTGCAACGTTTATTTCCCGATGGTCCGGCAAAACAGGCAACCAAGCTGGCGGGTCTCCCCAAACCGGTAAAATGTCTATAATTTTTTTGAGAATCATTTGCATTTAGAAATAGAAAGAATACAATGTCCGCGCAGCTTATGTTGTTATTTTATAGTTAAAAAGCCAAGGAGAAAAAATATGCGTTTATTAGCCACATTAGTCAGCGTACTTGCGTTAAGCTTAAGCAGTAGTCTTGCTTATGCGAAATTTAAAGTTGTCACTACGTTCACTATTATTCAAGATATGGCACAAAATGTCGCCGGTGATGCGGCAACCGTAGAGTCCATCACCAAGCCCGGAGCAGAAATCCATGATTACGAGCCAAC  >  1/1790459‑1790949                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                     cccATTGGGAAATTATCCATTTTGTGCGTGATTTTTATCAAGAATATAAAACCTCTCCTGCCATTCGTATGTTGGTGAAAGCCACCGCGCAAAAGCTGGGAGAAGATAAAGGCAATAGCCGTTATTTGCAACGTTTATTTCCCGATGGTCCGGCAAAACAGGCAACCAAGCTGGCGGGTCTCCCCAAACCGGTAAAATGTCTATAATTTTTTTGAGAATCATTTGCATTTAGAAATAGAAAGAATATAAt                                                                                                                                                                                                                                                   <  1:408336/250‑1 (MQ=255) cccATTGGGAAATTATCCATTTTGTGCGTGATTTTTATCAAGAATATAAAACCTCTCCTGCCATTCGTATGTTGGTGAAAGCCACCGCGCAAAAGCTGGGAGAAGATAAAGGCAATAGCCGTTATTTGCAACGTTTATTTCCCGATGGTCCGGCAAAACAGGCAACCAAGCTGGCGGGTCTCCCCAAACCGGTAAAATGTCTATAATTTTTTTGAGAATCATTTGCATTTAGAAATAGAAAGAATATAATg                                                                                                                                                                                                                                                  <  1:397093/251‑1 (MQ=255) cccATTGGGAAATTATCCATTTTGTGCGTGATTTTTATCAAGAATATAAAACCTCTCCTGCCATTCGTATGTTGGTGAAAGCCACCGCGCAAAAGCTGGGAGAAGATAAAGACAATAGCCGTTATTTGCAACGTTTATTTCCCGATGGTCCGGCAAAACAGGCAACCAAGCTGGCGGGTCTCCCCAAACCGGTAAAATGTCTATAATTTTTTTGAGAATCATTTGCATTTAGAAATAGAAAGAATATAATg                                                                                                                                                                                                                                                  <  1:88184/251‑1 (MQ=255)   cATTGGGAAATTATCCATTTTGTGCGTGATTTTTATCAAGAATATAAAACCTCTCCTGCCATTCGTATGTTGGTGAAAGCCACCGCGCAAAAGCTGGGAGAAGATAAAGGCAATAGCCGTTATTTGCAACGTTTATTTCCCGATGGTCCGGCAAAACAGGCAACCAAGCTGGCGGGTCTCCCCAAACCGGTAAAATGTCTATAATTTTTTTGAGAATCATTTGCATTTAGAAATAGAAAGAATATAATGTc                                                                                                                                                                                                                                                <  1:364958/251‑1 (MQ=255)                tCCATTTTGTGCGTGATTTTTATCAAGAATATAAAACCTCTCCTGCCATTCGTATGTTGGTGAAAGCCACCGCGCAAAAGCTGGGAGAAGATAAAGGCAATAGCCGTTATTTGCAACGTTTATTTCCCGATGGTCCGGCAAAACAGGCAACCAAGCTGGCGGGTCTCCCCAAACCGGTAAAATGTCTATCATTTTTTTGAGAATCATTTGCATTTAGAAATAGAAAGAATATAATg                                                                                                                                                                                                                                                  >  1:705205/1‑236 (MQ=255)                  cATTTTGTGCGTGATTTTTATCAAGAATATAAAACCTCTCCTGCCATTCGTATGTTGGTGAAAGCCACCGCGCAAAAGCTGGGAGAAGATAAAGGCAATAGCCGTTATTTGCAACGTTTATTTCCCGATGGTCCGGCAAAACAGGCAACCAAGCTGGCGGGTCTCCCCAAACCGGTAAAATGTCTATAATTTTTTTGAGAATCATTTGCATTTAGAAATAGAAAGAATATAATGTCCGCGCAGCTTAtgtt                                                                                                                                                                                                                                 <  1:123781/251‑1 (MQ=255)                          gCGTGATTTTTATCAAGAATATAAAACCTCTCCTGCCATTCGTATGTTGGTGAAAGCCACCGCGCAAAAGCTGGGAGAAGATAAAGGCAATAGCCGTTATTTGCAACGTTTATTTCCCGATGGTCCGGCAAAACAGGCAACCAAGCTGGCGGGTCTCCCCAAACCGGTAAAATGTCTATAATTTTTTTGAGAATCATTTGCATTTAGAAATAGAAAGAATATAAt                                                                                                                                                                                                                                                   >  1:155553/1‑225 (MQ=255)                           cGTGATTTTTATCAAGAATATAAAACCTCTCCTGCCATTCGTATGTTGGTGAAAGCCACCGCGCAAAAGCTGGGAGAAGATAAAGGCAATAGCCGTTATTTGCAACGTTTATTTCCCGATGGTCCGGCAAAACAGGCAACCAAGCTGGCGGGTCTCCCCAAACCGGTAAAATGTCTATAATTTTTTTGAGAATCATTTGCATTTAGAAATAGAAAGAATATAATGTCCGCGCAGCTTATGTTGTTAtttt                                                                                                                                                                                                                         <  1:365180/250‑1 (MQ=255)                               aTTTTTATCAAGAATATAAAACCTCTCCTGCCATTCGTATGTTGGTGAAAGCCACCGCGCAAAAGCTGGGAGAAGATAAAGGCAATAGCCGTTATTTGCAACGTTTATTTCCCGATGGTCCGGCAAAACAGGCAACCAAGCTGGCGGGTCTCCCCAAACCGGTAAAATGTCTATAATTTTTTTGAGAATCATTTGCATTTAGAAATAGAAAGAATATAAt                                                                                                                                                                                                                                                   <  1:780831/220‑1 (MQ=255)                                            atatAAAACCTCTCCTGCCATTCGTATGTTGGTGAAAGCCACCGCGCAAAAGCTGGGAGAAGATAAAGGCAATAGCCGTTATTTGCAACGTTTATTTCCCGATGGTCCGGCAAAACAGGCAACCAAGCTGGCGGGTCTCCCCAAACCGGTAAAATGTCTATAATTTTTTTGAGAATCATTTGCATTTAGAAATAGAAAGAATATAATGTCCGCGCAGCTTATGTTGTTATTTTATAGTTAAAAAGCCAAgg                                                                                                                                                                                                       >  1:836020/1‑251 (MQ=255)                                            atatAAAACCTCTCCTGCCATTCGTATGTTGGTGAAAGCCACCGCGCAAAAGCTGGGAGAAGATAAAGGCAATAGCCGTTATTTGCAACGTTTATTTCCCGATGGTCCGGCAAAACAGGCAACCAAGCTGGCGGGTCTCCCCAAACCGGTAAAATGTCTATAATTTTTTTGAGAATCATTTGCATTTAGAAATAGAAAGAATATAATGTCCGCGCAGCTTATGTTGTTATTTTATAGTTAAAAAGCCAAgg                                                                                                                                                                                                       >  1:341817/1‑251 (MQ=255)                                            atatAAAACCTCTCCTGCCATTCGTATGTTGGTGAAAGCCACCGCGCAAAAGCTGGGAGAAGATAAAGGCAATAGCCGTTATTTGCAACGTTTATTTCCCGATGGTCCGGCAAAACAGGCAACCAAGCTGGCGGGTCTCCCCAAACCGGTAAAATGTCTATAATTTTTTTGAGAATCATTTGCATTTAGAAATAGAAAGAATATAATGTCCGCGCAGCTTATGTTGTTATTTTATAGTTAAAAAGCCAAgg                                                                                                                                                                                                       >  1:35082/1‑251 (MQ=255)                                             tataAAACCTCTCCTGCCATTCGTATGTTGGTGAAAGCCACCGCGCAAAAGCTGGGAGAAGCTAAAGGCAATAGCCGTTATTTGCAACGTTTATTTCCCGATGGTCCGGCAAAACAGGCAACCAAGCTGGCGGGTCTCCCCAAACCGGTAAAATGTCTATAATTTTTTTGAGAATCATTTGCATTTAGAAATAGAAAGAATATAATGTCCGCGCAGCTTATGTTGTTATTTTATAGTTAAAAAGCCAAGga                                                                                                                                                                                                      >  1:568641/1‑251 (MQ=255)                                                 aaaCCTCTCCTGCCATTCGTATGTTGGTGAAAGCCACCGCGCAAAAGCTGGGAGAAGATAAAGGCAATAGCCGTTATTTGCAACGTTTATTTCCCGATGGTCCGGCAAAACAGGCAACCAAGCTGGCGGGTCTCCCCAAACCGGTAAAATGTCTATAATTTTTTTGAGAATCATTTGCATTTAGAAATAGAAAGAATATAATGTCCGCGCAGCTTATGTTGTTATTTTATAGTTAAAAAGc                                                                                                                                                                                                            <  1:302730/241‑1 (MQ=255)                                                    ccTCTCCTGCCATTCGTATGTTGGTGAAAGCCACCGCGCAAAAGCTGGGAGAAGATAAAGGCAATAGCCGTTATTTGCAACGTTTATTTCCCGATGGTCCGGCAAAACAGGCAACCAAGCTGGCGGGTCTCCCCAAACCGGTAAAATGTCTATAATTTTTTTGAGAATCATTTGCATTTAGAAATAGAAAGAATATAATGTCCGCGCAGCTTATGTTGTTATTTTATAGTTAAAAAGCCAAGGAGaaaaaa                                                                                                                                                                                               <  1:537832/251‑1 (MQ=255)                                                       ctcCTGCCATTCGTATGTTGGTGAAAGCCACCGCGCAAAAGCTGGGAGAAGATAAAGGCAATAGCCGTTATTTGCAACGTTTATTTCCCGATGGTCCGGCAAAACAGGCAACCAAGCTGGCGGGTCTCCCCAAACCGGTAAAATGTCTATAATTTTTTTGAGAATCATTTGCATTTAGAAATAGAAAGAATATAATGTCCGCGCAGCTTATGTTGTTATTTTATAGTTAAAAAGc                                                                                                                                                                                                            >  1:651478/1‑235 (MQ=255)                                                       ctcCTGCCATTCGTATGTTGGTGAAAGCCACCGCGCAAAAGCTGGGAGAAGATAAAAGCAATAGCCGTTATTTGCATCGTTTATTTCCCGATGGTCCGGCAAAACAGGCAACCAAGCTGGCGGGTCTCCCCAAACCGGTAAAATGTCTATAATTTTTTTGAGAATCATTTGCATTTAGAAATAGAAAGAATATAATg                                                                                                                                                                                                                                                  <  1:55102/197‑1 (MQ=255)                                                            gCCATTCGTATGTTGGTGAAAGCCACCGCGCAAAAGCTGGGAGAAGATAAAGGCAATAGCCGTTATTTGCAACGTTTATTTCCCGATGGTCCGGCAAAACAGGCAACCAAGCTGGCGGGTCTCCCCAAACCGGTAAAATGTCTATAATTTTTTTGAGAATCATTTGCATTTAGAAATAGAAAGAATATAATGTCCGCGCAGCTTATGTTGTTATTTTATAGTTAAAAAGCCAAGGAGaaaaaa                                                                                                                                                                                               >  1:381526/1‑243 (MQ=255)                                                              cATTCGTATGTTGGTGAAAGCCACCGCGCAAAAGCTGGGAGAAGATAAAGGCAATAGCCGTTATTTGCAACGTTTATTTCCCGATGGTCCGGCAAAACAGGCAACCAAGCTGGCGGGTCTCCCCAAACCGGTAAAATGTCTATAATTTTTTTGAGAATCATTTGCATTTAGAAATAGAAAGAATATAATGTCCGCGCAGCTTATGTTGTTATTTTATAGTTAAAAAGCCAAGGAGAAAAAATATGCGTtta                                                                                                                                                                                     <  1:850156/251‑1 (MQ=255)                                                               aTTCGTATGTTGGTGAAAGCCACCGCGCAAAAGCTGGGAGAAGATAAAGGCAATAGCCGTTATTTGCAACGTTTATTTCCCGATGGTCCGGCAAAACAGGCAACCAAGCTGGCGGGTCTCCCCAAACCGGTAAAATGTCTATAATTTTTTTGAGAATCATTTGCATTTAGAAATAGAAAGAATATAAt                                                                                                                                                                                                                                                   >  1:765918/1‑188 (MQ=255)                                                               aTTCGTATGTTGGTGAAAGCCACCGCGCAAAAGCTGGGAGAAGATAAAGGCAATAGCCGTTATTTGCAACGTTTATTTCCCGATGGTCCGGCAAAACAGGCAACCAAGCTGGCGGGTCTCCCCAAACCGGTAAAATGTCTATAATTTTTTTGAGAATCATTTGCATTTAGAAATAGAAAGAATATAATGTCCGCGCAGCTTATGTTGTTATTTTATAGTTAAAAAGCCAAGGAGAAAAAATATGCGTttat                                                                                                                                                                                    >  1:64543/1‑251 (MQ=255)                                                               aTTCGTATGTTGGTGAAAGCCACCGCGCAAAAGCTGGGAGAAGATAAAGGCAATAGCCGTTATTTGCAACGTTTATTTCCCGATGGTCCGGCAAAACAGGCAACCAAGCTGGCGGGTCTCCCCAAACCGGTAAAATGTCTATAATTTTTTTGAGAATCATTTGCATTTAGAAATAGAAAGAATATAATGTCCGCGCAGCTTATGTTGTTATTTTATAGTTAAAAAGCCAAGGAGAAAAAATATGCGTttat                                                                                                                                                                                    >  1:419654/1‑251 (MQ=255)                                                                 tCGTATGTTGGTGAAAGCCACCGCGCAAAAGCTGGGAGAAGATAAAGGCAATAGCCGTTATTTGCAACGTTTATTTCCCGATGGTCCGGCAAAACAGGCAACCAAGCTGGCGGGTCTCCCCAAACCGGTAAAATGTCTATAATTTTTTTGAGAATCATTTGCATTTAGAAATAGAAAGAATATAATGTCCGCGCAGCTTATGTTGTTATTTTATAGTTAAAAAGCCAAGGAGAAAAAATATGCGTttatta                                                                                                                                                                                  >  1:172854/1‑251 (MQ=255)                                                                     aTGTTGGTGAAAGCCACCGCGCAAAAGCTGGGAGAAGATAAAGGCAATAGCCGTTATTTGCAACGTTTATTTCCCGATGGTCCGGCAAAACAGGCAACCAAGCTGGCGGGTCTCCCCAAACCGGTAAAATGTCTATAATTTTTTTGAGAATCATTTGCATTTAGAAATAGAAAGAATATAATGTCCGCGCAGCTTATGTTGTTATTTTATAGTTAAAAAGCCAAGGAGAAAAAATATGCGTTTATTAGCca                                                                                                                                                                              >  1:747392/1‑251 (MQ=255)                                                                         tGGTGAAAGCCACCGCGCAAAAGCTGGGAGAAGATAAAGGCAATAGCCGTTATTTGCAACGTTTATTTCCCGATGGTCCGGCAAAACAGGCAACCAAGCTGGCGGGTCTCCCCAAACCGGTAAAATGTCTATAATTTTTTTGAGAATCATTTGCATTTAGAAATAGAAAGAATATAATGTCCgcgc                                                                                                                                                                                                                                           >  1:557563/1‑186 (MQ=255)                                                                                  ccACCGCGCAAAAGCTGGGAGAAGATAAAGGCAATAGCCGTTATTTGCAACGTTTATTTCCCGATGGTCCGGCAAAACAGGCAACCAAGCTGGCGGGTCTCCCCAAACCGGTAAAATGTCTATAATTTTTTTGAGAATCATTTGCATTTAGAAATAGAAAGAATATAATGTCCGCGCAGCTTATGTTGTTATTTTATAGTTAAAAAGCCAAGGAGAAAAAATATGCGTTTATTAGCCACATTAGTCAGCGt                                                                                                                                                                 >  1:266509/1‑251 (MQ=255)                                                                                                cTGGGAGAAGATAAAGGCAATAGCCGTTATTTGCAACGTTTATTTCCCGATGGTCCGGCAAAACAGGCAACCAAGCTGGCGGGTCTCCCCAAACCGGTAAAATGTCTATAATTTTTTTGAGAATCATTTGCATTTAGAAATAGAAAGAATATAATGTCCGCGCAGCTTATGTTGTTATTTTATAGTTAAAAAGCCAAGGAGAAAAAATATGCGTTTATTAGCCACATTAGTCAGCGTACTTGCGttaagct                                                                                                                                                   <  1:8668/251‑1 (MQ=255)                                                                                                cTGGGAGAAGATAAAGGCAATAGCCGTTATTTGCAACGTTTATTTCCCGATGGTCCGGCAAAACAGGCAACCAAGCTGGCGGGTCTCCCCAAACCGGTAAAATGTCTATAATTTTTTTGAGAATCATTTGCATTTAGAAATAGAAAGAATATAATGTCCGCGCAGCTTATGTTGTTATTTTATAGTTAAAAAGCCAAGGAGAAAAAATATGCGTTTATTAGCCACATTAGTCAGCGTACTTGCGttaagct                                                                                                                                                   <  1:263752/251‑1 (MQ=255)                                                                                                cTGGGAGAAGATAAAGGCAATAGCCGTTATTTGCAACGTTTATTTCCCGATGGTCCGGCAAAACAGGCAACCAAGCTGGCGGGTCTCCCCAAACCGGTAAAATGTCTATAATTTTTTTGAGAATCATTTGCATTTAGAAATAGAAAGAATATAATGTCCGCGCAGCTTATGTTGTTATTTTATAGTTAAAAAGCCAAGGAGAAAAAATATGCGTTTATTAGCCACATTAGTCAGCGTACTTGCGttaagc                                                                                                                                                    <  1:849960/250‑1 (MQ=255)                                                                                                cTGGGAGAAGATAAAGGCAATAGCCGTTATTTGCAACGTTTATTTCCCGATGGTCCGGCAAAACAGGCAACCAAGCTGGCGGGTCTCCCCAAACCGGTAAAATGTCTATAATTTTTTTGAGAATCATTTGCATTTAGAAATAGAAAGAATATAATGTCCGCGCAGCTTATGTTGTTATTTTATAGTTAAAAAGCCAAGGAGAAAAAATATGCGTTTATTAGCCACATTAGTCAGCGTACTTGCGttaagc                                                                                                                                                    <  1:582697/250‑1 (MQ=255)                                                                                                        agaTAAAGGCAATAGCCGTTATTTGCAACGTTTATTTCCCGATGGTCCGGCAAAACAGGCAACCAAGCTGGCGGGTCTCCCCAAACCGGTAAAATGTCTATAATTTTTTTGAGAATCATTTGCATTTAGAAATAGAAAGAATATAATGTCCGCGCAGCTTATGTTGTTATTTTATAGTTAAAAAGCCAAGGAGAAAAAATATGCGTTTATTAGCCACATTAGTCAGCGt                                                                                                                                                                 <  1:658978/229‑1 (MQ=255)                                                                                                        agaTAAAGGCAATAGCCGTTATTTGCAACGTTTATTTCCCGATGGTCCGGCAAAACAGGCAACCAAGCTGGCGGGTCTCCCCAAACCGGTAAAATGTCTATAATTTTTTTGAGAATCATTTGCATTTAGAAATAGAAAGAATATAATGTCCGCGCAGCTTATGTTGTTATTTTATAGTTAAAAAGCCAAGGAGAAAAAATATGCGTTTATTAGCCACATTAGTCAGCGTACTTGCGTTAAGCTTAAGCagt                                                                                                                                           >  1:366663/1‑251 (MQ=255)                                                                                                         gaTAAAGGCAATAGCCGTTATTTGCAACGTTTATTTCCCGATGGTCCGGCAAAACAGGCAACCAAGCTGGCGGGTCTCCCCAAACCGGTAAAATGTCTATAATTTTTTTGAGAATCATTTGCATTTAGAAATAGAAAGAATATAATGTCCGCGCAGCTTAtgt                                                                                                                                                                                                                                  >  1:688041/1‑163 (MQ=255)                                                                                                         gaTAAAGGCAATAGCCGTTATTTGCAACGTTTATTTCCCGATGGTCCGGCAAAACAGGCAACCAAGCTGGCGGGTCTCCCCAAACCGGTAAAATGTCTATAATTTTTTTGAGAATCATTTGCATTTAGAAATAGAAAGAATATAATGTCCGCGCAGCTTATGTTGTTATTTTATAGTTAAAAAGc                                                                                                                                                                                                            <  1:228287/185‑1 (MQ=255)                                                                                                         gaTAAAGGCAATAGCCGTTATTTGCAACGTTTATTTCCCGATGGTCCGGCAAAACAGGCAACCAAGCTGGCGGGTCTCCCCAAACCGGTAAAATGTCTATAATTTTTTTGAGAATCATTTGCATTTAGAAATAGAAAGAATATAATGTCCGCGCAGCTTATGTTGTTATTTTATAGTTAAAAAGCCAAGGAGAAAAAATATGCGTTTATTAGCcac                                                                                                                                                                             <  1:156061/216‑1 (MQ=255)                                                                                                          aTAAAGGCAATAGCCGTTATTTGCAACGTTTATTTCCCGATGGTCCGGCAAAACAGGCAACCAAGCTGGCGGGTCTCCCCAAACCGGTAAAATGTCTATAATTTTTTTGAGAATCATTTGCATTTAGAAATAGAAAGAATATAATGTCCGCGCAGCTTATGTTGTTATTTTATAGTTAAAAAGCCAAGGAGAAAAAATATGCGTTTATTAGCCACATTAGTCAGCGTACTTGCGTTAAGCTTAAGCagtag                                                                                                                                         <  1:206439/251‑1 (MQ=255)                                                                                                          aTAAAGGCAATAGCCGTTATTTGCAACGTTTATTTCCCGATGGTCCGGCAAAACAGGCAACCAAGCTGGCGGGTCTCCCCAAACCGGTAAAATGTCTATAATTTTTTTGAGAATCATTTGCATTTAGAAATAGAAAGAATATAATGTCCGCGCAGCTTATGTTGTTATTTTATAGTTAAAAAGCCAAGGAGAAAAAATATGCGTTTATTAGCCACATTAGTCAGCGTACTTGCGTTAAGCTTAAGCagtag                                                                                                                                         <  1:327095/251‑1 (MQ=255)                                                                                                               ggCAATAGCCGTTATTTGCAACGTTTATTTCCCGATGGTCCGGCAAAACAGGCAACCAAGCTGGCGGGTCTCCCCAAACCGGTAAAATGTCTATAATTTTTTTGAGAATCATTTGCATTTAGAAATAGAAAGAATATAATGTCCGCGCAGCTTATGTTGTTAtt                                                                                                                                                                                                                           <  1:106354/164‑1 (MQ=255)                                                                                                               ggCAATAGCCGTTATTTGCAACGTTTATTTCCCGATGGTCCGGCAAAACAGGCAACCAAGCTGGCGGGTCTCCCCAAACCGGTAAAATGTCTATAATTTTTTTGAGAATCATTTGCATTTAGAAATAGAAAGAATATAATGTCCGCGCAGCTTATGTTGTTATTTTATAGTTAAAAAGCCAAGGAGAAAAAATATGCGTTTATTAGCCACATTAGTc                                                                                                                                                                      <  1:533619/217‑1 (MQ=255)                                                                                                               ggCAATAGCCGTTATTTGCAACGTTTATTTCCCGATGGTCCGGCAAAACAGGCAACCAAGCTGGCGGGTCTCCCCAAACCGGTAAAATGTCTATAATTTTTTTGAGAATCATTTGCATTTAGAAATAGAAAGAATATAATGTCCGCGCAGCTTATGTTGTTATTTTATAGTTAAAAAGCCAAGGAGAAAAAATATGCGTTTATTAGCCACATTAGTCAGCGTACTTGCGTTAAGCTTAAGCAGTAGTCTTg                                                                                                                                    >  1:666792/1‑251 (MQ=255)                                                                                                                  aaTAGCCGTTATTTGCAACGTTTATTTCCCGATGGTCCGGCAAAACAGGCAACCAAGCTGGCGGGTCTCCCCAAACCGGTAAAATGTCTATAATTTTTTTGAGAATCATTTGCATTTAGAAATAGAAAGAATATAATGTCCGCGCAGCTTATGTTGTTATTTTATAGTTAAAAAGCCAAGGAGAAAAAATATGCGTTTATTAGCCACATTAGTCAGCGTACTTGCGTTAAGCTTAAGCAGTAGTCTTGCtt                                                                                                                                 <  1:705330/251‑1 (MQ=255)                                                                                                                       ccGTTATTTGCAACGTTTATTTCCCGATGGTCCGGCAAAACAGGCAACCAAGCTGGCGGGTCTCCCCAAACCGGTAAAATGTCTATAATTTTTTTGAGAATCATTTGCATTTAGAAATAGAAAGAATATAATGTCCGCGCAGCTTATGTTGTTATTTTATAGTTAAAAAGCCAAGGAGAAAAAATATGCGTTTATTAGCCACATTAGTCAGCGTACTTGCGTTAAGCTTAAGCAGTAGTCTTGCTTATGCg                                                                                                                            >  1:606845/1‑251 (MQ=255)                                                                                                                           tATTTGCAACGTTTATTTCCCGATGGTCCGGCAAAACAGGCAACCAAGCTGGCGGGTCTCCCCAAACCGGTAAAATGTCTATAATTTTTTTGAGAATCATTTGCATTTAGAAATAGAAAGAATATAATGTCCGCGCAGCTTATGTTGTTATTTtat                                                                                                                                                                                                                       <  1:357189/156‑1 (MQ=255)                                                                                                                           tATTTGCAACGTTTATTTCCCGATGGTCCGGCAAAACAGGCAACCAAGCTGGCGGGTCTCCCCAAACCGGTAAAATGTCTATAATTTTTTTGAGAATCATTTGCATTTAGAAATAGAAAGAATATAATGTCCGCGCAGCTTATGTTGTTATTTTATAGTTAAAAAGCCAAGGAGAAAAAATATGCGTTTATTAGCCACATTAGTCAGCGt                                                                                                                                                                 <  1:57236/210‑1 (MQ=255)                                                                                                                            aTTTGCAACGTTTATTTCCCGATGGTCCGGCAAAACAGGCAACCAAGCTGGCGGGTCTCCCCAAACCGGTAAAATGTCTATAATTTTTTTGAGAATCATTTGCATTTAGAAATAGAAAGAATATAAt                                                                                                                                                                                                                                                   >  1:788872/1‑127 (MQ=255)                                                                                                                                     gTTTATTTCCCGATGGTCCGGCAAAACAGGCAACCAAGCTGGCGGGTCTCCCCAAACCGGTAAAATGTCTATAATTTTTTTGAGAATCATTTGCATTTAGAAATAGAAAGAATATAATGTc                                                                                                                                                                                                                                                <  1:77343/121‑1 (MQ=255)                                                                                                                                            tCCCGATGGTCCGGCAAAACAGGCAACCAAGCTGGCGGGTCTCCCCAAACCGGTAAAATGTCTATAATTTTTTTGAGAATCATTTGCATTTAGAAATAGAAAGAATATAATGTCCGCGCAGCTTATGTTGTTATTTTATAGTTAAAAAGCCAAGGAGAAAAAATATGCGTTTATTAGCCACATTAGTCAGCGTACTTGCGTTAAGCTTAAGCAGTAGTCTTGCTTATGCGAAATTTAAAGTTGTCACTACg                                                                                                       <  1:822048/251‑1 (MQ=255)                                                                                                                                                          cAAAACAGGCAACCAAGCTGGCGGGTCTCCCCAAACCGGTAAAATGTCTATAATTTTTTTGAGAATCATTTGCATTTAGAAATAGAAAGAATATAATGTCCGCGCAGCTTATGTTGTTATTTTATAGTTAAAAAGCCAAGGAGAAAAAATATGCGTTTATTAGCCACATTAGTCAGCGTACTTGCGTTAAGCTTAAGCAGTAGTCTTGCTTATGCGAAATTTAAAGTTGTCACTACGTTCACTATTATTc                                                                                          <  1:617294/250‑1 (MQ=255)                                                                                                                                                               cAGGCAACCAAGCTGGCGGGTCTCCCCAAACCGGTAAAATGTCTATAATTTTTTTGAGAATCATTTGCATTTAGAAATAGAAAGAATATAATGTCCGCGCAGCTTATGTTGTTATTTTATAGTTAAAAAGCCAAGGAGAAAAAATATGCGTTTATTAGCCACATTAGTCAGCGTACTTGCGTTAAGCTTAAGCAGTAGTCTTGCTTATGCGAAATTTAAAg                                                                                                                  <  1:543361/221‑1 (MQ=255)                                                                                                                                                                                     tcCCCAAACCGGTAAAATGTCTATAATTTTTTTGAGAATCATTTGCATTTAGAAATAGAAAGAATATAATGTCCGCGCAGCTTATGTTGTTATTTTATAGTTAAAAAGCCAAGGAGAAAAAATATGCGTTTATTAGCCACATTAGTCAGCGTACTTGCGTTAAGCTTAAGCAGTAGTCTTGCTTATGCGAAATTTAAAGTTGTCACTACGTTCACtatta                                                                                             >  1:159691/1‑220 (MQ=255)                                                                                                                                                                                      ccccAAACCGGTAAAATGTCTATAATTTTTTTGAGAATCATTTGCATTTAGAAATAGAAAGAATATAATGTCCGCGCAGCTTATGTTGTTATTTTATAGTTAAAAAGCCAAGGAGAAAAAATATGCGTTTATTAGCCACATTAGTCAGCGTACTTGCGTTAAGCTTAAGCAGTAGTCTTGCTTATGCGAAATTTAAAGTTGTCACTACGTTCACTATTATTc                                                                                          >  1:184892/1‑222 (MQ=255)                                                                                                                                                                                      ccccAAACCGGTAAAATGTCTATAATTTTTTTGAGAATCATTTGCATTTAGAAATAGAAAGAATATAATGTCCGCGCAGCTTATGTTGTTATTTTATAGTTAAAAAGCCAAGGAGAAAAAATATGCGTTTATTAGCCACATTAGTCAGCGTACTTGCGTTAAGCTTAAGCAGTAGTCTTGCTTATGCGAAATTTAAAGTTGTCACTACGTTCACTATTATTCAAGATATGGCACAAAATGTCGCCGGTGAt                                                             >  1:411878/1‑251 (MQ=255)                                                                                                                                                                                      ccccAAACCGGTAAAATGTCTATAATTTTTTTGAGAATCATTTGCATTTAGAAATAGAAAGAATATAATGTCCGCGCAGCTTATGTTGTTATTTTATAGTTAAAAAGCCAAGGAGAAAAAATATGCGTTTATTAGCCACATTAGTCAGCGTACTTGCGTTAAGCTTAAGCAGTAGTCTTGCTTATGCGAAATTTAAAGTTGTCACTACGTTCACTATTATTCAAGATATGGCACAAAATGTCGCCGGTGAt                                                             >  1:394094/1‑251 (MQ=255)                                                                                                                                                                                             ccGGTAAAATGTCTATAATTTTTTTGAGAATCATTTGCATTTAGAAATAGAAAGAATATAATGTCCGCGCAGCTTATGTTGTTATTTTATAGTTAAAAAGCCAAGGAGAAAAAATATGCGTTTATTAGCCACATTAGTCAGCGTACTTg                                                                                                                                                            <  1:268376/149‑1 (MQ=255)                                                                                                                                                                                               ggTAAAATGTCTATAATTTTTTTGAGAATCATTTGCATTTAGAAATAGAAAGAATATAATGTCCGCGCAGCTTATGTTGTTATTTTATAGTTAAAAAGCCAAGGAGAAAAAATATGCGTTTATTAGCCACATTAGTCAGCGTACTTGCGTTAAGCTTAAGCAGTAGTCTTGCTTATGCGAAATTTAAAGTTGTCACTACGTTCACt                                                                                                 >  1:516975/1‑206 (MQ=255)                                                                                                                                                                                                       gTCTATAATTTTTTTGAGAATCATTTGCATTTAGAAATAGAAAGAATATAATGTCCGCGCAGCTTATGTTGTTATTTTATAGTTAAAAAGCCAAGGAGAAAAAATATGCGTTTATTAGcc                                                                                                                                                                               <  1:462425/120‑1 (MQ=255)                                                                                                                                                                                                       gTCTATAATTTTTTTGAGAATCATTTGCATTTAGAAATAGAAAGAATATAATGTCCGCGCAGCTTATGTTGTTATTTTATAGTTAAAAAGCCAAGGAGAAAAAATATGCGTTTATTAGCCACATTAGTCAGCGTACTTGCGTTAAGCTTAAGCAGTAGTCTTGCTTATGCGAAATTTAAAGTTGTCACTACGTTCACTATTATTCAAGATATGGCACAAAATGTCGCCGGTGATGCGGCAACCGTAGAGt                                             >  1:529695/1‑250 (MQ=255)                                                                                                                                                                                                            taATTTTTTTGAGAATCATTTGCATTTAGAAATAGAAAGAATATAATg                                                                                                                                                                                                                                                  >  1:399059/1‑48 (MQ=255)                                                                                                                                                                                                                           tCATTTGCATTTAGAAATAGAAAGAATATAATGTCCGCGCAGCTTATGTTGTTATTTTATAGTTAAAAAGCCAAGGAGAAAAAATATGCGTTTATTAGCCACATTAGTc                                                                                                                                                                      <  1:663624/109‑1 (MQ=255)                                                                                                                                                                                                                            cATTTGCATTTAGAAATAGAAAGAATATAATGTCCGCGCAGCTTATGTTGTTATTTTATAGTTAAAAAGCCAAGGAGAAAAAATATGCGTTTATTAGCCACATTAGTCAGCGTACTTGCGTTAAGCTTAAGCAGTAGTCTTGCt                                                                                                                                  >  1:495708/1‑144 (MQ=255)                                                                                                                                                                                                                            cATTTGCATTTAGAAATAGAAAGAATATAATGTCCGCGCAGCTTATGTTGTTATTTTATAGTTAAAAAGCCAAGGAGAAAAAATATGCGTTTATTAGCCACATTAGTCAGCGTACTTGCGTTAAGCTTAAGCAGTAGTCTTGCTTATGCGAAATTTAAAGTTGTCACTACGTTCACTATTATTCAAGATATGGCACAAAATGTCGCCGGTGATGCGGCAACCGTAGAGTCCATCACCAAGCCCGGAGCAGa                       >  1:128523/1‑251 (MQ=255)                                                                                                                                                                                                                            cATTTGCATTTAGAAATAGAAAGAATATAATGTCCGCGCAGCTTATGTTGTTATTTTATAGTTAAAAAGCCAAGGAGAAAAAATATGCGTTTATTAGCCACATTAGTCAGCGTACTTGCGTTAAGCTTAAGCAGTAGTCTTGCTTATGCGAAATTTAAAGTTGTCACTACGTTCACTATTATTCAAGATATGGCACAAAATGTCGCCGGTGATGCGGCAACCGTAGAGTCCATCACCAAGCCCGGAGCAGa                       >  1:486568/1‑251 (MQ=255)                                                                                                                                                                                                                                        gaaatagaaaGAATATAATGTCCGCGCAGCTTATGTTGTTATTTTATAGTTAAAAAGCCAAGgag                                                                                                                                                                                                     <  1:861040/65‑1 (MQ=255)                                                                                                                                                                                                                                         aaatagaaaGAATATAATGTCCGCGCAGCTTATGTTGTTATTTTATAGTTAAAAAGCCAAGGAGAAAAAATATGCGTTTATTAGCCACATTAGTCAGCGTACTTGCGTTAAGCTTAAGCAGTAGTCTTGCTTATGCGAAATTTAAAGTTGTCACTACGTTCACTATTATTCAAGATATGGCACAAAATGTCGCCGGTGATGCGGCAACCGTAGAGTCCATCACCAAGCCCGGAGCAGAAATCCATGATTa           <  1:141875/250‑1 (MQ=255)                                                                                                                                                                                                                                             agaaagaaTATAATGTCCGCGCAGCTTATGTTGTTATTTTATAGTTAAAAAGCCAAGGAGAAAAAATATGCGTTTATTAGCCACATTAGTCAGCGTACTTGCGTTAAGCTTAAGCAGTAGTCTTGCTTATGCGAAATTTAAAGt                                                                                                                 <  1:519217/144‑1 (MQ=255)                                                                                                                                                                                                                                                  gaaTATAATGTCCGCGCAGCTTATGTTGTTATTTTATAGTTAAAAAGCCAAGGAGAAAAAATATGCGTTTATTAGCCACATTAGTCAGCGTACTTg                                                                                                                                                            <  1:831325/96‑1 (MQ=255)                                                                                                                                                                                                                                                  gaaTATAATGTCCGCGCAGCTTATGTTGTTATTTTATAGTTAAAAAGCCAAGGAGAAAAAATATGCGTTTATTAGCCACATTAGTCAGCGTACTTGCGTTAAGCTTAAGCAGTAGTCTTGCTTATGCGAAATTTAAAGTTGTCACTACGTTCACTATTATTCAAGATATGGCACAAAATGTCGCCGGTGATGCGGCAACCGTag                                                >  1:574992/1‑204 (MQ=255)                                                                                                                                                                                                                                                  gaaTATAATGTCCGCGCAGCTTATGTTGTTATTTTATAGTTAAAAAGCCAAGGAGAAAAAATATGCGTTTATTAGCCACATTAGTCAGCGTACTTGCGTTAAGCTTAAGCAGTAGTCTTGCTTATGCGAAATTTAAAGTTGTCACTACGTTCACTATTATTCAAGATATGGCACAAAATGTCGCCGGTGATGCGGCAACCGTAGAGTCCATCACCAAGCCCGGAGCAGAAATCCATGATTACGAGCCAAc  >  1:789014/1‑250 (MQ=255)                                                                                                                                                                                                                                                  gaaTATAATGTCCGCGCAGCTTATGTTGTTATTTTATAGTTAAAAAGCCAAGGAGAAAAAATATGCGTTTATTAGCCACATTAGTCAGCGTACTTGCGTTAAGCTTAAGCAGTAGTCTTGCTTATGCGAAATTTAAAGTTGTCACTACGTTCACTATTATTCAAGATATGGCACAAAATGTCGCCGGTGATGCGGCAACCGTAGAGTCCATCACCAAGCCCGGAGCAGAAATCCATGATTACGAGCCAAc  >  1:738033/1‑250 (MQ=255)                                                                                                                                                                                                                                                   aaTATAATGTCCGCGCAGCTTATGTTGTTATTTTATAGTTAAAAAGCCAAGGAGAAAAAATATGCGTTTATTAGCCACATTAGTCAGCGTACTTGCGTTAAGCTTAAGCAGTAGTCTTGCTTATGCGAAATTTAAAGTTGTCACTACGTTCACTATTATTCAAGATATGGCACAAAATGTCGCCGGTGATGCGGCAACCGTAGAGTCCATCAc                                      <  1:821665/213‑1 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                     CCCATTGGGAAATTATCCATTTTGTGCGTGATTTTTATCAAGAATATAAAACCTCTCCTGCCATTCGTATGTTGGTGAAAGCCACCGCGCAAAAGCTGGGAGAAGATAAAGGCAATAGCCGTTATTTGCAACGTTTATTTCCCGATGGTCCGGCAAAACAGGCAACCAAGCTGGCGGGTCTCCCCAAACCGGTAAAATGTCTATAATTTTTTTGAGAATCATTTGCATTTAGAAATAGAAAGAATACAATGTCCGCGCAGCTTATGTTGTTATTTTATAGTTAAAAAGCCAAGGAGAAAAAATATGCGTTTATTAGCCACATTAGTCAGCGTACTTGCGTTAAGCTTAAGCAGTAGTCTTGCTTATGCGAAATTTAAAGTTGTCACTACGTTCACTATTATTCAAGATATGGCACAAAATGTCGCCGGTGATGCGGCAACCGTAGAGTCCATCACCAAGCCCGGAGCAGAAATCCATGATTACGAGCCAAC  >  1/1790459‑1790949

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