Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,627,027	T→C	S372S (AGT→AGC)	radA →	DNA recombination protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,627,027	0	T	C	96.4%	242.0 / NA	83	S372S (AGT→AGC)	radA	DNA recombination protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (31/49);  minor base A (2/1);  total (33/50)
Fisher's exact test for biased strand distribution p-value = 5.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGATGAAGTGACCTCCGGTAGCTCAGTGATGGTGGTATGGGAAGGAACGCGTCCACTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAACCCACGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTGTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGTGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGCGACAGACCGCTGCCGCAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATCTCTGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTGCCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGAC  >  NC_000913/4626783‑4627267                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                 gggaggAAGTGACCTCCGGCAGCGCCGTGATGGTGGTGTGGGAAGGAACGCGGCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACGCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGAc                                                                                                                                                                                                                                            <  2:124842/246‑1 (MQ=255) gagATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGGGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGAc                                                                                                                                                                                                                                            <  1:520473/249‑1 (MQ=255)   gATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACtt                                                                                                                                                                                                                                          <  2:134670/251‑1 (MQ=255)   gATGAAGTGACCTCCGGCAGCTCCGTGATGGTGGGGAGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACtt                                                                                                                                                                                                                                          >  1:19933/1‑251 (MQ=255)         gTGACCTCCGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGtt                                                                                                                                                                                                                                    <  1:79769/251‑1 (MQ=255)         gTGACCGCCGGCAGCTCCGGGAGGGTGGGGTGGGAAGGCACGCGTCCTTTGCTGGTGGAGATTCAGGCGCGGGTCGATCACTCGATGAGGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTGTGTGAACGTGGGCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGtt                                                                                                                                                                                                                                    <  2:67722/251‑1 (MQ=255)           gcccgccGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGGCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTa                                                                                                                                                                                                                                   <  2:116143/245‑1 (MQ=255)                ccggcagtgacgggATGGTGGGGTGGGAAGGAACGCGTGCGTTGCTGGTGGAGATTCGGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGACCAGGATGTGTTTGTGAACGGGGGCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTActgc                                                                                                                                                                                                                               <  2:250875/236‑1 (MQ=255)                ccGGCAGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTGGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTActgctg                                                                                                                                                                                                                             <  1:559294/251‑1 (MQ=255)                  ggcagctacgTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGACGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGACTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGCGGGCTGCAAAATGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGCCGTGAAGTTAACCGACACCAGCGCCGACTTAGCGTTACTGCTGGc                                                                                                                                                                                                                           >  2:614623/10‑251 (MQ=255)                    caGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGCGATCAAGGCGCTGGTCGATCACTCGCTGCTGGCGAATCCGCGCCGCGTGGCAGTGCGGCAGGAACCAAAACATCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGCTCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAACCCAGCGCCGACTTAGCGTTACTGCTGGCg                                                                                                                                                                                                                          >  1:369158/2‑250 (MQ=255)                    caGCTCCGTGATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCg                                                                                                                                                                                                                          >  1:377975/2‑250 (MQ=255)                             taTGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTGGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTcc                                                                                                                                                                                                                 <  2:427706/249‑1 (MQ=255)                             gATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCa                                                                                                                                                                                                                <  1:11480/251‑1 (MQ=255)                             gATGGTGGTGTGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGAGGGTTTCCa                                                                                                                                                                                                                <  2:107908/251‑1 (MQ=255)                              aTGGTGGTGTGGGAAGGATCGCGTCGGTTGCTGGTGGCGATGCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGGGGCAGTGGGGCTGGAACAAAACCTTCGGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACGTAGCGTTACTGCTGGCGATGGTTTCCAg                                                                                                                                                                                                               <  2:5273/251‑1 (MQ=255)                                ggtggtGTGGGAAGGAACGCGCCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGAGGGGCTGGAACAAAACCGTATGGCCATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCACATGGCCGATCAGGATGTGGTTGCGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGAGCCGACTTGGCGTTTCTGCTGGCGATGGTTTCCaacc                                                                                                                                                                                                             >  2:37388/1‑248 (MQ=255)                                   ggtgtGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCTCTGGTCGATCACGCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGGCTTGCACCGTCACGGGGGTCTGCAAATTGCCGATCAGGATGTGTTTGTGAACGTGGGCGGCGGCGTGAAGGTAGCCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTgc                                                                                                                                                                                                          <  1:325285/247‑1 (MQ=255)                                    gtttgtgaaggacCGCGTCCGTTGCTGGGGGAGAGATAGTCGCTGGTCGAGCACTCGATGCTGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCGGCTGGCGGTCTTGCACCGTCGCGGTGGGCTGCAAATGGCCGATCTGGATGTGTGTGTGAGCGTGGTCGGCGGCGTGAAGGGAACGGAGCCCAGCGCCGACTTAGCGTTACTGCTTGCGAGGGCTTCCATCCTgcg                                                                                                                                                                                                         <  2:545188/245‑1 (MQ=255)                                     tgtGGGAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTG‑AACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTgcgtg                                                                                                                                                                                                       >  2:270827/3‑249 (MQ=255)                                      gtgggaaggagcgggtccgtggggggtggAGATTCAGGCGCTTGTCTATCACTCGATGATGGCGAATCCGCGTCGCGGTGCAGTGTGGCTGGAACAAATCCGTCTGGCAATCCTGCTGGCTGTCGTGCTCCGGCACGGTGGGCTGCAAAGGGCCGATCGGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTCACCGAATCCAGCGCCGACTTAGCGGTTCGGCTGGCGAGGGTTTCCAGCCTgcgt                                                                                                                                                                                                        <  2:528169/226‑2 (MQ=255)                                        gggAAGGAACGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGAc                                                                                                                                                                                                     <  1:541860/251‑1 (MQ=255)                                            aggaaCGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGAc                                                                                                                                                                                                 <  1:295177/251‑1 (MQ=255)                                             ggaaCGCGTCCGTTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGAcc                                                                                                                                                                                                <  2:67028/251‑1 (MQ=255)                                                aCGCGTCCGTTGCTGGTGGGGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCt                                                                                                                                                                                             <  2:64639/251‑1 (MQ=255)                                                    gtccgttgctgGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGcc                                                                                                                                                                                          >  1:12899/7‑250 (MQ=255)                                                        gttgttgggggAGAGTCAGGCGCTGGTAGATCACGCGAGGATGGCGAATCCGCTCCGCGTGGCAGTGTGGCTGCAACAAAACCGTCTGGCAATCCTGCGGGCTGTCTTGCACCGTCACGGTGGTCTGCAAAGGGCCGATCAGGATGGGTTTGTGAACGTGGGCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACGTAGCGTTACTGCTGGCGATGGTTTCCAGCCGGCGTGACAGACCGCTGCCACAgg                                                                                                                                                                                     <  2:236166/242‑1 (MQ=255)                                                          tgctgGTGGATATGCAGGCGCTGGTCGATCTCTCGAGGATGGGGAATTCGCGCGGCGTGGCTGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGCCTTGGACCGTAACGGTGGTCTGCAAATGGCCGATCAGGATGTTTTTGTGAACGGTGTCGGCGTCGTGAAGGTAACCGAAACCAGCGCCGGCTTAGCGTTACTGCTGGCGATGGTTTCCAGACTGGGTGACGGACCGCTGCCTCAGGa                                                                                                                                                                                    <  2:296880/250‑1 (MQ=255)                                                           gctgGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATc                                                                                                                                                                                  <  1:358132/251‑1 (MQ=255)                                                               gtggAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACTAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGATACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCtggt                                                                                                                                                                              >  1:396074/1‑251 (MQ=255)                                                                  gagaTTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTggtggt                                                                                                                                                                           >  2:110846/1‑251 (MQ=255)                                                                    gaTTCAGGCGCTGGTCGATCCCTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGgtgt                                                                                                                                                                         <  2:260382/251‑1 (MQ=255)                                                                     aTTCAGGCGCTGGTCGATCACTCGATGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTGGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGtt                                                                                                                                                                        <  1:481648/251‑1 (MQ=255)                                                                         gggCGCTGGTCGATCACTCGATGAGGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCGGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGt                                                                                                                                                                    <  2:298548/250‑1 (MQ=255)                                                                               tGGTCGATCACTCGCTGATGGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTc                                                                                                                                                              >  1:487516/1‑251 (MQ=255)                                                                                ggTCGTTCACTCGATGAGGTCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAATCCGTCTGGCAAGCCTGCTGGCTGTCTTGCACCGTCACGGGGGTCTGCAAATGGCCGATCAGGATGTGTTTGGGAACGTGGTCGGCGGCGTGAAGGTAGCCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGt                                                                                                                                                               <  2:192886/249‑1 (MQ=255)                                                                                     aTCACTCGATGATGGCGAATCCGCGCCGCGAGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTg                                                                                                                                                        >  1:473823/1‑251 (MQ=255)                                                                                               gatggcgaatccggggcgcgTGGCAGGGGGGCTGGAACAAAACCGTCGGGCAATCCTGCTGGCTGTCTTGCTCCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCTTTACTGCTGGCGATGGTGTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGgag                                                                                                                                               <  2:418419/234‑1 (MQ=255)                                                                                               gatgGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGACCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAgggg                                                                                                                                                 >  1:450746/1‑248 (MQ=255)                                                                                                atgGCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGGGAAGGTAACCGAAACCAGCGCAGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGGGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGGGAAGTCGGGCGGGCAGGGgag                                                                                                                                               >  2:541907/1‑249 (MQ=255)                                                                                                  ggCGAATCCGCGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATcc                                                                                                                                           <  1:350882/251‑1 (MQ=255)                                                                                                    cgaatccgctgcgcgTGGAGGTGGTGCTGGAACAAAACCTTGTGGCAATCCTGCTGGCTGTCTTGCACCGTCCCGGTGGTCTGCAAATGGCCGATCAGGATGTGTTGGTGAACGTGGTCGGCGGCGTGCAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCGTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATACGc                                                                                                                                         <  2:354240/240‑1 (MQ=255)                                                                                                          cgcgccggGTTGCAGTGGGGCTGGAACAAAACCGTCTGGCAAGCCGGCTGGCGGTCTTGCATCGTCACGGTGGTCTGCAAAGGGCCGATCAGGATGTGTTTGTGAATGGGGTCGGCTGCGTGAAGGTATCCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCGGCCACAGGATCGGGGGGTGTGTGGTTAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTg                                                                                                                                   <  2:585259/249‑1 (MQ=255)                                                                                                                    ggCAGTGGGGCTGGACCCAAACCGTCTGGCAATCCTACTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGGCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACCGACCGCTGCCACAGGTTCTGGTTGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGCTCCGCCCGGTGCCCAGCgggc                                                                                                                         >  2:313269/1‑249 (MQ=255)                                                                                                                                  aaCAAAACCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTGGTGAACGTGGGCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAAtttca                                                                                                           <  2:384892/251‑5 (MQ=255)                                                                                                                                    cAAAACCGTCTGGCAATCCTGCGGGTTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAAtttcaga                                                                                                         <  2:362233/251‑7 (MQ=255)                                                                                                                                      aaaCCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAg                                                                                                       <  2:5724/251‑1 (MQ=255)                                                                                                                                       aaCCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAgc                                                                                                      >  1:487810/1‑251 (MQ=255)                                                                                                                                       aaCCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAgc                                                                                                      >  1:55617/1‑251 (MQ=255)                                                                                                                                        aCCGTCTGGCACTCCTGCTGGCTGTCGTGCCCCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAgcg                                                                                                     <  1:337828/251‑1 (MQ=255)                                                                                                                                        aCCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGACACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAgcg                                                                                                     >  1:53969/1‑251 (MQ=255)                                                                                                                                        aCCGTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACGGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAgcg                                                                                                     <  1:321767/251‑1 (MQ=255)                                                                                                                                           gTCTGGCAATCCTGCTGGCTGTCTTGCACCGTCACGGTGTTCTGCAAATAGCCGAGCAGGATGTGTTTGTTAACGTGTGCGGCGGCGTGAAGGTAACCGAAACCTGCGCCTTCTTAGCGTTACTGCTGGCGATGGTTTCCAGCTGGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAgcggcg                                                                                                  <  1:610333/251‑1 (MQ=255)                                                                                                                                                 cAATCCTGCTGGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGGCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAacac                                                                                            <  2:284588/251‑1 (MQ=255)                                                                                                                                                      ctgctgGCTGTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGttt                                                                                       >  1:161414/1‑251 (MQ=255)                                                                                                                                                        gctggctgTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTTCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGttt                                                                                       <  1:541907/249‑1 (MQ=255)                                                                                                                                                             ctgTCTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCgg                                                                                >  1:586620/1‑251 (MQ=255)                                                                                                                                                               tgctggCACCGTCACGGTGGGCTGCAAAGGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGGGTGAAGGTAACCGAAACCAGCGCCGACTCAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTTACAGACCGCTGCCACAGGATCTGGTGGTGTGTGGTGAAGTCGGGCTTGCAGTGGAGATCCGCCCGGTGCCCAGCTTTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCggg                                                                               <  1:302464/245‑1 (MQ=255)                                                                                                                                                                    gcgccGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCTGCGCCGACTGAGCGTTACTGCTGGCGATGGTTGCCAGCCTGCGTTACAGTCCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTTCCCAGCTGTCATGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTg                                                                         <  1:555816/248‑1 (MQ=255)                                                                                                                                                                    gCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGAtt                                                                          <  1:459261/250‑1 (MQ=255)                                                                                                                                                                         gTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGTGATTGTTCcg                                                                    >  2:364362/1‑251 (MQ=255)                                                                                                                                                                            aCGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCcggcg                                                                 <  2:487516/251‑1 (MQ=255)                                                                                                                                                                             cGGTGGTCGGCGAATGGTAGAACAGAATATGTTTCTGAACGTGGAAGGAGGCGAGAAGGTAACCGAAACAAGAGCCAAATTAGCGTTACTGCTCGCGATGGTTTCCAGCCTGAGTGACAGACCGCTGCAACAGGATATGGTGGTGTTTGGCGCAGTAGGGCTGGCAGGGGAGCTCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCtggcgg                                                                >  2:584517/1‑251 (MQ=255)                                                                                                                                                                                      gCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGACACACGGCTTTCGCCGGGCGATTGTTCCGGCGGCTAACGtac                                                       >  2:48942/1‑249 (MQ=255)                                                                                                                                                                                      gCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTGTGGTGAAGTCGGGCTGGCAGGGGAGATTGGCCCGGTTCCCAGCGGTCAGGAACGAATTTCAGATGCGGCTAAACACGGTTTTCGCCGGGCGATGGTTCCGGCGGCTAACGta                                                        <  1:242988/250‑2 (MQ=255)                                                                                                                                                                                         aaTGGCCGATCCGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCg                                                     >  1:251392/1‑250 (MQ=255)                                                                                                                                                                                          aTGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACggaccga                                                    >  2:280320/1‑244 (MQ=255)                                                                                                                                                                                            ggCCGATAAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGGGTTTGGTGAAGTCGGGCTGGCAGGGTAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGAGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTCACGTACCgcaa                                                  >  2:534463/1‑247 (MQ=255)                                                                                                                                                                                                gATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAgc                                             <  1:461587/251‑1 (MQ=255)                                                                                                                                                                                                 aTCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGAACCGAAAAAAgcg                                            >  2:218230/1‑251 (MQ=255)                                                                                                                                                                                                           gtTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGa                                  <  2:165790/251‑1 (MQ=255)                                                                                                                                                                                                            tttGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTATCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGAGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAACAAGCGCCGGAAGGGAt                                 >  1:51891/1‑251 (MQ=255)                                                                                                                                                                                                              tgtgAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGc                               >  2:39741/1‑251 (MQ=255)                                                                                                                                                                                                                   aCGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGAtt                          <  1:295407/251‑1 (MQ=255)                                                                                                                                                                                                                      tGGTCGGCGGCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGAttttt                       >  2:162447/1‑251 (MQ=255)                                                                                                                                                                                                                          cggcggcggGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCg                   <  2:141703/251‑1 (MQ=255)                                                                                                                                                                                                                          cggcggCGTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGACGGAAGGGATGCAGATTTTTGGCg                   >  2:202727/1‑251 (MQ=255)                                                                                                                                                                                                                          cggcggCGTGAAGGTAACCGAAACCAGAGCCGACTTAGTGTTACGGCTGGGGATGGTTTACAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGAGCCGCCCGGTGCCCAGCGGTCAGGAACGAATTGCAGAAGCGGCGAAACACGGTTTTCGCCGGGCTATTGTTCCGGCGGCTAACGTTCCGAAAAAAGCGCCGGAAGGGAGGCAGATTTGTGGCg                   <  2:579741/251‑1 (MQ=255)                                                                                                                                                                                                                               gcgTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTaaa              <  2:546089/251‑1 (MQ=255)                                                                                                                                                                                                                                 gTGAAGGTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGAGGATTTTTGGCGTTaaaa             >  2:616319/1‑250 (MQ=255)                                                                                                                                                                                                                                      ggTAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAActct       >  1:607675/1‑251 (MQ=255)                                                                                                                                                                                                                                        tAACCGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTcc     >  1:331568/1‑251 (MQ=255)                                                                                                                                                                                                                                           ccGAAACCAGCGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGAc  <  1:39741/251‑1 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                 GCGATGAAGTGACCTCCGGTAGCTCAGTGATGGTGGTATGGGAAGGAACGCGTCCACTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAACCCACGCCGCGTGGCAGTGGGGCTGGAACAAAACCGTCTGGCAATCCTGCTGGCTGTGTTGCACCGTCACGGTGGTCTGCAAATGGCCGATCAGGATGTGTTTGTGAACGTGGTCGGCGGCGTGAAGGTAACCGAAACCAGTGCCGACTTAGCGTTACTGCTGGCGATGGTTTCCAGCCTGCGCGACAGACCGCTGCCGCAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATCTCTGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTGCCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGAC  >  NC_000913/4626783‑4627267

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