Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	Exported	522,457	4 bp→CC	coding (2536‑2539/2592 nt) coding (385‑388/462 nt)	BCHEFKBB_00548 → BCHEFKBB_00549 ←	restriction endonuclease EcoPI subunit R[NGO_0546] hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	522,457	0	T	.	95.6%	310.3 / 5.8	68	coding (2536/2592 nt) coding (388/462 nt)	BCHEFKBB_00548 BCHEFKBB_00549	restriction endonuclease EcoPI subunit R[NGO_0546] hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (3/0);  new base . (32/33);  total (35/33)
Fisher's exact test for biased strand distribution p-value = 2.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*	Exported	522,458	0	A	.	98.5%	333.6 / ‑1.2	68	coding (2537/2592 nt) coding (387/462 nt)	BCHEFKBB_00548 BCHEFKBB_00549	restriction endonuclease EcoPI subunit R[NGO_0546] hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base . (34/33);  minor base G (1/0);  total (35/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.82e-02
*	Exported	522,459	0	T	C	100.0%	265.1 / NA	68	Y846Y (TAT→TAC)  Y129R (TAT→CGT)	BCHEFKBB_00548 BCHEFKBB_00549	restriction endonuclease EcoPI subunit R[NGO_0546] hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (35/33);  total (35/33)
*	Exported	522,460	0	A	C	100.0%	264.5 / NA	68	M847L (ATG→CTG)  Y129R (TAT→CGT)	BCHEFKBB_00548 BCHEFKBB_00549	restriction endonuclease EcoPI subunit R[NGO_0546] hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (35/33);  total (35/33)

