Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,269,856	G→T	8.3%	intergenic (‑46/‑2)	ldrB ← / → rdlB	toxic polypeptide, small/sRNA antisense regulator affects LdrB translation; proposed addiction module in LDR‑B repeat, with toxic peptide LdrB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,269,856	0	G	T	8.3%	222.3 / 6.3	75	intergenic (‑46/‑2)	ldrB/rdlB	toxic polypeptide, small/sRNA antisense regulator affects LdrB translation; proposed addiction module in LDR‑B repeat, with toxic peptide LdrB
Reads supporting (aligned to +/- strand):  ref base G (29/38);  new base T (3/4);  total (32/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

GAGGCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGTGATAAAGACAAAAATGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTA‑CCTCTCAACGTGCGGGGGTTT‑TCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGGTGTGGATACTACGTCTCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATAAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTCGTTTCATGAGTGCTCTGGATGATGCTTCCAGCTCGGGTTGCC  >  NC_000913/1269572‑1270128                                                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                                                    gAGGCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGTGATAAAGACAAAAATGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAgtg                                                                                                                                                                                                                                                                                   >  1:297021/1‑287 (MQ=35) gAGGCTTCCAGCTCGGGTTGCAAATATTTACTTGTGGAAGTGATAAAGCCAAAAATGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTG‑Tc                                                                                                                                                                                                                                                                                <  1:809497/289‑1 (MQ=25)   ggCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGTGATAAAGACAAAAATGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGCAGTG‑TCTg                                                                                                                                                                                                                                                                              >  1:28405/1‑289 (MQ=21)    gCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGTGATAAAGACAAAAATGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTG‑TCTgg                                                                                                                                                                                                                                                                             >  1:279554/1‑289 (MQ=34)        ccAGCTCGGGTTGCCAATATTTACTTGTGGAAGTGATAAAGACAAAAATGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTG‑TCTGGtt                                                                                                                                                                                                                                                                           >  2:67439/1‑287 (MQ=34)        ccAGCTCGGGTTGCCAATATTTACTTGTGGAAGTGATAAAGACAAAAATGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTG‑TCTGGTTTc                                                                                                                                                                                                                                                                         <  2:679471/289‑1 (MQ=34)        ccAGCTCGGGTTGCCAATATTTACTTGTGGAAGTGATAAAGACAAAAATGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTG‑TCTGGTTTc                                                                                                                                                                                                                                                                         >  1:202500/1‑289 (MQ=34)                               ctgttgGAAGTGATAAAGACAAAAAGGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCtgttgt                                                                                                                                                                                                                                                  <  2:417053/285‑1 (MQ=25)                               cTTGTGGAAGTGATAAAGACAAAAATGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCtgttgt                                                                                                                                                                                                                                                  >  2:100509/1‑289 (MQ=35)                                   tgGAAGTGATAAAGACAAAAATGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAATTGTTCTGG‑TTCAAGATTAGCCCCCGTTCtgttgt                                                                                                                                                                                                                                                  >  2:820292/1‑285 (MQ=21)                                   tgGAAGTGATAAAGACAAAAATGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAATTGTTCTGG‑TTCAAGATTAGCCCCCGTTCtgttgt                                                                                                                                                                                                                                                  <  1:820292/285‑1 (MQ=21)                                   tgGAAGTGATAAAGACAAAAATGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAgg                                                                                                                                                                                                                                              >  1:629066/1‑289 (MQ=35)                                   tgGAAGTGATAAAGACAAAAATGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTCTCACCTTACTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGCAATAATTCCCGCTAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGACATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTG‑TCTGGTTTCATGATTAGACCCCGTTCTGTTGTCAgg                                                                                                                                                                                                                                              >  1:276512/1‑289 (MQ=255)                                       aGTGATAAAGACAAAAATGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGtttt                                                                                                                                                                                                                                          <  1:100509/289‑1 (MQ=35)                                             aaaGACAAAAATGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTA‑Cctc                                                                                                                                                                                                                                    >  1:137800/1‑288 (MQ=34)                                                     aaaTGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATTACAACATTTTGCGGCTATTCTTGAATTG‑TCTGGTTTCAAGTTTTGCCCCCGTTCTGTTGTCAGGTTTTA‑CTTCCCAACGttg                                                                                                                                                                                                                           >  2:499193/1‑287 (MQ=255)                                                                   tGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTA‑CCTCTCAACGTGCGGGGGTTT‑Tctct                                                                                                                                                                                                             <  1:151473/288‑1 (MQ=34)                                                                       aCCCGTGGGGGCGGTGTCGGGCGCATATTGCCATCACGGAAGCCTGACGCCCGTTTTCACCTTACTTGCGGTTACGCCACCAGCCGACAATCGCTGCGGTGATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCTGCTTTTTTTTAATGACAACATATCGCGGCTATTCTTGAGCTG‑TCTGGTTTCAGGATTAGACGTCTATAAGTTGCAAGGTTTTA‑CCTCTCAACGTGCGGGGGTTT‑TCTCTTTCCa                                                                                                                                                                                                        <  2:540610/289‑1 (MQ=11)                                                                         cccTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTGTTCAGGGTTTA‑ACTCGCAAGGTGTGGGGGTTTATCT‑TTTCCAg                                                                                                                                                                                                       >  2:552097/1‑288 (MQ=11)                                                                             gCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATCTAGCCGCCTTGTTGTAATGACAACATTTTTCGGCTATTCCTTAAGTG‑TCTGGTTTCAAGGTTAGCCCCCGTTCTGTTGTCAGGTTTTA‑CCTATCAACGTTGCGGGGTCT‑TCTTTTTACAggcacc                                                                                                                                                                                                  >  2:504979/1‑284 (MQ=11)                                                                                   ggTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTTCGGCTTTTCTTGAAGTT‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTA‑CCTCCCAACGTGCGGGGGTTT‑TCCCTTTCCAGCACCAAATGc                                                                                                                                                                                             >  2:382357/1‑288 (MQ=255)                                                                                      ttCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTA‑CCTCTCAACGTGCGGGGGTTT‑TCTCTTTCCAGCAACCAATGccac                                                                                                                                                                                          <  2:286860/288‑1 (MQ=34)                                                                                       tCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTA‑CCTCTCAACGTGCGGGGGTTT‑TCTCTTTCCAGCAACCAATGccacca                                                                                                                                                                                        <  1:67439/289‑1 (MQ=34)                                                                                              cATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTA‑CCTCTCAACGTGCGGGGGTTT‑TCTCTTTCCAGCAACCAATGCCACCAGGGATaa                                                                                                                                                                                 >  2:334201/1‑289 (MQ=34)                                                                                               atatTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTA‑CCTCTCAACGTGCGGGGGTTT‑TCTCTTTCCAGCAACCAATGCCACCAGGGATaaa                                                                                                                                                                                <  1:100433/289‑1 (MQ=34)                                                                                                tatTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTA‑CCTCTCAACGTGCGGGGGTTT‑TCTCTTTCCAGCAACCAAt                                                                                                                                                                                               >  1:515486/1‑273 (MQ=34)                                                                                                tatTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTA‑CCTCTCAACGTGCGGGGGTTT‑TCTCTTTCCAGCAACCAAt                                                                                                                                                                                               <  2:515486/273‑1 (MQ=34)                                                                                                 atTGCCATCACGGCAGCCTGCCGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTA‑CCTCTCAACGTGCGGGGGTTT‑TCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGc                                                                                                                                                                              <  1:504979/289‑1 (MQ=21)                                                                                                         cACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGACGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTA‑CCTCTCAACGTGCGGGGGTTT‑TCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCa                                                                                                                                                                       >  1:462916/1‑288 (MQ=25)                                                                                                         cACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTA‑CCTCTCAACGTGCGGGGGTTT‑TCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCa                                                                                                                                                                       >  1:211110/1‑288 (MQ=35)                                                                                                           cGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTA‑CCTCTCAACGTGCGGGGGTTT‑TCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACa                                                                                                                                                                    >  1:504245/1‑289 (MQ=37)                                                                                                               aGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTA‑CCTCTCAACGTGCGGGGGTTT‑TCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTgc                                                                                                                                                                >  1:657108/1‑289 (MQ=37)                                                                                                                 ccTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTA‑CCTCTCAACGTGCGTGGGTTT‑TCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTgcg                                                                                                                                                               >  1:601978/1‑288 (MQ=37)                                                                                                                  cTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTA‑CCTCTCAACGTGCGGGGGTTT‑TCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTgcg                                                                                                                                                               <  1:847469/287‑1 (MQ=37)                                                                                                                      cGCCCGTTTTCACCTTACTTCCGTGTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTA‑CCTCTCAACGTGCGGGGGTTT‑TCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCTCAc                                                                                                                                                         <  1:552097/289‑1 (MQ=37)                                                                                                                      cGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGCGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTA‑CCTCTCAACGTGCGGGGGTTT‑TCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCTCAc                                                                                                                                                         <  1:240608/289‑1 (MQ=37)                                                                                                                         ccGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTA‑CCTCTCAACGTGCGGGGGTTT‑TCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCgg                                                                                                                                                      <  1:499193/289‑1 (MQ=37)                                                                                                                                    ttACTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTA‑CCTCTCAACGTGCGGGGGTTT‑TCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGc                                                                                                                                           <  2:297021/289‑1 (MQ=39)                                                                                                                                    ttACTTCCGGTTAATCCACCAGCAGACAAGCGCTGCGGTAATAAGTCCCGCCAGGATCGGTGCGGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTCTGA‑CCTCTCAACGTGCGACGCTCT‑GCATATTGCAGCCACCGTCGCAACTGTAGATGAAGCCGCCGCAACATTGCGCCTCACCGGATAATTCCGGc                                                                                                                                           <  2:206889/289‑1 (MQ=11)                                                                                                                                      ccTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTA‑CCTCTCAACGTGCGGGGGTTT‑TCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCtt                                                                                                                                         <  2:504245/288‑1 (MQ=39)                                                                                                                                      aCTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTA‑CCTCTCAACGTGCGGGGGTTT‑TCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCtt                                                                                                                                         <  2:202500/289‑1 (MQ=39)                                                                                                                                         tCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTA‑CCTCTCAACGTGCGGGGGTTT‑TCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGgt                                                                                                                                      <  2:137800/289‑1 (MQ=39)                                                                                                                                         tCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTA‑CCTCTCAACGTGCGGGGGTTT‑TCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGgt                                                                                                                                      <  2:279554/289‑1 (MQ=39)                                                                                                                                         tCCGATTACGACACCAGCCGACAATCTCTGCCGTAATAAGTCCCGCCAGGAACGGTGCCGCCAGGTCGTGCCAGAAAATCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTTAATTG‑TCTTGGGTCAAGAGTAGCCCCCGTTCTGTTGTCATGTTTTA‑CCTCTCAACGTGCGGGGGTTT‑TCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCTCCCCGGATAATTCCGGCTTGgt                                                                                                                                      <  2:40701/289‑1 (MQ=11)                                                                                                                                                 ccccaccaTCCTACAATGGGTGCGGTAATAATTCCCGCCAGGAGCGGGGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGCATTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTA‑CCTCTCAACGTGCGGGGGTTT‑TCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGGTGTGGAtac                                                                                                                              <  2:276512/287‑1 (MQ=255)                                                                                                                                                           cGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTA‑CCTCTCAACGTGCGGGGGTTT‑TCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGGTGTGGATACTACGTCTCaa                                                                                                                    <  2:318145/289‑1 (MQ=39)                                                                                                                                                                                  cccGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTA‑CCTCTCAACGTGCGGGGGTTT‑TCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGGTGTGGATACTACGTCTCAATTCATCTTCACTTCATCCCTGaa                                                                                             >  1:806371/1‑289 (MQ=39)                                                                                                                                                                                    cGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTA‑CCTCTCAACGTg                                                                                                                                                                                                                            <  2:4834/161‑1 (MQ=35)                                                                                                                                                                                    cGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTA‑CCTCTCAACGTg                                                                                                                                                                                                                            >  1:4834/1‑161 (MQ=35)                                                                                                                                                                                       cAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTA‑CCTCTCAACGTGCGGGGGTTT‑TCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGGTGTGGATACTACGTCTCAATTCATCTTCACTTCATCCCTGAAATGt                                                                                         <  1:382357/288‑1 (MQ=39)                                                                                                                                                                                       cAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTA‑CCTCTCAACGTGCGGGGGTTT‑TCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGGTGTGGATACTACGTCTCAATTCATCTTCACTCCATCCCTGAAATGt                                                                                         <  2:657108/288‑1 (MQ=39)                                                                                                                                                                                            tCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTA‑CCTCTCAACGTGCGGGGGTTT‑TCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGGTGTGGATACTACGTCTCAATTCATCTTCACTTCATCCCTGAAATGTTTGCa                                                                                    <  1:764352/288‑1 (MQ=39)                                                                                                                                                                                            tCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTA‑CCTCTCAACGTGCGGGGGTTT‑TCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGGTGTGGATACTACGTCTCAATTCATCTTCACTTCATCCCTGAAATGTTTGCa                                                                                    <  2:806371/288‑1 (MQ=39)                                                                                                                                                                                            tCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTA‑CCTCTCAACGTGCGGGGGTTT‑TCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGGTGTGGATACTACGTCTCAATTCATCTTCACTTCATCCCTGAAATGTTTGCa                                                                                    <  1:22935/288‑1 (MQ=39)                                                                                                                                                                                                   cgccAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTA‑CCTCTCAACGTGCGGGGGTTT‑TCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGGTGTGGATACTACGTCTCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATAAaga                                                                             <  2:445159/288‑1 (MQ=39)                                                                                                                                                                                                     ccAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTA‑CCTCTCAACGTGCGGGGGTTT‑TCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGgt                                                                                                                                      <  1:697272/229‑1 (MQ=39)                                                                                                                                                                                                     ccAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTA‑CCTCTCAACGTGCGGGGGTTT‑TCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGgt                                                                                                                                      >  2:697272/1‑229 (MQ=39)                                                                                                                                                                                                          tCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTA‑CCTCTCAACGTGCGGGGGTTT‑TCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCACCATTGCGCCTCACCGGATAATTCCGCCTTGGTGTGGATACTACGTCTCAATTCATCTTCACTTCATCCCTGAAACGTTTCCAATAAAGAGTACAtt                                                                      >  1:75644/1‑288 (MQ=255)                                                                                                                                                                                                                    aaGTAATGGCAGACTGCGCAAAGGTAATATAGCCGCCTTGTTGTAATGACAACATTTTGCGCGTAGTCTTGAAGTG‑TGTGGATGCACGATCAGACCCCGTACTGGTGTCACGTTTTA‑GCTCTCAACGAGCAGTGTTTC‑TTTCTCGGCAGCAACAAATGACCCCAGGGATAAGGCCCACGAAACATTGCGCGTCACACGATAATTCCGCCTTGGATAGCATAGTCCGCGTCAGCTGTGTATAAGTGCCTGCCTGAAATGTGTGCAATAAAGAGTACATTCCGGCgtt                                                              <  2:194669/286‑1 (MQ=2)                                                                                                                                                                                                                         aTGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTA‑CCTCTCAACGTGCGGGGGTTT‑TCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGGTGTGGATACTACGTCTCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATAAAGAGTACATTCCGGc                                                                 >  2:265696/1‑278 (MQ=38)                                                                                                                                                                                                                         aTGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTA‑CCTCTCAACGTGCGGGGGTTT‑TCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGGTGTGGATACTACGTCTCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATAAAGAGTACATTCCGGc                                                                 <  1:265696/278‑1 (MQ=38)                                                                                                                                                                                                                             cAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTA‑CCTCTCAACGTGCGGGGGTTT‑TCTCTTTCCAGCAACCaa                                                                                                                                                                                                >  2:622749/1‑147 (MQ=32)                                                                                                                                                                                                                             cAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTA‑CCTCTCAACGTGCGGGGGTTT‑TCTCTTTCCAGCAACCaa                                                                                                                                                                                                <  1:622749/147‑1 (MQ=32)                                                                                                                                                                                                                                aCTGCGCGAGCGTCCTATAGCAGCGTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTA‑CCTCTCAACGTGCGGGGGTTT‑TCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGGTGTGGATACTACGTCTCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATAAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCa                                               <  2:629066/289‑1 (MQ=255)                                                                                                                                                                                                                                   gcgcGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAATTGTTCTGG‑TTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTA‑CCTCTCAACGTGCGGGGGTTT‑TCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGGTGTGGATACTACGTCTCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATAAAGAGTACATTCCGGCTTTTCAACAg                                                       >  1:340853/1‑278 (MQ=12)                                                                                                                                                                                                                                   gcgcGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAATTGTTCTGG‑TTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTA‑CCTCTCAACGTGCGGGGGTTT‑TCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGGTGTGGATACTACGTCTCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATAAAGAGTACATTCCGGCTTTTCAACAg                                                       <  2:340853/278‑1 (MQ=12)                                                                                                                                                                                                                                             cATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTA‑CCTCTCAACGTGCGGGGGTTT‑TCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGGTGTGGATACTACGTCGCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATAAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTCGTTTCATGa                                   <  2:211110/288‑1 (MQ=38)                                                                                                                                                                                                                                                       ccTTGTTGTAATGACAACATTTTGCGGCTATTCTTTAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTA‑CCTCTCAACGTGCGGGGGTTT‑TCTCTTTCCAGCAACCACTGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGGTGTGGATACTACGTCTCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATAAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTCGTTTCCTGAGTGCCCTGgat                        >  1:552543/1‑289 (MQ=255)                                                                                                                                                                                                                                                       ccTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTTTGTTGTCAGGTTTTA‑CCTCTCAACGTGCGGGGGTTT‑TCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGGTGTGGATACTACGTCTCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATAAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTCGTTTCATGAGTGCTCTGgat                        >  1:228467/1‑289 (MQ=255)                                                                                                                                                                                                                                                       ccTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTA‑CCTCTCAACGTGCGGGGGTTT‑TCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGGTGTGGATACTACGTCTCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATAAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTCGTTTCATGAGTGCTCTGgat                        >  1:252970/1‑289 (MQ=255)                                                                                                                                                                                                                                                              gtAATGACAACATTTTGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTA‑CCTCTCAACGTGCGGGGGTTT‑TCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGGTGTGGATACTACGTCTCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATAAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTCGTTTCATGAGTGCTCTGGATGATGCt                  >  1:492182/1‑288 (MQ=255)                                                                                                                                                                                                                                                                      aacaTTTTGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTA‑CCTCTCAACGTGCGGGGGTTT‑TCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGGTGTGGATACTACGTCTCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATAAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTCGTTTCATGAGTGCTCTGGATGATGCTTCCAGCTCg         >  2:183235/1‑289 (MQ=255)                                                                                                                                                                                                                                                                       acaTTTTGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTA‑CCTCTCAACGTGCGGGGGTTT‑TCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGGTGTGGATACTACGTCTCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATAAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTCGTTTCATGAGTGCTCTGGATGATGCTTCCAGCTCgg        <  2:601978/289‑1 (MQ=255)                                                                                                                                                                                                                                                                        caTTTTGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTGTACCCTCTCAACGTGCGGGGGTTT‑TCTCTTTCCAGCAACCAATGCCACCAGGGATAAATCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGGTGTGGATACTACGTCTCAATTCATCTTCACTTCATCCCTGAAATGATTGCAATAAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTCGTTTCATGAGTGCTCTGGATGATGCTTCCAGCTCgg        <  1:858277/289‑1 (MQ=255)                                                                                                                                                                                                                                                                            tttcGGCTATTCTTAAAGTG‑TCTGGTTTCAAGATTACCCCCCGTTCTGTTGTCAGGTTTTA‑CCTCTCAACGTGCGGGGGTTT‑TCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGGTGTGGATACTACGTCTCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATAAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTCGTTTCATGAGTGCTCTGGATGATGCTTCCAGCTCGGGTTGc   <  2:75644/286‑1 (MQ=255)                                                                                                                                                                                                                                                                             tGCGGCTATGCTTGAAGTG‑TCTGGTTTCAAGTTTATCCCCCATTCTGTAGTCAGGTTTTA‑CCTCTCAACGTGCGGGGGTTT‑TCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGGTGTGGATACTACGTCTCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATAAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTCGTTTCATGAGTGCTCTGGATGATGCTTCCAGCTCGGGTTGcc  <  2:552543/289‑1 (MQ=255)                                                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                                                    GAGGCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGTGATAAAGACAAAAATGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTA‑CCTCTCAACGTGCGGGGGTTT‑TCTCTTTCCAGCAACCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGGTGTGGATACTACGTCTCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATAAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTCGTTTCATGAGTGCTCTGGATGATGCTTCCAGCTCGGGTTGCC  >  NC_000913/1269572‑1270128

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