Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	1,486,949	T→C	C34C (TGT→TGC)	unknown →	Flp pilus assembly membrane protein TadE

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	1,486,949	0	T	C	100.0%	112.5 / NA	42	C34C (TGT→TGC)	unknown	Flp pilus assembly membrane protein TadE
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (14/28);  total (14/28)

TCAAAAATAATGACATTGATTACGCAAAAGATATTATTGTGAAGGAAAACATTAATACATCTCCTGATGATTTAATTAATGCTCTCAAAAAAACTGAACGCGTATCTAATAATATCGCCAGATAAAAAAAGGTTCTAT‑AATGAAAAAATTCCTATCGAATAAGAAGGGAGCATCAACCGTTGAGTTTGCCTTGACGGTTGTATTTTATTTCTTTGTCGTATTCCTAATTCTAGAATTCTGTCGTATATCTATTACAACTGCGTACTGGGACTTGGCTATCACTGAAAGTGCAAGAATTGCTAAAAATAGAACTGCAGAGGGTAATGATTATGCGGTAGAATTTGAGAAAGCTCTAAAACAACAGCTTGTTTATCAAGAATCTTCAACTTTAGGTTATCTCGCCCGTTTGGATAAAAATGGCGGTTACAAGATTGATGTTAAATATGTTGATTGTGGTAGTGAGAGTTCTTGTGTTAAGTC  >  1/1486709‑1487188                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                tCAAGAATAACGAAATTGATTACGCAAAAGATATTATTGTGAAGGAAAACATTAACACATCCCCTGATGATTTAGTTAATGCCCTCAAGAAAACTGAACGCGTGTCTAGTAATATCACCAGAT‑‑AAAAAGGTTCTGTAAATGGAAAAATTCCTATCAAATAAAAAGGGAGCGTCAGCCGTTGAGTTTGCCTTGACGGTTGTATTCTATTTCTTTATCGTATTCTTAATTCTAGAATTCTGCCGTATATCTa                                                                                                                                                                                                                                       <  1:167299/250‑1 (MQ=255)      aataacGAAATTGATTACGCAAAAGATATTATTGTGAAGGAAAACATTAACACATCCCCTGATGATTTAGTTAATGCCCTCAAGAAAACTGAACGCGTGTCTAGTAATATCACCAGAT‑‑AAAAAGGTTCTGTAAATGGAAAAATTCCTATCAAATAAAAAGGGAGCGTCAGCCGTTGAGTTTGCCTTGACGGTTGTATTCTATTTCTTTATCGTATTCTTAATTCTAGAATTCTGCCGTATATCTATTAcc                                                                                                                                                                                                                                  <  1:87235/250‑2 (MQ=255)      aataacGAAATTGATTACGCAAAAGATATTATTGTGAAGGAAAACATTAACACATCCCCTGATGATTTAGTTAATGCCCTCAAGAAAACTGAACGCGTGTCTAGTAATATCACCAGAT‑‑AAAAAGGTTCTGTAAATGGAAAAATTCCTATCAAATAAAAAGGGAGCGTCAGCCGTTGAGTTTGCCTTGACGGTTGTATTCTATTTCTTTATCGTATTCTTAATTCTAGAATTCTGCCGTATATCTATTAcc                                                                                                                                                                                                                                  <  1:209028/250‑2 (MQ=255)      aataacGAAATTGATTACGCAAAAGATATTATTGTGAAGGAAAACATTAACACATCCCCTGATGATTTAGTTAATGCCCTCAAGAAAACTGAACGCGTGTCTAGTAATATCACCAGAT‑‑AAAAAGGTTCTGTAAATGGAAAAATTCCTATCAAATAAAAAGGGAGCGTCAGCCGTTGAGTTTGCCTTGACGGTTGTATTCTATTTCTTTATCGTATTCTTAATTCTAGAATTCTGCCGTATATCTATTAcc                                                                                                                                                                                                                                  <  1:59003/250‑2 (MQ=255)      aataacGAAATTGATTACGCAAAAGATATTATTGTGAAGGAAAACATTAACACATCCCCTGATGATTTAGTTAATGCCCTCAAGAAAACTGAACGCGTGTCTAGTAATATCACCAGAT‑‑AAAAAGGTTCTGTAAATGGAAAAATTCCTATCAAATAAAAAGGGAGCGTCAGCCGTTGAGTTTGCCTTGACGGTTGTATTCTATTTCTTTATCGTATTCTTAATTCTAGAATTCTGCCGTATATCTATTAc                                                                                                                                                                                                                                   <  1:199123/249‑1 (MQ=255)             aaaTTGATTACGCAAAAGATATTATTGTGAAGGAAAACATTAACACATCCCCTGATGATTTAGTTAATGCCCTCAAGAAAACTGAACGCGTGTCTAGTAATATCACCAGAT‑‑AAAAAGGTTCTGTAAATGGAAAAATTCCTATCAAATAAAAAGGGAGCGTCAGCCGTTGAGTTTGCCTTGACGGTTGTATTCTATTTCTTTATCGTATTCTTAATTCTAGAATTCTGCCGTATATCTATTACCACTGcat                                                                                                                                                                                                                           <  1:786540/248‑3 (MQ=255)                 tGATTACGCAAAAGATATTATTGTGAAGGAAAACATTAACACATCCCCTGATGATTTAGTTAATGCCCTCAAGAAAACTGAACGCGTGTCTAGTAATATCACCAGAT‑‑AAAAAGGTTCTGTAAATGGAAAAATTCCTATCAAATAAAAAGGGAGCGTCAGCCGTTGAGTTTGCCTTGACGGTTGTATTCTATTTCTTTATCGTATTCTTAATTCTAGAATTCTGCCGTATATCTATTACCACTGcatattgg                                                                                                                                                                                                                      <  1:658225/251‑8 (MQ=255)                              gATATTATTGTGAAGGAAAACATTAACACATCCCCTGATGATTTAGTTAATGCCCTCAAGAAAACTGAACGCGTGTCTAGTAATATCACCAGAT‑‑AAAAAGGTTCTGTAAATGGAAAAATTCCTATCAAATAAAAAGGGAGCGTCAGCCGTTGAGTTTGCCTTGACGGTTGTATTCTATTTCTTTATCGTATTCTTAATTCTAGAATTCTGCCGTATATCTATTACCACTGCATATTGGGatttagccataa                                                                                                                                                                                                         >  1:17258/1‑240 (MQ=255)                                          aaGGAAAACATTAACACATCCCCTGATGATTTAGTTAATGCCCTCAAGAAAACTGAACGCGTGTCTAGTAATATCACCAGAT‑‑AAAAAGGTTCTGTAAATGGAAAAATTCCTATCAAATAAAAAGGGAGCGTCAGCCGTTGAGTTTGCCTTGACGGTTGTATTCTATTTCTTTATCGTATTCTTAATTCTAGAATTCTGCCGTATATCTATTACCACTGCATATTGGGATTTAGCCATAACTGAAAGTGCaa                                                                                                                                                                                             >  1:447683/1‑251 (MQ=255)                                             gAAAACATTAACACATCCCCTGATGATTTAGTTAATGCCCTCAAGAAAACTGAACGCGTGTCTAGTAATATCACCAGAT‑‑AAAAAGGTTCTGTAAATGGAAAAATTCCTATCAAATAAAAAGGGAGCGTCAGCCGTTGAGTTTGCCTTGACGGTTGTATTCTATTTCTTTATCGTATTCTTAATTCTAGAATTCTGCCGTATATCTATTACCACTGCATATTGGGATTTAGCCATAACTGAAAGTGCAAg                                                                                                                                                                                            <  1:887285/249‑1 (MQ=255)                                              aaaaCATTAACACATCCCCTGATGATTTAGTTAATGCCCTCAAGAAAACTGAACGCGTGTCTAGTAATATCACCAGAT‑‑AAAAAGGTTCTGTAAATGGAAAAATTCCTATCAAATAAAAAGGGAGCGTCAGCCGTTGAGTTTGCCTTGACGGTTGTATTCTATTTCTTTATCGTATTCTTAATTCTAGAATTCTGCCGTATATCTATTACCACTGCATATTGGGATTTAGCCATAACTGAAAGTGCAAgga                                                                                                                                                                                          <  1:698914/250‑3 (MQ=255)                                              aaaaCATTAACACATCCCCTGATGATTTAGTTAATGCCCTCAAGAAAACTGAACGCGTGTCTAGTAATATCACCAGAT‑‑AAAAAGGTTCTGTAAATGGAAAAATTCCTATCAAATAAAAAGGGAGCGTCAGCCGTTGAGTTTGCCTTGACGGTTGTATTCTATTTCTTTATCGTATTCTTAATTCTAGAATTCTGCCGTATATCTATTACCACTGCATATTGGGATTTAGCCATAACTGAAAGTGCAAgga                                                                                                                                                                                          <  1:586392/250‑3 (MQ=255)                                                        caCATCCCCTGATGATTTAGTTAATGCCCTCAAGAAAACTGAACGCGTGTCTAGTAATATCACCAGAT‑‑AAAAAGGTTCTGTAAATGGAAAAATTCCTATCAAATAAAAAGGGAGCGTCAGCCGTTGAGTTTGCCTTGACGGTTGTATTCTATTTCTTTATCGTATTCTTAATTCTAGAATTCTGCCGTATATCTATTACCACTGCATATTGGGATTTAGCCATAACTGAAAGTGCAAGGATCGCTaaaaac                                                                                                                                                                               <  1:494047/250‑2 (MQ=255)                                                            tccccTGATGATTTAGTTAATGCCCTCAAGAAAACTGAACGCGTGTCTAGTAATATCACCAGAT‑‑AAAAAGGTTCTGTAAATGGAAAAATTCCTATCAAATAAAAAGGGAGCGTCAGCCGTTGAGTTTGCCTTGACGGTTGTATTCTATTTCTTTATCGTATTCTTAATTCTAGAATTCTGCCGTATATCTATTACCACTGCATATTGGGATTTAGCCATAACTGAAAGTGCAAGGATCGCTAAAAACAGaa                                                                                                                                                                           >  1:35572/4‑251 (MQ=255)                                                            tccccTGATGATTTAATTAATGCCCTCAAGAAAACTGAACGCGTGTCTAGTAATATCACCAGAT‑‑AAAAAGGTTCTGTAAATGGAAAAATTCCTATCAAATAAAAAGGGAGCGTCAGCCGTTGAGTTTGCCTTGACGGTTGTATTCTATTTCTTTATCGTATTCTTAATTCTAGAATTCTGCCGTATATCTATTAcc                                                                                                                                                                                                                                  <  1:14055/193‑2 (MQ=255)                                                                  gatgatTTAGTTAATGCCCTCAAGAAAACTGAACGCGTGTCTAGTAATATCACCAGAT‑‑AAAAAGGTTCTGTAAATGGAAAAATTCCTATCAAATAAAAAGGGAGCGTCAGCCGTTGAGTTTGCCTTGACGGTTGTATTCTATTTCTTTATCGTATTCTTAATTCTAGAATTCTGCCGTATATCTATTACCACTGCATATTGGGATTTAGCCATAACTGAAAGTGCAAgga                                                                                                                                                                                          <  1:462452/230‑3 (MQ=255)                                                                             taatGCCCTCAAGAAAACTGAACGCGTGTCTAGTAATATCACCAGAT‑‑AAAAAGGTTCTGTAAATGGAAAAATTCCTATCAAATAAAAAGGGAGCGTCAGCCGTTGAGTTTGCCTTGACGGTTGTATTCTATTTCTTTATCGTATTCTTAATTCTAGAATTCTGCCGTATATCTATTACCACTGCATATTGGGATTTAGCCATAACTGAAAGTGCAAGGATCGCTAAAAACAGAACGGCGGAAGGTAACGAt                                                                                                                                                          <  1:482743/251‑1 (MQ=255)                                                                                                 aaCGCGTGTCTAGTAATATCACCAGAT‑‑AAAAAGGTTCTGTAAATGGAAAAATTCCTATCAAATAAAAAGGGAGCGTCAGCCGTTGAGTTTGCCTTGACGGTTGTATTCTATTTCTTTATCGTATTCTTAATTCTAGAATTCTGCCGTATATCTATTACCACTGCATATTGGGATTTAGCCATAACTGAAAGTGCAAGGATCGCTAAAAACAGAACGGCGGAAGGTAACGATTATGCTGCAGAATTTgaaaa                                                                                                                                      <  1:423956/251‑4 (MQ=255)                                                                                                        gtCTAGTAATATCACCAGAT‑‑AAAAAGGTTCTGTAAATGGAAAAATTCCTATCAAATAAAAAGGGAGCGTCAGCCGTTGAGTTTGCCTTGACGGTTGTATTCTATTTCTTTATCGTATTCTTAATTCTAGAATTCTGCCGTATATCTATTAccac                                                                                                                                                                                                                                <  1:81610/153‑4 (MQ=255)                                                                                                        gtCTAGTAATATCACCAGAT‑‑AAAAAGGTTCTGTAAATGGAAAAATTCCTATCAAATAAAAAGGGAGCGTCAGCCGTTGAGTTTGCCTTGACGGTTGTATTCTATTTCTTTATCGTATTCTTAATTCTAGAATTCTGCCGTATATCTATTACCACTGCATATTGGGatttagcca                                                                                                                                                                                                            <  1:857768/173‑9 (MQ=255)                                                                                                             gtaataTCACCAGAT‑‑AAAAAGGTTCTGTAAATGGAAAAATTCCTATCAAATAAAAAGGGAGCGTCAGCCGTTGAGTTTGCCTTGACGGTTGTATTCTATTTCTTTATCGTATTCTTAATTCTAGAATTCTGCCGTATATCTATTACCACTGcatatt                                                                                                                                                                                                                        <  1:520550/156‑6 (MQ=255)                                                                                                             gtaataTCACCAGAT‑‑AAAAAGGTTCTGTAAATGGAAAAATTCCTATCAAATAAAAAGGGAGCGTCAGCCGTTGAGTTTGCCTTGACGGTTGTATTCTATTTCTTTATCGTATTCTTAATTCTAGAATTCTGCCGTATATCTATTACCACTGCATATTGGGATTTAGCCATAACTGAAAGTGCAAgga                                                                                                                                                                                          >  1:632849/2‑185 (MQ=255)                                                                                                             gtaataTCACCAGAT‑‑AAAAAGGTTCTGTAAATGGAAAAATTCCTATCAAATAAAAAGGGAGCGTCAGCCGTTGAGTTTGCCTTGACGGTTGTATTCTATTTCTTTATCGTATTCTTAATTCTAGAATTCTGCCGTATATCTATTACCACTGCATATTGGGATTTAGCCATAACTGAAAGTGCAAGGATCGCTAAAAACAGAACGGCGGAAGGTAACGATTATGCTGCAGAATTTGAAAAAGCTCTAAaaca                                                                                                                          >  1:179337/2‑251 (MQ=255)                                                                                                                     tcaccagataaaaaGGTTCTGTAAATGGAAAAATTCCTATCAAATAAAAAGGGAGCGTCAGCCGTTGAGTTTGCCTTGACGGTTGTATTCTATTTCTTTATCGTATTCTTAATTCTAGAATTCTGCCGTATATCTa                                                                                                                                                                                                                                       >  1:288878/10‑136 (MQ=255)                                                                                                                          agataaaaaGGTTCTGTAAATGGAAAAATTCCTATCAAATAAAAAGGGAGCGTCAGCCGTTGAGTTTGCCTTGACGGTTGTATTCTATTTCTTTATCGTATTCTTAATTCTAGAATTCTGCCGTATATCTATTACCACTGCATATTGGGATTTAGCCATAACTGAAAGTg                                                                                                                                                                                                >  1:594033/5‑170 (MQ=255)                                                                                                                           gataaaaaGGTTCTGTAAATGGAAAAATTCCTATCAAATAAAAAGGGAGCGTCAGCCGTTGAGTTTGCCTTGACGGTTGTATTCTATTTCTTTATCGTATTCTTAATTCTAGAATTCTGCCGTATATCTATTAccac                                                                                                                                                                                                                                <  1:355721/134‑4 (MQ=255)                                                                                                                           gataaaaaGGTTCTGTAAATGGAAAAATTCCTATCAAATAAAAAGGGAGCGTCAGCCGTTGAGTTTGCCTTGACGGTTGTATTCTATTTCTTTATCGTATTCTTAATTCTAGAATTCTGCCGTATATCTATTAc                                                                                                                                                                                                                                   >  1:340139/4‑134 (MQ=255)                                                                                                                           gataaaaaGGTTCTGTAAATGGAAAAATTCCTATCAAATAAAAAGGGAGCGTCAGCCGTTGAGTTTGCCTTGACGGTTGTATTCTATTTCTTTATCGTATTCTTAATTCTAGAATTCTGCCGTATATCTATTACCACTGCATATTGGGATTTAGCCATAACTGAAAGTGCAAGGATCGCTAAAAACAGAACGGCGGAAGGTAACGATTATGCTGCAGAATTTgaaaa                                                                                                                                      >  1:570457/4‑224 (MQ=255)                                                                                                                           gataaaaaGGTTCTGTAAATGGAAAAATTCCTATCAAATAAAAAGGGAGCGTCAGCCGTTGAGTTTGCCTTGACGGTTGTATTCTATTTCTTTATCGTATTCTTAATTCTAGAATTCTGCCGTATATCTATTACCACTGCATATTGGGATTTAGCCATAACTGAAAGTGCAAGGATCGCTAAAAACAGAACGGCGGAAGGTAACGATTATGCTGCAGAATTTGAAAAAGCTCTAAAACAACAACTTATTTa                                                                                                              >  1:596994/4‑251 (MQ=255)                                                                                                                            ataaaaaGGTTCTGTAAATGGAAAAATTCCTATCAAATAAAAAGGGAGCGTCAGCCGTTGAGTTTGCCTTGACGGTTGTATTCTATTTCTTTATCGTATTCTTAATTCTAGAATTCTGCCGTATATCTATTACCACTGCATATTGGGATTTAGCCATAACTGAAAGTGCAAGGATCGCTAAAAACAGAACGGCGGAAGGTAACGATTATGCTGCAGAATTTgaaaa                                                                                                                                      >  1:663824/3‑223 (MQ=255)                                                                                                                            ataaaaaGGTTCTGTAAATGGAAAAATTCCTATCAAATAAAAAGGGAGCGTCAGCCGTTGAGTTTGCCTTGACGGTTGTATTCTATTTCTTTATCGTATTCTTAATTCTAGAATTCTGCCGTATATCTATTACCACTGCATATTGGGATTTAGCCATAACTGAAAGTGCAAGGATCGCTAAAAACAGAACGGCGGAAGGTAACGATTATGCTGCAGAATTTGAAAAAGCTCTAAAACAACAACttatttat                                                                                                             >  1:138808/3‑251 (MQ=255)                                                                                                                            ataaaaaGGTTCTGTAAATGGAAAAATTCCTATCAAATAAAAAGGGAGCGTCAGCCGTTGAGTTTGCCTTGACGGTTGTATTCTATTTCTTTATCGTATTCTTAATTCTAGAATTCTGCCGTATATCTATTACCACTGCATATTGGGATTTAGCCATAACTGAAAGTGCAAGGATCGCTAAAAACAGAACGGCGGAAGGTAACGATTATGCTGCAGAATTTGAAAAAGCTCTAAAACAACAACttatttat                                                                                                             >  1:304702/3‑251 (MQ=255)                                                                                                                                                 aaaaaTTCCTATCAAATAAAAAGGGAGCGTCAGCCGTTGAGTTTGCCTTGACGGTTGTATTCTATTTCTTTATCGTATTCTTAATTCTAGAATTCTGCCGTATATCTATTACCACTGCATATTGGGATTTAGCCATAACTGAAAGTGCAAGGATCGCTAAAAACAGAACGGCGGAAGGTAACGATTATGCTGCAGAATTTGAAAAAGCTCTAAAACAACAACTTATTTATCAAGAAACCTCAACTATAggc                                                                                        <  1:335046/251‑2 (MQ=255)                                                                                                                                                               aaTAAAAAGGGAGCGTCAGCCGTTGAGTTTGCCTTGACGGTTGTATTCTATTTCTTTATCGTATTCTTAATTCTAGAATTCTGCCGTATATCTATTACCACTGCATATTGGGATTTAGCCATAACTGAAAGTGCAAGGATCGCTAAAAACAGAACGGCGGAAGGTAACGATTATGCTGCAGAATTTGAAAAAGCTCTAAAACAACAACTTATTTATCAAGAAACCTCAACtat                                                                                            <  1:371139/233‑3 (MQ=255)                                                                                                                                                               aaTAAAAAGGGAGCGTCAGCCGTTGAGTTTGCCTTGACGGTTGTATTCTATTTCTTTATCGTATTCTTAATTCTAGAATTCTGCCGTATATCTATTACCACTGCATATTGGGATTTAGCCATAACTGAAAGTGCAAGGATCGCTAAAAACAGAACGGCGGAAGGTAACGATTATGCTGCAGAATTTGAAAAAGCTCTAAAACAACAACTTATTTATCAAGAAACCTCAACTATAGGCTATCTTGCCCGttt                                                                          >  1:111173/1‑251 (MQ=255)                                                                                                                                                                   aaaagGGAGCGTCAGCCGTTGAGTTTGCCTTGACGGTTGTATTCTATTTCTTTATCGTATTCTTAATTCTAGAATTCTGCCGTATATCTATTACCACTGCATATTGGGATTTAGCCATAACTGAAAGTGCAAgga                                                                                                                                                                                          <  1:446240/133‑3 (MQ=255)                                                                                                                                                                                          tttGCCTTGACGGTTGTATTCTATTTCTTTATCGTATTCTTAATTCTAGAATTCTGCCGTATATCTATTACCACTGCATATTGGGATTTAGCCATAACTGAAAGTGCAAGGATCGCTAAAAACAGAACGGCGGAAGGTAACGATTATGCTGCAGAATTTGAAAAAGCTCTAAAACAACAACTTATTTATCAAGAAACCTCAACTATAGGCTATCTTGCCCGTTTGGATAAAAATGGCGGTTatgaaattg                                                <  1:274106/250‑9 (MQ=255)                                                                                                                                                                                                                      ttatCGTATTCTTAATTCTAGAATTCTGCCGTATATCTATTACCACTGCATATTGGGATTTAGCCATAACTGAAAGTGCAAGGATCGCTAAAAACAGAACGGCGGAAGGTAACGATTATGCTGCAGAATTTGAAAAAGCTCTAAAACAACAACTTATTTATCAAGAAACCTCAACTATAGGCTATCTTGCCCGTTTGGATAAAAATGGCGGTTATGAAATTGATGTTAAATATGTTGATTGTGGTAGTgag                   <  1:690141/248‑1 (MQ=255)                                                                                                                                                                                                                      ttatCGTATTCTTAATTCTAGAATTCTGCCGTATATCTATTACCACTGCATATTGGGATTTAGCCATAACTGAAAGTGCAAGGATCGCTAAAAACAGAACGGCGGAAGGTAACGATTATGCTGCAGAATTTGAAAAAGCTCTAAAACAACAACTTATTTATCAAGAAACCTCAACTATAGGCTATCTTGCCCGTTTGGATAAAAATGGCGGTTATGAAATTGATGTTAAATATGTTGATTGTGGTAGTgag                   <  1:726941/248‑1 (MQ=255)                                                                                                                                                                                                                        atCGTATTCTTAATTCTAGAATTCTGCCGTATATCTATTACCACTGCATATTGGGATTTAGCCATAACTGAAAGTGCAAGGATCGCTAAAAACAGAACGGCGGAAGGTAACGATTATGCTGCAGAATTTGAAAAAGCTCTAAAACAACAACTTATTTATCAAGAAACCTCAACTATAGGCTATCTTGCCCGTTTGGATAAAAATGGCGGTTATGAAATTGATGTTAAATATGTTGATTGTGGt                         <  1:292311/242‑1 (MQ=255)                                                                                                                                                                                                                              ttcttAATTCTAGAATTCTGCCGTATATCTATTACCACTGCATATTGGGATTTAGCCATAACTGAAAGTGCAAGGATCGCTAAAAACAGAACGGCGGAAGGTAACGATTATGCTGCAGAATTTGAAAAAGCTCTAAAACAACAACTTATTTATCAAGAAACCTCAACTATAGGCTATCTTGCCCGTTTGGATAAAAATGGCGGTTATGAAATTGATGTTAAATATGTTGATTGTGGTAGTGAGAGCTCTtg           <  1:620439/247‑1 (MQ=255)                                                                                                                                                                                                                                       cTAGAATTCTGCCGTATATCTATTACCACTGCATATTGGGATTTAGCCATAACTGAAAGTGCAAGGATCGCTAAAAACAGAACGGCGGAAGGTAACGATTATGCTGCAGAATTTGAAAAAGCTCTAAAACAACAACTTATTTATCAAGAAACCTCAACTATAGGCTATCTTGCCCGTTTGGATAAAAATGGCGGTTATGAAATTGATGTTAAATATGTTGATTGTGGTAGTGAGAGCTCTTGTATTaaatc  <  1:360761/251‑4 (MQ=255)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                TCAAAAATAATGACATTGATTACGCAAAAGATATTATTGTGAAGGAAAACATTAATACATCTCCTGATGATTTAATTAATGCTCTCAAAAAAACTGAACGCGTATCTAATAATATCGCCAGATAAAAAAAGGTTCTAT‑AATGAAAAAATTCCTATCGAATAAGAAGGGAGCATCAACCGTTGAGTTTGCCTTGACGGTTGTATTTTATTTCTTTGTCGTATTCCTAATTCTAGAATTCTGTCGTATATCTATTACAACTGCGTACTGGGACTTGGCTATCACTGAAAGTGCAAGAATTGCTAAAAATAGAACTGCAGAGGGTAATGATTATGCGGTAGAATTTGAGAAAGCTCTAAAACAACAGCTTGTTTATCAAGAATCTTCAACTTTAGGTTATCTCGCCCGTTTGGATAAAAATGGCGGTTACAAGATTGATGTTAAATATGTTGATTGTGGTAGTGAGAGTTCTTGTGTTAAGTC  >  1/1486709‑1487188

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