CTGGATTACCGGCTTGTCGGCAGCAAAATCCATCCGACCAAATTTACCAATAAAGACGGCAAACCGCGCGCGGTGAAAAAAGCAGATTTGGGCAGATTTGAAGATACGGAGCACCAGTATATGCCGC‑TACTATCTCTTCGGCGAGATTTTCTACGATTCGGATATAGAACATGAAAATGTCGCCAACAACCAAATTGAAGGCGTAATCGTATTTACCAAAATACCGAGAAACTCCATCAAAA  >  Exported/522340‑522581                                                                                                                      ||                                                                                                                             cTGGATTACCGGCTTGTCGGCAGCAAAATCCATCCGACCAAATTTACCAATAAAGACGGCAAACCGCGCGCGGTGAAAAAAGCAGATTTGGGCAGATTTGAAGATACGGAGCACCGG‑‑CCTGCCGCCGGCTATCTCTTCGGCGAGAtttt                                                                                              >  2:43121/1‑149 (MQ=255)   ggATTACCGGCTTGTCGGCAGCAAAATCCATCCGACCAAATTTACCCATAAAGACGGCAAACCGCGCGCGGTGAAAAAAGCAGATTTGGGCAGATTTGAAGATACGGAGCACCGG‑‑CCTGCCGCCGGCTATCTCTTCGGCGAGATTTTCTa                                                                                           >  1:730256/1‑150 (MQ=255)   ggATTACCGGCTTGTCGGCAGCAAAATCCATCCGACCAAATTTACCAATAAAGACGGCAAACCGCGCGCGGTGAAAAAAGCAGATTTGGGCAGATTTGAAGATACGGAGCACCGG‑‑CCTGCCGCCGGCTATCTCTTCGGCGAGATTTTCTa                                                                                           >  1:750573/1‑150 (MQ=255)   ggATTACCGGCTTGTCGGCAGCAAAATCCATCCGACCAAATTTACCAATAAAGACGGCAAACCGCGCGCGGTGAAAAAAGCAGATTTGGGCAGATTTGAAGATACGGAGCACCGG‑‑CCTGCCGCCGGCTATCTCTTCGGCGAGATTTTCTa                                                                                           >  2:579802/1‑150 (MQ=255)   ggATTACCGGCTTGTCGGCAGCAAAATCCATCCGACCAAATTTACCAATAAAGACGGCAAACCGCGCGCGGTGAAAAAAGCAGATTTGGGCAGATTTGAAGATACGGAGCACCGG‑‑CCTGCCGCCGGCTATCTCTTCGGCGAGATTTTCTa                                                                                           >  2:423233/1‑150 (MQ=255)   ggATTACCGGCTTGTCGGCAGCAAAATCCATCCGACCAAATTTACCAATAAAGACGGCAAACCGCGCGCGGTGAAAAAAGCAGATTTGGGCAGATTTGAAGATACGGAGCACCGG‑‑CCTGCCGCCGGCTATCTCTTCGGCGAGATTTTCTa                                                                                           >  1:374039/1‑150 (MQ=255)    gATTACCGGCTTGTCGGCAGCAAAATCCATCCGACCAAATTTACCACGAAAGACGGCAAACCGCGCGCGGTGAAAAAACCAGATTTGGGCAGATTGGAAGATACGGAGCACC‑‑TGCCTGCCGCCGGCTATCTCTTCGGCGAGATTTTCTACg                                                                                         >  2:1013002/1‑151 (MQ=255)    gATTACCGGCTTGTCGGCAGCAAAATCCATCCGACCAAATTTACCAATAAAGACGGCAAACCGCGCGCGGTGAAAAAAGCAGATTTGGGCAGATTTGAAGATACGGAGCACCGG‑‑CCTGCCGCCGGCTATCTCTTCGGCGAGATTTTCTACg                                                                                         >  2:12109/1‑151 (MQ=255)              ttGTCGGCAGCAAAATCCATCCGACCAAATTTACCAATAAAGACGGCAAACCGCGCGCGGTGAAAAAAGCAGATTTGGGCAGATTTGAAGATACGGAGCACCGG‑‑CCTGCCGCCGGCTATCTCTTCGGCGAGATTTTCTACGATTCGGatat                                                                               <  1:421446/151‑1 (MQ=255)                 tCGGCAGCAAAATCCATCCGACCAAATTTACCAATAAAGACGGCAAACCGCGCGCGGTGAAAAAAGCAGATTTGGGCAGATTTGAAGATACGGAGCACCGG‑‑CCTGCCGCCGGCTATCTCTTCGGCGAGATTTTCTACGATTCGGATATAg                                                                             <  1:253643/150‑1 (MQ=255)                  cGGCAGCAAAATCCATCCGACCAAATTTACCAATAAAGACGGCAAACCGCGCGCGGTGAAAAAAGCAGATTTGGGCAGATTTGAAGATACGGAGCACCGG‑‑CCTGCCGCCGGCTATCTCTTCGGCGAGATTTTCTACGATTCGGATATAGaa                                                                           <  2:1124075/151‑1 (MQ=255)                  cGGCAGCAAAATCCATCCGACCAAATTTACCAATAAAGACGGCAAACCGCGCGCGGTGAAAAAAGCAGATTTGGGCAGATTTGAAGATACGGAGCACCGG‑‑CCTGCCGCCGGCTATCTCTTCGGCGAGATTTTCTACGATTCGGATATAGaa                                                                           <  1:150737/151‑1 (MQ=255)                  cGGCAGCAAAATCCATCCGACCAAATTTACCAATAAAGACGGCAAACCGCGCGCGGTGAAAAAAGCAGATTTGGGCAGATTTGAAGATACGGAGCACCGG‑‑CCTGCCGCCGGCTATCTCTTCGGCGAGATTTTCTACGATTCGGATATAGaa                                                                           <  1:796045/151‑1 (MQ=255)                    gcagcaAAATCCATCCGACCATATTTACCAATAAAGACGGCAAACCGCGCGCGGTGAAAAAAGCAGATTTGGGCAGATTTGAAGATACGGAGCACCGG‑‑CCTGCCGCCGGCTATCTCTTCGGGGAGATTTTCTAGGATTCGGATATAGAAc                                                                          <  2:13933/150‑1 (MQ=255)                       gcaAAATCCATCCGACCAAATTTACCAATAAAGACGGCAAACCGCGCGCGGTGAAAAAAGCAGATTTGGGCAGATTTGAAGATACGGAGCACCGG‑‑CCTGCCGCCGGCTATCTCTTCGGCGAGATTTTCTACGATTCGGATATAGAACATGa                                                                      >  1:338372/1‑151 (MQ=255)                       gcaAAATCCATCCGACCAAATTTACCAATAAAGACGGCAAACCGCGCGCGGTGAAAAAAGCAGATTTGGGCAGATTTGAAGATACGGAGCACCGG‑‑CCTGCCGCCGGCTATCTCTTCGGCGAGATTTTCTACGATTCGGATATAGAACATGa                                                                      >  2:648763/1‑151 (MQ=255)                       gcaAAATCCATCCGACCAAATTTACCAATAAAGACGGCAAACCGCGCGCGGTGAAAAAAGCAGATTTGGGCAGATGTGAAGATACGGAGCACCGG‑‑CCTGCCGCCGGCTATCTCTTCGGCGAGATTTTCTACGATTCGGatcgagaccatgt                                                                      >  1:238790/1‑140 (MQ=255)                        caAAATCCATCCGACCAAATTTACCAATAACGACGGCAAACCGCGCGCGGTGAAAAAAGCAGATTTGGGCAGATTTGAAGATACGGCGCACCGG‑‑CCTGCCGCCGGCTATCTCTTCGGCGAGATTTTCTACGATTCGGATATAGACCATgca                                                                     <  2:1086813/151‑3 (MQ=255)                        caAAATCCATCCGACCAAATTTACCAATAAAGACGGCAAACCGCGCGCGGTGAAAAAAGCAGATTTGGGCAGATTTGAAGATACGGAGCACCGG‑‑CCTGCCGCCGGCTATCTCTTCGGCGAGATTTTCTACGATTCGGATATAGAACATGaa                                                                     <  1:944321/151‑1 (MQ=255)                          aaaaccatccGACCAAATTTACCAATAAAGACGTCAAACCGCGCGCGGTGAAAAAATCAGAATTGGGAAGATTTGAAGATACGGAGCACCGG‑‑CCTGCCGCCGGCTATCTCTTCGGAGAGATTTTCTACGATTCGGATATAGAACAAGaaaa                                                                   >  2:934107/5‑151 (MQ=255)                          aaaTCCATCCGACCAAATTTACCAATAAAGACGGCAAACCGCGCGCGGTGAAAAAAGCAGATTTGGGCAGATTTGAAGATACGGAGCACCGG‑‑CCTGCCGCCGGCTATCTCTTCGGCGAGATTTTCTACGATTCGGATATAGAACATGaaaa                                                                   <  2:998984/151‑1 (MQ=255)                          aaaTCCATCCGACCAAATTTACCAATAAAGACGGCAAACCGCGCGCGGTGAAAAAAGCAGATTTGGGCAGATTTGAAGATACGGAGCACCGG‑‑CCTGCCGCCGGCTATCTCTTCGGCGAGATTTTCTACGATTCGGATATAGAACATGaaaa                                                                   >  1:1206589/1‑151 (MQ=255)                          aaaTCCATCCGACCAAATTTACCAATAAAGACGGCAAACCGCGCGCGGTGAAAAAAGCAGATTTGGGCAGATTTGAAGATACGGAGCACCGG‑‑CCTGCCGCCGGCTATCTCTTCGGCGAGATATTCTACGATTCGGATATAGAACATGaaaa                                                                   >  2:934086/1‑151 (MQ=255)                               catccGACCAAATTTACCAATAAAGACGGCAAACCGCGCGCGGAGAAAAAAGCAGATTTGGGCAGATTTGAAGATACGGAGCACCGG‑‑CCTGCAGCCGGCTATATCTTCGGCGAGATTTTCTACGATTCGGATATAGAACCTGAAAATGTCg                                                              >  2:1165624/1‑151 (MQ=255)                                 tccGACCAAATTTACCAAGACAGACGGCAAACCGCGCGCGGTGAAAAAAGCAGATTTGGGCAGATTTGAAGATACGGAGCACCGG‑‑CCTGCCGCCGGCTATCTCTTCGGCGAGATTTTCTACGATTCGGATATAGAACATGAAAATGTCGcc                                                            <  1:411986/151‑1 (MQ=255)                                    gACCAAATTTACCAATAAAGACGGCAAACCGCGCGCGGTGAAAAAAGCAGATTTGGGCAGATTTGAAGATACGGAGCTCC‑GT‑CCTGCCGCCGGCTATCTCTTCGGCGAGATTTTCTACTATTCGGATATAGAACATGAAAATGTCGCcaa                                                          >  2:1343751/1‑150 (MQ=255)                                    gACCAAATTTACCAATAAAGACGGCAAACCGCGCGCGGTGAAAAAAGCAGATTTGGGCAGAATTGAAGATACGGAGCACC‑GT‑CCTGCCGCCGGCTATATCTTCGGCGAGATATTCTACGATTCGGATATAGAACATGAAAATGTCGCcaa                                                          >  2:1343750/1‑150 (MQ=255)                                        aaaTTTACCAATAAAGACGGCAAACCGCGCGCGGTGAAAAAAGCAGATTTGGGCAGATTTGAAGATACGGAGCACCGG‑‑CCTGCCGCCGGCTATCTCTTCGGCGAGATTTTCTACGATTCGGATATAGAACATGAAAATGTCGCCaacaac                                                      >  1:658788/1‑150 (MQ=255)                                        aaaTTTACCAATAAAGACGGCAAACCGCGCGCGGTGAAAAAAGCAGATTTGGGCAGATTTGAAGATACGGAGCACCGG‑‑CCTGCCGCCGGCTATCTCTTCGGCGAGATTTTCTACGATTCGGATATAGAACATGAAAATGTCGCCaacaac                                                      >  2:1041347/1‑150 (MQ=255)                                         aaTTTACCAATAAAGACGGCAAACCGCGCGCGGTGAAAAAAGCAGATTTGGGCAGATTTGAAGATACGGAGCACCGG‑‑CCTGCCGCCGGCTATCTCTTCGGCGAGATTTTCTACGATTCGGATATAGAACATGAAAATGTCGCCAACAACCa                                                    >  1:627941/1‑151 (MQ=255)                                            ttACCAATAAAGACGGCAAACCGCGCGCGGTGAAAAAAGCAGATTTGGGCAGATTTGAAGATACGGAGCACCGG‑‑CCTGCCGCCGGCTATCTCTTCGGCGAGATTTTCTACGATTCGGATATAGAACATGAAAATGTCGCCAACAACCAAAt                                                 <  2:268728/151‑1 (MQ=255)                                            ttACCAATAAAGACGGCAAACCGCGCGCGGTGAAAAAAGCAGATTTGGGCAGATTTGAAGATACGGAGCACCGG‑‑CCTGCCGCCGGCTATCTCTTCGGCGAGATTTTCTACGATTCGGAGATAGAACATGAAAATGTCGCCAACAACCAAAt                                                 <  2:985501/151‑1 (MQ=255)                                            tgaCCAATAAAGACGGCAAACCGCGCGCGGTGAAAAAAGCAGATTTGGGCAGATTTGAAGATACGGAGCACCGG‑‑CCTGCCGCCGGCTATCTCTTCGGCGAGATTTTCTACGATTCGGATATAGAACATGAAAATGTCGCCAACAACCAAAt                                                 <  2:177422/149‑1 (MQ=255)                                               ccAATAAAGACGGCAAACCGCGCGCGGTGAAAAAAGCAGATTTGGGCAGATTTGAAGATACGGAGCACCGG‑‑CCTGCCGCCGGCTATCTCTTCGGCGAGATTTTCTACGATTCGGATATAGAACATGAAAATGTCGCCAACAACCAAATTGa                                              >  1:1289302/1‑151 (MQ=255)                                                 aagaaaGACGGCAAACCGCGCGCGGTGAAAAAAGGAGATGTGGGCCGATTTGAAGATACGGAGCACCGG‑‑CCTGCCGCCGGCTATCGCTTCGGCGAGATTTTCTACGTGTCGGATATAGAACTTGGCTATGTCGCCAGCAACCAAATTGAAg                                            >  2:884107/4‑151 (MQ=255)                                                 aaTAAAGACGGCAAACCGCGCGCGGTGAAAAAAGCAGATTTGGGCAGATTTGAAGATACGGAGCACCGG‑‑CCTGCCGCCGGCTATCTCTTCGGCGAGATTTTCTACGATTCGGATATAGAACATGAAAATGTCGCCAACAACCAAATTGAAg                                            >  1:856423/1‑151 (MQ=255)                                                  aTAAAGACGGCAGACCGCGCGCGGTGAAAAAAGCAGATTTGGGCAGATTTGAAGATACGGAGCACCGG‑‑CCTGCCGCCGGCTATCTCTTCGGCGAGATTTTCTACGATTCGGATATAGAACATGAAAATGTCGCCAACAACCAAATTGAAgg                                           <  2:449062/151‑1 (MQ=255)                                                   tAAAGACGGCAAACCGCGCGCGGTGAAAAAAGCAGATTTGGGCAGATTTGAAGATACGGAGCACCGG‑‑CCTGCCGCCGGCTATCTCTTCGGCGAGATTTTCTACGATTCGGATATAGAACATGAAAATGTCGCCAACAACCAAATTGAAgg                                           >  2:1189666/1‑150 (MQ=255)                                                     aaGACGGCAAACCGCGCGCGGTGAAAAAAGCAGATTTGGGCAGATTTGAAGATACGGAGCACCGG‑‑CCTGCCGCCGGCTATCTCTTCGGCGAGATTTTCTACGATTCGGATATAGAACATGAAAATGTCGCCAACAACCAAATTGAAGGCGt                                        >  1:249594/1‑151 (MQ=255)                                                       gACGGCAAACCGCGCGCGGTGAAAAAAGCAGATTTGGGCAGATTTGAAGATACGGAGCACCGG‑‑CCTGCCGCCGGCTATCTCTTCGGCGAGATTTTCTACGATTCGGATATAGAACATGAAAATGTCGCCAACAACCAAATTGAAGGCGTaa                                      >  1:598295/1‑151 (MQ=255)                                                          ggCAAACCGCGCGCGGTGAAAAAAGCAGATTTGGGCAGATTTGAAGATTCGGAGCACCGG‑‑CCTGCCGCCGGCTATCTCTTCGGCGAGATTTTCTACGATTCGGATATAGAACATGAAAATGTCGCCATCAACCAAATTGAAGGCGTAATCg                                   >  2:158273/1‑151 (MQ=255)                                                           gCAAACCGCGCGCGGTGAAAACAGCAGATTTGGGCAGATTTGAAGATACGGAGCACCGG‑‑CCTGCCGCCGGCTATCTCTTCGGCGAGATTTTCTACGATTCGGATATAGAACATGAAAATGTCGCCAACAACCAAATTGAAGGCGTAATCg                                   <  2:1101629/150‑1 (MQ=255)                                                           gCAAACCGCGCGCGGTGAAAAAAGCAGATTTGGGCAGATTTGAAGATACGGAGCACCGG‑‑CCTGCCGCCGGCTATCTCTTCGGCGAGATTTTCTACGATTCGGATATAGAACATGAAAATGTCTCCAACAACCAAATTGAAGGCGTAATCGt                                  >  1:1209900/1‑151 (MQ=255)                                                           gCAAACCGCGCGCGGTGAAAAAAGCAGATTTGGGCAGATTTGAAGATACGGAGCACCGG‑‑CCTGCCGCCGGCTATCTCTTCGGCGAGATTTTCTACGATTCGGATATAGAACATGAAAATGTCGCCAACAACCAAATTGAAGGCGTAATCg                                   <  2:373927/150‑1 (MQ=255)                                                           gCAAACCGCGCGCGGTGAAAAAAGCAGATTTGGGCAGATTTGAAGATACGGAGCACCGG‑‑CCTGCCGCCGGCTATCTCTTCGGCGAGATTTTCTACGATTCGGATATAGAACATGAAAATGTCGCCAACAACCAAATTGAAGGCGTAATCg                                   <  1:1306433/150‑1 (MQ=255)                                                           gCAAACCGCGCGCGGTGAAAAAAGCAGATTTGGGCAGATTTGAAGATACGGAGCACCGG‑‑CCTGCCGCCGGCTATCTCTTCGGCGAGATTTTCTACGATTCGGATATAGAACATGAAAATGTCGCCAACAACCAAATTGAAGGCGTAATCg                                   <  2:760151/150‑1 (MQ=255)                                                           gCAAACCGCGCGCGGTGAAAAAAGCAGATTTGGGCAGATTTGAAGATACGGAGCACCGG‑‑CCTGCCGCCGGCTATCTCTTCGGCGAGATTTTCTACGATTCGGATATAGAACATGAAAATGTCGCCAACAACCAAATTGAAGGCGTAATCGt                                  >  2:732405/1‑151 (MQ=255)                                                             aaaCCGCGCGCGGTGAAAAAAGCAGATTTGGGCAGATTTGAAGATACGGAGCACCGG‑‑CCTGCCGCCGGCTATCTCTTCGGCGAGATTTTCTACGATTCGGATATAGAACATGAAAATGTCGCCAACAACCAAATTGAAGGCGTAATCGTAt                                <  1:1115915/151‑1 (MQ=255)                                                             aaaCCGCGCGCGGTGAAAAAAGCAGATTTGGGCAGATTTGAAGATACGGAGCACCGG‑‑CCTGCCGCCGGCTATCTCTTCGGCGAGATTTTCTACGATTCGGATATAGAACATGAAAATGTCGCCAACAACCAAATTGAAGGCGTAATCGTAt                                <  1:826542/151‑1 (MQ=255)                                                                       cgGTGAAAAAAGCAGATTTGGGCAGATTTGAAGATACGGAGCACCGG‑‑CCTGCCGCCGGCTATCTCTTCGGCGAGATTTTCTACGATTCGGATATAGAACATGAAAATGTCGCCAACAACCAAATTGAAGGCGTAATCGTATTTACCAAAAt                      <  2:534038/151‑1 (MQ=255)                                                                       cgGTGAAAAAAGCAGATTTGGGCAGATTTGAAGATACGGAGCACCGG‑‑CCTGCCGCCGGCTATCTCTTCGGCGAGATTTTCTACGATTCGGATATAGAACATGAAAATGTCGCCAACAACCAAATTGAAGGCGTAATCGTATTTACCAAAAt                      <  2:534065/151‑1 (MQ=255)                                                                       cgGTGAAAAAAGCAGATGTGGGCAGATTTGAAGATACGGAGCACCGG‑‑CCTGCCGCCGGCTATCTCTTCGGCGAGATTTTCTACGATTCGGATATAGAACATGAAAATGTCGCCAACAACCAAATTGAAGGCGTAATCGTATTTACCAAAAt                      <  2:1112439/151‑1 (MQ=255)                                                                          tGAATAAAGCAGATTTGGGCTGATTTGAAGATACGGAGCACCGG‑‑CCTGCCGCCGGCTATCTCTTCGGCGAGATTTTCTACTATTCGGATATAGAACATGAAAATGTCGCCAACAACCAAATTGAAGGCGTAAGCGTATTTACCAAAATAc                    <  2:91090/150‑1 (MQ=255)                                                                          tGAAAAAAGCAGATTTGGGCAGATTTGAAGATACGGAGCACCGG‑‑CCTGCCGCCGGCTATCTCTTCGGCGAGATTTTCTACGATTCGGATATAGAACATGAAAATGTCGCCAACAACCAAATTGAAGGCGTAATCGTATTTACCAAAATAc                    <  1:1181436/150‑1 (MQ=255)                                                                          tGAAAAAAGCAGATTTGGGCAGATTTGAAGATACGGAGCACCGG‑‑CCTGCCGCCGGCTATCTCTTCGGCGAGATTTTCTACGATTCGGATATAGAACATGAAAATGTCGCCAACAACCAAATTGAAGGCGTAATCGTATTTACCAAAATAc                    <  1:417839/150‑1 (MQ=255)                                                                          tGAAAAAAGCAGATTTGGGCAGATTTGAAGATACGGAGCACCGG‑‑CCTGCCGCCGGCTATCTCTTCGGCGAGATTTTCTACGATTCGGATATAGAACATGAAAATGTCGCCAACAACCAAATTGAAGGCGTAATCGTATTTACCAAAATAc                    <  2:91086/150‑1 (MQ=255)                                                                          tGAAAAAAGCAGATTTGGGCAGATTTGAAGATACGGAGCACCGG‑‑CCTGCCGCCGGCTATCTCTTCGGCGAGATTTTCTACGATTCGGATATAGAACATGAAAATGTCGCCAACAACCAAATTGAAGGCGTAATCGTATTTACCAAAATAc                    <  2:91083/150‑1 (MQ=255)                                                                          tGAAAAAAGCAGATTTGGGCAGATTTGAAGATACGGAGCACCGG‑‑CCTGCCGCCGGCTATCTCTTCGGCGAGATTTTCTACGATTCGGATATAGAACATGAAAATGTCGCCAACAACCAAATTGAAGGCGTAATCGTATTTACCAAAATAc                    <  1:352557/150‑1 (MQ=255)                                                                                   cAGATTTGGGCAGATTTGAAGATACGGAGCACCGG‑‑CCTGCCGCCGGCTATCTCTTCGGCGAGATTTTCTACGATTCGGATATAGAACATGAAAATGTCGCCAACAACCAAATTGAAGGCGTAATCGTATTTACCAAAATACCGAGAAAc            <  2:566022/149‑1 (MQ=255)                                                                                     gagttGGGCAGATTTGAAGATACGGAGCACCGG‑‑CCTGCCGCCGGCGATCTCTGCGGCGAGAATTTCTACGATTCGGATATAGAACATGAAAATGTCGCCAACCACCAAATTGAAGGCTTAATCGTATTTACCAAAATACCGAGAAACTCCa        >  2:537924/4‑151 (MQ=255)                                                                                     gATTTGGGCAGATTTGAAGATACGGAGCACCGG‑‑CCTGCCGCCGGCTATCTCTTCGGCGAGATTTTCTACGATTCGTATATAGACCATGAAAATGTCGCCAACAACCAAATTGAAGGCGTAATCTTATTTACAAAAATACCGAGAAACTCCa        >  2:964244/1‑151 (MQ=255)                                                                                     gATTTGGGCAGATTTGAAGATACGGAGCACCGG‑‑CCTGCCGCCGGCTATCTCTTCGGCGAGATTTTCTACGATTCGGATATAGAACATGAAAATGTCGCCAACAACCAAATTGAAGGCGTAATCGTATTTACCAAAATACCGAGAAACTCCa        >  1:520557/1‑151 (MQ=255)                                                                                     gATTTGGGCAGATTTGAAGATACGGAGCACCGG‑‑CCTGCCGCCGGCTATCTCTTCGGCGAGATTTTCTACGATTCGGATATAGAACATGAAAATGTCGCCAACAACCAAATTGAAGGCGTAATCGTATTTACCAAAATACCGAGAAACTCCa        >  1:520551/1‑151 (MQ=255)                                                                                      aTTTGGGCAGATTTGAAGATACGGAGCACCGG‑‑CCTGCCGCCGGCTATCTCTTCGGCGAGATTTTCTACGATTCGGATATAGAACATGAAAATGTCGCCAACAACCAAATTGAAGGCGTAATCGTATTTACCAAAATACCGAGAAACTc          <  2:1271864/148‑1 (MQ=255)                                                                                      aTTTGGGCAGATTTGAAGATACGGAGCACCGG‑‑CCTGCCGCCGGCTATCTCTTCGGCGAGATTTTCTACGATTCGGATATAGAACATGAAAATGTCGCCAACAACCAAATTGAAGGCGTAATCGTATTTACCAAAATACCGAGAAACTc          <  2:338273/148‑1 (MQ=255)                                                                                         tGGGCAGATTTGAAGATACGGAGCACCGG‑‑CCTGCCGCCGGCTATCTCTTCGGCGAGATTTTCTACGATTCGGATATAGAACATGAAAATGTCGCCAACAACCAAATTGAAGGCGTAATCGTATTTACCAAAATACCGAGAAACTCCATCaa    >  1:1155869/1‑151 (MQ=255)                                                                                          gggCAGATTTGAAGATACGGAGCACCGG‑‑CCTGCCGCCGGCTATCTCTTCGGCGAGATTTTCTACGATTCGGATATAGAACATGAAAATGTCGCCAACAACCAAATTGAAGGCGTAATCGTATTTACCAAAATACCGAGAAACTCCATCaaa   >  1:860414/1‑151 (MQ=255)                                                                                            gCAGATTTGAAGATACGGAGCACCGG‑‑CCTGCCGCCGGCTATCTCTTCGGCGAGATTTTCTACGATTCGGATATAGAACATGAAAATGTCGCCAACAACCAAATTGAAGGCGTAATCGTATTTACCAAAATACCGAGAAACTCCATCaaaa  <  2:982795/150‑1 (MQ=255)                                                                                                                      ||                                                                                                                             CTGGATTACCGGCTTGTCGGCAGCAAAATCCATCCGACCAAATTTACCAATAAAGACGGCAAACCGCGCGCGGTGAAAAAAGCAGATTTGGGCAGATTTGAAGATACGGAGCACCAGTATATGCCGC‑TACTATCTCTTCGGCGAGATTTTCTACGATTCGGATATAGAACATGAAAATGTCGCCAACAACCAAATTGAAGGCGTAATCGTATTTACCAAAATACCGAGAAACTCCATCAAAA  >  Exported/522340‑522581

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