Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	173,223	A→G	9.3%	S588G (AGT→GGT)	fhuB →	iron(3+)‑hydroxamate import ABC transporter permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	173,223	0	A	G	9.3%	248.8 / 4.0	97	S588G (AGT→GGT)	fhuB	iron(3+)‑hydroxamate import ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base A (40/43);  new base G (9/2);  total (52/45)
Fisher's exact test for biased strand distribution p-value = 5.27e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GACCTGGATTTCCGGTTCGACCTACAACGCGACCGATGCGCAGGTCTGGCGCACCGGAATTGTGATGGTGATTTTGCTGGCGATTACCCCG‑CTGTGCCGCCGCT‑GGCTGACCATTTTACCGCTGGGTGGTGATACCGCCCGAGCCGTA‑GGAATGGCGCTGACGCCGACGCGAATTGCGCTGCTGCTGTTAGCGGCTTGCCTGACGGCGACCGCGACGATGACTATTGGACCGTTGAGTTTTGTTGGTTTAATGGCACCGCATATTGCGCGGATGATGGGCTTTCGACGGACGATGCCACACATCGTAATTTCGGCGCTGGTGGGTGGTTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGG  >  NC_000913/172988‑173460                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                              gACCTGGATTTCCGGTTCGACCTACAACGCGACCGATGCGCAGGTCTGGCGCACCGGAATTGTGATGGTGATTTTGCTGGCGATTACCCCG‑CTGTGCCGCCGCT‑GGCTGACCATTTTACCGCTGGGTGGTGATACCGCCCGAGCCGTA‑GGAATGGCGCTGACGCCGACGCGAATTGCGCTGCTGCTGTTAGCGGCTTGCCTGACGGCGACCGCGACGATGACTATTGGACGGTTGCGtt                                                                                                                                                                                                                                            >  1:805359/1‑239 (MQ=255) gACCTGGATTTCCGGTTCGACCTACAACGCGACCGATGCGCAGGTCTGGCGCACCGGAATTGTGATGGTGATTTTGCTGGCGATTACCCCG‑CTGTGCCGCCGCT‑GGCTGACCATTTTACCGCTGGGTGGTGATACCGCCCGAGCCGTA‑GGAATGGCGCTGACGCCGACCCGAATTGCGCTGCTGCTGTTAGCGGCTTGCCTGACGGCGACCGCGACGATGACCATTGGGCCGTGGAGtt                                                                                                                                                                                                                                            >  1:1319852/1‑239 (MQ=255)  aCCTGGATTTCCGGTTCGACCTACAACGCGACCGATGCGCAGGTCTGGCGCACCGGAATTGTGATGGTGATTTTGCTGGCGATTACCCCG‑CTGTGCCGCCGCT‑GGCTGACCATTTTACCGCTGGGTGGTGATACCGCCCGAGCCGTA‑GGAATGGCGCTGACGCCGGCGCGAATTGCGCTGCTGCTGTTAGCGGCTTGCCTGACGGAGACCGCGACGATGAAGAATGGAACGTTGAGttt                                                                                                                                                                                                                                           >  2:382402/1‑239 (MQ=255)            ccGGTTCGACCTACAACGCGACCGATGCGCAGGTCTGGCGCACCGGAATTGTGATGGTGATTTTGCTGGCGATTACCCCG‑CTGTGCCGCCGCT‑GGCTGACCATTTTACCGCTGGGTGGTGATACCGCCCGAGCCGTA‑GGAATGGCGCTGACGCCGACGCGAATTGCGCTGCTGCTGTTAGCGGCTTGCCTGACGGCGACCGCGACGATGGCTATTGGACCCTTGGGTTTTGTTGGTTTa                                                                                                                                                                                                                                 >  1:9669/1‑239 (MQ=255)                 tCGCCCTACAACGCGACCGATGCGCAGGTCGGGCGCACCGGAATTGTGATGGTGATTTTGCTGGCGATTCCCCCG‑CTGTGCCGCCGCT‑GGCTGACCATTTTACCGCTGGGTGGTGATACCGCCCGAGCCGTA‑GGAATGGCGCTGACGCCGACGCGAATTGCGCTGCTGCTGTTAGCGGCTTGCCTGACGGCGACCGCGACGATGACTATTGGACCGTTGAGTTTTGTTGGTTTAATGGc                                                                                                                                                                                                                            <  1:262508/239‑1 (MQ=255)                 tCGACCTACAACGCGCCCGATGCGCAGGTCTGGCGCACCGGCATTGTGCTGGGGCTTTTGCTGGCGATTACCCCG‑CTGTCCCGCCGCT‑GGCTGCCCATTTTACCGCTGGGTGGGGATACCGCCCGAGCCGTA‑GGAATGGCGCTGCTGCCGACGCGAATTGCGCTCCTGCTGCTAGCGCACTGCCTGACGGCCTCCGCGCAGCCGTCTATTGGACCCGTGGGTTTTGTATGtttcaggg                                                                                                                                                                                                                             >  2:1243481/1‑233 (MQ=255)                 tCGACCTACAACGCGACCGATGCGCAGGTCTGGCGCACCGGAATTGTGATGGTGATTTTGCTGGCGATTACCCCG‑CTGTGCCGCCGCT‑GGCTGACCATTTTACCGCTGGGTGGTGATACCGCCCGAGCCGTA‑GGAATGGCGCTGACGCCGACGCGAATTGCGCTGCTGCTGTTAGCGGCTTGCCTGACGGCGCCCGCGACGATGACCATTGGACCGTTGAGTTTTGTTGGTTTAATgg                                                                                                                                                                                                                             >  1:182451/1‑238 (MQ=255)                    aCCTACAACGCGACCGATGCGCAGGTCTGGCGCACCGGAATTGTGATGGTGATTTTGCTGGCGATTACCCCG‑CTGTGCCGCCGCT‑GGCTGACCATTTTACCGCTGGGTGGGGATACCGCCCGAGCCGTA‑GGCATGGCGCTGACGCCGACGCGAATTGCGCTGCTGCTGTTAGCGGCTTGCCTTGCGGCGGCCGGGCCGAGGCAGCTTTGGCCGGTGGGTTTTTTTGGGTTAATGGCAcc                                                                                                                                                                                                                         >  2:1065474/1‑239 (MQ=255)                               gACCGATGCGCAGGTCTGGCGCACCGGAATTGTGATGGTGATTTTGCTGGCGATTACCCCG‑CTGTGCCGCCGCT‑GGCTGACCATTTTACCGCTGGGTGGTGATACCGCCCGAGCCGTA‑GGAATGGCGCTGACGCCGACGCGAATTGCGCTGCTGCTGTTAGCGGCTTGCCTGACGGCGACCGCGACGATGACTATTGGGCCGTTGGGTTTTTTTGGTTTAATGGCACCGGACATTgggg                                                                                                                                                                                                              >  2:219741/1‑236 (MQ=255)                               gACCGATGCGCAGGTCTGGCGCACCGGAATTGTGATGGTGATTTTGCTGGCGATTACCCCG‑CTGTGCCGCCGCT‑GGCTGACCATTTTACCGCTGGGTGGTGATACCGCCCGAGCCGTA‑GGAATGGCGCTGACGCCGACGCGAATTGCGCTGCTGCTGTTAGCGGCTTGCCTGACGGCGACCGCGACGATGACTATTGGACCGTTGGGTTTTGTTGGTTTATTGGCACCGCATATTgcgc                                                                                                                                                                                                              >  1:420120/1‑239 (MQ=255)                               gACCGATGCGCAGGTCTGGCGCACCGGAATTGTGATGGTGATTTTGCTGGCGATTACCCCG‑CTGTGCCGCCGCT‑GGCTGACCATTTTACCGCTGGGTGGTGATACCGCCCGAGCCGTA‑GGAATGGCGCTGACGCCGACGCGAATTGCGCTGCTGCTGTTAGCGGCTTGCCTGACGGCGACCCCGACGATGACTATTGGACCGTTGAGTTTTGTTGGTTTAATGGCACCGCATATTggcc                                                                                                                                                                                                              >  1:369401/1‑236 (MQ=255)                                                 gcgcACCGGAATTGTGATGGTGATTTTGCTGGCGATTACCCCG‑CTGTGCCGCCGCT‑GGCTGACCATTTTACCGCTGGGTGGTGATACCGCCCGAGCCGTA‑GGAATGGCGCTGACGCCGACGCGAATTGCGCTGCTGCTGTTAGCGGCTTTCCTGACGGCGACCGCGGCGCTGACTATTGGACCGTTGGGTTTTTTTGGGTTATGGGCCCCCCCTCTTTCGCGGGTGTGGGGtttttgtc                                                                                                                                                                                            >  2:779824/1‑235 (MQ=255)                                                 gcgcACCGGAATTGTGATGGTGATTTTGCTGGCGATTACCCCG‑CTGTGCCGCCGCT‑GGCTGACCATTTTACCGCTGGGTGGTGATACCGCCCGAGCCGTA‑GGAATGGCGCTGACGCCGACGCGAATTGCGCTGCTGCTGTTAGCGGCTTGCCTGACGGCGACCGCGACGATGACTATTGGACCGTTGAGTTTTGTTGGTTTAATGGCACCGCATATTGCGCGGATAATGGGCttttgcc                                                                                                                                                                                            >  1:355702/1‑235 (MQ=255)                                                    acccGGAATTGTGATGGTGATTTTGCTGGCGATTACCCCG‑CTGTGCCGCCGCT‑GGCTGACCATTTTACCGCTGGGTGGTGATACCGCCCGAGCCGTA‑GGAATGGCGCTGACGCCGACGCGAATTGCGCTGCTGCTGTTAGCGGCTTGCCTGACGGCGACCGCGACGATGACTATTGGACCGTTGAGTTTTGTTGGTTTAATGGCACCGCATATTGCGCGGATGATGGGCTTTCgacgga                                                                                                                                                                                         <  1:42580/237‑1 (MQ=255)                                                     acgggATTTGGGATGGTGATTTGGCGGGCGATTCCCCCG‑GTGGGCCGCCGCT‑GGCTGACCATTTTCCCGGTGGGTGGTGATACCGCCCGAGCCGTA‑GGAATGGCGCTGACGCCGACGCGAATTGCGCTGCTGCTGTTAGCGGCTTGCCTGACGGCGACCGCGACGATGACTATTGGACCGTTGAGTTTTGTTGGTTTAATGGCACCGCATATTGCGCGGATGATGGGCTTTCgacgga                                                                                                                                                                                         <  1:1014413/235‑1 (MQ=255)                                                      ccGGAATTGTGATGGTGATTTTGCTGGCGATTACCCCG‑CTGTGCCGCCGCT‑GGCTGACCATTTTACCGCTGGGTGGTGATACCGCCCGAGCCGTA‑GGAATGGCGCTGACGCCGACGCGAATTGCGCTGCTGCTGTTAGCGGCTTGCCTGACGGCGACCGCGACGATGACTATTGGGCCGTTGAGTTTTGTTGGTTTAATGGCACCGCATATTGCGCGGATGATGGGCTTTCgacggacg                                                                                                                                                                                       >  1:439148/1‑239 (MQ=255)                                                      ccGGAATTGTGATGGTGATTTTGCTGGCGATTACCCCG‑CTGTGCCGCCGCT‑GGCTGACCATTTTACCGCTGGGTGGTGATACCGCCCGAGCCGTA‑GGAATGGCGCTGACGCCGACGCGAATTGCGCTGCTGCTGTTAGCGGCTTGCCTGACGGCGACCGCGACGATGACTATTGGGCCGTTGAGTTTTGTTGGTTTAATGGCACCGCATATTGCGCGGAGGATGGGCTTTCgccggcgg                                                                                                                                                                                       >  1:777870/1‑236 (MQ=255)                                                      ccGGAATTGTGATGGTGATTTTGCTGGCGATTACCCCG‑CTGTGCCGCCGCT‑GGCTGACCATTTTACCGCTGGGTGGTGATACCGCCCGAGCCGTA‑GGAATGGCGCTGACGCCGACGCGAATTGCGCTGCTGCTGTTAGCGGCTTGCCTGACGGCGACCGCGACGATGACTATTGGACCGTTGAGTTTTGTTGGTTTAATGGCACCGCCTATTGCGCGGATGATGGGCTTTCgacggacg                                                                                                                                                                                       >  1:22800/1‑239 (MQ=255)                                                      ccGGAATTGTGATGGTGATTTTGCTGGCGATTACCCCG‑CTGTGCCGCCGCT‑GGCTGACCATTTTACCGCTGGGTGGTGATACCGCCCGAGCCGTA‑GGAATGGCGCTGACGCCGACGCGAATTGCGCTGCTGCTGTTAGCGGCTTGCCTGACGGCGACCGCGACGATGACTATTGGACCGTTGAGTTTTGTTGGTTTAATGGCACCGCATATTGCGCGGATGATGGGCTTTCgccggcgg                                                                                                                                                                                       >  1:1040947/1‑236 (MQ=255)                                                            ttGTGATGGTGATTTCGGTGGCGATTCCCCCG‑TTGGCCCGCCGCT‑GGCTGACCATTTTACCGCGGGGTGGTGATACCGCCCGAGCCGTA‑GGAATGGCGCTGACGCCGACGCGAATTGCGCTGCTGCTGTTAGCGGCTTGCCTGACGGCGACCGCGACGATGACTATTGGACCGTTGAGTTTTGTTGGTTTAATGGCACCGCATATTGCGCGGATGATGGGCTTTCGACGGACGATGcc                                                                                                                                                                                  <  1:867385/238‑1 (MQ=255)                                                            ttGGGATGGTGATTTTGCTGGCGATTCCCCCG‑CTGTGCCGCCGCT‑GGCTGCCCATTTGACCGCTGGGTGGTGATACCGCCCGAGCCGTA‑GGAATGGCGCTGACGCCGACGCGAATTGCGCTGCTGCTGTTAGCGGCTTGCCTGACGGCGACCGCGACGATGACTATTGGACCGTTGAGTTTTGTTGGTTTAATGGCACCGCATATTGCGCGGATGATGGGCTTTCGACGGACGATGCca                                                                                                                                                                                 <  1:823781/239‑1 (MQ=255)                                                                 aTGGTGATTTTGCTGGCGATTACCCCG‑CTGTGCCGCCGCT‑GGCTGACCATTTTACCGCTGGGTGGTGATACCGCCCGAGCCGTA‑GGAATGGCGCTGACGCCGACGCGAATTGCGCTGCTGCTGTTAGCGGCTTGCCTGACGGCGACCGCGACGATGACTATTGGACCGTTGAGTTTTGTTGGTTTAATGGCACGGCATATTGCGCGGATGATGGGCTTTCTACGGCGGAGGCcacaca                                                                                                                                                                             >  1:294626/1‑238 (MQ=255)                                                                 aTGGTGATTTTGCTGGCGATTACCCCG‑CTGTGCCGCCGCT‑GGCTGACCATTTTACCGCTGGGTGGTGATACCGCCCGAGCCGTA‑GGAATGGCGCTGACGCCGACGCGAATTGCGCTGCTGCTGTTAGCGGCTTGCCTGACGGCGACCGCGACGATGACTATTGGACCGTTGAGTTTTGTTGGTTTAATGGCACCGCATATTGCGCGGATGATGGGCTTTCGACGGACGATGCcacaca                                                                                                                                                                             >  1:113255/1‑238 (MQ=255)                                                                    gTGATTTTGCTGGCGATTACCCCG‑CTGTGCCGCCGCT‑GGCTGACCATTTTACCGCTGGGTGGTGATACCGCCCGAGCCGTA‑GGAATGGCGCTGACGCCGACGCGAATTGCGCTGCTGCTGTTAGCGGCTTGCCTGACGGCGACCGCGACGCTGACTATTGGACCGTTGCGTTTTGTTGGTTTAATGGCACCCCATATTTCGCGGATGAGGGGCTTTCGACGGACGATGCGACCAATCGt                                                                                                                                                                         >  1:276151/1‑239 (MQ=255)                                                                    gTGATTTTGCTGGCGATTACCCCG‑CTGTGCCGCCGCT‑GGCTGACCATTTTACCGCTGGGTGGTGATACCGCCCGAGCCGTA‑GGAATGGCGCTGACGCCGACGCGAATTGCGCTGCTGCTGTTAGCGGCTTGCCTGACGGCGACCGCGACGATGACTATTGGACCGTTGAGTTTTGTTGGTTTAATGGCCCCGCATATTGCGCGGCTGATGGGCTTTCTACGGACGATGCCACACATCGt                                                                                                                                                                         >  1:1002269/1‑239 (MQ=255)                                                                              gggcggtccccccG‑CTGTGGCGCGGGT‑GGTGGCCCAGTTTGCCGTGGGGTGGGGACCCCGCCGGGGCCGGA‑GGAAGGGCGCTGACGCCGGCGCGAATTGCGCTGCTGCTGTTAGCGGCTTGCCTGACGGCGACCGCGACGATGACTATTGGACCGTTGAGTTTTGTTGGTTTAATGGCACCGCATATTGCGCGGATGATGGGCTTTCGACGGACGATGCCACACATCGTAATTTCGgcg                                                                                                                                                               <  2:113255/230‑1 (MQ=255)                                                                                gCGATTACCCCG‑CTGTGCCGCCGCT‑GGCTGACCATTTTACCGCTGGGTGGTGATACCGCCCGAGCCGTA‑GGAATGGCGCTGACGCCGACGCGAATTGCGCTGCTGCTGTTAGCGGCTTGCCTGACGGCGACCGCGCAGATGCTAATTGGGACGGTGGGTTTGGTTGGTTTAATGGCACCGCGTATTGCCGGGTGGGATGGCTTTCGACGGACGTTGCCACACCTCGGAATTACGgcggg                                                                                                                                                             >  2:777632/1‑237 (MQ=255)                                                                                  aaTTCCCGCG‑CTGGGCCG‑CGCTGGGCGGCCCATTGTACCGTGGGGTGGTGATACCGCCCGAGCCGTA‑GGAATGGCGCTGACGCCGACGCGAATTGCGCTGCTGCTGTTAGCGGCTTGCCTGACGGCGACCGCGACGATGACTATTGGACCGTTGAGTTTTGTTGGTTTAATGGCACCGCATATTGCGCGGATGATGGGCTTTCGACGGACGATGCCACACATCGTAATTTCGGCGCtg                                                                                                                                                            <  1:724709/237‑1 (MQ=255)                                                                                   gatcccccGCCTGGGCCGCCGCG‑GGCGGCCCATTTCACCGGGGGGTGGTGATCCCGCCCGAGCCGTA‑GGAATGGCGCTGACGCCGACGCGAATTGCGCTGCTGCTGTTAGCGGCTTGCCTGACGGCGACCGCGACGATGACTATTGGACCGTTGAGTTTTGTTGGTTTAATGGCACCGCATATTGCGCGGATGATGGGCTTTCGACGGACGATGCCACACATCGTAATTTCGGCGCtg                                                                                                                                                            <  1:1355748/234‑1 (MQ=255)                                                                                   aTTACCCCG‑CTGTGCCGCCGCT‑GGCTGACCATTTTACCGCTGGGTGGTGATACCGCCCGAGCCGTA‑GGAATGGCGCTGACGCCGACGCGAATTGCGCTGCTGCTGTTAGCGGCTTGCCTGACGGCAACCGCGACGATGACTATTGGACCGTTGGGGTTTGTTGGGTTAATGGCACCGCATATTGCGCGGATGATGGGCTTTCTACGGGCGATGCCACACATCGTCATTTCGGCGCtggt                                                                                                                                                          >  1:273047/1‑239 (MQ=255)                                                                                    ttcccccG‑TTGTCCCGCCGCT‑GGCTGACCATTTTACCGCGGGGTGGTGATCCCGCCCGAGCCGTA‑GGAATGGCGCTGACGCCGACGCGAATTGCGCTGCTGCTGTTAGCGGCTTGCCTGACGGCGACCGCGACGATGACTATTGGACCGTTGAGTTTTGTTGGTTTAATGGCACCGCATATTGCGCGGATGATGGGCTTTCGACGGACGATGCCACACATCGTAATTTCGGCGCtggt                                                                                                                                                          <  1:758139/235‑1 (MQ=255)                                                                                    ttACCCCG‑CGTGGCCGCCGCG‑GGCGGCCCATTTTCCCGCGGGGTGGTGATACCGCCCGAGCCGTA‑GGAATGGCGCTGACGCCGACGCGAATTGCGCTGCTGCTGTTAGCGGCTTGCCTGACGGCGACCGCGACGATGACTATTGGACCGTTGAGTTTTGTTGGTTTAATGGCACCGCATATTGCGCGGATGATGGGCTTTCGACGGACGATGCCACACATCGTAATTTCGGCGCtggt                                                                                                                                                          <  1:766798/238‑1 (MQ=255)                                                                                       cccctccTGTGCCGGCGTT‑GGTTGACCATTTTGCCGGTGGGTTGTGATCCCGCCCGAGCCGTA‑GGAATGGCGCTGACGCCGACGCGAATTGCGCTGCTGCTGTTAGCGGCTTGCCTGACGGCGACCGCGACGATGACTATTGGACCGTTGAGTTTTGTTGGTTTAATGGCACCGCATATTGCGCGGATGATGGGCTTTCGACGGACGATGCCACACATCGTAATTTCGGCGCTggtggg                                                                                                                                                       <  1:222757/233‑1 (MQ=255)                                                                                               gtgcgcggggt‑gggtgTCCTTTTTCCCCCTGGGGGGTAATCGCCCCCGAGCCGGA‑GGGATTGCGCTGACGCCGACGCGAATTGCGCGGCTGCTGTTTGCGGCTTGCCTGACGGCGACCGCGACGATGACTATTGGACCGTTGAGTTTTGTTGGTTTAATGGCACCGCATATTGCGCGGATGATGGGCTTTCGACGGACGATGCCACACATCGTAATTTCGGCGCTGGTGGGTGGTTTa                                                                                                                                                <  2:1040947/228‑1 (MQ=255)                                                                                                tgCCGCCGCT‑GGCTGACCATTTTACCGCTGGGTGGTGATACCGCCCGAGCCGTA‑GGAATGGCGCTGACGCCGACGCGAATTGCGCTGCTGCTGTTAGCGGCTTGCCTGACGGCCACCGCGACGATGACTATTGGACCGTTGAGTTTTGTTGGTTTAATGGCCCCGCATATTGCGCGGATGATGGGCTTTCGACGGTCGATGCCACACAACCGAATTTCGGCGCGGGGGGGGGGtttta                                                                                                                                               >  2:258936/1‑236 (MQ=255)                                                                                                       cgtcggcggACCATTTGGCCGGTGGGTGGTGATACCGCCCGAGCCGTA‑GGAATGGCGCTGACGCCGACGCGAATTGCGCTGCTGCTGTTAGCGGCTTGCCTGACGGCGACCGCGACGATGACTATTGGACCGTTGAGTTTTGTTGGTTTAATGGCACCGCATATTGCGCGGATGATGGGCTTTCGACGGACGATGCCACACATCGTAATTTCGGCGCTGGTGGGTGGTTTACTGCTGgt                                                                                                                                        <  1:607216/235‑1 (MQ=255)                                                                                                           ggctgACCATTTTCCCGCTGGGTTGTGATCCCGCCCGAGCCGTA‑GGAATGGCGCTGACGCCGACGCGAATTGCGCTGCTGCTGTTAGCGGCTTGCCTGACGGCGACCGCGACGATGACTATTGGACCGTTGAGTTTTGTTGGTTTAATGGCACCGCATATTGCGCGGATGATGGGCTTTCGACGGACGATGCCACACATCGTAATTTCGGCGCTGGTGGGTGGTTTACTGCTGGTGTTc                                                                                                                                    <  1:954085/239‑1 (MQ=255)                                                                                                                     tttcccggggggTGGTGGTCCCCGCCGGGCCGTA‑GGAATGGCCCTGACGCCGACGCGAATTGCGCTGCTTCTGTTAGCGGCTTGCCTGACGGCGACCGCGACGATGACTATTGGACCGTTGAGTTTTGTTGGTTTAATGGCACCGCATATTGCGCGGATGATGGGCTTTCGACGGACGATGCCACACATCGTAATTTCGGCGCTGGTGGGTGGTTTACTGCTGGTGTTCGCTGACTGgt                                                                                                                          <  1:10582/230‑1 (MQ=255)                                                                                                                      ttACCGCTGGGTGGTGATACCGCCCGAGCCGTA‑GGAATGGCGCTGACGCCGACGCGAATTGCGCTGCTGCTGTTAGCGGCTTGCCTGACGGCGACCGCGACGATGACCATTGGTCCGTTGCGTTTTGTTGGTTTAATGGCACCGCCTATTgc                                                                                                                                                                                                                >  1:1091126/1‑152 (MQ=255)                                                                                                                         ccGGTGGGGGGTGAACCCGCCCGAGCCGTA‑GGAATGGCGCTGACGCCGACGCGAATTGCGCTGCTGCTGTTAGCGGCTTGCCTGACGGCGACCGCGACGATGACTATTGGACCGTTGAGTTTTGTTGGTTTAATGGCACCGCATATTGCGCGGATGATGGGCTTTCGACGGACGATGCCACACATCGTAATTTCGGCGCTGGTGGGTGGTTTACTGCTGGTGTTCGCTGACTGGTGTgg                                                                                                                      <  1:869489/239‑1 (MQ=255)                                                                                                                            cTGGGTGGTGATACCGCCCGAGCCGTA‑GGAATGGCGCTGACGCCGACGCGAATTGCGCTGCTGCTGTTAGCGGCTTGCCTGACGGCGACCGCGACGATGACTATTGGACCGTTGGGTTTTGTTTGTTTAATGGCACCGCATATTGCGCGGATGATGGGCTTTCGACGGACGATGCCACACATCGTAATTTCGGCGCTGGGGGGTGGTTTACTGCTGGTGTTCGCTGACTGGTGGGGGCg                                                                                                                   >  1:534399/1‑239 (MQ=255)                                                                                                                            cTGGGTGGTGATACCGCCCGAGCCGTA‑GGAATGGCGCTGACGCCGACGCGAATTGCGCTGCTGCTGTTAGCGGCTTGCCTGACGGCGACCGCGACGATGACTATTGGACCGTTGAGTTTTGTTGGTTTAATGGCCCCGCATATTGCGCGGATGATGGGCTTTCGACGGACGATGCCCCACCTCGTACTTTCGGCGCTGGTGGGTGGTTTACTGCTGGTGTTCGCTGACTGTGGTGGGCg                                                                                                                   >  1:1271446/1‑239 (MQ=255)                                                                                                                               ggtggtggtCCCGCCCGCGTCGTATGGATTGGCG‑TGACGCCGACGGGAATTGCGCTGCTGCTGTTAGCGGCTTGCCTGACGGCGACCGCGACGATGACTATTGGACCGTTGAGTTTTGTTGGTTTAATGGCACCGCATATTGCGCGGATGATGGGCTTTCGACGGACGATGCCACACATCGTAATTTCGGCGCTGGTGGGTGGTTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGa                                                                                                                 <  2:151076/238‑1 (MQ=255)                                                                                                                               ggtggtGATACCGCCCGAGCCGTA‑GGAATGGCGCTGACGCCGACGCGAATTGCGCTGCTGCTGTTAGCGGCTTGCCTGACGGCGACCGCGACGATGACTATTGGACCGTTGAGTTTTGTTGGTTTAATGGCACCGCATATTGCGCGGATGATGGGCTTTCGACGGACGAGGCCCCCCATCGTAATTTCGGCGCGGGGGGGGGGTTTACTGCTGGGGTTCGCTGACTGGTGTgggggggt                                                                                                                >  1:759402/1‑237 (MQ=255)                                                                                                                               ggtgggGATCCCGCCCGGGCCGTA‑GGAATGGCGCTGACGCCGACGCGAATTGCGCTGCTGCTGTTAGCGGCTTGCCTGACGGCGACCGCGACGATGACTATTGGACCGTTGAGTTTTGTTGGTTTAATGGCACCGCATATTGCGCGGATGATGGGCTTTCGACGGACGATGCCACACATCGTAATTTCGGCGCTGGTGGGTGGTTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGa                                                                                                                 <  1:1041992/238‑1 (MQ=255)                                                                                                                               ggtgggCATCCCGCCCGAGCCGTA‑GGAATGGCGCTGACGCCGACGCGAATTGCGTTGCTGCTGTTAGCGGCTTGCCTGACGGCGACCGCGACGATGACTATTGGACCGTTGAGTTTTGTTGGTTTAATGGCACCGCATATTGCGCGGATGATGGGCTTTCGACGGACGATGCCACACATCGTAATTTCGGCGCTGGTGGGTGGTTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGa                                                                                                                 <  1:878126/238‑1 (MQ=255)                                                                                                                                ggggttaacccgccgGAGCCGTA‑GGAATGGCGCTGGCGCCGACGCGAATTGCGCTGCTGCTGTTAGCGGCTTGCCTGACGGCGACCGCGACGATGACTATTGGACCGTTGAGTTTTGTTGGTTTAATGGCACCGCATATTGCGCGGATGATGGGCTTTCGACGGACGATGCCACACATCGTAATTTCGGCGCTGGTGGGTGGTTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGAt                                                                                                                <  1:919652/229‑1 (MQ=255)                                                                                                                                        gccGCCCGGGCCGTA‑GGAATGGCGCTGACGGCGACGCAAATTGCGCTGCTGCTGTTAGCGGCTTGCCTGACGGCGACCGCGACGATGACTATTGGACCGTTGAGTTTTGTTGGTTTAATGGCACCGCATATTGCGCGGATGATGGGCTTTCGACGGACGATGCCACACATCGTAATTTCGGCGCTGGTGGGTGGTTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGtt                                                                                                       <  1:75391/238‑1 (MQ=255)                                                                                                                                         ccGCCCGGGCCTGA‑AGAAGGGCGCTGACGCCGACGCGAATTGCGCTGGTGCTGTTAGCGGGTTGCCTGACGGCGACCGCGACGATGACTATTGGACCGTTGAGTTTTGTTGGTTTAATGGCACCGCATATTGCGCGGATGATGGGCTTTCGACGGACGATGCCACACATCGTAATTTCGGCGCTGGTGGGTGGTTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGttt                                                                                                      <  1:1136213/239‑1 (MQ=255)                                                                                                                                         ccGCCCGAGCCGTA‑GGAATGGCGCTGACGCCGACGCGAATTGCGCTGCTGCTGTTAGCGGCTTGCCTGACGGCGACCGCGACGATGACTATTGGACCGTTGAGTTTTGTTGGTTTAATGGCAcc                                                                                                                                                                                                                         <  2:321865/124‑1 (MQ=255)                                                                                                                                         ccGCCCGAGCCGTA‑GGAATGGCGCTGACGCCGACGCGAATTGCGCTGCTGCTGTTAGCGGCTTGCCTGACGGCGACCGCGACGATGACTATTGGACCGTTGAGTTTTGTTGGTTTAATGGCAcc                                                                                                                                                                                                                         >  1:321865/1‑124 (MQ=255)                                                                                                                                         ccGCCCGAGCCGTA‑GGAATGGCGCTGACGCCGACGCGAATTGCGCTGCTGCTGTTAGCGGCTTGCCTGACGGCGACCGCGACGATGACTATTGGACCGTTGAGTTTTGTTGGTTTAATGGCACCGCATATTGCGCGGATGATGGGCTTTCGACGGACGATGCCACACATCGTAATTTCGGCGCTGGTGGGTGGTTTACTGCTGGTGTTCGATGACTGGTGGGGGCGGATGGGGCTGttt                                                                                                      >  1:151076/1‑239 (MQ=255)                                                                                                                                            cccGGGCCCTA‑GGAATGGCGCGTACGCCGCAGCGAATTGCGCTTGTGCTGTTAGCGGCTTGCCTGACGGCGACCGCGACGATGACTATTGGACCGTTGAGTTTTGTTGGTTTAATGGCACCGCATATTGCGCGGATGATGGGCTTTCGACGGACGATGCCACACATCGTAATTTCGGCGCTGGTGGGTGGTTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTttcca                                                                                                   <  2:420120/239‑1 (MQ=255)                                                                                                                                            cccGAGCCGTA‑GGAATGGCGCTGACGCCGACGCGAATTGCGCTGCTGCTGTTAGCGGCTTGCCTGACGGCGACCGCGACGATGACTATTGGACCGTTGAGTTTTGTTGGTTTAATGGCACCGCATATTGCGCGGATGATGGGCTTTCGACGGACGATGCCACACATCGTAATTTCGGCGCTGGTGGGTGGTTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTttcca                                                                                                   >  1:417982/1‑239 (MQ=255)                                                                                                                                                    gTA‑GGAATGGCGCTGACGCCGACGCGAATTGCGCTGCTGCTGTTAGCGGCTTGCCTGACGGCGACCGCGACGATGACTATTGGACCGTTGAGTTTTGTTGGTTTAATGGCACCGCATATTGCGCGGATGATGGGCTTTCGACGGACGATGCCACACATCGTAATTTCGGCGCTGGTGGGGGGTTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGa                                                                                            >  1:1194370/1‑238 (MQ=255)                                                                                                                                                     tA‑GGAATGGCGCTGACGCCGACGCGAATTGCGCTGCTGCTGTTAGCGGCTTGCCTGACGGCGACCGCGACGATGACTATTGGCCCGTTGAGTTTTGTTGGTTTAATGGCACCGCATATTGCGCGGGTGATGGGCTTTCGACGGACCCTGCCCCACTTCGTAATTTCGCCGCTGGGGGGGTGTTTGCTGGTGGGGTGTGCTGACTGGTTTGGGCGGGTGGGGGTGGTTCCCTTCCCGATc                                                                                          >  2:901829/1‑239 (MQ=255)                                                                                                                                                      a‑GGAATGGCGCTGCCGCCGACGCGAATTGCGCTGCTTCTGTTAGCGGCTTGCCTGACGGCGACCGCGACGTTGACTATTGGACCGTTGAGTTTTGTTGGTTTAATGGCACCGCATATTGCGCGGATGATGGGCTTTCGACGGACGATGCCACACATCGTAATTTCGGCGCTGGTGGGTGGTTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATc                                                                                          <  1:1158788/238‑1 (MQ=255)                                                                                                                                                             ggCGCTGACGCCGACGCGAATTGCGCTGCTGCTGTTAGCGGCTTGCCTGACGGCGACCGCGACGATGACTATTGGACCGTTGAGTTTTGTTGGTTTAATGGCACCGCATATTGCGCGGATGATGGGCTTTCGACGGACGATGCCACACATCGTAATTTCGGCGCTGGTGGGTGGTTTACTGCTGGTGTTCGCTGACTGGGGTGGGCGGATGGGGCTGTTTCCATTCCAGATCCcggcgg                                                                                   >  1:184014/1‑239 (MQ=255)                                                                                                                                                                 cTGACGCCGACGCGAATTGCGCTGCTGCTGTTAGCGGCTTGCCTGACGGCGACCGCGACGATGACTATTGGACCGTTGAGTTTTGTTGGTTTAATGGCACCGCATATTGCGCGGATGATGGGCTTTCGACGGACGATGCCACACATCGTAATTTCGGCGCTGGGGGGTGGTTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCTTTCCGAATCCCGGCggg                                                                                  >  1:1291082/1‑236 (MQ=255)                                                                                                                                                                    aCGCCGACGCGAATTGCGCTGCTGCTGTTAGCGGCTTGCCTGACGGCGACCGCGACGATGACTATTGGACCGTTGAGTTTTGTTGGTTTAATGGCACCGCATATTGCGCGGATGATGGGCTTTCGACGGACGATGCCACACATCGTAATTTCGGCGCTGGTGGGTGGtt                                                                                                                                                  <  2:890651/169‑1 (MQ=255)                                                                                                                                                                    aCGCCGACGCGAATTGCGCTGCTGCTGTTAGCGGCTTGCCTGACGGCGACCGCGACGATGACTATTGGACCGTTGAGTTTTGTTGGTTTAATGGCACCGCATATTGCGCGGATGATGGGCTTTCGACGGACGATGCCACACATCGTAATTTCGGCGCTGGTGGGTGGtt                                                                                                                                                  >  1:890651/1‑169 (MQ=255)                                                                                                                                                                    aCGCCGACGCGAATTGCGCTGCTGCTGTTAGCGGCTTGCCTGACGGCGACCGCGACGATGACTATTGGACCGTTGAGGTTTGTTGGTTTAATGGCACCGCATATTGCGCGGATGATGGGCTTTCGACGGACGATGCCCACCATCGTAATTTCGGCGCTGGTGGGTGGTTTACTGCTGGGGTTCGCTGACTGGTGTGGGCGGATGGGGCTGTTTCCCTTCCAGATCCCGGCGGGgctgct                                                                            >  1:830191/1‑239 (MQ=255)                                                                                                                                                                      gCCGACGGCAATTGCCCTTCTGCTGTTAGCGGCTTGCCTGACGGCGACCGCGACGATGACTATTGGACCGTTGAGTTTTGTTGGTTTAATGGCACCGCATATTGCGCGGATGATGGGCTTTCGACGGACGATGCCACACATCGTAATTTCGGCGCTGGTGGGTGGTTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGt                                                                          <  1:777632/239‑1 (MQ=255)                                                                                                                                                                      cccgtcccgATTTTCGTTGTTGCTTTTTGGGGTTTGCTTGACGGCGACCGCGACGATGGCTATTGGGCCGTTGGGTTTTGTTTGTTTAATGGCACCGCATATTGCGCGGATGATGGGCTTTCGACGGACGATGCCACACATCGTAATTTCGGCGCTGGTGGGTGGTTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGt                                                                          <  2:534399/232‑1 (MQ=255)                                                                                                                                                                        cGACGGGAATTGCGCGGCTGTTTTTAGCGGTTTGCCTGGCGGGACCCGGGACGATGACTTTTGGACCGTTGAGTTTTGTTGGTTTAATGGCACCGCATATTGCGCGGATGATGGGCTTTCGACGGACGATGCCACACATCGTAATTTCGGCGCTGGTGGGTGGTTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCa                                                                        <  2:1291082/239‑1 (MQ=255)                                                                                                                                                                                       tgctgctgTTAGCGGCTTGCCTGACGGCGACCGCGACGATGACTATTGGACCGTTGAGTTTTGTTGGTTTAATGGCACCGCATATTGCGCGGATGATGGGCTTTCGACGGACGATGCCACACATCGTAATTTCGGCGCTGGTGGGTGGTTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGgcgcg                                                         <  1:382402/239‑1 (MQ=255)                                                                                                                                                                                        ttttgtgttgtcGGTTTGCGTGACGGGGACCGCGCAGCTGACTATTGGCCGGTTGAGTTTGTTGGGTTTAATGGCACGGCATATTGGCGGGGTGATGGGCTTTGGACGGACGATGCCACACATCGTAATTTCGGCGCTGGTGGGTGGTTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGgcgcgc                                                        <  2:1319852/228‑1 (MQ=255)                                                                                                                                                                                        gctgctgTTAGCGGCTTGCCTGACGGCGACCGCGACGATGACTATTGGACCTTTGAGTTTTGTTGGTTTAATGGCACCGCATATTGCGCGGATGATGGGCTTTCGACGGACGATGCCACACATCGTAATTTCGGCGCTGGTGGGTGGTTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGgcgcgc                                                        <  1:353321/239‑1 (MQ=255)                                                                                                                                                                                        gcgtcgtttAGGGGTTTGCCTTGCGGGGGCCGGGACGATGACGTTTGGGCCGGTGAGTTTTGTTGGTTTAATGGCACCGCATATTGCGGGGATGATGGGCTTTCGACGGACGATGCCACACATCGTAATTTCGGCGCTGGTGGGTGGTTTACTGCGGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGgcgcgc                                                        <  2:1194370/232‑1 (MQ=255)                                                                                                                                                                                            ctgTTAGCGGCTTGCCTGACGGCGACCGCGACGATGACTATTGGACCGTTGAGTTTTGTTGGTTTAATGGCACCGCATATTGCGCGGCTGATGGGCTTTCGACGGACGATGCCACACATCGTAATTTCGGCGCTGGTGGGTGGTTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGCCGGGGCTGCTGTCAACCTTTTTCGGGGCGCCata                                                    >  1:896579/1‑239 (MQ=255)                                                                                                                                                                                              gtttgtgggttGCCTGTCGGCGACCGGGGGGATGACTTTTGGACCGTTGATTTTTGTTGGTTTAAGGGCACCGCATATTGCGCGGATGATGGGCTTTCGACGGACGATGCCACACATCGTAATTTCGGCGCTGGTGGGTGGTTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATAtt                                                  <  2:654538/230‑1 (MQ=255)                                                                                                                                                                                                 aGCGGCTTGGCTGACGGTGACCGCGACGATTGCTATTGGACCGTTGAGTTTTGGTGGTTTAATGGCACCGCATATTGCGCGGATGATGGGCTTTCGACGGACGATGCCACACATCGTAATTTCGGCGCTGGTGGGTGGTTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTa                                               <  1:253840/239‑1 (MQ=255)                                                                                                                                                                                                 aGCGGCTTGCCTGACGGCGACCGCGACGATGACTATTGGACCGTTGAGTTTTGTTGGTTTAATGGCACCGCATATTGCGCGGATGATGGGCTTTCGACGGACGATGCCACACATCGTAATTTCGGCGCTGGTGGGTGGTTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCGTGTGTCAACCTTTATTGGCGCGCCCTATTTTa                                               >  1:1121696/1‑239 (MQ=255)                                                                                                                                                                                                   cGGCTTGCCTGACGGCGACCGCGACGATGACTATTGGACCGTTGAGTTTTGTTGGTTTAATGGCACCGCATATTGCGCGGATGATGGGCTTTCGACGGACGATGCCACACATCGTAATTTCGGCGCTGGTGGGTGGTTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATTGGGGCGCCATAtttctt                                              >  1:695650/1‑235 (MQ=255)                                                                                                                                                                                                       ttGCCTGACGGCGACCGCGACGATGACTATTGGACCGTTGAGTTTTGTTGGTTTAATGGCACCGCATATTGCGCGGATGATGGGCTTTCGACGGACGATGCCACACATCGTAATTTCGGCGCTGGTGGGTGGTTTa                                                                                                                                                <  2:894572/136‑1 (MQ=255)                                                                                                                                                                                                       ttGCCTGACGGCGACCGCGACGATGACTATTGGACCGTTGAGTTTTGTTGGTTTAATGGCACCGCATATTGCGCGGATGATGGGCTTTCGACGGACGATGCCACACATCGTAATTTCGGCGCTGGTGGGTGGTTTa                                                                                                                                                >  1:894572/1‑136 (MQ=255)                                                                                                                                                                                                         gCCTGACGGCGACCGCGACGATGACTATTGGACCGTTGAGTTTTGTTGGTTTAATGGCACCGCATATTGCGCGGATGATGGGCTTTCGACGGACGATGCCACACATCGTAATTTCgg                                                                                                                                                                 <  1:656890/117‑1 (MQ=255)                                                                                                                                                                                                         gCCTGACGGCGACCGCGACGATGACTATTGGACCGTTGAGTTTTGTTGGTTTAATGGCACCGCATATTGCGCGGATGATGGGCTTTCGACGGACGATGCCACACATCGTAATTTCgg                                                                                                                                                                 >  2:656890/1‑117 (MQ=255)                                                                                                                                                                                                         gCCTGACGGCGACCGCGACGATGACTATTGGACCGTTGAGTTTTGTTGGTTTAATGGCACCGCATATTGCGCGGATGATGGGCTTTCGACGGACGATGCCACACATCGTAATTTCGGCGCTGGGGGGGGGTTTACCGCTGGGGTTCGCTGCCTGGTGTGGGCGGATGGGGCTGTTTCCTTTCCGAATCCCGGCGGGGGGGGGTTCAACCTTTTTCAGCGGGCCCCATTTTTTTTAtttt                                       >  2:1032631/1‑238 (MQ=255)                                                                                                                                                                                                                     cgggtgcggtGACTGTTGGGCCGTGGAGGTTTGTTGGTTTAAGGGCCCCGGATGTTGGCGGGAGGGGGGGTTTTCGACGGACGATGCCACACATCGTAATTTCGGCGCTGGTGGGTGGTTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCa                            <  2:182451/229‑1 (MQ=255)                                                                                                                                                                                                                     cggggtcgaTGACTATTGGACCGTTGAGTTTTGTTGGTTTAATGGCACCGCATATTGCGCGGATGATGGGCTTTCGACGGACGATGCCACACATCGTAATTTCGGCGCTGGTGGGTGGTTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCa                            <  1:383327/232‑1 (MQ=255)                                                                                                                                                                                                                                     tGGACCGTTGAGTTTTGTTGGTTTAATGGCACCGCATATTGCGCGGATGATGGGCTTTCGACGGACGATGCCACACATCGTAATTTCGGCGCTGGTGGGTGGTTTACTGCTGGGGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGTTCCCGGGGGGGCTGCTGTCTACCTTTATTGGGCGGCCATATTTTTTTTTTTTTTTTGGAAAGCAGACCCCGTAAATTTTcc           >  2:439685/1‑239 (MQ=255)                                                                                                                                                                                                                                       gACCGTTGAGTTTTGTTGGTTTAATGGCACCGCATATTGCGCGGATGATGGGCTTTCGACGGACGATGCCACACATCGTAATTTCGGCGCTGGTGGGTGGTTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGGTGGTAACCTTTATCGGCGGCCCCTTTTTTTTCTTTTTGTTGGGAAAGGAGAGCCCTTAATTTTTTCat         >  2:1303162/1‑239 (MQ=255)                                                                                                                                                                                                                                       gACCGTTGAGTTTTGTTGGTTTAATGGCACCGCATATTGCGCGGATGATGGGCTTTCGACGGACGATGCCACACATCGTAATTTCGGCGCTGGTGGGTGGTTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCCCCATCTTTTATCTATTTGTTTGGAAAGCAGAGCCGGTAAATTTTCCat         >  1:1121654/1‑239 (MQ=255)                                                                                                                                                                                                                                           gTTGGGTTTTGTTGGTTTAATGGCACCGCATATTGCGCGGATGATGGGCTTTCGACGGACGATGCCACACATCGTAATTTCGGCGCTGGTGGGTGGTTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATaa     <  2:417982/239‑1 (MQ=255)                                                                                                                                                                                                                                           gTTGAGTTTTGTTTGTTTAATGGCACCGCATATTGCGCGGATGATGGGCTTTCGACGGACGATGCCACACATCGTAATTTCGGCGCTGGTGGGTGGTTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATaa     <  1:258936/239‑1 (MQ=255)                                                                                                                                                                                                                                           gTTGAGTTTTGTTGGTTTAATGGCACCGCATATTGCGCGGATGATGGGCTTTCGACGGACGATGCCACACATCGTAATTTCGGCGCTGGTGGGTGGTTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATaa     <  1:439685/239‑1 (MQ=255)                                                                                                                                                                                                                                           gTTGAGTTTTGTTGGTTTAATGGCACCGCATATTGCGCGGATGATGGGCTTTCGACGGACGATGCCACACATCGTAATTTCGGCGCTGGTGGGTGGTTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATaa     <  1:326948/239‑1 (MQ=255)                                                                                                                                                                                                                                            ttGAGTTTTGTTGGTTTAATGGCACCGCATATTGCGCGGATGATGGGCTTTCGACGGACGATGCCACACATCGTAATTTCGGCGCTGGTGGGTGGTTTACTGCTGGTGTTCGCTGACTGGTGTgg                                                                                                                      >  2:796464/1‑125 (MQ=255)                                                                                                                                                                                                                                            ttGAGTTTTGTTGGTTTAATGGCACCGCATATTGCGCGGATGATGGGCTTTCGACGGACGATGCCACACATCGTAATTTCGGCGCTGGTGGGTGGTTTACTGCTGGTGTTCGCTGACTGGTGTgg                                                                                                                      <  1:796464/125‑1 (MQ=255)                                                                                                                                                                                                                                            ttGAGTTTTGTTGGTTTAATGGCACCGCATATTGCGCGGATGATGGGCTTTCGACGGACGATGCCACACATCGTAATTTCGGCGCTGGTGGGTGGTTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTGAACCTTTATCGGCGCGCCCATTTTTATCTATTTGTTGAGAAAACAGAGCCGTTAATTTTTCCCTATAAc    >  1:631213/1‑239 (MQ=255)                                                                                                                                                                                                                                            ttGAGTTTTGTTGGTTTAATGGCACCGCATATTGCGCGGATGATGGGCTTTCGACGGACGATGCCACACATCGTAATTTCGGCGCTGGTGGGTGGTTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGGCCGCCACATTTTTTCTTTTTGTTTGGAAAGCAGAGCCGTTAATTTTTCCCTTTAAc    >  1:821598/1‑239 (MQ=255)                                                                                                                                                                                                                                            ttGAGTTTTGTTGGTTTAATGGCACCGCATATTGCGCGGATGATGGGCTTTCGACGGACGATGCCACACATCGTAATTTCGGCGCTGGTGGGTGGTTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAAc    >  1:60359/1‑239 (MQ=255)                                                                                                                                                                                                                                            ttGAGTTTTGTTGGTTTAATGGCACCGCATATTGCGCGGATGATGGGCTTTCGACGGACGATGCCACACATCGTAATTTCGGCGCTGGTGGGTGGTTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCGGTCACCCTTTTTTCGCGCGCCCAATTTTTTTTTTTTGTTGGGAAAACAAAGCCGTTAATTTTTTCaaaaaaa    >  2:1184492/1‑233 (MQ=255)                                                                                                                                                                                                                                            ttGAGTTTTGTTGGTTTAATGGCACCGCATATTGCGCGGATGATGGGCTTTCGACGGACGATGCCACACATCGTAATTTCGGCGCTGGTGGGTGGTTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCGGTCAACCTTTATCGGCGCGCCATATTTTTTCTATTTGTTGAGAAAGCAGAGCCGTTCATTTTTCCCTTTAAc    >  1:886525/1‑239 (MQ=255)                                                                                                                                                                                                                                            ttGAGTTTTGTTGGTTTAATGGCACCGCATATTGCGCGGATGATGGGCTTTCGACGGACGATGCCACACATCGTAATTTCGGCGCTGGTGGGTGGTTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCCTGCAACCTTTATCGGTGCGCCCTATTTTTCCTATTTTTTGAGAAAAAAGAACCGTTAAATTTTCCCAATaaa    >  2:657689/1‑238 (MQ=255)                                                                                                                                                                                                                                              gAGTTTTGGTGGTTTAATGGCACCGCATGTTGCGCGGATGATGGGCTTTCGACGGACGATGCCACACATCGTAATTTCGGCGCTGGTGGGTGGTTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACgg  <  1:901829/239‑1 (MQ=255)                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                              GACCTGGATTTCCGGTTCGACCTACAACGCGACCGATGCGCAGGTCTGGCGCACCGGAATTGTGATGGTGATTTTGCTGGCGATTACCCCG‑CTGTGCCGCCGCT‑GGCTGACCATTTTACCGCTGGGTGGTGATACCGCCCGAGCCGTA‑GGAATGGCGCTGACGCCGACGCGAATTGCGCTGCTGCTGTTAGCGGCTTGCCTGACGGCGACCGCGACGATGACTATTGGACCGTTGAGTTTTGTTGGTTTAATGGCACCGCATATTGCGCGGATGATGGGCTTTCGACGGACGATGCCACACATCGTAATTTCGGCGCTGGTGGGTGGTTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGG  >  NC_000913/172988‑173460

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