BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	66,688	G→A	intergenic (‑138/+147)	araD ← / ← araA	L‑ribulose‑5‑phosphate 4‑epimerase AraD/L‑arabinose isomerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	66,688	0	G	A	97.1%	274.3 / NA	107	intergenic (‑138/+147)	araD/araA	L‑ribulose‑5‑phosphate 4‑epimerase AraD/L‑arabinose isomerase
Reads supporting (aligned to +/- strand): ref base G (3/0); major base A (56/44); minor base C (3/0); total (63/44)
Fisher's exact test for biased strand distribution p-value = 2.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 107 total aligned reads displayed.

GCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAAT‑CGCAGCGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCA‑GGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCAGATAAGACGCGTCAGCGTCGCATCAGGCGTTACATACCGGATGCGGCTACTTAGCGACGAAACCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGTGTGTCGTTATCAATCACCGTGATTTCAA  >  NC_000913/66445‑66926

  tCGCGATAAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAAT‑CGGAACGCGTAGGc                                                                                                                                                                                                                                        >  2:137351/1‑251 (MQ=255)

    gcgATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAAT‑CGGAACGCGTAGGCCt                                                                                                                                                                                                                                      >  1:512221/1‑251 (MQ=255)

     cgATCAACGGCGCTGACGTTGCCCCATGTGAGCGTTACCAGGGTATGTGTTGGCAGCGCCAGGTGGGCCTCTAATACCTGGCGTTTGAGATATTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAAT‑CGGAACGCGTAGGCCTg                                                                                                                                                                                                                                     <  2:235722/251‑1 (MQ=21)

     cgATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCTCCAGGTTGGCCTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAAT‑CGGAACGCGTAGGCCTg                                                                                                                                                                                                                                     <  1:50976/251‑1 (MQ=255)

     cgATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAAT‑CGGAACGCGTAGGCCTg                                                                                                                                                                                                                                     >  1:104429/1‑251 (MQ=255)

      gATCAACGGCGCTGCCGTTGCCCCAGGTGAGGGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCGCTAATACCGGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAAT‑CGGAACGCGTAGGCCTGa                                                                                                                                                                                                                                    <  2:434521/251‑1 (MQ=21)

          aaCGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAAT‑CGGAACGCGTAGGCCTGATAAg                                                                                                                                                                                                                                <  1:534189/251‑1 (MQ=255)

             ggCGCTGACTTTGACCCATGTGAGCGGGACCAGGTTGTGTTTTGGCAGGGCCAGGTTGGCCTCTAATACCTGGCTGTTGAGATCTTCTAACATGTTGACGCCTTCTTGCCGGATGCGCGTTGCTTATCAGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAAT‑CGGAACGCGTAGGCCTGATAAGAcg                                                                                                                                                                                                                             <  2:418529/251‑1 (MQ=21)

                gcTGACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGTTTGCGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAAT‑CGGAACGCGTAGGCCTGATAAGACgcgc                                                                                                                                                                                                                          >  2:575701/1‑251 (MQ=255)

                   gccgttgacagtagtgAGCGGGCCCAGGTTATGTGTTGGCAGCGCCAGGTTGGCAGAGAAGGCCGGGAGTTTGTGATCTGCTAATAGGTTGACTCTCTCGTTTCGGAGTCGCGTTGCTGACCCGGTCTACAACATCGCGGCGTGGAGGCCTGAGAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAAT‑CGGAACGCGTAGGCCTGAGAAGACGCGCCAg                                                                                                                                                                                                                       <  2:246490/238‑1 (MQ=9)

                   gACGTTGCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGATTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCGGTGCATGCCGCATGCGCTTTGCATATCCGGCCTACAAAAT‑CGGAACACGTAGGCCTGATAAGACGCGCCa                                                                                                                                                                                                                        >  2:652328/1‑250 (MQ=25)

                   gACGTAGCCCCATGTGGGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAAT‑CGGAACGCGTAGGCCTGATAAGACGCGCCAg                                                                                                                                                                                                                       <  2:192334/251‑1 (MQ=255)

                     cGTTGCCCCATGTTAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACACAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGTCGTGGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAAT‑CGGACCGCGTAGGCCTGATACGAcgctcacgat                                                                                                                                                                                                                     >  2:135576/1‑244 (MQ=25)

                      gTTGCCCCATGTGAGCGTGACCAGGTTCTGTTTTGGCTGCGCCAGGTTGGCCTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAAT‑CGGAACGCGTAGGCCTGATAAGACGCGCCAGCGt                                                                                                                                                                                                                    >  1:771602/1‑251 (MQ=255)

                         gCCCCATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGGTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAAT‑CGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCg                                                                                                                                                                                                                  <  2:422055/250‑1 (MQ=255)

                          cccaaTGTGAGCGTGCCCAGGTTATGTTTTGGCAGAGCCAGGTTGGCCTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGCCTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGACAGCGTCGCATCAGGCGTTGAATGCCGGAGGCGCTTTGCTTATCCGGCCTACAACAT‑CGGAACGCGTCGGCCTGATAAGACGCGCCAGCTTCGc                                                                                                                                                                                                                 >  2:737221/5‑250 (MQ=21)

                           cccATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGACAGGGTGGCCTCTAATACCTGGCGTTTGCGATCTTCTAACATGGTGCCTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGTTAAGACGCGCCAGCGTCGCATCAGGCGTTGCATGCCGGATGCGCTTTGATTATCCGGCCTACAACATAGGCA‑CGCGACGGCCTGATAAGACGCGCCAGAGACGCAt                                                                                                                                                                                                               >  2:408356/1‑251 (MQ=18)

                            ccATGTGAGCGTGACCAGGTTAGGTTTTGGCAGCGCCCGGTTGGCCTCTAATTCTGGGGGGTTGATATCTTCTAACATGTTGTCTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCGGCGGGTAGGCCTGATAAGACGCTCCATCGTCGCATCAGGGGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAAT‑CGGAACGCGTAGGCCTGATAAGACGCGCCAGCGGCGCATc                                                                                                                                                                                                              <  2:469732/251‑1 (MQ=18)

                             cATGTGAGCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAAT‑CGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCa                                                                                                                                                                                                             <  1:642175/251‑1 (MQ=255)

                                   ggCGTGACCAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAAT‑CGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCA‑GGCGtt                                                                                                                                                                                                      <  1:746966/250‑1 (MQ=255)

                                           cAGGTTATGTTTTGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAAT‑CGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCA‑GGCGTTGAATGCCg                                                                                                                                                                                              >  1:296435/1‑251 (MQ=255)

                                              gttatgtTTTGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGCATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAAT‑CGGACCGCGTAGGCCTGATAAGACGCTCCAGCGACGCATCA‑CGCGTTGAATGCCgg                                                                                                                                                                                             >  2:662150/5‑249 (MQ=21)

                                               ttGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCTGGCGTTTGAGATCTTCTAACGTGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAA‑T‑CGCAGCGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCA‑GGCGATGAATGCCGGATgc                                                                                                                                                                                         >  2:496030/1‑251 (MQ=255)

                                                tatgtTTTGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAAT‑CGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCA‑GGCGTTGAATGCCGGATgc                                                                                                                                                                                         >  1:387439/3‑251 (MQ=255)

                                                 atgtTTTGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAAT‑CGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATAA‑GGCGTTGAATGCCGGATgcg                                                                                                                                                                                        >  1:988079/2‑251 (MQ=255)

                                                    ttttGGCAGCGCCAGGTTGGCCTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAAT‑CGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCA‑GGCGTTGAATGCCGGATGCGCtt                                                                                                                                                                                     >  1:120580/1‑251 (MQ=255)

                                                        ggCAGCGCCAGGTTGGCCTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGACTACAAAAT‑CGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCA‑GGCGTTGAATGCAGGATGCGCTTTGCt                                                                                                                                                                                 >  1:571900/1‑251 (MQ=255)

                                                         acAGCGCCAGTTTGGCCTCTAATACCTGGCGGTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGCATGCGCTTTGCTTATCCGGCATACAAAATAGCAGAGTGGAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGATTTGCTTATCCGGCCTACAAACT‑CGGGACGCGTAGGCCTGATAAGACGCGACAGCCACACATCA‑GGCTGTGGACGCCGGATGCGCCTTGCtt                                                                                                                                                                                >  2:141894/2‑251 (MQ=18)

                                                           aGCGCCAGGTTGGCCTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAAT‑CGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCA‑GGCGTTGAATGCCGGATGCGCTTTGCTTAt                                                                                                                                                                              >  1:513033/1‑251 (MQ=255)

                                                              gcCAGGTTGGCCTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAAT‑CGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCA‑GGCGTTGAATGCCGGATGCGCTTTGCTTATCCg                                                                                                                                                                           >  1:342012/1‑251 (MQ=255)

                                                              gcCAGGTTGGCCTCTAAGTCCTGGCGTTTGAGATCTTCTAACATGTTGACCCCTTCGTGCCGGACGCGCTTTGCTTACCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAAT‑CGGAACGCGTAGTCCTGATAAGACGCGCCAGCGTCGCATCA‑GGCGTTGAATGCCGGATGCGCTTTGCTTATCCg                                                                                                                                                                           <  1:620176/251‑1 (MQ=21)

                                                                cAGGTTGGCCTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCTGGATGCGCTTTGCTTATCCGGCCTACAAAAT‑CGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCA‑GGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGc                                                                                                                                                                         >  1:632705/1‑251 (MQ=255)

                                                                     tGGCCTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAAT‑CGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCA‑GGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTAc                                                                                                                                                                     >  1:454029/1‑250 (MQ=255)

                                                                      ggcctCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGCCTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAAT‑CGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCA‑GGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACa                                                                                                                                                                    <  1:814474/246‑1 (MQ=255)

                                                                         ctCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTACTTCGTGCCGGATGCGCTTTGCTTATCCGTCCTACAAAATCGCCGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTCCAAAAT‑CGGAGAGCGTAGCCCTGGTAAAACGCGCCAGCGTCGCATCA‑GGCGTTGCATGCCGGATGCGCGTTGCTTCTACGGcccaacaaa                                                                                                                                                                 >  2:717794/2‑243 (MQ=18)

                                                                            acaTACCGGGCGTTTGAGATCTTCTAACCTGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCAGCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAAT‑CGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCTCATCA‑GGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGc                                                                                                                                                             <  2:971348/249‑1 (MQ=21)

                                                                                   tGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAAT‑CGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCA‑GGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCgtgt                                                                                                                                                      >  1:708176/1‑251 (MQ=255)

                                                                                           gagaTCTTCTACCATGTTGACTACTTCGAGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCCGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTCCAAAAT‑CGGAACGCGTCAGCCTGATAAGACGCGCCAGCGTCGCATCA‑GGCGTTGCATGCACGATGCGCTTTGCGTATCGGGCCTACCAAATCGCAGGGTGTAGGcctg                                                                                                                                               >  2:489726/1‑248 (MQ=17)

                                                                                             gaTCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGCCGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAAT‑CGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCA‑GGCGTTGCATGCCGGATGCGCTTTGCTTATCCGGCCTCCAAAATCGCAGCGTGTCGGCCTGAt                                                                                                                                             >  2:404430/1‑250 (MQ=18)

                                                                                                 ttctAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAAT‑CGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCA‑GGCGTTGAATGCCGGATGCACTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAggc                                                                                                                                        >  2:414524/1‑249 (MQ=255)

                                                                                                   ctAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGTCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCTCATCAGGCGTTGAATGACGGATGCGCTTTGCTTATCCGGCCTACCAAAT‑CGGAACGCGTAGGCCTGATAAGACTCGCCAGCGTCGCATAC‑GACGTTGCATGCCGGATGCGCTTTGCTTATCCGGCCTACGAAATCGCAGCGTGTAGGCCTGATATGAcgc                                                                                                                                      >  2:173395/1‑251 (MQ=255)

                                                                                                     aataTGTTGACTCCTTCGTGCCAGATGCGCTTTGCTTATCCGGCCTACAAAAGCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGAGGCGCTTTGCTTATCCGGCCTACAAAAT‑CGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCA‑GGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGAcgc                                                                                                                                      <  1:438982/246‑1 (MQ=21)

                                                                                                      aCATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAAT‑CGGAACGCGTACGCCTGCTAAGACGCGCCAGCGTCGCATCA‑GGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGCTAAGACGCGTc                                                                                                                                   >  2:241907/1‑251 (MQ=255)

                                                                                                          gTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAAT‑CGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCA‑GGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCg                                                                                                                               >  1:495778/1‑251 (MQ=255)

                                                                                                            tGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAAT‑CGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCA‑GGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTc                                                                                                                             >  1:233832/1‑251 (MQ=255)

                                                                                                                 ccTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAAT‑CGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCA‑GGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATc                                                                                                                        >  1:442209/1‑251 (MQ=255)

                                                                                                                  cTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAGTCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAAT‑CGGAACTCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCC‑GGCGTTAAATGCCGGATGCGCTTTGCTTATCCGGCCTACCAAATCGCAGCGCGTAGGCCTGATAAGACGCGGCAGCGgcgccac                                                                                                                        >  2:708766/1‑247 (MQ=18)

                                                                                                                    tCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCCGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACCAAAT‑CGGAACGAGTAGGCCTGATAAGCAGCGCCAGCGTCGCA‑CACGGCGTTGAATGCCGGATGCGCTTTGCTTATCCTGCCGACAACCTCGAAGCGCGTAGGCCTGATGAGACGCGTCAGCGTCGTATCAgg                                                                                                                     >  2:908744/1‑251 (MQ=17)

                                                                                                                    tCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAAT‑CGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCA‑GGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGTCCTGATAAGACGCGTCAGCGTCGCATCa                                                                                                                       <  1:883853/249‑1 (MQ=255)

                                                                                                                     cGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAAT‑CGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCA‑GGCGTTGAATGCCGGATGCGCTTTGATTATCCGGCCTACAAAATCGCCGCGTGTGGGCCTGTTCAGACGCGTCAGCGTCGCATCAggg                                                                                                                    >  2:929328/1‑250 (MQ=18)

                                                                                                                       tGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCCACAAAAT‑CGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCA‑GGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGt                                                                                                                  >  1:369577/1‑251 (MQ=255)

                                                                                                                        gCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAAT‑CGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCA‑GGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGtt                                                                                                                 >  1:135556/1‑251 (MQ=255)

                                                                                                                        gCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAAT‑CGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCA‑GGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGtt                                                                                                                 >  1:306699/1‑251 (MQ=255)

                                                                                                                          cGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAAT‑CGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCA‑GGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTAc                                                                                                               <  1:8746/251‑1 (MQ=255)

                                                                                                                            gATGCGCTGTGCTCGTCCGGCCGACAAATGCGCAGTGTTGAGGCCTGATAAGACGCGCCCGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCGTATCCGGCCTACAAAAT‑CGGAACGCGTAGGCCTGATAAGGCGCGCCAGTGTCGCATCA‑GGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACac                                                                                                             <  2:524771/251‑2 (MQ=17)

                                                                                                                                 gcTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAAT‑CGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCA‑GGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACAcctg                                                                                                        <  2:413114/251‑3 (MQ=255)

                                                                                                                                    ttGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAAT‑CGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCA‑GGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGAt                                                                                                      <  1:791142/250‑1 (MQ=255)

                                                                                                                                    ttGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAAT‑CGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCA‑GGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGAt                                                                                                      <  1:176282/250‑1 (MQ=255)

                                                                                                                                     tGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAAT‑CGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCA‑GGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGc                                                                                                    >  1:313327/1‑251 (MQ=255)

                                                                                                                                          aTCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAAT‑CGGAACGCGTAGGCCTGATAAGACTCGCCAGCGTCGCATCA‑GGCGTTGAATGCCGGATGCGCTTTGCTTATCAGACCTACAAAATCGCAGCGTGTAGGCCTGATCAGAAGCGTCCGCGTCGCATCAGGCGGTACACAcctgcttcgaccc                                                                                               >  2:734619/1‑240 (MQ=18)

                                                                                                                                          aTCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAAT‑CGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCA‑GGCGTTGAATGCCGGATGCGCTTTGCTTATCCTGCCTACAAAATCGCAGCGTGTAGGCCTGCTAAGACGCGTCGGGGTCGCATCAGGCGTTACACACCTGATGCgact                                                                                                >  2:47930/1‑247 (MQ=21)

                                                                                                                                           tCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAAT‑CGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCA‑GGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTAc                                                                                              >  1:95149/1‑251 (MQ=255)

                                                                                                                                              ggCCTACAAATGCGCAGCGTGTAGGCCTAAGAAGACGCGACAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCGTAGCCGGCCTACAAAAT‑CGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCA‑GGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTa                                                                                           <  2:512221/251‑1 (MQ=18)

                                                                                                                                               gCCTACAAAATCGCCGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAAT‑CGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCA‑GGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATGGAAACGGGTCGGCCTGATAAGACGCGTCACAGTCGCATCAcgcggccccaacctgaacccaccacttac                                                                                          >  2:412330/1‑226 (MQ=18)

                                                                                                                                                       aaTCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAAT‑CGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCA‑GGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGaaa                                                                                  <  1:655068/251‑1 (MQ=255)

                                                                                                                                                        aTCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAAT‑CGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCA‑GGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGaaat                                                                                 >  1:511647/1‑250 (MQ=255)

                                                                                                                                                              gCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCGGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAAT‑CGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCA‑GGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTaa                                                                           <  2:442348/251‑1 (MQ=18)

                                                                                                                                                              gCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAAT‑CGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCA‑GGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTa                                                                            >  1:271955/1‑250 (MQ=18)

                                                                                                                                                                  gtAGGCCGGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGGCGGATGCGCTTTGCTTATCCGGCCTACAAAAT‑CGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCA‑GGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATac                                                                        <  2:420781/250‑1 (MQ=18)

                                                                                                                                                                           aTAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAAT‑CGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCA‑GGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGCCGAAATCCGTAATACACTTCGTTcc                                                              >  1:612003/1‑251 (MQ=18)

                                                                                                                                                                            tAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAAT‑CGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCA‑GGCGTTGAATGCCGGCTGCGCTTTGCTTATCCGGCCTACAAAATCGCACCGTGTATGCCTGCTAAGACGCGTCAGCGTCGCACCAGGCGTTACACACCTGATGCGCCTACTTAGCGACGAAATCCGTAATCCCCTTCTTTccc                                                             >  2:693596/1‑250 (MQ=18)

                                                                                                                                                                             aaGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAAT‑CGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCA‑GGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCca                                                             <  1:970517/250‑1 (MQ=18)

                                                                                                                                                                             aaGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAAT‑CGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCA‑GGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCca                                                             <  1:524073/250‑1 (MQ=18)

                                                                                                                                                                                    gcccgCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCGTATCCGGCCTACAAAAT‑CGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCA‑GGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCcagcgcagcg                                                     <  2:341194/247‑1 (MQ=18)

                                                                                                                                                                                      cAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAAT‑CGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCA‑GGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGcgcg                                                   >  1:551457/1‑251 (MQ=18)

                                                                                                                                                                                         cGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAAT‑CGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCA‑GGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCt                                                >  1:987664/1‑251 (MQ=18)

                                                                                                                                                                                             gCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAAT‑CGGAACGCGTAGGCCTGATAAGACGCGCCCGCGTCGCATCA‑GGCGTTGAATGCCGGATGCGCGTTGCTTATCCGGCCTACCCAACCGCAGCGTGTAGGCGGGATGAGCCGCGTCAGCGTCGCATCAGGCGTTAACCACATGCTTCGACTACTTAGCGCAGAAATCCGTAATACACTTCGGTCCCGCGCAGcgcgggttta                                             >  2:5616/1‑244 (MQ=9)

                                                                                                                                                                                             gCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAAT‑CGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCA‑GGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTa                                             >  1:475848/1‑250 (MQ=18)

                                                                                                                                                                                                 cAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAAT‑CGCAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCA‑GGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCTTGTAGGCCTGCTCAGACTCTTCAGCGTCGCATCAGGCGTTACACACCTGCTGCGACTACATATCGACGACATCCGTAATACCCTTCGTTCCAGCGCAGCGCGGCTTTACACGc                                        >  2:6759/1‑251 (MQ=17)

                                                                                                                                                                                                  aGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAAT‑CGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCA‑GGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGCTAAGACGAGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATACGTAATACACTTCGTTCCCGCGCAGCGCGTCTTTAAACgat                                       >  2:575496/1‑249 (MQ=18)

                                                                                                                                                                                                  aGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAAT‑CGCAGCGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCA‑GGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCAGATAAGACGCGTCAGCGTCGCATCAGGCGTTACATACCGGATGCGGCTACTTAGCGACGAAACCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCt                                       >  1:439771/1‑251 (MQ=255)

                                                                                                                                                                                                   ggCGTTGAATGTCGGCTGCGCTTTGCTTACCCGGCCTACAAAAT‑CGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCTCATCA‑GGCGTTGAAGGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTg                                      <  2:135323/251‑1 (MQ=18)

                                                                                                                                                                                                     cGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAAT‑CGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCA‑GGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAACCGCGGGc                                    >  1:954197/1‑251 (MQ=18)

                                                                                                                                                                                                         gAATGCCGGTAGCTCTTTGCTTATCCGGCCAACCAAAA‑CGGAACGCGGAGGCCTGAAAAGACGCGCCGGCGTCGCATCA‑GGTGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAAGCCGTAATACACTGCGTTCCAGCGCAGCGCGTTTTTAAACGCTGGCAGGc                                <  2:571900/251‑1 (MQ=18)

                                                                                                                                                                                                         gAATGCCGGATTCGCTTTGCTTATCCGGCCTACAAAAT‑CGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCA‑GGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGGAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAgg                                 <  2:799294/250‑1 (MQ=18)

                                                                                                                                                                                                                ggATGCGCTTTGCTTATCCGGCCTACAAGAT‑CGGACCGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCA‑GGCGTTGCATGCCGGATGCGCTTTGCTTATCCGGCCTACAAACTCGCAGCGTGTAGGCCTGATAAGACGCGTCAGGGTCGCATCAGGCATTACACCCCTTATGCGACTCCTTAGCGACACAATCCGAAATACGCTTGGTTCCAGCGCCTGGCGTCTTTAAACGCTAGCAGGCgtggtgc                         >  2:475929/1‑247 (MQ=9)

                                                                                                                                                                                                                 gATGCGCTTTGCTTATCCGGCCTACAAAAT‑CGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCA‑GGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACCCCCCTGATGCGACTACTTAGCGTCGAAATCCGTAATAAACTTCGTTCCAGCGCACCGCGTCTTTAACCGCTGGCAGGCggttatcg                        >  2:223981/1‑244 (MQ=17)

                                                                                                                                                                                                                       cTTTGCTTATCCGGCCTACAAAAT‑CGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCA‑GGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGCCGAAATCCGTAATACACTTCGTTCCAGAGCAGCGCGTCTTTAAACGCTGGCAGGCgggtattcttatc                   >  1:537598/1‑238 (MQ=18)

                                                                                                                                                                                                                       cTTTGCTTATCCGGCCTACAAAAT‑CGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCA‑GGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGGATCGTTatc                   >  1:372730/1‑250 (MQ=18)

                                                                                                                                                                                                                         ttGCTTATCCGGCCTACAAAAT‑CGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCA‑GGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGTCGCGTCAGCGTCGCATCAGACGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGCTCCAGCGCAGCGCGTCTTTAAACGCTGACAGGCggatatcgtttacaat                >  2:13493/1‑236 (MQ=18)

                                                                                                                                                                                                                         ttGCTGATCCGGCCTACTAAAT‑CGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCA‑GGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAAGCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTGCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTatca                  <  2:143570/249‑1 (MQ=18)

                                                                                                                                                                                                                           gCTTATCCGGCCTACAAAAT‑CGGAACGCGTAGGCCTGATAAGACGCGCCAGCGGCGCATCA‑GGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTatcaatca              <  1:177175/251‑1 (MQ=21)

                                                                                                                                                                                                                           gCTTAGCCGGCCTACAAAAT‑CGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCA‑GGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTatcaatca              <  1:600919/251‑1 (MQ=21)

                                                                                                                                                                                                                            cTTATCCGGCCTACAAAAT‑CGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCA‑GGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTCCACACCTGATGCGCCTACTTAGCGACGACATCCGTAATACCCTTCGTTCCAGCGCAGCGCGACTTTAAACGCTGGCCGGCggggatcgttgtcactcac             >  2:691235/1‑233 (MQ=17)

                                                                                                                                                                                                                             ttATCCGGCCTACAAAAT‑CGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCA‑GGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGGAACGTTATCAATCAcc            >  1:158268/1‑251 (MQ=18)

                                                                                                                                                                                                                              tcaccgtccagcAAAAT‑CGGAACGCGTAGGCCTGCTTAGACGCGCCAGCGTCGCATCA‑GGCGTTGAATGCCGGATTCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTGAGCGACGAAATCCGTAATTCACTTCGTTCCAGCGCCGCGCGTCTTTAAACGCTGTCAGGCGGGTATCGTTATCAATCAc             <  2:67473/238‑1 (MQ=18)

                                                                                                                                                                                                                                cccGGCCTACAAAAT‑CGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCA‑GGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTg         <  1:737221/250‑1 (MQ=21)

                                                                                                                                                                                                                                   ggCCTACAAAAT‑CGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCA‑GGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGCCTAATTAGCGACGAAATCCGTAATAAACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCggggatcggttttaatcaacctggtt      >  2:878655/1‑226 (MQ=18)

                                                                                                                                                                                                                                    cccTACAAAAT‑CGTTACGCGTAGGCCTGATAATGCGCGGCGGCGTCGCAGCA‑GGCGTTGAATGCCGGATGCGATTTGCTTACACGGGCTACAAAATCGCATCGTGTAGGCCTGATAAGACGCGTCAGCGTAGCATCAGACGTTACACACCTTAGGCGTCTATTTAGCGACGTAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTTGCAGGCGGGTATCGTTATCAATCTCCGTTAttt     <  2:835880/250‑1 (MQ=255)

                                                                                                                                                                                                                                        aCAAAAT‑CGGAACGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCA‑GGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCaa  <  1:575701/250‑1 (MQ=21)

GCTCGCGATCAACGGCGCTGACGTTGCCCCATGTGAGCGTGACCAGGTTGTGTTTTGGCAGCGCCAGGTTGGCTTCTAATACCTGGCGTTTGAGATCTTCTAACATGTTGACTCCTTCGTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAAT‑CGCAGCGCGTAGGCCTGATAAGACGCGCCAGCGTCGCATCA‑GGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCAGATAAGACGCGTCAGCGTCGCATCAGGCGTTACATACCGGATGCGGCTACTTAGCGACGAAACCCGTAATACACTTCGTTCCAGCGCAGCGCGTCTTTAAACGCTGGCAGGCGTGTGTCGTTATCAATCACCGTGATTTCAA  >  NC_000913/66445‑66926

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 12 ≤` `ATCG/ATCG` `< 14 ≤` `ATCG/ATCG` `< 31 ≤` `ATCG/ATCG` `< 39 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`