Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	1,483,195	5 bp→TAACA	coding (144‑148/1281 nt)	unknown →	Type II/IV secretion system ATP hydrolase TadA/VirB11/CpaF, TadA subfamily

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	1,483,195	0	C	T	100.0%	123.5 / NA	47	S48S (AGC→AGT)	unknown	Type II/IV secretion system ATP hydrolase TadA/VirB11/CpaF, TadA subfamily
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (20/27);  total (20/27)
*	1	1,483,196	0	G	A	100.0%	124.9 / NA	47	G49N (GGT→AAC)	unknown	Type II/IV secretion system ATP hydrolase TadA/VirB11/CpaF, TadA subfamily
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (20/27);  total (20/27)
*	1	1,483,197	0	G	A	100.0%	123.9 / NA	47	G49N (GGT→AAC)	unknown	Type II/IV secretion system ATP hydrolase TadA/VirB11/CpaF, TadA subfamily
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (20/27);  total (20/27)
*	1	1,483,198	0	T	C	100.0%	133.9 / NA	47	G49N (GGT→AAC)	unknown	Type II/IV secretion system ATP hydrolase TadA/VirB11/CpaF, TadA subfamily
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (20/27);  total (20/27)
*	1	1,483,199	0	G	A	97.9%	119.2 / ‑5.7	47	A50T (GCT→ACT)	unknown	Type II/IV secretion system ATP hydrolase TadA/VirB11/CpaF, TadA subfamily
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (20/26);  minor base T (0/1);  total (20/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

AAAAACAGAGTTTATTTTCTTCATTTTTTTCAAAACTCATTAAACAGTAAATAGTGGCTAAACATGCTGACCAAACAACAAAAAATCTTACTTAGAAGTGAAGTACTAAGTAATCTTGATATTGAAAAGATTGATGAACTTCAAAGTGAAAGAAGCTCACTGGTAAACGAACTGGTTCAGATTGTTAACCGTGTTGCCAATAAAAGCGGTGCTTATTTGACCTCAGCCGATACTCTAGTTATGGCTGAAATTGTTGCTGATGAGATTGAGGGTTATGGACCATTACGAGATCTCATGGCTGACGATACTATAAATGATATTTTAGTTAATGGCCCTAATGACATCTGGGTTGAACGTGCTGGTATTCTGGAAAAAACTGATAAAGAATTTGTTAGTAATGAACAATTAACTGATATTGCA  >  1/1482989‑1483408                                                                                                                                                                                                               |||||                                                                                                                                                                                                                  aaaaaCAAAGTTTATTTTCCTCATTTTTTACAAAACTCATTAAACAATAAGTAGTGGTTAAACATGCTGACCAAACAACAAAAAATCTTACTTAGAAGTGAAGTATTGAGCAATCTTGATGTTGAAAAGATTGATGAACTTCAAAGTGAAAGAAGCTCATTAGTAAACGAACTGGTTCAGATTGTTAATCGTATTGCCAATAAAAGTAACTCTTATTTGACGTCATCCGATGCACTGGTTATGGCTGAAAt                                                                                                                                                                           <  1:73968/251‑1 (MQ=255)          gTTTATTTTCCTCATTTTTTACAAAACTCATTAAACAATAAGTAGTGGTTAAACATGCTGACCAAACAACAAAAAATCTTACTTAGAAGTGAAGTATTGAGCAATCTTGATGTTGAAAAGATTGATGAACTTCAAAGTGAAAGAAGCTCATTAGTAAACGAACTGGTTCAGATTGTTAATCGTATTGCCAATAAAAGTAACACTTATTTGACGTCATCCGATGCACTGGTTATGGCTGAAATTGTTGCtga                                                                                                                                                                  >  1:62286/1‑251 (MQ=255)                            ttacAAAACTCATTAAACAATAAGTAGTGGTTAAACATGCTGACCAAACAACAAAAAATCTTACTTAGAAGTGAAGTATTGAGCAATCTTGATGTTGAAAAGATTGATGAACTTCAAAGTGAAAGAAGCTCATTAGTAAACGAACTGGTTCAGATTGTTAATCGTATTGCCAATAAAAGTAACACTTATTTGACGTCATCCGATGCACTGGTTATGGCTGAAATTGTTGCTGATGAGATTGAGGGTTATgg                                                                                                                                                <  1:61077/248‑1 (MQ=255)                                     tCATTAAACAATAAGTAGTGGTTAAACATGCTGACCAAACAACAAAAAATCTTACTTAGAAGTGAAGTATTGAGCAATCTTGATGTTGAAAAGATTGATGAACTTCAAAGTGAAAGAAGCTCATTAGTAAACGAACTGGTTCAGATTGTTAATCGTATTGCCAATAAAAGTAACACTTATTTGACGTCATCCGATGCACTGGTTATGGCTGAAATTGTTGCTGATGAGATTGAg                                                                                                                                                        <  1:391690/234‑1 (MQ=255)                                     tCATTAAACAATAAGTAGTGGTTAAACATGCTGACCAAACAACAAAAAATCTTACTTAGAAGTGAAGTATTGAGCAATCTTGATGTTGAAAAGATTGATGAACTTCAAAGTGAAAGAAGCTCATTAGTAAACGAACTGGTTCAGATTGTTAATCGTATTGCCAATAAAAGTAACACTTATTTGACGTCATCCGATGCACTGGTTATGGCTGAAATTGTTGCTGATGAGATTGAGGGTTATGGTCCGTTACg                                                                                                                                       >  1:646085/1‑251 (MQ=255)                                             caataagtAGTGGTTAAACATGCTGACCAAACAACAAAAAATCTTACTTAGAAGTGAAGTATTGAGCAATCTTGATGTTGAAAAGATTGATGAACTTCAAAGTGAAAGAAGCTCATTAGTAAACGAACTGGTTCAGATTGTTAATCGTATTGCCAATAAAAGTAACACTTATTTGACGTCATCCGATGCACTGGTTATGGCTGAAATTGTTGCTGATGAGATTGAGGGTTATGGTCCGTTACGGGAtctcg                                                                                                                               >  1:386819/8‑250 (MQ=255)                                             caataagtAGTGGTTAAACATGCTGACCAAACAACAAAAAATCTTACTTAGAAGTGAAGTATTGAGCAATCTTGATGTTGAAAAGATTGATGAACTTCAAAGTGAAAGAAGCTCATTAGTAAACGAACTGGTTCAGATTGTTAATCGTATTGCCAATAAAAGTAACACTTATTTGACGTCATCCGATGCACTGGTTATGGCTGAAATTGTTGCTGATGAGATTGAGGGTTATGGTCCGTTACGGGAtctcg                                                                                                                               >  1:123341/8‑250 (MQ=255)                                                    tAGTGGTTAAACATGCTGACCAAACAACAAAAAATCTTACTTAGAAGTGAAGTATTGAGCAATCTTGATGTTGAAAAGATTGATGAACTTCAAAGTGAAAGAAGCTCATTAGTAAACGAACTGGTTCAGATTGTTAATCGTATTGCCAATAAAAGTAACACTTATTTGACGTCATCCGATGCACTGGTTATGGCTGAAATTGTTGCTGATGAGATTGAGGGTTATGGTCCGTTACGGGATCTCGTGGCTg                                                                                                                         <  1:551862/250‑1 (MQ=255)                                                       tggttAAACATGCTGACCAAACAACAAAAAATCTTACTTAGAAGTGAAGTATTGAGCAATCTTGATGTTGAAAAGATTGATGAACTTCAAAGTGAAAGAAGCTCATTAGTAAACGAACTGGTTCAGATTGTTAATCGTATTGCCAATAAAAGTAACACTTATTTGACGTCATCCGATGCACTGGTTATGGCTGAAATTGTTGCTGATGAGATTg                                                                                                                                                          >  1:330389/5‑214 (MQ=255)                                                        ggttAAACATGCTGACCAAACAACAAAAAATCTTACTTAGAAGTGAAGTATTGAGCAATCTTGATGTTGAAAAGATTGATGAACTTCAAAGTGAAAGAAGCTCATTAGTAAACGAACTGGTTCAGATTGTTAATCGTATTGCCAATAAAAGTAACACTTATTTGACGTCATCCGATGCACTGGTTATGGCTGAAATTGTTGCTGATGAGAt                                                                                                                                                            <  1:71867/208‑1 (MQ=255)                                                        ggttAAACATGCTGACCAAACAACAAAAAATCTTACTTAGAAGTGAAGTATTGAGCAATCTTGATGTTGAAAAGATTGATGAACTTCAAAGTGAAAGAAGCTCATTAGTAAACGAACTGGTTCAGATTGTTAATCGTATTGCCAATAAAAGTAACACTTATTTGACGTCATCCGATGCACTGGTTATGGCTGAAATTGTTGCTGATGAGATTGAGGGt                                                                                                                                                     >  1:612660/4‑218 (MQ=255)                                                        ggttAAACATGCTGACCAAACAACAAAAAATCTTACTTAGAAGTGAAGTATTGAGCAATCTTGATGTTGAAAAGATTGATGAACTTCAAAGTGAAAGAAGCTCATTAGTAAACGAACTGGTTCAGATTGTTAATCGTATTGCCAATAAAAGTAACACTTATTTGACGTCATCCGATGCACTGGTTATGGCTGAAATTGTTGCTGATGAGATTGAGGGTTATgg                                                                                                                                                >  1:879319/4‑223 (MQ=255)                                                        ggttAAACATGCTGACCAAACAACAAAAAATCTTACTTAGAAGTGAAGTATTGAGCAATCTTGATGTTGAAAAGATTGATGAACTTCAAAGTGAAAGAAGCTCATTAGTAAACGAACTGGTTCAGATTGTTAATCGTATTGCCAATAAAAGTAACACTTATTTGACGTCATCCGATGCACTGGTTATGGCTGAAATTGTTGCTGATGAGATTGAGGGTTATGGTCCGTTACGGGATCTCGTGGCTGACGAt                                                                                                                    >  1:42600/4‑251 (MQ=255)                                                        ggttAAACATGCTGACCAAACAACAAAAAATCTTACTTAGAAGTGAAGTATTGAGCAATCTTGATGTTGAAAAGATTGATGAACTTCAAAGTGAAAGAAGCTCATTAGTAAACGAACTGGTTCAGATTGTTAATCGTATTGCCAATAAAAGTAACACTTATTTGACGTCATCCGATGCACTGGTTATGGCTGAAATTGTTGCTGATGAGATTGAGGGTTATGGTCCGTTACGGGATCTCGTGGCTGACGAt                                                                                                                    >  1:367179/4‑251 (MQ=255)                                                         gttAAACATGCTGACCAAACAACAAAAAATCTTACTTAGAAGTGAAGTATTGAGCAATCTTGATGTTGAAAAGATTGATGAACTTCAAAGTGAAAGAAGCTCATTAGTAAACGAACTGGTTCAGATTGTTAATCGTATTGCCAATAAAAGTAACACTTATTTGACGTCATCCGATGCACTGGTTATGGCTGAAATTGTTGCTGATGAGATTGAGGGTTATGGTCCGTTACGGGATCTCGTGGCTGACGAt                                                                                                                    <  1:167461/248‑1 (MQ=255)                                                                aTGCTGACCAAACAACAAAAAATCTTACTTAGAAGTGAAGTATTGAGCAATCTTGATGTTGAAAAGATTGATGAACTTCAAAGTGAAAGAAGCTCATTAGTAAACGAACTGGTTCAGATTGTTAATCGTATTGCCAATAAAAGTAACACTTATTTGACGTCATCCGATGCACTGGTTATGGCTGAAATTGTTGCTGATGAGATTGAGGGTTATGGTCCGTTACGGGATCTCGTGGCTGACGATACTATaaa                                                                                                            <  1:722872/251‑1 (MQ=255)                                                                   cTGACCAAACAACAAAAAATCTTACTTAGAAGTGAAGTATTGAGCAATCTTGATGTTGAAAAGATTGATGAACTTCAAAGTGAAAGAAGCTCATTAGTAAACGAACTGGTTCAGATTGTTAATCGTATTGCCAATAAAAGTAACACTTATTTGACGTCATCCGATGCACTGGTTATGGCTGAAATTGTTGCTGATGAGATTGAGGGTTATGGTCCGTTACGGGATCTCGTGGCTGACGATACTATAAATg                                                                                                          <  1:328925/250‑1 (MQ=255)                                                                   cTGACCAAACAACAAAAAATCTTACTTAGAAGTGAAGTATTGAGCAATCTTGATGTTGAAAAGATTGATGAACTTCAAAGTGAAAGAAGCTCATTAGTAAACGAACTGGTTCAGATTGTTAATCGTATTGCCAATAAAAGTAACACTTATTTGACGTCATCCGATGCACTGGTTATGGCTGAAATTGTTGCTGATGAGATTGAGGGTTATGGTCCGTTACGGGATCTCGTGGCTGACGATACTATAAATg                                                                                                          <  1:764746/250‑1 (MQ=255)                                                                   cTGACCAAACAACAAAAAATCTTACTTAGAAGTGAAGTATTGAGCAATCTTGATGTTGAAAAGATTGATGAACTTCAAAGTGAAAGAAGCTCAATAGTAAACGAACTGGTTCAGATTGTTAATCGTATTGCCAATAAAAGTAACACTTATTTGACGTCATCCGATGCACTGGTTATGGCTGAAATTGTTGCTGATGAGATTGAGGGTTATGGTCCGTTACGGGATCTCGTGGCTGACGATACTATAAATg                                                                                                          <  1:184408/250‑1 (MQ=255)                                                                      aCCAAACAACAAAAAATCTTACTTAGAAGTGAAGTATTGAGCAATCTTGATGTTGAAAAGATTGATGAACTTCAAAGTGAAAGAAGCTCATTAGTAAACGAACTGGTTCAGATTGTTAATCGTATTGCCAATAAAAGTAACACTTATTTGACGTCATCCGATGCACTGGTTATGGCTGAAATTGTTGCTGATGAGATTGAGGGTTATGGTCCGTTACGGGATCTCGTGGCTGACGATACTATAAATGatat                                                                                                      <  1:8576/251‑1 (MQ=255)                                                                             aacaaAAAATCTTACTTAGAAGTGAAGTATTGAGCAATCTTGATGTTGAAAAGATTGATGAACTTCAAAGTGAAAGAAGCTCATTAGTAAACGAACTGGTTCAGATTGTTAATCGTATTGCCAATAAAAGTAACACTTATTTGACGTCATCCGATGCACTGGTTATGGCTGAAATTGTTGCTGATGAGATTGAGGGTTATGGTCCGTTACGGGATCTCGTGGCTGACGATACTATAAATGATATCTTAGtt                                                                                               <  1:536723/251‑1 (MQ=255)                                                                                           cttaGAAGTGAAGTATTGAGCAATCTTGATGTTGAAAAGATTGATGAACTTCAAAGTGAAAGAAGCTCATTAGTAAACGAACTGGTTCAGATTGTTAATCGTATTGCCAATAAAAGTAACACTTATTTGACGTCATCCGATGCACTGGTTATGGCTGAAATTGTTGCTGATGAGATTGAGGGTTATGGTCCGTTACGGGATCTCGTGGCTGACGATACTATAAATGATATCTTAGTTAATGGCCCTAATGa                                                                                 <  1:840605/251‑1 (MQ=255)                                                                                                  gtgaagtATTGAGCAATCTTGATGTTGAAAAGATTGATGAACTTCAAAGTGAAAGAAGCTCATTAGTAAACGAACTGGTTCAGATTGTTAATCGTATTGCCAATAAAAGTAACACTTATTTGACGTCATCCGATGCACTGGTTATGGCTGAAATTGTTGCTGATGAGATTGAGGGTTATGGTCCGTTACGGGATCTCGTGGCTGACGATACTATAAATGATATCTTAGTTAATGGCCCTAATGa                                                                                 <  1:152822/244‑1 (MQ=255)                                                                                                       gtattgagcaaTCTTGATGTTGAAAAGATTGATGAACTTCAAAGTGAAAGAAGCTCATTAGTAAACGAACTGGTTCAGATTGTTAATCGTATTGCCAATAAAAGTAACACTTATTTGACGTCATCCGATGCACTGGTTATGGCTGAAATTGTTGCTGATGAGATTGAGGGTTATGGTCCGTTACGGGATCTCGTGGCTGACGATACTATAAATGATATCTTAGTTAATGGCCCTAATGATATCTGGGTTGa                                                                     >  1:590920/10‑251 (MQ=255)                                                                                                       gtattgagcaaTCTTGATGTTGAAAAGATTGATGAACTTCAAAGTGAAAGAAGCTCATTAGTAAACGAACTGGTTCAGATTGTTAATCGTATTGCCAATAAAAGTAACACTTATTTGACGTCATCCGATGCACTGGTTATGGCTGAAATTGTTGCTGATGAGATTGAGGGTTATGGTCCGTTACGGGATCTCGTGGCTGACGATACTATAAATGATATCTTAGTTAATGGCCCTAATGATATCTGGGTTGa                                                                     >  1:676058/10‑251 (MQ=255)                                                                                                         attgagcaaTCTTGATGTTGAAAAGATTGATGAACTTCAAAGTGAAAGAAGCTCATTAGTAAACGAACTGGTTCAGATTGTTAATCGTATTGCCAATAAAAGTAACACTTATTTGACGTCATCCGATGCACTGGTTATGGCTGAAATTGTTGCTGATGAGATTGAGGGTTATGGTCCGTTACGGGATCTCGTGGCTGGCGATACTATAAATGATATCTTAGTTAATGGCCCTAATGATATCTGGGTTGAAc                                                                   <  1:36229/244‑1 (MQ=255)                                                                                                         attgagcaaTCTTGATGTTGAAAAGATTGATGAACTTCAAAGTGAAAGAAGCTCATTAGTAAACGAACTGGTTCAGATTGTTAATCGTATTGCCAATAAAAGTAACACTTATTTGACGTCATCCGATGCACTGGTTATGGCTGAAATTGTTGCTGATGAGATTGAGGGTTATGGTCCGTTACGGGATCTCGTGGCTGACGATACTATAAATGATATCTTAGTTAATGGCCCTAATGATATCTGGGTTGAAc                                                                   <  1:860441/244‑1 (MQ=255)                                                                                                         attgagcaaTCTTGATGTTGAAAAGATTGATGAACTTCAAAGTGAAAGAAGCTCATTAGTAAACGAACTGGTTCAGATTGTTAATCGTATTGCCAATAAAAGTAACACTTATTTGACGTCATCCGATGCACTGGTTATGGCTGAAATTGTTGCTGATGAGATTGAGGGTTATGGTCCGTTACGGGATCTCGTGGCTGACGATACTATAAATGATATCTTAGTTAATGGCCCTAATGATATCTGGGTTGAAc                                                                   <  1:177570/244‑1 (MQ=255)                                                                                                            gagcaaTCTTGATGTTGAAAAGATTGATGAACTTCAAAGTGAAAGAAGCTCATTAGTAAACGAACTGGTTCAGATTGTTAATCGTATTGCCAATAAAAGTAACACTTATTTGACGTCATCCGATGCACTGGTTATGGCTGAAATTGTTGCTGATGAGATTGAGGGTTATGGTCCGTTACGGGATCTCGTGGCTGACGATACTATAAATGATATCTTAGTTAATGGCCCTAATGATATCTGGGt                                                                        >  1:188727/5‑243 (MQ=255)                                                                                                                 aTCTTGATGTTGAAAAGATTGATGAACTTCAAAGTGAAAGAAGCTCATTAGTAAACGAACTGGTTCAGATTGTTAATCGTATTGCCAATAAAAGTAACACTTATTTGACGTCATCCGATGCACTGGTTATGGCTGAAATTGTTGCTGATGAGATTGAGGGTTATGGTCCGTTACggga                                                                                                                                    <  1:198914/178‑4 (MQ=255)                                                                                                                 aTCTTGATGTTGAAAAGATTGATGAACTTCAAAGTGAAAGAAGCTCATTAGTAAACGAACTGGTTCAGATTGTTAATCGTATTGCCAATAAAAGTAACACTTATTTGACGTCATCCGATGCACTGGTTATGGCTGAAATTGTTGCTGATGAGATTGAGGGTTATGGTCCGTTACGGGATCTCGTGGCTGACGATACTATAAATg                                                                                                          <  1:221242/204‑1 (MQ=255)                                                                                                                 aTCTTGATGTTGAAAAGATTGATGAACTTCAAAGTGAAAGAAGCTCATTAGTAAACGAACTGGTTCAGATTGTTAATCGTATTGCCAATAAAAGTAACACTTATTTGACGTCATCCGATGCACTGGTTATGGCTGAAATTGTTGCTGATGAGATTGAGGGTTATGGTCCGTTACGGGATCTCGTGGCTGACGATACTATAAATGat                                                                                                        <  1:496865/206‑1 (MQ=255)                                                                                                                 aTCTTGATGTTGAAAAGATTGATGAACTTCAAAGTGAAAGAAGCTCATTAGTAAACGAACTGGTTCAGATTGTTAATCGTATTGCCAATAAAAGTAACACTTATTTGACGTCATCCGATGCACTGGTTATGGCTGAAATTGTTGCTGATGAGATTGAGGGTTATGGTCCGTTACGGGATCTCGTGGCTGACGATACTATAAATGat                                                                                                        >  1:176770/1‑206 (MQ=255)                                                                                                                 aTCTTGATGTTGAAAAGATTGATGAACTTCAAAGTGAAAGAAGCTCATTAGTAAACGAACTGGTTCAGATTGTTAATCGTATTGCCAATAAAAGTAACACTTATTTGACGTCATCCGATGCACTGGTTATGGCTGAAATTGTTGCTGATGAGATTGAGGGTTATGGTCCGTTACGGGATCTCGTGGCTGACGATACTATAAATGATATCTTAGTTAATGGCCCTAATg                                                                                  >  1:816649/1‑228 (MQ=255)                                                                                                                     tGATGTTGAAAAGATTGATGAACTTCAAAGTGAAAGAAGCTCATTAGTAAACGAACTGGTTCAGATTGTTAATCGTATTGCCAATAAAAGTAACACTTATTTGACGTCATCCGATGCACTGGTTATGGCTGAAATTGTTGCTGATGAGATTGAGGGTTATGGTCCGTTACGGGATCTCGTGGCTGACGATACTATAAATGATATCTTAGTTAATGGCCCTAATGATATCTGGGTTGAACGTGCCGGTATTc                                                       <  1:211039/251‑1 (MQ=255)                                                                                                                                                aaGTGAAAGAAGCTCATTAGTAAACGAACTGGTTCAGATTGTTAATCGTATTGCCAATAAAAGTAACACTTATTTGACGTCATCCGATGCACTGGTTATGGCTGAAATTGTTGCTGATGAGATTGAGGGTTATGGTCCGTTACGGGATCTCGTGGCTGACGATACTATAAATGATATCTTAGtt                                                                                               >  1:308875/1‑184 (MQ=255)                                                                                                                                                  gTGAAAGAAGCTCATTAGTAAACGAACTGGTTCAGATTGTTAATCGTATTGCCAATAAAAGTAACACTTATTTGACGTCATCCGATGCACTGGTTATGGCTGAAATTGTTGCTGATGAGATTGAGGGTTATGGTCCGTTACGGGATCTCGTGGCTGACGATACTATAAATGATATCTTAGTTAATGGCCCTAATGATATc                                                                             >  1:264397/1‑200 (MQ=255)                                                                                                                                                  gTGAAAGAAGCTCATTAGTAAACGAACTGGTTCAGATTGTTAATCGTATTGCCAATAAAAGTAACACTTATTTGACGTCATCCGATGCACTGGTTATGGCTGAAATTGTTGCTGATGAGATTGAGGGTTATGGTCCGTTACGGGATCTCGTGGCTGACGATACTATAAATGATATCTTAGTTAATGGCCCTAATGATATCTGGGTTGAACGTGCCGGTATTc                                                       <  1:23339/222‑1 (MQ=255)                                                                                                                                                           gctcattagTAAACGAACTGGTTCAGATTGTTAATCGTATTGCCAATAAAAGTAACACTTATTTGACGTCATCCGATGCACTGGTTATGGCTGAAATTGTTGCTGATGAGATTGAGGGTTATGGTCCGTTACGGGATCTCGTGGCTGACGATACTATAAATGATATCTTAGTTAATGGCCCTAATGATATCTg                                                                           <  1:645237/185‑1 (MQ=255)                                                                                                                                                           gctcattagTAAACGAACTGGTTCAGATTGTTAATCGTATTGCCAATAAAAGTAACACTTATTTGACGTCATCCGATGCACTGGTTATGGCTGAAATTGTTGCTGATGAGATTGAGGGTTATGGTCCGTTACGGGATCTCGTGGCTGACGATACTATAAATGATATCTTAGTTAATGGCCCTAATGATATCTGGGTTGAACGTGCCGGTATTCTGGAAAAAACTGATAAAGAATTCGTTAGTAATGAAc                   >  1:378392/9‑249 (MQ=255)                                                                                                                                                                ttagTAAACGAACTGGTTCAGATTGTTAATCGTATTGCCAATAAAAGTAACACTTATTTGACGTCATCCGATGCACTGGTTATGGCTGAAATTGTTGCTGATGAGATTGAGGGTTATGGTCCGTTACGGGATCTCGTGGCTGACGATACTATAAATGATATCTTAGTTAATGGCCCTAATGATATCTGGGTTGAACGTGCCGGt                                                           <  1:473317/201‑1 (MQ=255)                                                                                                                                                                ttagTAAACGAACTGGTTCAGATTGTTAATCGTATTGCCAATAAAAGTAACACTTATTTGACGTCATCCGATGCACTGGTTATGGCTGAAATTGTTGCTGATGAGATTGAGGGTTATGGTCCGTTACGGGATCTCGTGGCTGACGATACTATAAATGATATCTTAGTTAATGGCCCTAATGATATCTGGGTTGAACGTGCCGGTATTCTGGAAAAAACTGATAAAGAATTCGTTAGTAATGAACAATTaa             >  1:719901/4‑250 (MQ=255)                                                                                                                                                                    tAAACGAACTGGTTCAGATTGTTAATCGTATTGCCAATAAAAGTAACACTTATTTGACGTCATCCGATGCACTGGTTATGGCTGAAGTTGTTGCTGATGAGATTGAGGGTTATGGTCCGTTACAGGATCTCGTGGCTGACGATACTATAAATGATATCTTAGTTAATGGCCCTAATGATATCTGGGTTGAACGTGCCGGTATTCTGGAAAAAACTGATAAAGAATTCGTTAGTAATGAACAATTAACTGa         <  1:91403/250‑1 (MQ=255)                                                                                                                                                                    tAAACGAACTGGTTCAGATTGTTAATCGTATTGCCAATAAAAGTAACACTTATTTGACGTCATCCGATGCACTGGTTATGGCTGAAATTGTTGCTGATGAGATTGAGGGTTATGGTCCGTTACGGGATCTCGTGGCTGACGATACTATAAATGATATCTTAGTTAATGGCCCTAATGATATCTGGGTTGAACGTGCCGGTATTCTGGAAAAAACTGATAAAGAATTCGTTAGTAATGAACAATTAACTGat        >  1:203429/1‑251 (MQ=255)                                                                                                                                                                          aaCTGGTTCAGATTGTTAATCGTATTGCCAATAAAAGTAACACTTATTTGACGTCATCCGATGCACTGGTTATGGCTGAAATTGTTGCTGATGAGATTGAGGGTTATGGTCCGTTACGGGATCTCGTGGCTGACGATACTATAAATGATATCTTAGTTAATGGCCCTAATGATATCTGGGTTGAACGTGCCGGTATTCTGGAAAAAACTGATAAAGAATTCGTTAGTAATGAACAATTAACTGATATTGCa  >  1:865880/1‑251 (MQ=255)                                                                                                                                                                                    gATTGTTAATCGTATTGCCAATAAAAGTAACACTTATTTGACGTCATCCGATGCACTGGTTATGGCTGAAATTGTTGCTGATGAGATTGAGGGTTATGGTCCGTTACGGGATCTCGTGGCTGACGATACTATAAATGATATCTTAGTTAATGGc                                                                                         <  1:547307/154‑1 (MQ=255)                                                                                                                                                                                    gATTGTTAATCGTATTGCCAATAAAAGTAACACTTATTTGACGTCATCCGATGCACTGGTTATGGCTGAAATTGTTGCTGATGAGATTGAGGGTTATGGTCCGTTACGGGATCTCGTGGCTGACGATACTATAAATGATATCTTAGTTAATGGCCCTAATGATATCTGGGTTGAACGTGCCGGTATTCTGGAAAAAACTGATAAAGAATTCGTTAGTAATGAACAATTAACt           <  1:338096/232‑1 (MQ=255)                                                                                                                                                                                                               |||||                                                                                                                                                                                                                  AAAAACAGAGTTTATTTTCTTCATTTTTTTCAAAACTCATTAAACAGTAAATAGTGGCTAAACATGCTGACCAAACAACAAAAAATCTTACTTAGAAGTGAAGTACTAAGTAATCTTGATATTGAAAAGATTGATGAACTTCAAAGTGAAAGAAGCTCACTGGTAAACGAACTGGTTCAGATTGTTAACCGTGTTGCCAATAAAAGCGGTGCTTATTTGACCTCAGCCGATACTCTAGTTATGGCTGAAATTGTTGCTGATGAGATTGAGGGTTATGGACCATTACGAGATCTCATGGCTGACGATACTATAAATGATATTTTAGTTAATGGCCCTAATGACATCTGGGTTGAACGTGCTGGTATTCTGGAAAAAACTGATAAAGAATTTGTTAGTAATGAACAATTAACTGATATTGCA  >  1/1482989‑1483408

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