Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	1,453,130	T→C	L111L (CTT→CTC)	unknown →	Coupling protein VirD4, ATPase required for T‑DNA transfer

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	1,453,130	0	T	C	100.0%	161.3 / NA	49	L111L (CTT→CTC)	unknown	Coupling protein VirD4, ATPase required for T‑DNA transfer
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (23/26);  total (23/26)

AATTCGTTTAATCGAAGCATTTAAGCGCTACTACCCTACTCTAACGCTGGATGCCATTTGGGGTAATTATCCGAAAATATGGCAGTCTCTCCTTATTAGTTTTGTAGTGAGCTTACTAATCATGCTGAAACTATTCTATCAACCTAAAAAATCCCTTTATGGTAACGCTGAATTTGCCACAAAATCAGAAGTACAAAAAATGGGTTTATTGGAAGAAAAAAACGGGATTATTGTCGGCAAACTTGGTAATAAACTTCTTCGCTTTATCGGTGCACAATTCGTTTCGTTGGCCGCACCGACCCGAGGTGGGAAAGGAATAGGGATTGTCATACCGAATTTATTAGAGTGGCAACAATCAGTTGTTGTGAAAGATGTTAAAAACGAATGTTATCGCATTACCAGTAAATACCGTCAAAAAATATTGGGGCAAAAAGTCTATCG  >  1/1452887‑1453327                                                                                                                                                                                                                                                    |                                                                                                                                                                                                      aaTTCGTTTAATCGAAGCATTTAAGCGCTACTACCCTACTCTAACGCTGGATGCCATTTGGGGTAATTATCCGAAAATATGGCAGTCTCTCCTTATTAGTTTTGTAGTGAGCTTACTAATCATGCTGAAACTATTCTATCAACCTAAAAAATCCCTTTATGGTAACGCTGAATTTGCCACAAAATCAGAAGTACAAAAAATGGGTTTATTGGAAGAAAAGAACGGGATTATTGTCGGCAAACTCGGtaata                                                                                                                                                                                                >  1:233116/1‑251 (MQ=255)    tCGTTTAATCGAAGCATTTAAGCGCTACTACCCTACTCTAACGCTGGATGCCATTTGGGGTAATTATCCGAAAATATGGCAGTCTCTCCTTATTAGTTTTGTAGTGAGCTTACTAATCATGCTGAAACTATTCTATCAACCTAAAAAATCCCTTTATGGTAACGCTGAATTTGCCACAAAATCAGAAGTACAAAAAATGGGTTTATTGGAAGAAAAAAACGGGATTATTGTCGGCAAACTCGGTAATAAAc                                                                                                                                                                                             <  1:127528/251‑1 (MQ=255)    tCGTTTAATCGAAGCATTTAAGCGCTACTACCCTACTCTAACGCTGGATGCCATTTGGGGTAATTATCCGAAAATATGGCAGTCTCTCCTTATTAGTTTTGTAGTGAGCTTACTAATCATGCTGAAACTATTCTATCAACCTAAAAAATCCCTTTATGGTAACGCTGAATTTGCCACAAAATCAGAAGTACAAAAAATGGGTTTATTGGAAGAAAAAAACGGGATTATTGTCGGCAAACTCGGTAATAAAc                                                                                                                                                                                             <  1:434275/251‑1 (MQ=255)    tCGTTTAATCGAAGCATTTAAGCGCTACTACCCTACTCTAACGCTGGATGCCATTTGGGGTAATTATCCGAAAATATGGCAGTCTCTCCTTATTAGTTTTGTAGTGAGCTTACTAATCATGCTGAAACTATTCTATCAACCTAAAAAATCCCTTTATGGTAACGCTGAATTTGCCACAAAATCAGAAGTACAAAAAATGGGTTTATTGGAAGAAAAAAACGGGATTATTGTCGGCAAACTCGGTAATAAAc                                                                                                                                                                                             <  1:367819/251‑1 (MQ=255)                 gCATTTAAGCGCTACTACCCTACTCTAACGCTGGATGCCATTTGGGGTAATTATCCGAAAATATGGCAGTCTCTCCTTATTAGTTTTGTAGTGAGCTTACTAATCATGCTGAAACTATTCTATCAACCTAAAAAATCCCTTTATGGTAACGCTGAATTTGCCACAAAATCAGAAGTACAAAAAATGGGTTTATTGGAAGAAAAAAACGGGATTATTGTCGGCAAACTCGGtaat                                                                                                                                                                                                 >  1:13636/1‑234 (MQ=255)                  cATTTAAGCGCTACTACCCTACTCTAACGCTGGATGCCATTTGGGGTAATTATCCGAAAATATGGCAGTCTCTCCTTATTAGTTTTGTAGTGAGCTTACTAATCATGCTGAAACTATTCTATCAACCTAAAAAATCCCTTTATGGTAACGCTGAATTTGCCACAAAATCAGAAGTACAAAAAATGGGTTTATTGGAAGAAAAAAACGGGATTATTGTCGGCAAACTCGGTAATAAActgctccgtttta                                                                                                                                                                                 <  1:565339/249‑12 (MQ=255)                   aTTTAAGCGCTACTACCCTACTCTAACGCTGGATGCCATTTGGGGTAATTATCCGAAAATATGGCAGTCTCTCCTTATTAGTTTTGTAGTGAGCTTACTAATCATGCTGAAACTATTCTATCAACCTAAAAAATCCCTTTATGGTAACGCTGAATTTGCCACAAAATCAGAAGTACAAAAAATGGGTTTATTGGAAGAAAAAAACGGGATTATTGTCGGCAAACTCGGTAATAAACTGCTCCGTTTTATCg                                                                                                                                                                              <  1:866902/251‑1 (MQ=255)                            ctactaCCCTACTCTAACGCTGGATGCCATTTGGGGTAATTATCCGAAAATATGGCAGTCTCTCCTTATTAGTTTTGTAGTGAGCTTACTAATCATGCTGAAACTATTCTATCAACCTAAAAAATCCCTTTATGGTAACGCTGAATTTGCCACAAAATCAGAAGTACAAGAAATGGGTTTATTGGAAGAAAAAAACGGGATTATTGTCGGCAAACTCGGTAATAAACTGCTCCGTTTTATCGGTGCGCAAt                                                                                                                                                                     >  1:446355/1‑251 (MQ=255)                            ctactaCCCTACTCTAACGCTGGATGCCATTTGGGGTAATTATCCGAAAATATGGCAGTCTCTCCTTATTAGTTTTGTAGTGAGCTTACTAATCATGCTGAAACTATTCTATCAACCTAAAAAATCCCTTTATGGTAACGCTGAATTTGCCACAAAATCAGAAGTACAAAAAATGGGTTTATTGGAAGAAAAAAACGGGATTATTGTCGGCAAACTCGGTAATAAACTGCTCCGTTTTATCGGTGCGCAAt                                                                                                                                                                     >  1:79817/1‑251 (MQ=255)                             tactacCCTACTCTAACGCTGGATGCCATTTGGGGTAATTATCCGAAAATATGGCAGTCTCTCCTTATTAGTTTTGTAGTGAGCTTACTAATCATGCTGAAACTATTCTATCAACCTAAAAAATCCCTTTATGGTAACGCTGAATTTGCCACAAAATCAGAAGTACAAAAAATGGGTTTATTGGAAGAAAAAAACGGGATTATTGTCGGCAAACTCGGTAATAAACTGCTCCGTTTTATCGGTGCGCAAtt                                                                                                                                                                    <  1:451936/251‑1 (MQ=255)                                 acCCTACTCTAACGCTGGATGCCATTTGGGGTAATTATCCGAAAATATGGCAGTCTCTCCTTATTAGTTTTGTAGTGAGCTTACTAATCATGCTGAAACTATTCTATCAACCTAAAAAATCCCTTTATGGTAACGCTGAATTTGCCACAAAATCAGAAGTACAAAAAATGGGTTTATTGGAAGAAAAAAACGGGATTATTGTCGGCAAACTCGGTAATAAACTGCTCCGTTTTATCGGTGCGCAAtttgtt                                                                                                                                                                >  1:750024/1‑251 (MQ=255)                                 acCCTACTCTAACGCTGGATGCCATTTGGGGTAATTATCCGAAAATATGGCAGTCTCTCCTTATTAGTTTTGTAGTGAGCTTACTAATCATGCTGAAACTATTCTATCAACCTAAAAAATCCCTTTATGGTAACGCTGAATTTGCCACAAAATCAGAAGTACAAAAAATGGGTTTATTGGAAGAAAAAAACGGGATTATTGTCGGCAAACTCGGTAATAAACTGCTCCGTTTTATCGGTGCGCAAtttgtt                                                                                                                                                                >  1:71594/1‑251 (MQ=255)                                      aCTCTAACGCTGGATGCCATTTGGGGTAATTATCCGAAAATATGGCAGTCTCTCCTTATTAGTTTTGTAGTGAGCTTACTAATCATGCTGAAACTATTCTATCAACCTAAAAAATCCCTTTATGGTAACGCTGAATTTGCCACAAAATCAGAAGTACAAAAAATGGGTTTATTGGAAGAAAAAAACGGGATTATTGTCGGCAAACTCGGTAATAAACTGCTCCGTTTTATCGGTGCGCAAtttgtttca                                                                                                                                                             >  1:296252/1‑248 (MQ=255)                                      aCTCTAACGCTGGATGCCATTTGGGGTAATTATCCGAAAATATGGCAGTCTCTCCTTATTAGTTTTGTAGTGAGCTTACTAATCATGCTGAAACTATTCTATCAACCTAAAAAATCCCTTTATGGTAACGCTGAATTTGCCACAAAATCAGAAGTACAAAAAATGGGTTTATTGGAAGAAAAAAACGGGATTATTGTCGGCAAACTCGGTAATAAACTGCTCCGTTTTATCGGTGCGCAAtttgtttca                                                                                                                                                             >  1:13814/1‑248 (MQ=255)                                      aCTCTAACGCTGGATGCCATTTGGGGTAATTATCCGAAAATATGGCAGTCTCTCCTTATTAGTTTTGTAGTGAGCTTACTAATCATGCTGAAACTATTCTATCAACCTAAAAAATCCCTTTATGGTAACGCTGAATTTGCCACAAAATCAGAAGTACAAAAAATGGGTTTATTGGAAGAAAAAAACGGGATTATTGTCGGCAAACTCGGTAATAAACTGCTCCGTTTTATCGGTGCGCAAtttgtttca                                                                                                                                                             >  1:184116/1‑248 (MQ=255)                                      aCTCTAACGCTGGATGCCATTTGGGGTAATTATCCGAAAATATGGCAGTCTCTCCTTATTAGTTTTGTAGTGAGCTTACTAATCATGCTGAAACTATTCTATCAACCTAAAAAATCCCTTTATGGTAACGCTGAATTTGCCACAAAATCAGAAGTACAAAAAATGGGTTTATTGGAAGAAAAAAACGGGATTATTGTCGGCAAACTCGGTAATAAACTGCTCCGTTTTATCGGTGCGCAAt                                                                                                                                                                     >  1:616356/1‑241 (MQ=255)                                       ctctAACGCTGGATGCCATTTGGGGTAATTATCCGAAAATATGGCAGTCTCTCCTTATTAGTTTTGTAGTGAGCTTACTAATCATGCTGAAACTATTCTATCAACCTAAAAAATCCCTTTATGGTAACGCTGAATTTGCCACAAAATCAGAAGTACAAAAAATGGGTTTATTGGAAGAAAAAAACGGGATTATTGTCGGCAAACTCGGTAATAAACTGCTCCGTTTTAt                                                                                                                                                                                <  1:665709/229‑1 (MQ=255)                                       ctctAACGCTGGATGCCATTTGGGGTAATTATCCGAAAATATGGCAGTCTCTCCTTATTAGTTTTGTAGTGAGCTTACTAATCATGCTGAAACTATTCTATCAACCTAAAAAATCCCTTTATGGTAACGCTGAATTTGCCACAAAATCAGAAGTACAAAAAATGGGTTTATTGGAAGAAAAAAACGGGATTATTGTCGGCAAACTCGGTAATAAACTGCTCCGTTTTATCGGTGCGCAAtttgtttca                                                                                                                                                             >  1:354455/1‑247 (MQ=255)                                       ctctAACGCTGGATGCCATTTGGGGTAATTATCCGAAAATATGGCAGTCTCTCCTTATTAGTTTTGTAGTGAGCTTACTAATCATGCTGAAACTATTCTATCAACCTAAAAAATCCCTTTATGGTAACGCTGAATTTGCCACAAAATCAGAAGTACAAAAAATGGGTTTATTGGAAGAAAAAAACGGGATTATTGTCGGCAAACTCGGTAATAAACTGCTCCGTTTTATCGGTGCGCAAtttgtttca                                                                                                                                                             >  1:294198/1‑247 (MQ=255)                                       ctctAACGCTGGATGCCATTTGGGGTAATTATCCGAAAATATGGCAGTCTCTCCTTATTAGTTTTGTAGTGAGCTTACTAATCATGCTGAAACTATTCTATCAACCTAAAAAATCCCTTTATGGTAACGCTGAATTTGCCACAAAATCAGAAGTACAAAAAATGGGTTTATTGGAAGAAAAAAACGGGATTATTGTCGGCAAACTCGGTAATAAACTGCTCCGTTTTATCGGTGCGCAAtttgtttca                                                                                                                                                             >  1:128522/1‑247 (MQ=255)                                       ctctAACGCTGGATGCCATTTGGGGTAATTATCCGAAAATATGGCAGTCTCTCCTTATTAGTTTTGTAGTGAGCTTACTAATCATGCTGAAACTATTCTATCAACCTAAAAAATCCCTTTATGGTAACGCTGAATTTGCCACAAAATCAGAAGTACAAAAAATGGGTTTATTGGAAGAAAAAAACGGGATTATTGTCGGCAAACTCGGTAATAAACTGCTCCGTTTTATCGGTGCGCAAtttgtttca                                                                                                                                                             >  1:54421/1‑247 (MQ=255)                                       ctctAACGCTGGATGCCATTTGGGGTAATTATCCGAAAATATGGCAGTCTCTCCTTATTAGATTTGTAGTGAGCTTACTAATCATGCTGAAACTATTCTATCAACCTAAAAAATCCCTTTATGGTAACGCTGAATTTGCCACAAAATCAGAAGTACAAAAAATGGGTTTATTGGAAGAAAAAAACGGGATTATTGTCGGCAAACTCGGTAATAAACTGCTCCGTTTTATCGGTGCGCAAtttgtttca                                                                                                                                                             >  1:523495/1‑247 (MQ=255)                                         ctAACGCTGGATGCCATTTGGGGTAATTATCCGAAAATATGGCAGTCTCTCCTTATTAGTTTTGTAGTGAGCTTACTAATCATGCTGAAACTATTCTATCAACCTAAAAAATCCCTTTATGGTAACGCTGAATTTGCCACAAAATCAGAAGTACAAAAAATGGGTTTATTGGAAGAAAAAAACGGGATTATTGTCGGCAAACTCGGTAATAAACTGCTCCGTTTTATCGGTGCGCAAtttgtttcactgg                                                                                                                                                         >  1:127311/1‑245 (MQ=255)                                            aCGCTGGATGCCATTTGGGGTAATTATCCGAAAATATGGCAGTCTCTCCTTATTAGTTTTGTAGTGAGCTTACTAATCATGCTGAAACTATTCTATCAACCTAAAAAATCCCTTTATGGTAACGCTGAATTTGCCACAAAATCAGAAGTACAAAAAATGGGTTTATTGGAAGAAAAAAACGGGATTATTGTCGGCAAACTCGGTAATAAACTGCTCCGTTTTATCgg                                                                                                                                                                             >  1:18733/1‑227 (MQ=255)                                             cGCTGGATGCCATTTGGGGTAATTATCCGAAAATATGGCAGTCTCTCCTTATTAGTTTTGTAGTGAGCTTACTAATCATGCTGAAACTATTCTATCAACCTAAAAAATCCCTTTATGGTAACGCTGAATTTGCCACAAAATCAGAAGTACAAAAAATGGGTTTATTGGAAGAAAAAAACGGGATTATTGTCGGCAAACTCGGTAATAAACTGCTCCGTTTTATCGGTGCGCAAtttgtttcactggcagcg                                                                                                                                                    >  1:738979/1‑241 (MQ=255)                                                              ggTAATTATCCGAAAATATGGCAGTCTCTCCTTATTAGTTTTGTAGTGAGCTTACTAATCATGCTGAAACTATTCTATCAACCTAAAAAATCCCTTTATGGTAACGCTGAATTTGCCACAAAATCAGAAGTACAAAAAATGGGTTTATTGGAAGAAAAAAACGGGATTATTGTCGGCAAACTCGGtaataa                                                                                                                                                                                               <  1:743348/191‑1 (MQ=255)                                                                                         ctcCTTATTAGTTTTGTAGTGAGCTTACTAATCATGCTGAAACTATTCTATCAACCTAAAAAATCCCTTTATGGTAACGCTGAATTTGCCACAAAATCAGAAGTACAAAAAATGGGTTTATTGGAAGAAAAAAACGGGATTATTGTCGGCAAACTCGGTAATAAAc                                                                                                                                                                                             <  1:99469/166‑1 (MQ=255)                                                                                           ccTTATTAGTTTTGTAGTGAGCTTACTAATCATGCTGAAACTATTCTATCAACCTAAAAAATCCCTTTATGGTAACGCTGAATTTGCCACAAAATCAGAAGTACAAAAAATGGGTTTATTGGAAGAAAAAAACGGGATTATTGTCGGCAAACTCGGTAATAAACTGCTCCGTTTTATCGGTGCGCAAtttgtttcactg                                                                                                                                                          >  1:600522/1‑195 (MQ=255)                                                                                                   gTTTTGTAGTGAGCTTACTAATCATGCTGAAACTATTCTATCAACCTAAAAAATCCCTTTATGGTAACGCTGAATTTGCCACAAAATCAGAAGTACAAAAAATGGGTTTATTGGAAGAAAAAAACGGGATTATTGTCGGCAAACTCGGTAATAAACTGCTCCGTTTTATCgg                                                                                                                                                                             <  1:348220/172‑1 (MQ=255)                                                                                                        gTAGTGAGCTTACTAATCATGCTGAAACTATTCTATCAACCTAAAAAATCCCTTTATGGTAACGCTGAATTTGCCACAAAATCAGAAGTACAAAAAATGGGTTTATTGGAAGAAAAAAACGGGATTATTGTCGGCAAACTCGGTAATAAAc                                                                                                                                                                                             >  1:675303/1‑151 (MQ=255)                                                                                                             gAGCTTACTAATCATGCTGAAACTATTCTATCAACCTAAAAAATCCCTTTATGGTAACGCTGAATTTGCCACAAAATCAGAAGTACAAAAAATGGGTTTATTGGAAGAAAAAAACGGGATTATTGTCGGCAAACTCGGtaat                                                                                                                                                                                                 <  1:771551/142‑1 (MQ=255)                                                                                                               gCTTACTAATCATGCTGAAACTATTCTATCAACCTAAAAAATCCCTTTATGGTAACGCTGAATTTGCCACAAAATCAGAAGTACAAAAAATGGGTTTATTGGAAGAAAAAAACGGGATTATTGTCGGCAAACTCGGTAATAAACTGCTCCGTTTTATCGGTGCGCAAtt                                                                                                                                                                    <  1:62394/169‑1 (MQ=255)                                                                                                                                        cTATCAACCTAAAAAATCCCTTTATGGTAACGCTGAATTTGCCACAAAATCAGAAGTACAAAAAATGGGTTTATTGGAAGAAAAAAACGGGATTATTGTCGGCAAACTCGGTAATAAActgctccgttt                                                                                                                                                                                   <  1:365133/129‑10 (MQ=255)                                                                                                                                               ccTAAAAAATCCCTTTATGGTAACGCTGAATTTGCCACAAAATCAGAAGTACAAAAAATGGGTTTATTGGAAGAAAAAAACGGGATTATTGTCGGCAAACTCGGTAATAAACTGCTCCGTTTTATCGGTGCGCAATTTGTTTCACTGGCAGCGCCGACACGGGGAGGTAAAAGCGTAGGGATTGTTATTCCAAATTTGCTGGCATGGGAACAATCTGTGGTTGTGAAGGACATTAAAAACGAATGCTATCg                                                  <  1:251483/251‑1 (MQ=255)                                                                                                                                                      aaTCCCTTTATGGTAACGCTGAATTTGCCACAAAATCAGAAGTACAAAAAATGGGTTTATTGGAAGAAAAAAACGGGATTATTGTCGGCAAACTCGGTAATAAACTGCTCCGTTTTATCGGTGCGCAATTTGTTTCACTGGCAGCGCCGACACGGGGAGGTAAAAGCGTAGGGATTGTTATTCCAAATTTGCTGGCATGGGAACAATCTGTGGTTGTGAAGGACATTAAAAACGAATGCTATCgtatcacc                                           <  1:673060/251‑8 (MQ=255)                                                                                                                                                      aaTCCCTTTATGGTAACGCTGAATTTGCCACAAAATCAGAAGTACAAAAAATGGGTTTATTGGAAGAAAAAAACGGGATTATTGTCGGCAAACTCGGTAATAAACTGCTCCGTTTTATCGGTGCGCAATTTGTTTCACTGGCAGCGCCGACACGGGGAGGTAAAAGCGTAGGGATTGTTATTCCAAATTTGCTGGCATGGGAACAATCTGTGGTTGTGAAGGACATTAAAAACGAATGCTATCgtatcacc                                           <  1:690776/251‑8 (MQ=255)                                                                                                                                                      aaTCCCTTTATGGTAACGCTGAATTTGCCACAAAATCAGAAGTACAAAAAATGGGTTTATTGGAAGAAAAAAACGGGATTATTGTCGGCAAACTCGGTAATAAACTGCTCCGTTTTATCGGTGCGCAATTTGTTTCACTGGCAGCGCCGACACGGGGAGGTAAAAGCGTAGGGATTGTTATTCCAAATTTGCTGGCATGGGAACAATCTGTGGTTGTGAAGGACATTAAAAACGAATGCTATCgtatcacc                                           <  1:819439/251‑8 (MQ=255)                                                                                                                                                           cTTTATGGTAACGCTGAATTTGCCACAAAATCAGAAGTACAAAAAATGGGTTTATTGGAAGAAAAAAACGGGATTATTGTCGGCAAACTCGGTAATAAACTGCTCCGTTTTATCGGTGCGCAATTTGTTTCACTGGCAGCGCCGACACGGGGAGGTAAAAGCGTAGGGATTGTTATTCCAAATTTGCTGGCATGGGAACAATCTGTGGTTGTGAAGGACATTAAAAACGAATgc                                                       <  1:681836/234‑2 (MQ=255)                                                                                                                                                           cTTTATGGTAACGCTGAATTTGCCACAAAATCAGAAGTACAAAAAATGGGTTTATTGGAAGAAAAAAACGGGATTATTGTCGGCAAACTCGGTAATAAACTGCTCCGTTTTATCGGTGCGCAATTTGTTTCACTGGCAGCGCCGACACGGGGAGGTAAAAGCGTAGGGATTGTTATTCCAAATTTGCTGGCATGGGAACAATCTGTGGTTGTGAAGGACATTAAAAACGAATGCTATCGTATCACCAGTaa                                      >  1:431897/1‑251 (MQ=255)                                                                                                                                                           cTTTATGGTAACGCTGAATTTGCCACAAAATCAGAAGTACAAAAAATGGGTTTATTGGAAGAAAAAAACGGGATTATTGTCGGCAAACTCGGTAATAAACTGCTCCGTTTTATCGGTGCGCAATTTGTTTCACTGGCAGCGCCGACACGGGGAGGTAAAAGCGTAGGGATTGTTATTCCAAATTTGCTGGCATGGGAACAATCTGTGGTTGTGAAGGACATTAAAAACGAATGCTATCGTATCACCAGTaa                                      <  1:81259/251‑1 (MQ=255)                                                                                                                                                           cTTTATGGTAACGCTGAATTTGCCACAAAATCAGAAGTACAAAAAATGGGTTTATTGGAAGAAAAAAACGGGATTATTGTCGGCAAACTCGGTAATAAACTGCTCCGTTTTATCGGTGCGCAATTTGTTTCACTGGCAGCGCCGACACGGGGAGGTAAAAGCGTAGGGATTGTTATTCCAAATTTGCTGGCATGGGAACAATCTGTGGTTGTGAAGGACATTAAAAACGAATGCTATCGTATCACCAGTaa                                      <  1:207452/251‑1 (MQ=255)                                                                                                                                                               aTGGTAACGCTGAATTTGCCACAAAATCAGAAGTACAAAAAATGGGTTTATTGGAAGAAAAAAACGGGATTATTGTCGGCAAACTCGGTAATAAACTGCTCCGTTTTATCGGTGCGCAATTTGTTTCACTGGCAGCGCCGACACGGGGAGGTAAAAGCGTAGGGATTGTTATTCCAAATTTGCTGGCATGGGAACAATCTGTGGTTGTGAAGGACATTAAAAACGAATGCTATCGTATCACCAGTAAAt                                    <  1:40726/249‑1 (MQ=255)                                                                                                                                                               aTGGTAACGCTGAATTTGCCACAAAATCAGAAGTACAAAAAATGGGTTTATTGGAAGAAAAAAACGGGATTATTGTCGGCAAACTCGGTAATAAACTGCTCCGTTTTATCGGTGCGCAATTTGTTTCACTGGCAGCGCCGACACGGGGAGGTAAAAGCGTAGGGATTGTTATTCCAAATTTGCTGGCATGGGAACAATCTGTGGTTGTGAAGGACATTAAAAACGAATGCTATCGTATCACCAGTAAATat                                  <  1:86244/251‑2 (MQ=255)                                                                                                                                                                   tAACGCTGAATTTGCCACAAAATCAGAAGTACAAAAAATGGGTTTATTGGAAGAAAAAAACGGGATTATTGTCGGCAAACTCGGTAATAAACTGCTCCGTTTTATCGGTGCGCAATTTGTTTCACTGGCAGCGCCGACACGGGGAGGTAAAAGCGTAGGGATTGTTATTCCAAATTTGCTGGCATGGGAACAATCTGTGGTTGTGAAGGACATTAAAAACGAATGCTATCGTATCACCAGTAAATatcgac                              <  1:201076/251‑6 (MQ=255)                                                                                                                                                                                 ccACAAAATCAGAAGTACAAAAAATGGGTTTATTGGAAGAAAAAAACGGGATTATTGTCGGCAAACTCGGTAATAAACTGCTCCGTTTTATCGGTGCGCAATTTGTTTCACTGGCAGCGCCGACACGGGGAGGTAAAAGCGTAGGGATTGTTATTCCAAATTTGCTGGCATGGGAACAATCTGTGGTTGTGAAGGACATTAAAAACGAATGCTATCGTATCACCAGTAAAt                                    <  1:761182/231‑1 (MQ=255)                                                                                                                                                                                        aTCAGAAGTACAAAAAATGGGTTTATTGGAAGAAAAAAACGGGATTATTGTCGGCAAACTCGGTAATAAACTGCTCCGTTTTATCGGTGCGCAATTTGTTTCACTGGCAGCGCCGACACGGGGAGGTAAAAGCGTAGGGATTGTTATTCCAAATTTGCTGGCATGGGAACAATCTGTGGTTGTGAAGGACATTAAAAACGAATGCTATCGTATCACCAGTAAATATCGACAAAAAATATTAGGGCaaaagg         >  1:163853/1‑249 (MQ=255)                                                                                                                                                                                             aaGTACAAAAAATGGGTTTATTGGAAGAAAAAAACGGGATTATTGTCGGCAAACTCGGTAATAAACTGCTCCGTTTTATCGGTGCGCAATTTGTTTCACTGGCAGCGCCGACACGGGGAGGTAAAAGCGTAGGGATTGTTATTCCAAATTTGCTGGCATGGGAACAATCTGTGGTTGTGAAGGACATTAAAAACGAATGCTATCGTATCACCAGTAAATATCGACAAAAAATATTAGGGCaaaaggtatat    <  1:877799/251‑8 (MQ=255)                                                                                                                                                                                               gTACAAAAAATGGGTTTATTGGAAGAAAAAAACGGGATTATTGTCGGCAAACTCGGTAATAAACTGCTCCGTTTTATCGGTGCGCAATTTGTTTCACTGGCAGCGCCGACACGGGGAGGTAAAAGCGTAGGGATTGTTATTCCAAATTTGCTGGCATGGGAACAATCTGTGGTTGTGAAGGACATTAAAAACGAATGCTATCGTATCACCAGTAAATATCGACAAAAAATATTAGGGCAAAAGGTATATCg  >  1:410405/1‑251 (MQ=255)                                                                                                                                                                                                tACAAAAAATGGGTTTATTGGAAGAAAAAAACGGGATTATTGTCGGCAAACTCGGTAATAAACTGCTCCGTTTTATCGGTGCGCAATTTGTTTCACTGGCAGCGCCGACACGGGGAGGTAAAAGCGTAGGGATTGTTATTCCAAATTTGCTGGCATGGGAACAATCTGTGGTTGTGAAGGACATTAAAAACGAATGCTATCGTATCACCAGTAAATATCGACAAAAAATATTAGGGCAAAAGGTATATCg  <  1:71111/250‑1 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                      AATTCGTTTAATCGAAGCATTTAAGCGCTACTACCCTACTCTAACGCTGGATGCCATTTGGGGTAATTATCCGAAAATATGGCAGTCTCTCCTTATTAGTTTTGTAGTGAGCTTACTAATCATGCTGAAACTATTCTATCAACCTAAAAAATCCCTTTATGGTAACGCTGAATTTGCCACAAAATCAGAAGTACAAAAAATGGGTTTATTGGAAGAAAAAAACGGGATTATTGTCGGCAAACTTGGTAATAAACTTCTTCGCTTTATCGGTGCACAATTCGTTTCGTTGGCCGCACCGACCCGAGGTGGGAAAGGAATAGGGATTGTCATACCGAATTTATTAGAGTGGCAACAATCAGTTGTTGTGAAAGATGTTAAAAACGAATGTTATCGCATTACCAGTAAATACCGTCAAAAAATATTGGGGCAAAAAGTCTATCG  >  1/1452887‑1453327

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