Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	1,047,106	0	G	A	38.0%	102.0 / 28.5	71	S61F (TCC→TTC)	unknown	Sulfate transporter, CysZ‑type
Reads supporting (aligned to +/- strand):  ref base G (33/11);  new base A (0/27);  total (33/38)
Fisher's exact test for biased strand distribution p-value = 4.65e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GCTGTAAACCAGCTTTTGCCATTCCCGTACCAGCATCCGCGGCACATCTTTCATAAAATCCCACAGGGACATATCCTGCAAATTTTCTCCGGTCAGCATTTTTTCCACTTTTTCCGCCAGCAAGCC‑ATTAAAAGGCGCGGCGATAAAGCCCGACAGCGTGGTAAAGACGAAGTAAAACAGAACGAGAATCATAAGAATGGAAAGCAGCAGCAAAA‑CGCCGCT‑CAACCAA‑CTCAGCCAATCGGGAATGTAGGAAATAATGCCGTC‑AATATAACCGCCGATTTGGTTCACAAAAAGCCAAAACAGACCGCCTAATAACAGCATATTTAACAGGATCGGCATGATCACAAAACGGCGCAGCCCTTTCTGTAAAATCAAATGCCAACCCATCACAAAATATTCAAAACCGGATTTCATTTCCGCGTGTTTCAACATAGCGTTATCCCTAAGAAAAAAGAGGTTCAGGATACAAAGTGCGGTGAAAAAAATCAAACTTTGTC  >  1/1046856‑1047352                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                        gCTGTAAACCAGCTTTTGCCATTCCCGTACCAGCATCCGCGGCACATCTTTCATAAAATCCCACAGGGACATATCCTGCAAATTTTCTCCGGTCAGCATTTTTTCCACTTTTTCCGCCAGCAAGCC‑ATTAAAAGGCGCGGCGATAAAGCCCGACAGCGTGGTAAAGACGAAGTAAAACAGAACGAGAATCATAAGAATGGAAAGCAGCAGCAAAA‑CGCCGCT‑CAACCAA‑CTCAGCCAATCGGGAATGTAgg                                                                                                                                                                                                                                                         >  1:455639/1‑251 (MQ=255)     tAAACCAGCTTTTGCCATTCCCGTACCAGCATCCGCGGCACATCTTTCATAAAATCCCACAGGGACATATCCTGCAAATTTTCTCCGGTCAGCATTTTTTCCACTTTTTCCGCCAGCAAGCC‑ATTAAAAGGCGCGGCGATAAAGCCCGACAGCGTGGTAAAGACGAAGTAAAACAGAACGAGAATCATAAGAATGGAAAGCAGCAGCAAAA‑CGCCGCT‑CAACCAA‑CTCAGCCAATCGGGAATGTAGGAaat                                                                                                                                                                                                                                                     <  1:476619/251‑1 (MQ=255)                      ttCCCGTACCAGCATCCGCGGCACATCTTTCATAAAATCCCACAGGGACATATCCTGCAAATTTTCTCCGGTCAGCATTTTTTCCACTTTTTCCGCCAGCAAGCC‑ATTAAAAGGCGCGGCGATAAAGCCCGACAGCGTGGTAAAGACGAAGTAAAACAGAACGAGAATCATAAGAATGGAAAGCAGCAGCAAAA‑CGCCGCT‑CAACCAA‑CTCAGCCTATCGGGAATGTAGGAAATAATGCCGTC‑AATATAAc                                                                                                                                                                                                                                   <  1:551997/251‑1 (MQ=255)                      ttCCCGTACCAGCATCCGCGGCACATCTTTCATAAAATCCCACAGGGACATATCCTGCAAATTTTCTCCGGTCAGCATTTTTTCCACTTTTTCCGCCAGCAAGCC‑ATTAAAAGGCGCGGCGATAAAGCCCGACAGCGTGGTAAAGACGAAGTAAAACAGAACGAGAATCATAAGAATGGAAAGCAGCAGCAAAA‑CGCCGCT‑CAACCAA‑CTCAGCCAATCGGGAATGTAGGAAATAATGCCGTC‑AATATAAc                                                                                                                                                                                                                                   <  1:294421/251‑1 (MQ=255)                       tCCCGTACCAGCATCCGCGGCACCTCTTTCATAAAATCCCACAGGGACATATCCTGCAAATTTTCTCCGGTCAGCATTTTTTCCACTTTTTCCGCCAGCAAGCC‑ATTAAAAGGCGCGGCGATAAAGCCCGACAGCGTGGTAAAGACGAAGTAAAACAGAACGAGAATCATAAGAATGGAAAGCAGCAGCAAAA‑CGCCGCT‑CAACCAA‑CTCAGCCAATCGGGAATGTAGGAAATAATGCCGTC‑AATATAAc                                                                                                                                                                                                                                   >  1:636802/1‑250 (MQ=255)                       tCCCGTACCAGCATCCGCGGCACATCTTTCATAAAATCCCACAGGGACATATCCTGCAAATTTTCTCCGGTCAGCATTTTTTCCACTTTTTCCGCCAGCAAGCC‑ATTAAAAGGCGCGGCGATAAAGCCCGACAGCGTGGTAAAGACGAAGTAAAACAGAACGAGAATCATAAGAATGGAAAGCAGCAGCAAAA‑CGCCGCT‑CAACCAA‑CTCAGCCAATCGGGAATGTAGGAAATAATGCCGTC‑AATATAAc                                                                                                                                                                                                                                   >  1:79384/1‑250 (MQ=255)                        cccGTACCAGCATCCGCGGCACATCTTTCATAAAATCCCACAGGGACATATCCTGCAAATTTTCTCCGGTCAGCATTTTTTCCACTTTTTCCGCCAGCAAGCC‑ATTAAAAGGCGCGGCGATAAAGCCCGACAGCGTGGTAAAGACGAAGTAAAACAGAACGAGAATCATAAGAATGGAAAGCAGCAGCAAAA‑CGCCGCT‑CAACCAA‑CTCAGCCAATCGGGAATGTAGGAAATAATGCCGTC‑AATATAAccg                                                                                                                                                                                                                                 >  1:74002/1‑251 (MQ=255)                           gTACCAGCATCCGCGGCACATCTTTCATAAAATCCCACAGGGACATATCCTGCAAATTTTCTCCGGTCAGCATTTTTTCCACTTTTTCCGCCAGCAAGCC‑ATTAAAAGGCGCGGCGATAAAGCCCGACAGCGTGGTAAAGACGAAGTAAAACAGAACGAGAATCATAAGAATGGAAAGCAGCAGCAAAA‑CGCCGCT‑CAACCAA‑CTCAGCCAATCGGGAATGTAGGAAATAATGCCGTC‑AATATAAccgccg                                                                                                                                                                                                                              >  1:590579/1‑251 (MQ=255)                                   aTCCGCGGCACATCTTTCATAAAATCCCACAGGGACATATCCTGCAAATTTTCTCCGGTCAGCATTTTTTCCACTTTTTCCGCCAGCAAGCC‑ATTAAAAGGCGCGGCGATAAAGCCCGACAGCGTGGTAAAGACGAAGTAAAACAGAACGAGAATCATAAGAATGGAAAGCAGCAGCAAAA‑CGCCGCT‑CAACCAA‑CTCAGCCAATCGGGAATGTAGGAAATAATGCCGTC‑AATATAACCGCCGATTTGGtt                                                                                                                                                                                                                      >  1:294713/1‑251 (MQ=255)                                     ccGCGGCACATCTTTCATAAAATCCCACAGGGACATATCCTGCAAATTTTCTCCGGTCAGCATTTTTTCCACTTTTTCCGCCAGCAAGCC‑ATTAAAAGGCGCGGCGATAAAGCCCGACAGCGTGGTAAAGACGAAGTAAAACAGAACGAGAATCATAAGAATGGAAAGCAGCAGCAAAA‑CGCCGCT‑CAACCAA‑CTCAGCCAATCGGGAATGTAGGAAATAATGCCGTC‑AATATAACCGCCGATTTGGTTca                                                                                                                                                                                                                    >  1:457592/1‑251 (MQ=255)                                                              cacaGGGACATATCCTGCAAATTTTCTCCGGTCAGCATTTTTTCCACTTTTTCCGCCAGCAAGCC‑ATTAAAAGGCGCGGCGATAAAGCCCGACAGCGTGGTAAAGACGAAGTAAAACAGAACGAGAATCATAAGAATGGAAAGCAGCAGCAAAA‑CGCCGCT‑CAACCAA‑CTCAGCCAATCGGGAATGTAGGAAATAATGCCGTC‑AATATAAc                                                                                                                                                                                                                                   <  1:2167/211‑1 (MQ=255)                                                                caGGGACATATCCTGCAAATTTTCTCCGGTCAGCATTTTTTCCACTTTTTCCGCCAGCAAGCC‑ATTAAAAGGCGCGGCGATAAAGCCCGACAGCGTGGTAAAGACGAAGTAAAACAGAACGAGAATCATAAGAATGGAAAGCAGCAGCAAAA‑CGCCGCT‑CAACCAA‑CTCAGCCAATCGGGAATGTAGGAAATAATGCCGTC‑AATATAACCGCCGATTTGGTTCACAAAAAGCCAAAACAGACCGCCtaata                                                                                                                                                                                         >  1:178895/1‑251 (MQ=255)                                                                  gggACATATCCTGCAAATTTTCTCCGGTCAGCATTTTTTCCACTTTTTCCGCCAGCAAGCC‑ATTAAAAGGCGCGGCGATAAAGCCCGACAGCGTGGTAAAGACGAAGTAAAACAGAACGAGAATCATAAGAATGGAAAGCAGCAGCAAAA‑CGCCGCT‑CAACCAA‑CTCAGCCAATCGGGAATGTAGGAAATAATGCCGTC‑AATATAACCGCCGATTTGGTTCACAAAAAGCCAAAACAGACCGCCtaataa                                                                                                                                                                                        >  1:483094/1‑250 (MQ=255)                                                                   ggACATATCCTGCAAATTTTCTCCGGTCAGCATTTTTTCCACTTTTTCCGCCGGAAAGCC‑ATTAAAAGGCGCGGCGATAAAGCCCGAAAGCGGGGAAAAGACGAAGAAAAACAAAACGAGAATCATAAGAATGGAAAGCAGCAGCAAAA‑CGCCGCT‑CAACCAA‑CTAAGCCAACCGGAATTGTAGGaa                                                                                                                                                                                                                                                       >  1:434929/1‑187 (MQ=255)                                                                   ggACATATCCTGCAAATTTTCTCCGGACAGCATTTTTTCCACTTTTTCCGCCAGCAAGCC‑ATTAAAAGGCGCGGCGATAAAGCCCGACAGCGTGGTAAAGACGAAGTAAAACAGAACGAGAATCATAAGAATGGAAAGCAGCAGCAAAA‑CGCCGCT‑CAACCAA‑CTCAGCCAATCGGGAATGTAGGAAATAATGCCGTC‑AATATAACCGCCGATTTGGTTCACAAAAAGCCAAAACAGACCGCCTAATAACa                                                                                                                                                                                      >  1:56712/1‑251 (MQ=255)                                                                   ggACATATCCTACAAATTTTCCCCGGTCAGCATTTTTTCCACTTTTCCCGCCAGCAAGCC‑ATAAAAAGGCGCGGCGATAAAGCCCGACAGCGTGGTAAAGACGAAGTAAAACAGAACGAGAATCATAAGAATGGAAAGCAGCAGAAAAC‑CGCCGCT‑CAACCAA‑CTCAGCCAATCGGGAATGTAGAAAATAATGCCGTC‑AATATAACCGCCGATTTGGTTCACAAAAAGCCAAAACAGACCGCCTAATAAc                                                                                                                                                                                       <  1:218411/250‑1 (MQ=255)                                                                   gaacatacccggaaaaTTTTCTCCGGCCAGCATTTTTCCCACTTTTCCCGCAAGCAAGCA‑ATAAAAAGGCGCGGCGAAAAACCCCGACAGCGTGGTAAAGACGAAGTAAAACAGAACAAGAATCAAAAGAAGGAAAAGCAGCAACAAAC‑CGCCGCT‑CAACCAA‑CCCACCCAATCGGGAATGTAGAAAATAATGCCGTA‑AATATAACCGCCGATTTGGTTCACAAAAAGCCAAAACAGACCGCCTAATAAc                                                                                                                                                                                       <  1:234376/237‑1 (MQ=255)                                                                   gaaCATATCCTGCAAATTTTCTCCGGTCAGCATTTTTCCCATTTTTCCCGCCAGCAAGCC‑ATTAAAAGGCGCGGCAAAAAAGCCCGACAGCGGGGTAAAGACAAAGTAAAACAGAACGAAAATCATAAGAATGGAAAGCAGCAGCAAAC‑CGCCGCT‑AAACCAA‑CTCAGCAAATGGGAAATGTAGGAAATAATGCCGTA‑AATATAACCGCCGATTTGGTTCACAAAAAGCCAAAACAGACCGCCTAATAAc                                                                                                                                                                                       <  1:627352/248‑1 (MQ=255)                                                                    aaCATACCCTGAAAATTTTCCCGGGTCAGCTTTTTTTCCCCTTTTCCCGCCAGCAAGCC‑ATAAAAAGGGGCGGCAAAAAAGCCCGACAGCGTGGAAAAGACGAAGAAAAACAGAACGAGAATCATAAGAAGGAAAAGCAGCAGAAAAA‑CGCCGCT‑AAACCAA‑CTCAGCCAATCGGAAATGTAGAAAATAATGCCGTC‑AATAAAACCGCCGATTTGGTTCACAAAAAGCCAAAACAGACCGCCTAATAACAg                                                                                                                                                                                     <  1:100265/250‑1 (MQ=255)                                                                     aCATATCCTGAAAATTTTCCCCGGTCAGCATTTTTCCCACTTTTCCCGCCAAAAAGCC‑ATAAAAAGGCGCGGCAAAAAAGCCCGACAGCGTGGTAAACACGAAGAAAAACAGAACGAAAATCAAAAGAATGGAAAGCAGCAGCAAAA‑CGCCGCT‑CACCCAA‑CTCAGCCAATCGGGAATGTAGGAAATAATGCCGTC‑AATATAACCGCCGATTTGGTTCACAAAAAGCCAAAACAGACCGCCTAATAACAg                                                                                                                                                                                     <  1:288491/250‑1 (MQ=255)                                                                       atacccTGCAAATTTCCCCCGGTAAGCATTTTTCCCACTTTTCCCCCCAGC‑AGCCAATAAAAAGGCGCGGAGAAAAAGCCCAACAGCTTGGAAAAAAGGAAATAAAACAAAACTAAAATCATAAGAGAGGAAACAAGCAACAAAA‑CCCCCCA‑CACCCAC‑CTCACCCAATCGGTAATGTAGAAAAAAATGCCGAC‑AATATAACCGCCGATTTGGTTCACAAAAAGCCAAAACAGACCGCCTAATAACAGCat                                                                                                                                                                                  <  1:446170/247‑1 (MQ=255)                                                                              caaaaTTTCCCCCGGTAATCATTTTTCCCACTTTTCCCGCCAGCAAGAC‑ATAAAAAGGCGCGGCGATAAACCCCGAGAGCGGGATAAAGAGGAAAAAAAACAAAACGAAAATTATAAAAATGGAAAGCACCAGCAAAG‑CACC‑CTCCACCCAA‑CTCAACAAATCGGTAATTAAGAAAAAAATGCCGTC‑AATATAACCGCCGATTTGGTTCACAAAAAGCCAAAACAGACCGCCTAATAACAGCATATTTAAc                                                                                                                                                                           <  1:137201/249‑1 (MQ=255)                                                                              aaaaaTTTCCCCCGGACAGCATTTTTCCCCCTTTTCCCGCCAGAAAGCA‑ATAAAAAGGCGGGGCAATAAACCCGGACGGCGTGGTAAAAACAAAGTAAAACAAAACGAGAATCAAAAGAAGGGAAAGCAGCAGCAAAA‑CGCCGCT‑CAACCAA‑CTCACCCATCCGGAAAGGAAGAAAAAAATGCCTTC‑AATATAACCGCCGATTTGGTTCACAAAAAGCCAAAACAGACCGCCTAATAACAGCATATTTAAc                                                                                                                                                                           <  1:269064/249‑1 (MQ=255)                                                                                 aaTTTTCTCCGGTCAGCATTTTTTCCACTTTTTCCGCCAGCAAGCC‑ATTAAAAGGCGCGGCGATAAAGCCCGACAGCGTGGTAAAGACGAAGTAAAACAGAACGAGAATCATAAGAATGGAAAGCAGCAGCAAAA‑CGCCGCT‑CAACCAA‑CTCAGCCAATCGGGAATGTAGGAAATAATGCCGTC‑AATATAACCGCCGATTTGGCTCACAAAAAGCCAAAACAGACCGCCTAATAACAGCATATTTAACAgg                                                                                                                                                                        >  1:590924/1‑251 (MQ=255)                                                                                 aaTTTTCCCCGGTCAGCTTTTTTCCCCCTTTTCCCGCCAGCAAGCA‑ATAAAAAGGCGCGGCAAAAAAGCCCGACGGCGTGGAAAAGAAGAAATAAAAAAGAACGAGAATAAAAAGAAGGGAAAGCAGCAGCAAAA‑CGCCGCT‑CACCCAC‑CCCAGCCAAGCGGGAATGTAGGAAATAATGCCGTA‑AATATAACCGCCGATTTGGTTCACAAAAAGCCAAAACAGACCGCCTAATAACAGCATATTTAACAgg                                                                                                                                                                        <  1:177552/251‑1 (MQ=255)                                                                                        cccGGTCAGCTTTTTTCCCACTTTTCCCGCCAGCAAGCC‑ATAAAAAGGCGCGGCAATAAAGCCCGACAGCGTGGTAAAGACGAAGAAAAACAGAACGAGAATCATAAGAATGGAAAGCAGCAGCAAAA‑CGCCGCT‑CAACCAC‑CTCAGCCAATCGGGAATGTAGAAAATAATGCCGTA‑AATATAACCGCCGATTTGGTTCACAAAAAGCCAAAACAGACCGCCTAATAACAGCATATTTAACAGGATCGGCa                                                                                                                                                                 <  1:460909/250‑1 (MQ=255)                                                                                                  ttttttcCCCCTTTTCCCGCCAGAAAGCC‑ATTAAAAGGCGCGGCGATAAAGCCCGACAGCGTGGAAAAGACGAAGTAAAACAGAACGAAAATCATAAGAATGGAAAGCAGCAGCAAAA‑CGCCCCT‑CAACCAA‑CCCAGCAAATCGGGAATGTAGGAAATAATGCCGTC‑AATATAACCGCCGATTTGGTTCACAAAAAGCCAAAACAGACCGCCTAATAACAGCATATTTAACAGGATCGGCATGATCACaaa                                                                                                                                                       <  1:10209/250‑1 (MQ=255)                                                                                                    ttttccccttttcccccccaccagtcc‑ATAAAAAGGCGCGTCGAAAAACCCCACCAGCGTGAAAAAGACAAAGAAAAACAGAACGAAAACCATAAAAAGGAAAAACCGCAACAAAA‑CGCCCCT‑CAACCAA‑CCCACCAAAGCGGGAATGAAGAAAATAAGGCCGTA‑AATATAACCGCCGATTTGGTTCACAAAAAGCCAAAACAGACCGCCTAATAACAGCATATTTAACAGGATCGGCATGATCACAAAAc                                                                                                                                                     <  1:537265/226‑1 (MQ=255)                                                                                                    ttttcccccTTTTCCCCCCAGCAAGCC‑ATAAAAAGGCGCGGCGAAAAACCCCGACAGCGTGGAAAAGACAAAGAAAAACAGAACGAGAATCAAAAAAATGAAAAGCAGCAGAAAAC‑CGCCGCT‑AAACCAA‑CTCAGCCAAGGGGGAATGTAGAAAATAATGCCGTA‑AATATAACCGCCGATTTGGTTCACAAAAAGCCAAAACAGACCGCCTAATAACAGCATATTTAACAGGATCGGCATGATCACAAAAc                                                                                                                                                     <  1:505939/243‑1 (MQ=255)                                                                                                                  ccGCAAGCAAGCA‑ATAAAAAGGCGCGGCAAAAAAGCCCGACAGCGTGGAAAAGACAAAGAAAAACAGAACGAAAATAAAAAAAAGGAAAAGCAGCAGCAAAC‑CGCCGCT‑AAACCAA‑CCCAGCCAATCGGGAATGTAGAAAATAATGCCGTC‑AATATAACCGCCGATTTGGTTCACAAAAAGCCAAAACAGACCGCCTAATAACAGCATATTTAACAGGATCGGCATGATCACAAAACGGCGCAGCCCTTTc                                                                                                                                       <  1:341676/251‑1 (MQ=255)                                                                                                                    cccAGCAAGCC‑ATAAAAAGGCGCGGCGATAAAGCCCGACAGCGTGGTAAAGACGAAGTAAAACAGAACGAGAATCATAAGAATGAAAAGCAGCAGCAAAA‑CGCCGCT‑CAACCAA‑CTCAGCCAATCGGGAATGAAGGAAATAATGCCGTA‑AATATAACCGCCGATTTGGTTCACAAAAAGCCAAAACAGACCGCCTAATAACAGCATATTTAACAGGATCGGCATGATCACAAAACGGCGCAGCCCTTTCTg                                                                                                                                     <  1:170965/250‑1 (MQ=255)                                                                                                                      cAGCAAGCC‑ATTAAAAGGCGCGGCGATAAAGCCCGACAGCGTGGTAAAGACGAAGTAAAACAGAACGAGAATCATAAGAATGGAAAGCAGCAGCAAAA‑CGCCGCT‑CAACCAA‑CTCAGCCAATCGGGAATGTAGGAAATAATGCCGTC‑AATATAACCGCCGATTTGGTTCACAAAAAGCCAAAACAGAc                                                                                                                                                                                                  >  1:76971/1‑188 (MQ=255)                                                                                                                                ataaaaaGGCGCGGCAATAAACCCCGACAGCGTGGTAAAGAAGAAGTAAAACAGAACGAGAATCATAAGAAGGAAAAGCAGCAGAAAAA‑CGCCGCT‑CAACCAC‑CTCAGCCAATGGGAAATGTAGAAAATAATGCCGTC‑AATATAACCGCCGATTTGGTTCACAAAAAGCCAAAACAGACCGCCTAATAACAGCATATTTAACAGGATCGGCATGATCACAAAACGGCGCAGCCCTTTc                                                                                                                                       <  1:281111/235‑1 (MQ=255)                                                                                                                                aTTAAAAGGCGCGGCGATAAAGCCCGACAGCGTGGTAAAGACGAAGTAAAACAGAACGAGAATCATAAGAATGGAAAGCAGCAGCAAAA‑CGCCGCT‑CAACCAA‑CTCAGCCAATCGGGAATGTAGGAAATAATGCCGTC‑AATATAACCGCCGATTTGGTTCACAAAAAGCCAAAACAGACCGCCTAATAACAGCATATTTAAc                                                                                                                                                                           >  1:464390/1‑202 (MQ=255)                                                                                                                                aTTAAAAGGCGCGGCGATAAAGCCCGACAGCGTGGTAAAGACGAAGTAAAACAGAACGAGAATCATAAGAATGGAAAGCAGCAGCAAAA‑CGCCGCT‑CAACCAA‑CTCAGCCAATCGGGAATGTAGGAAATAATGCCGTC‑AATATAACCGCCGATTTGGTTCACAAAAAGCCAAAACAGACCGCCTAATAACAGCATATTTAACAGGATCGGCATGATCACAAAACGGCGCAGCCCTTTCTGTAAAATCAAAt                                                                                                                          >  1:385957/1‑251 (MQ=255)                                                                                                                                aTTAAAAGGCGCGGCGATAAAGCCCGACAGCGTGGTAAAGACGAAGTAAAACAGAACGAGAATCATAAGAATGGAAAGCAGCAGCAAAA‑CGCCGCT‑CAACCAA‑CTCAGCCAATCGGGAATGTAGGAAATAATGCCGTC‑AATATAACCGCCGATTTGGTTCACAAAAAGCCAAAACAGACCGCCTAATAACAGCATATTTAACAGGATCGGCATGATCACAAAACGGCGCAGCCCTTTCTGTAAAATCAAAt                                                                                                                          >  1:390634/1‑251 (MQ=255)                                                                                                                                aTTAAAAGGCGCGGCGATAAAGCCCGACAGCGTGGAAAAGACGAAGAAAAACAGAACGAAAATCATAAGAAGGGAAAGCAGCAGAAAAA‑CGCCGCT‑CACCCAA‑CTCAGCCAATCGGGAATGTAGGAAATAATGCCGTC‑AATATAACCGCCGATTTGGTTCACAAAAAGCCAAAACAGACCGCCTAATAAc                                                                                                                                                                                       <  1:203158/190‑1 (MQ=255)                                                                                                                                aTTAAAAGGCGCGGCGATAAAGCCAGACAGCGTGGTAAAGACGAAGTAAAACAGAACGAGAATCATAAGAATGGAAAGCAGCAGCAAAA‑CGCCGCT‑CAACCAA‑CTCAGCCAATCGGGAATGTAGGAAATAATGCCGTC‑AATATAACCGCCGATTTGGTTCACAAAAAGCCAAAACAGACCGCCTAATAACAGCATATTTAACAGGATCGGCATGATCACAAAACGGCGCAGCCCTTTCTGTAAAATCAAAt                                                                                                                          >  1:427923/1‑251 (MQ=255)                                                                                                                                             gaaaaaaaccccGACAGCGTGGAAAAGACGAAGAAAAACAAAACGAGAATAAAAAGAAGGAAAAGCAGAAGAAAAC‑CGCCGCC‑CAACCAA‑CCCAGCCAATGGGGAATGTAGAAAAAAAGGCCGTA‑AATAAAACCGCCGATTTGGTTCACAAAAAGCCAAAACAGACCGCCTAATAACAGCATATTTAACAGGATCGGCATGATcac                                                                                                                                                          <  1:133187/197‑1 (MQ=255)                                                                                                                                                   aaGCCCGACAGCGTGGTAAAGACGAAGTAAAACAGAACGAGAATCATAAGAATGGAAAGCAGCAGCAAAA‑CGCCGCT‑CAACCAA‑CTCAGCCAATCGGGAATGTAGGAAATAATGCCGTC‑AATATAACCGCCGATTTGGTTCACAAAAAGCCAAAACAGACCGCCTAATAACAGCATATTTAACAGGATCGGCATGATCACAAAACGGCGCAGCCCTTTCTGTAAAATCAAATGCCAACCCATCACaaaa                                                                                                         >  1:506847/1‑249 (MQ=255)                                                                                                                                                      cccaaCAGCGTGTAAAAAACAAAGAAAAACAGAACGAAAACAAAAAGAAGGAAAAGCAGCAAAAAAA‑CCCCCCT‑CACCCAA‑CCCCGCCAATCGGAAATGTAGAAAATAAGGCCTTA‑AATAAAACCGCCGATTTGGTTCACAAAAAGCCAAAACAGACCGCCTAATAACAGCATATTTAACAGGATCGGCATGATCACAAAACGGCGCAGCCCTTTCTGTAAAATCAAATGCCAACCCATCACAAAATATTc                                                                                                    <  1:592988/247‑1 (MQ=255)                                                                                                                                                                 ggaaaaGACGAAGAAAAACAGAACAAAAATAAAAAGAAGAGAAAGCAGCACAAAAA‑CGCCGCC‑CAACAAA‑CTCAGCCAATCGGGAATGTAGAAAATAATGCCTTC‑AATAAAACCGCCGATTTGGTTCACAAAAAGCCAAAACAGACCGCCTAATAACAGCATATTTAACAGGATCGGCATGATCACAAAACGGCGCAGCCCTTTCTGTAAAATCAAATGCCAACCCATCACAAAat                                                                                                        <  1:424695/233‑1 (MQ=255)                                                                                                                                                                 ggTAAAGACGAAGTAAAACAGAACGAGAATCATAAGAATGGAAAGCAGCAGCAAAA‑CGCCGCT‑CAACCAA‑CTCAGCCAATCGGGAATGTAGGAAATAATGCCGTC‑AATATAACCGCCGATTTGGTTCACAAAAAGCCAAAACAGACCGCCTAATAACAGCATATTTAACAGGATCGGCATGATCACAAAACGGCGCAGCCCTTTCTGt                                                                                                                                    >  1:346448/1‑208 (MQ=255)                                                                                                                                                                 ggTAAAGACGAAGTAAAACAGAACGAGAATCATAAGAATGGAAAGCAGCAGCAAAA‑CGCCGCT‑CAACCAA‑CTCAGCCAATCGGGAATGTAGGAAATAATGCCGTC‑AATATAACCGCCGATTTGGTTCACAAAAAGCCAAAACAGACCGCCTAATAACAGCATATTTAACAGGATCGGCATGATCACAAAACGGCGCAGCCCTTTCTGTAAAATCAAATGCCAACCCATCACAAAATATTCAAAACCGGatt                                                                                         >  1:93653/1‑251 (MQ=255)                                                                                                                                                                 ggTAAAGACGAAGTAAAACAGAACGAGAATCATAAGAATGGAAAGCAGCAGCAAAA‑CGCCGCT‑CAACCAA‑CTCAGCCAATCGGGAATGTAGGAAATAATGCCGTC‑AATATAACCGCCGATTTGGTTCACAAAAAGCCAAAACAGACCGCCTAATAACAGCATATTTAACAGGATCGGCATGATCACAAAACGGCGCAGCCCTTTCTGTAAAATCAAATGCCAACCCATCACAAAATATTCAAAACCGGatt                                                                                         >  1:182138/1‑251 (MQ=255)                                                                                                                                                                 ggTAAAGACGAAGTAAAACAGAACGAGAATCATAAGAATGGAAAGCAGCAGCAAAA‑CGCCGCT‑CAACCAA‑CTCAGCCAATCGGGAATGTAGGAAATAATGCCGTC‑AATATAACCGCCGATTTGGTTCACAAAAAGCCAAAACAGACCGCCTAATAACAGCATATTTAACAGGATCGGCATGATCACAAAACGGCGCAGCCCTTTCTGTAAAATCAAATGCCAACCCATCACAAAATATTCAAAACCGGatt                                                                                         >  1:207836/1‑251 (MQ=255)                                                                                                                                                                 ggTAAAGACGAAGTAAAACAGAACGAGAATCATAAGAATGGAAAGCAGCAGCAAAA‑CGCCGCT‑CAACCAA‑CTCAGCCAATCGGGAATGTAGGAAATAATGCCGTC‑AATATAACCGCCGATTTGGTTCACAAAAAGCCAAAACAGACCGCCTAATAACAGCATATTTAACAGGATCGGCATGATCACAAAACGGCGCAGCCCTTTCTGTAAAATCAAATGCCAACCCATCACAAAATATTCAAAACCGGatt                                                                                         >  1:33471/1‑251 (MQ=255)                                                                                                                                                                  gTAAAGACGAAGTAAAACAGAACGAGAATCATAAGAATGGAAAGCAGCAGCAAAA‑CGCCGCT‑CAACCAA‑CTCAGCCAATCGGGAATGTAGGAAATAATGCCGTC‑AATATAACCGCCGATTTGGTTCACAAAAAGCCAAAACAGACCGCCTAATAACAGCATATTTAACAGGATCGGCa                                                                                                                                                                 >  1:263461/1‑178 (MQ=255)                                                                                                                                                                      aGACGAAGTAAAACAGAACGAGAATCATAAGAATGGAAAGCAGCAGCAAAA‑CGCCGCT‑CAACCAA‑CTCAGCCAATCGGGAATGTAGGAAATAATGCCGTC‑AATATAACCGCCGATTTGGTTCACAAAAAGCCAAAACAGACCGCCTAATAACAGCat                                                                                                                                                                                  >  1:310800/1‑157 (MQ=255)                                                                                                                                                                        aCGAAGAAAAAAAGAACAAAAATCATAAAAATGAAAAGCAGCAAAAAAC‑CGCCGCC‑CAACCAC‑CCCACCAAAGCGGAAATGTAGAAAATAATGCCGTA‑AATAAAACCGCCGATTTGGTTCACAAAAAGCCAAAACAGACCGCCTAATAACAGCATATTTAACAGGATCGGCATGATCACAAAACGGCGCAGCCCTTTCTGTAAAATCAAATGCCAACCCATCACAAAATATTCAAAACCGGatttcatttc                                                                                  <  1:569921/251‑1 (MQ=255)                                                                                                                                                                           aaGTAAAACAGAACGAGAATCATAAGAATGGAAAGCAGCAGCAAAA‑CGCCGCT‑CAACCAA‑CTCAGCCAATCGGGAATGTAGGAAATAATGCCGTC‑AATATAACCGCCGATTTGGTTCACAAAAAGCCAAAACAGACCGCCTAATAACAGCATATTTAACAGGATCGGCATGATCACAAAACGGCGCAGCCCTTTCTGTAAAATCAAATGCCAACCCATcaca                                                                                                            >  1:46613/1‑222 (MQ=255)                                                                                                                                                                             gTAAAACAGAACGAGAATCATAAGAATGGAAAGCAGCAGCAAAA‑CGCCGCT‑CAACCAA‑CTCAGCCAATCGGGAATGTAGGAAATAATGCCGTC‑AATATAACCGCCGATTTGGTTCACAAAAAGCCAAAACAGACCGCCTAATAACAGCATATTTAACAGGATCGGCATGATCACAAAACGGCGCAGCCCTTTCTGTAAAATCAAATGCCAACCCATCACAAAATATTCAAAAc                                                                                               >  1:211572/1‑233 (MQ=255)                                                                                                                                                                             gTAAAACAGAACGAGAATCATAAGAATGGAAAGCAGCAGCAAAA‑CGCCGCT‑CAACCAA‑CTCAGCCAATCGGGAATGTAGGAAATAATGCCGTC‑AATATAACCGCCGATTTGGTTCACAAAAAGCCAAAACAGACCGCCTAATAACAGCATATTTAACAGGATCGGCATGATCACAAAACGGCGCAGCCCTTTCTGTAAAATCAAATGCCAACCCATCACAAAATATTCAAAACCGGATTTCATTTCCGCgt                                                                             >  1:302325/1‑251 (MQ=255)                                                                                                                                                                                            aaTCATAAGAATGGAAAGCAGCA‑CAAAACCCCCGCT‑CAACCAA‑CCCAGCAAATCGGAAATGTAGAAAAAAATGCCGTA‑AATATAACCGCCGATTTGGTTCACAAAAAGCCAAAACAGACCGCCTAATAACAGCATATTTAACAGGATCGGCATGATCACAAAACGGCGCAGCCCTTTCTGTAAAATCAAATGCCAACCCATCACAAAATATTCAAAACCGGATTTCATTTCCGCGTGTTTCAACATAGCGt                                                              <  1:50135/251‑1 (MQ=255)                                                                                                                                                                                                 aaaaaagggAAAGCAGCAGAAAAC‑CGCCCCT‑CAACAAA‑CTCAGCAAAGCGGAAATGTAGAAAAAAACGCCGTA‑AATATAACCGCCGATTTGGTTCACAAAAAGCCAAAACAGACCGCCTAATAACAGCATATTTAACAGGATCGGCATGATCACAAAACGGCGCAGCCCTTTCTGTAAAATCAAATGCCAACCCATCACAAAATATTCAAAACCGGATTTCATTTCCGCGTGTTTCAACATAGCGTTATcc                                                         <  1:612392/244‑1 (MQ=255)                                                                                                                                                                                                 aaaGAAAGGAAAGCAGCAGCAAAA‑C‑CCCCC‑CAAACAACCCCACCCAATCGGGAATGAAGAAAATAATGCCGTA‑AATATAACCGCCGATTTGGTTCACAAAAAGCCAAAACAGACCGCCTAATAACAGCATATTTAACAGGATCGGCATGATCACAAAACGGCGCAGCCCTTTCTGTAAAATCAAATGCCAACCCATCACAAAATATTCAAAACCGGATTTCATTTCCGCGTGTTTCAACATAGCGTTATcc                                                         <  1:435733/250‑1 (MQ=255)                                                                                                                                                                                                  aagaaggaaaaGCAGCAGAAAAC‑CGCCGCT‑CAACCAC‑CCCAGCCAACCGGGAATGTAGAAAATAATGCCGTA‑AATAAAACCGCCGATTTGGTTCACAAAAAGCCAAAACAGACCGCCTAATAACAGCATATTTAACAGGATCGGCATGATCACAAAACGGCGCAGCCCTTTCTGTAAAATCAAATGCCAACCCATCACAAAATATTCAAAACCGGATTTCATTTCCGCGTGTTTCAACATAGCGTTATccc                                                        <  1:157370/243‑1 (MQ=255)                                                                                                                                                                                                   aagaaggaaaaccagaagaaaacccccgcT‑CAACCAA‑CTCAGCCAAGCGGAAATGTAGAAAATAATGCCGTA‑AATATAACCGCCGATTTGGTTCACAAAAAGCCAAAACAGACCGCCTAATAACAGCATATTTAACAGGATCGGCATGATCACAAAACGGCGCAGCCCTTTCTGTAAAATCAAATGCCAACCCATCACAAAATATTCAAAACCGGATTTCATTTCCGCGTGTTTCAACATAGCGTTATccc                                                        <  1:161102/226‑1 (MQ=255)                                                                                                                                                                                                         gAAAGCAGCAGCAAAA‑CGCCGCT‑CAACCAA‑CTCAGCCAATCGGGAATGTAGGAAATAATGCCGTC‑AATATAACCGCCGATTTGGTTCACAAAAAGCCAAAACAGACCGCCTAATAACAGCATATTTAACAGGATCGGCATGATCACAAAACGGCGCAGCCCTTTCTGTAAAATCAAATGCCAACCCATCACAAAATATTCAAAACCGGATTTCATTTCCGCGTGTTTCAAc                                                                     >  1:626487/1‑231 (MQ=255)                                                                                                                                                                                                                   acaAAA‑CGCCGCC‑CAACCAA‑CTCAGCAAATCGGAAATGTAGAAAAAAATGCCGTC‑AATATAACCGCCGATTTGGTTCACAAAAAGCCAAAACAGACCGCCTAGTAACAGCATATTTAACAGGATCGGCATGATCACAAAACGGCGCAGCCCTTTCTGTAAAATCAAATGCCAACCCATCACAAAATATTCAAAACCGGATTTCATTTCCGCGTGTTTCAACATAGCGTTATCCCTAAGAAAAAAGAGGTTc                                       <  1:501917/250‑1 (MQ=255)                                                                                                                                                                                                                           cccgcT‑CAACCAC‑CCCACCCAAGCGGAAATGTAGAAAATAAT‑CCGGCAAATATAACCGCCGATTTGGTTCACAAAAAGCCAAAACAGACCGCCTAATAACAGCATATTTAACAGGATCGGCATGATCACAAAACGGCGCAGCCCTTTCTGTAAAATCAAATGCCAACCCATCACAAAATATTCAAAACCGGATTTCATTTCCGCGTGTTTCAACATAGCGTTATCCCTAAGAAAAAAGAGGTTCAGGATAc                                <  1:467052/250‑1 (MQ=255)                                                                                                                                                                                                                              ccT‑CACCCAA‑CTCTGCTAATGGGTAATGTAGAAAATAATGCCTTT‑ATTATAACCGCCGATTTGGTTCAC‑AAAAGCCAAAACAGACCGCCTAATAACAGCATATTTAACAGGATCGGCATGATCACAAAACGGCGCAGCCCTTTCTGTAAAATCAAATACCAACCCATCCCAAAATATTCAAAACCGGATTTCATTTCCGCGTGTTTCAACATAGCGTTATCCCTAAGAAAAAAGAGGTTCAGGATACAAAg                            <  1:629694/250‑1 (MQ=255)                                                                                                                                                                                                                               cccaaaccaa‑ccCAGCCAACCGGAAATTAAGAAAATAATGCCGAC‑AAAAAAACCGCCGATTTGGTTCACAAAAAGCCAAAACAGACCGCCTAATAACAGCATATTTAACAGGATCGGCATGATCACAAAACGGCGCAGCCCTTTCTGTAAAATCAAATGCCAACCCATCACAAAATATTCAAAACCGGATTTCATTTCCGCGTGTTTCAACATAGCGTTATCCCTAAGAAAAAAGAGGTTCAGGATACaaa                             <  1:386347/247‑1 (MQ=255)                                                                                                                                                                                                                                 ccaaccaacCTCAGCAAATGGGAAAGGTAGAAAATAATGCCTAA‑AATATAACCGCCGATTTGGTTCACAAAAAGCCAAAACAGACCGCCTAATAACAGCATATTTAACAGGATCGGCATGATCACAAAACGGCGCAGCCCTTTCTGTAAAATCAAATGCCAACCCATCACAAAATATTCAAAACCGGATTTCATTTCCGCGTGTTTCAACATAGCGTTATCCCTAAGAAAAAAGAGGTTCAGGATACAAAg                            <  1:347820/250‑1 (MQ=255)                                                                                                                                                                                                                                  caaccaa‑cTCACCCAATCGGGAATGTAGAAAATAATGCCGTC‑AATATAACCGCCGATTTGGTTCACAAAAAGCCAAAACAGACCGCCTAATAACAGCATATTTAACAGGATCGGCATGATCACAAAACGGCGCAGCCCTTTCTGTAAAATCAAATGCCAACCCATCACAAAATATTCAAAACCGGATTTCATTTCCGCGTGTTTCAACAt                                                                   <  1:581282/210‑1 (MQ=255)                                                                                                                                                                                                                                     ccaa‑cTAAACCAATCGGGAATGTAGAAAATAATGCCTTC‑AATATAACCGCCGATTTGGTTCACAAAAAGCCAAAACAGACCGCCTAATAACAGCATATTTAACAGGATCGGCATGATCACAAAACGGCGCAGCCCTTTCTGTAAAATCAAATGCCAACCCATCACAAAATATTCAAAACCGGATTTCATTTCCGCGTGTTTCAACATAGCGTTATCCCTAAGAAAAAAGAGGTTCAGGATACAAAGTGCgg                       <  1:322581/251‑1 (MQ=255)                                                                                                                                                                                                                                      caa‑cTCAGCCAATCGGGAATGTAGGAAATAATGCCGTC‑AATATAACCGCCGATTTGGTTCACAAAAAGCCAAAACAGACCGCCTAATAACAGCATATTTAACAGGATCGGCATGATCACAAAACGGCGCAGCCCTTTCTGTAAAATCAAATGCCAACCCATCACAAAATATTCAAAACCGGATTTCATTTCCGCGTGTTTCAACATAGCGt                                                              >  1:26799/1‑211 (MQ=255)                                                                                                                                                                                                                                        a‑cTCAGCCAATCGCGAATGGAGGAAATAATGCCGTC‑AATATAACCGCCGATTTGGTTCACAAAAAGCCAAAACCGACCGCATAATAACAGCATATTTAACAGGATCGGCATGATCACAAAACGGCGCAGCCCTTTCTGTAAAATCAAATGCCAACCCATCACAAAATATTCAAAACCGGATTTCATTTCCGCGTGTTTCAACATAGCGTTATCCCTAAGAAAAAAGAGGTTAAGGCTACAAAGAGCGGTGa                    >  1:358916/1‑251 (MQ=255)                                                                                                                                                                                                                                             aGCCAATCGGGAATGTAGGAAATAATGCCGTC‑AATATAACCGCCGATTTGGTTCACAAAAAGCCAAAACAGACCGCCTAATAACAGCATATTTAACAGGATCGGCATGATCACAAAACGGCGCAGCCCTTTCTGTAAAATCAAATGCCAACCCATCACAAAATATTCAAAACCGGATTTCATTTCCCCGGGTTTAACCATAGCGTAATCCCAAAAAAAAAAGAGTTTAAGATTACAAAGTGCGGGGaaaaa                >  1:512345/1‑251 (MQ=255)                                                                                                                                                                                                                                                          tGTAGGAAATAATGCCGTC‑AATATAACCGCCGATTTGGTTCACAAAAAGCCAAAACAGACCGCCTAATAACAGCATATTTAACAGGATCGGCATGATCACAAAACGGCGCAGCCCTTTCTGTAAAATCAAATGCCAACCCATCACAAAATATTCAAAACCGGATTTCATTTCCGCGTGTTTCAACATAGCGTTATCCCTAAGAAAAAAGAGGTTCAGGATACAAAGCGCGGTGAAAAAAATCAAACTTTGt   >  1:211720/1‑251 (MQ=255)                                                                                                                                                                                                                                                           gTAGGAAATAATGCCGTC‑AATATAACCGCCGATTTGGTTCACAAAAAGCCAGAACAGACCGCCTAATAACAGCATATTTAACAGGATCGGCATGATCACAAAACGGCGCAGCCCTTTCTGTAAAATCAAATGCCAACCCATCACAAAATATTCAAAACCGGATTTCATTTCCGCGTGTTTCAACATAGCGTTATCCCTAAGAAAAAAGAGGTTCAGGATACAAAGTGCGGTGAAAAAAATCAAACTTTGt   >  1:64151/1‑250 (MQ=255)                                                                                                                                                                                                                                                            tAGGAAATAATGCCGTC‑AATATAACCGCCGATTTGGTTCACAAAAAGCCAAAACAGACCGCCTAATAACAGCATATTTAACAGGATCGGCATGATCACAAAACGGCGCAGCCCTTTCTGTAAAATCAAATGCCAACCCATCACAAAATATTCAAAACCGGATTTCATTTCCGCGTGTTTCAACATAGCGTTATCCCTAAGAAAAAAGAGGTTCAGGATACAAAGTGCGGTGAAAAAAATCAAACTTTGTc  <  1:399096/250‑1 (MQ=255)                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                        GCTGTAAACCAGCTTTTGCCATTCCCGTACCAGCATCCGCGGCACATCTTTCATAAAATCCCACAGGGACATATCCTGCAAATTTTCTCCGGTCAGCATTTTTTCCACTTTTTCCGCCAGCAAGCC‑ATTAAAAGGCGCGGCGATAAAGCCCGACAGCGTGGTAAAGACGAAGTAAAACAGAACGAGAATCATAAGAATGGAAAGCAGCAGCAAAA‑CGCCGCT‑CAACCAA‑CTCAGCCAATCGGGAATGTAGGAAATAATGCCGTC‑AATATAACCGCCGATTTGGTTCACAAAAAGCCAAAACAGACCGCCTAATAACAGCATATTTAACAGGATCGGCATGATCACAAAACGGCGCAGCCCTTTCTGTAAAATCAAATGCCAACCCATCACAAAATATTCAAAACCGGATTTCATTTCCGCGTGTTTCAACATAGCGTTATCCCTAAGAAAAAAGAGGTTCAGGATACAAAGTGCGGTGAAAAAAATCAAACTTTGTC  >  1/1046856‑1047352

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