Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,351	A→G	L206L (CTA→CTG)	thrC →	threonine synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,351	0	A	G	93.7%	222.8 / 2.3	79	L206L (CTA→CTG)	thrC	threonine synthase
Reads supporting (aligned to +/- strand):  ref base A (3/2);  new base G (32/42);  total (35/44)
Fisher's exact test for biased strand distribution p-value = 6.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.20e-01

GACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTTACCGAATGTGAAAGTGGTTATCCTCTATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAACTGAAAGTGGCGCTAGGGTTAAACTCGGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCTGTTGCGCAGCTGCCGCAGGAGACGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTCTGCCGGTGAAACGTTTTATTGCTGCGACCAACGTGAACGATACCGTGCCACGTTTCCTGCACGACGGTCAGTGG  >  NC_000913/4114‑4597                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                       gacaaTTCTGACCGCGACCTCCGGGGATAACGGAGTGGCAGTGGCTTATGCTTTCTACGGTTTACCGAATGTGAAAGTGGTTATCCGCTAGCCACGAGGCAAAAACAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTca                                                                                                                                                                                                                                            <  2:547929/246‑2 (MQ=255) gACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTTACCGAATGTGAAAGTGGTTATCCTCTATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTcag                                                                                                                                                                                                                                           <  1:19743/251‑3 (MQ=255) gACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTTACCGAATGTGAAAGTGGTTATCCTCTATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTca                                                                                                                                                                                                                                            <  1:491467/250‑2 (MQ=255)     aTTCTGACCGCGACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTTACCGAATGTGAAAGTGGTTATCCTCTATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTaa                                                                                                                                                                                                                                       <  1:259999/251‑1 (MQ=255)          gACCGCGACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTTACCGAATGTGAAAGTGGTTATCCTCTATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTTAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGa                                                                                                                                                                                                                                  <  1:467670/251‑1 (MQ=255)             cgcgACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTTACCGAATGTGAAAGTGGTTATCCTCTATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTa                                                                                                                                                                                                                               <  1:190235/251‑1 (MQ=255)                       ggTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTTACCGAATGTGAAAGTGGTTATCCTCTATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCg                                                                                                                                                                                                                     >  1:316306/1‑251 (MQ=255)                        gtgcgaCCGGAGAGACAGTGGCTAATGCTTTCTACGGTTTACCGAATGGGAAAGAGGTTATCCTCTATACACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGTCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCgg                                                                                                                                                                                                                    <  2:120943/246‑2 (MQ=255)                        gTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTTACCGAATGTGAAAGTGGTTATCCTCTATCCACGAGGACATATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCTGTTGATGATGA‑AGGGCTGACAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGt                                                                                                                                                                                                                    >  2:364255/1‑251 (MQ=255)                                ggAGCCGCAGGGGCTCATGCTTTCTACCGTTTACCGAAGGTGAAAGTGGTTATCCTGTATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGGACATTGGGCGGCAATATCGAAACTTTGGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGTTGAGGCAGGCGTTTGATGATGA‑AGAGCTTAAAGTGTCGCTGGGGTTAATCACAGCTAACTCGCTTAACATGAGTCGTTTGCTGgc                                                                                                                                                                                                            <  2:655131/251‑1 (MQ=255)                                    cGGCAGAGGCTCATGCTATCTACGGTTTACCGAATGTGCAAGTGGTTATACTCTATCCACGAGGCAAAATCAGTCCCCTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAAATGTTGCCATAGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCAGAAAGGGGCGCTGTGGTTAAACTCAACTAACTCGATTAACATCAGCTGTTTGCTGgcgac                                                                                                                                                                                                         >  2:560991/1‑248 (MQ=255)                                     ggCAGTGGCTCATTATTTCTCCGGATTACCGAATGTGAAAGGGGTTATCCTCTATCCAAGCGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCAAAACTGTTGCAATAGACGGCGATTTAGATGACTGTCAGGCGCTGGAGAAGCAGGAGTTTAATGATGA‑CGAGCTGAAAGGGGCGCTGGGGGTAAAATCCGCTAACTCGTTTAACAACAGCAGTTTGCGGGCGCAgg                                                                                                                                                                                                       >  2:468910/1‑250 (MQ=255)                                        aGTGGCTCATGCTTTCTACGGTTTACCGAATGTGAAAGTGGTTATCCTCTATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGAttt                                                                                                                                                                                                    <  1:325572/251‑1 (MQ=255)                                                        tACGGTTTACCGCATGTGAACGAGGTTATCCTCTATCCACGCGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGCTTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGCGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGCTTTGCTCCTACTTTGAAGc                                                                                                                                                                                    >  2:503454/1‑251 (MQ=255)                                                         aCGGTTTACCGAATGTGAAAGTGGTTATCCTCTATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAg                                                                                                                                                                                     <  1:262686/249‑1 (MQ=255)                                                         aCGGTTTACAAAAAGTGAAAGTGTTTAGCCTGTATCCACGAGGCAAATTCAATCCACTGCAAGCAATACTGTTCTGTACATTGGGCGGCAATATCGAACGTGATGCCACCGACGTCGATTTCGATGCCTGTCAGGCGCTGGTGGTGCAGGCGTTTGGTGATGA‑AGAACGAAAAGTGGCGCTAGGGTTCAACTCGGATAACTCGATTAAAATCAGCCTATCGCTGGCGCAGATGTGCTAAGACTTTGAAGcc                                                                                                                                                                                   >  2:268457/1‑250 (MQ=255)                                                            gTTTACCGAATGTGAAAGTGGTTATCCTCTATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGtt                                                                                                                                                                                >  1:4511/1‑251 (MQ=255)                                                             tttACCGAATGTGAAAGTGGTTATCCGCTATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGtt                                                                                                                                                                                <  2:456319/250‑1 (MQ=255)                                                                   gAATGTGAAAGTGGTAATACTCTATCCACGAGGCAAAATCAGTCCAATGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGCAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGAAGGCGTTTGCTGATGA‑AGACCTGAAAGAGGCGCTAGGGTTAAACACGGCTAAATCGATTAAAAACAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCTGTTgcgc                                                                                                                                                                            >  2:593325/1‑248 (MQ=255)                                                                      tgtgAAAGTGGTTATCCTCTATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCCGCTAACTCGATTACCATCAGCCGTTTGCTGGCGCGGATTTGCTACTACTTTGAAGcag                                                                                                                                                                                  >  2:418725/1‑237 (MQ=255)                                                                      tgtgAAAGTGGTTATCCTCTATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGcag                                                                                                                                                                                  <  1:418725/239‑3 (MQ=255)                                                                                            ccACGAGGCAAAATCATGCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCAGCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACAGCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCt                                                                                                                                                <  1:114863/251‑1 (MQ=255)                                                                                                    cAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAACGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGtct                                                                                                                                        >  2:279846/1‑251 (MQ=255)                                                                                                    cAAAATCAGGCCCCTGAAAGAAAAACTGGTCTGTACATTGGGAGGCAATATCGAAACTGTTGCCACCGACGGCGATGTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAACTGAAAGTGTCACTAGGGTTAAAATCGGCTAACTCGATGAACATCAGCCGTTTGCTGGCGCAGATTTGCTCCTACTTGGAAGCTGTTGCGCGGCTGCCGCAGGAGACGCGCAACCAGCTGGTTGtct                                                                                                                                        <  2:51029/251‑1 (MQ=255)                                                                                                        aTCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGt                                                                                                                                    <  1:537250/251‑1 (MQ=255)                                                                                                           aGTCCACTGCAAGAAAAACTGTTCTGGACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGACAGTGGCGCTGGGGTTAAACTCAGCTCACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCATTTGCGCAGCGGCCGCAGAAAGCGCGCACCCACCTGGTTGtcgcgtgtcc                                                                                                                                 >  2:581461/1‑243 (MQ=255)                                                                                                                    cAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGaa                                                                                                                         >  1:539105/1‑250 (MQ=255)                                                                                                                    cAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCAGCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGaaa                                                                                                                        <  1:286116/251‑1 (MQ=255)                                                                                                                     aaGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAAATCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAAc                                                                                                                       >  1:379782/1‑251 (MQ=255)                                                                                                                      aGAACAACAGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCTATCAACGGCGCTTTCGATGTCTGACAAGCCCTGGTGAAGCAGACGCTTGATGATGC‑GGAGCTGAAAGGGGCGGTGGGGGTAAACTCAGCTAACTCGATTACCATCAGCGGTGTGCTGGCGCAGGTTGGATACTACTTTGAAGCAGTTGCGTCGCTCCCGCAGGAAGCGCACAACCAGCTGGATGTACCGGTGcccacaggaaact                                                                                                                      >  2:13155/1‑240 (MQ=255)                                                                                                                      aGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACt                                                                                                                      >  2:536611/1‑251 (MQ=255)                                                                                                                             cTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGCCGGCGATTTCGATGCCTGTAAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCgg                                                                                                               >  2:34376/1‑250 (MQ=255)                                                                                                                                     tACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGTGTTGAcgg                                                                                                       >  2:372779/1‑251 (MQ=255)                                                                                                                                      aCATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGtc                                                                                                                                         >  1:137070/1‑216 (MQ=255)                                                                                                                                      aCATTGGGCGGCAATATCGAAACTGTGGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGtc                                                                                                                                         <  2:137070/216‑1 (MQ=255)                                                                                                                                       cATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGAcggcg                                                                                                     >  1:397026/1‑251 (MQ=255)                                                                                                                                          tGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGt                                                                                                                                    >  1:487639/1‑217 (MQ=255)                                                                                                                                           gggCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGt                                                                                                  <  1:313348/250‑1 (MQ=255)                                                                                                                                           gggCGGCAATAACGAACCTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGTAGCAGGCGTTTGATGATGA‑AGAACTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACAGCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTc                                                                                                 <  2:122206/251‑1 (MQ=255)                                                                                                                                               ggcAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGACTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCACACGGAAACTTCGGCGATTTGCCGGCgggggtgtt                                                                                             >  2:319302/1‑245 (MQ=255)                                                                                                                                                 cAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCCGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGCAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCTGt                                                                                                                                                                                 >  1:645176/1‑165 (MQ=255)                                                                                                                                                 cAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATAAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGtacgctgg                                                                                           >  1:380436/1‑244 (MQ=255)                                                                                                                                                    tatCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTggggg                                                                                        >  1:393309/1‑248 (MQ=255)                                                                                                                                                              gTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCg                                                                               <  1:104744/250‑1 (MQ=255)                                                                                                                                                                 gCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTCt                                                                           <  1:536611/251‑1 (MQ=255)                                                                                                                                                                  ccatggccggagagttcggggCCTGTCAGGCGCTGTTGAAGCGTGCGTTTGAGGA‑GAGAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACGCGAGTAACTTCAGCCGTTGGCTGGCGCAGATTTGCTACTACTTTGAAGCTGTTGCGCAGCTGCCTCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGGCACTCGGTCTg                                                                          <  2:588811/231‑1 (MQ=255)                                                                                                                                                                  ccATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGGAGGAGTTTGATGATGA‑AGAGCTGAAAGTGGCGATGGGGTTAAACTCAGCTAACTCGATTAAAATCAGCCGTTTGCTGGCGCAGAAGTGCTACTACTTTAAAGAAGTTGCGAAGATGCAGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAACATTCGGCGATTTGACGGCGGGActgtgggcgggagggacaggggtg                                                                          >  2:659677/1‑230 (MQ=255)                                                                                                                                                                  ccATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCCGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCCGGtgggctgggggaggaccgtggatgg                                                                          >  2:514265/1‑227 (MQ=255)                                                                                                                                                                  ccATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCg                                                                                                                            <  2:452215/201‑1 (MQ=255)                                                                                                                                                                  ccATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCg                                                                                                                            >  1:452215/1‑201 (MQ=255)                                                                                                                                                                    aTCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACAGCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTCTGcc                                                                        <  1:245566/251‑1 (MQ=255)                                                                                                                                                                    aTCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAACTGAAAGTGGCGCTAGGGTTAAACTCGGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCTGTTGCGCAGCTGCCGCAGGAGACGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGGCGGCTGGCggcgacaacccgcgcggc                                                                        >  2:312113/1‑234 (MQ=255)                                                                                                                                                                       gACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGtt                                                                                                                                                                                >  1:146957/1‑144 (MQ=255)                                                                                                                                                                       gACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGtt                                                                                                                                                                                <  2:146957/144‑1 (MQ=255)                                                                                                                                                                        aCGGCGATTTCGATTCCTGTCAGGCGCTGGTGAAGCAGGCGGTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTCTGCCGGTg                                                                    <  1:330595/251‑1 (MQ=255)                                                                                                                                                                         gggCGATTTGAATGCCTGTCAGGCGTGGGTTAAGCAGGCGTTGGAGGAGGA‑AGAGCGTAAAGTGGCGCTGGGTTTAAACTCAGCTAACGCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTCTGCCGGTGa                                                                   <  2:314301/250‑1 (MQ=255)                                                                                                                                                                         cGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTg                                                                                                                                   >  1:401063/1‑187 (MQ=255)                                                                                                                                                                         cGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTGGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTg                                                                                                                                   <  2:401063/187‑1 (MQ=255)                                                                                                                                                                          ggCGATTTCGATGCCTGTCAGGCGCTGGCGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGTGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAAACAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTCTGCCGGTGaa                                                                  <  1:259989/251‑1 (MQ=255)                                                                                                                                                                            cGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGGATGCTGGCggaggtccaccgtctgccggggaga                                                                 >  1:232509/1‑226 (MQ=255)                                                                                                                                                                              aTTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTCTGCCGGTGAAACGt                                                              <  1:617040/251‑1 (MQ=255)                                                                                                                                                                                 tCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTCTGCCGGTGAAACGttt                                                            <  1:431307/250‑1 (MQ=255)                                                                                                                                                                                       ccTGGCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACGCAGCTAAGCCGATGAACATCAGCCGGTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCTATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTCTGCCGGTGAAACGTTTTATtgct                                                     <  2:397026/251‑1 (MQ=255)                                                                                                                                                                                            cAGGCGCTGGTGAAGCAGGCGGTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTCTGCCGGTGAAACGTTTTATTGCTGCGAc                                                <  2:316306/251‑1 (MQ=255)                                                                                                                                                                                                      ggagaccggCGTTTGATGATGA‑AGAGCTGAAAGGGGCGCTGGGGTAAAACTAAGCTAACTCGATTAACAACAGCGGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTCTGCCGGTGAAACGTTTTATTGCTGCGACCAGCGTGAAc                                      <  2:4511/244‑1 (MQ=255)                                                                                                                                                                                                            aGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGTTTAAACTCAGCTAACTCGCTTAACATCGGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCACCCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGGCTGCTGGCGAAGTCACTCGGTCTGCCGGTGAAACGTTTTa                                                          <  2:542452/225‑1 (MQ=255)                                                                                                                                                                                                             ggCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTCTGCCGGTGAAACGTTTTATTGCTGCGACCAACGTGAACGATACCg                               <  1:268960/251‑1 (MQ=255)                                                                                                                                                                                                                   tgatgatga‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTCTGCCGGTGAAACGTTTTTTTTCTGCGACCAACGTGAACGATACCGTGcccc                         >  1:193019/1‑249 (MQ=255)                                                                                                                                                                                                                      ggatga‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACAGCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTct                                                                                                <  2:19485/176‑1 (MQ=255)                                                                                                                                                                                                                       gatga‑AGAGCTGAAAGTGGCGCTGGGGTTAAACCCAGCTAACTCGATTAACAGCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCGGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTCTGCCGGTGAAACGTTTTATTGCTGCGACCAACGTGAACGATACCGTGCCACGttt                     <  2:193019/251‑1 (MQ=255)                                                                                                                                                                                                                             agaaCTGAAAGTGGCGCTAGGGTTAAACTCGGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCTGTTGCGCAGCTGCCGCAGGAGACGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTCTGCCGGTGAAACGTTTTATTGCTGCGACCAACGTGAACGATACCGTGCCACGTTTCCTGc                <  1:312113/251‑1 (MQ=255)                                                                                                                                                                                                                              gagcTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTgg                                                                                                                                              <  1:156575/121‑1 (MQ=255)                                                                                                                                                                                                                              gagcTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTgg                                                                                                                                              >  2:156575/4‑124 (MQ=255)                                                                                                                                                                                                                              gagcTGAAAGTGGCGCGGGGGTTAACCTCAGCGAACTCGATGAACAGCAGCCGTTTGCTGGCGCAGATTTGGTACCACTTTGAAGCAGTTGCGCAGCTGCCGTAGGAAGCGCGCAAACAGATGGTTGTCTCGGTGCCAAGCGGAAACGTCGGCGATTTTACGGCGGGTCTGCTGGCGAAGTCACTCGGTCTGCCGGTGAAACGTTTTATTGCTGCGACCAACGTGAACGATACCGTGCCACGTTTCCTGCa               <  2:384063/248‑1 (MQ=255)                                                                                                                                                                                                                                    aaaGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTCTGCCGGTGAAACGTTTTATTGCTGCGACCAACGTGAACGATACCGTGCCACGTTTCCTGCACGACgg         <  1:36973/251‑1 (MQ=255)                                                                                                                                                                                                                                        tGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTAACTCGGTCTGCCGGTGGAAAGTTTTTTTGTTGCGACCCACGTGAAAGATAACGTGCCCCGTTTTCTGcccgccgggccg     >  1:428663/1‑240 (MQ=255)                                                                                                                                                                                                                                         ggCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCTGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGGCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCATGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTCTGCCGGTGAAACGTTTTATTGCTGCGACCAACGTGAACGATACCGTGCCACGTTTCCTGCACGACGGTCAg     <  1:222900/250‑1 (MQ=255)                                                                                                                                                                                                                                           cgcTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTCTGCCGGTGAAACGTTTTATTGCTGCGACCAACGTGAACGATACCGTGCCACGTTTCCTGCACGACGGTCAGTgg  <  1:610749/251‑1 (MQ=255)                                                                                                                                                                                                                                           cgcTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGGCTGCTGGCGAAGGCACTCGGGCTGCCGGTGGAACGGTTTATTGCTGGGACCAAAGTGAAAGACACCGGGGCACGCTTCCTGCCCGACGGTCAGTgg  >  1:104366/1‑251 (MQ=255)                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                       GACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTTACCGAATGTGAAAGTGGTTATCCTCTATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAACTGAAAGTGGCGCTAGGGTTAAACTCGGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCTGTTGCGCAGCTGCCGCAGGAGACGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTCTGCCGGTGAAACGTTTTATTGCTGCGACCAACGTGAACGATACCGTGCCACGTTTCCTGCACGACGGTCAGTGG  >  NC_000913/4114‑4597

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