Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,880,375	G→T	5.6%	L283M (CTG→ATG)	recF ←	gap repair protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,880,375	0	G	T	5.6%	216.9 / 2.5	71	L283M (CTG→ATG)	recF	gap repair protein
Reads supporting (aligned to +/- strand):  ref base G (34/33);  new base T (2/2);  total (36/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.85e-01

CTTTGATACTGGAGGAGTCATAAGAATTCGACATCAACGTTTCTCGCTCATTTATACTTGGGTTAATCCGTTATTTTACCCTTTTCCACGGTAAACATCTTCGAATTTTCGTCCGACATGTCTATAACGTGTTCAGCACTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCGATCGCGTTCAAAATTACGTTCCAGCACCTCAGCATATTCTGTCTCTTTCTCCCAGCCGCGCTGGAAAGAGAAAGTCAGAGAAAACTCAGGGAGAAATTGCTTACAGGTATCAGCCATATCGGCCGCGATACCGGCGCTATACTCCGCGCGCCAGGTGCTGATTTGCTC  >  NC_000913/3880086‑3880666                                                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                                                                    ttttGATATTGTAGGTGTCATAAGAATTGGACATCATCGTTTTTCGCTCTTTTATACTTGGGTTAATCCGTTATTTTACCCTTTTCCACGGTAAACATCTTCGAATTTTCGTCCGACATGTCTATAACGTGTTCAGCACTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGgcgc                                                                                                                                                                                                                                                                                          <  1:573370/300‑1 (MQ=255)    tgttACTGGAGGGGTCATAAGAATTCGACATCAACGTTTCTCGTTCATTTATACTTGGGTTAATCCGTTATTTTACCCTTTTCCACGGTAAACATCTTCGAATTTTCGTCCGACATGTCTATAACGTGTTCAGCACTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCTTACGTCTGGCGcaag                                                                                                                                                                                                                                                                                       <  1:223161/298‑3 (MQ=255)       tACTGGAGGAGTCATAAGAATTCGACATCAACGTTTCTCGCTCATTTATACTTGGGTTAATCCGTTATTTTACCCTTTTCCACGGTAAACATCTTCGAATTTTCGTCCGACATGTCTATAACGTGTTCAGCACTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATAATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTAGCGCCAGACGTAAGGCCcacagcc                                                                                                                                                                                                                                                                                    >  1:142336/1‑298 (MQ=255)           ggaggaGTCATAAGAATTCGACATCAACGTTTCTCGCTCATTTATACTTGGGTTAATCCGTTATTTTACCCTTTTCCACGGTAAACATCTTCGAATTTTCGTCCGACATGTCTATAACGTGTTCAGCACTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACac                                                                                                                                                                                                                                                                                >  1:224013/1‑300 (MQ=255)              ggagTCATAAGAATTCGACATCAACGTTTCTCGCTCATTTATACTTGGGTTAATCCGTTATTTTACCCTTTTCCACGGTAAACATCTTCGAATTTTCGTCCGACATGTCTATAACGTGTTCAGCACTGATCGCGCTGACAACGACCTGTCATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAACATAGTCAACCACGTAGAGACACCACCGCCACCTTTCCAGGATGAGAAACACCACTTGACCCAGAcacacaccacaacacaacagctt                                                                                                                                                                                                                                                                             >  1:6232/1‑279 (MQ=255)                    aTAAGAATTCGACATCAACGTTTCTCGCTCATTTATACTTGGGTTAATCCGTTATTTTACCCTTTTCCACGGTAAACATCTTCGAATTTTCGTCCGACATGTCTATAACGTGTTCAGCACTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCAGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCACTCAGGTAGAGACACCGCCGCCCGCTTTAACGGGTGAGGACCTCTCCTTGCGCAAGACGTAAGGCAAACATACAAAGCATAAgaag                                                                                                                                                                                                                                                                       >  2:627652/1‑298 (MQ=255)                            tCGACATCAACGTTTCTCGTTCATTTATACTTGGGTTAATCCGTTATTTTTCCCTTTTCCACGGTAAACATCTTCGAATTTGCGTCCGACATGTCTATAACGTGTTCAGCACTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCAcgcg                                                                                                                                                                                                                                                               <  1:169893/301‑1 (MQ=255)                                  tCAACGTTTCTCGCTCATTTATACTTGGGTTAATCCGTTATTTTACCCTTTTCCACGGTAAACATCTTCGAATTTTCGTCCGACATGTCTATAACGTGTTCAGCACTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCCACAGCTAAAGCTGCCCACACGCAAAgg                                                                                                                                                                                                                                                         >  1:623337/1‑301 (MQ=255)                                   cAACGTTTCTCGCTCATTTATACTTGGGTTAATCCGTTATTTTACCCTTTTCCACGGTAAACATCTTCGAATTTTCGTCCGACATGTCTATAACGTGTTCAGCACTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCATTCACGGGTGAGGAACTCTCCTTGAGCCAGACGGAAGGCGCACATCAAAAGCTTAAGCTGCCCACCCAATAAggg                                                                                                                                                                                                                                                        >  1:244687/1‑300 (MQ=255)                                                   ttttTACTTTGTTTAATCCTTTATTTTTCCCTTTTCCACGGTATACATCTTCGAATTTTCGGCCGACATGTCTATAACGTGTTCAGCACTGATCGCGCTGAAAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCGCGCGACGCTCATCATCAAGCTCAGTGGCAAAATCATCTATCGGGTAGAGACACCGCCGCCCGCTATCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATGTTCCACCGGCGCa                                                                                                                                                                                                                                        <  2:204539/301‑1 (MQ=255)                                                    ttttgtgggggTTGATCGGTTATTTTACCTGTTTCTACGGGAATCATCTTCGAATTTTCGTCCGACATGTATATAACGTGTTCAGCACTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTAGCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCAc                                                                                                                                                                                                                                       <  1:281295/293‑1 (MQ=255)                                                      ataCTTGGGTTAATCCGTTATTTTACTCTTTTCCACGGTAAACATCTTCGAATTTTCGTCCGACATGTCTATAACGTGTCCAGCACTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCAccc                                                                                                                                                                                                                                     >  2:263416/1‑300 (MQ=255)                                                         cTTGGGTTAATCCGTTATTTTACCCTTTTCCACGGTAAACATCTTCGAATTTTCGTCCGACATGTCTATAACGTGTTCAGCACTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTAAAGCTGCCCACGCGATAAGGTATCTTACACCGGCGCACCGtag                                                                                                                                                                                                                                  >  2:143981/1‑299 (MQ=255)                                                           cagacgtgtgctcttccgatctcTTTTCCACGGTAAACATCTTCGAATTTTCGTCCGACATGTCTATAACGTGTTCAGCACTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGgc                                                                                                                                                                                                                                <  2:312050/279‑1 (MQ=255)                                                             ggTTAATCCGTTATTTTACCCTTTTCCACGGTAAACATCTTCGAATTTTCGTCCGACATGTCTATAACGTGTTCAGCACTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGACAAAAACATCAAACAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGCGGAACTCCCATTGAGCCAGACGTAAGGCGCCACTCCACAGCTTAAGCTGCCCACGCGAAAAGGTATCTACCACCGGCCCACCgaccgccc                                                                                                                                                                                                                              >  1:478309/1‑294 (MQ=255)                                                                    tctttATTTTTCCCTTTTCCACGGTAAACATCTTCGAATTTTCGTCCGACATGTCTATAACGTGTTCAGCACTGCTCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTACGCGGCTGGCAAGCAGACCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTGAAGCTGCCCACGCGATAAGGTATCGTCCACCGGCGCACCGTCGGCGCGAATGCg                                                                                                                                                                                                                       <  1:156055/298‑1 (MQ=255)                                                                     ggTTTTTTTACCCTTTTCTCGGGTAAACTTCTTCGAATTTTCGTCCGACATGTCTATAACGTGTTCAGCACTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGt                                                                                                                                                                                                                      <  1:138453/300‑1 (MQ=255)                                                                     cGTTATTTTACCCTTTTCCACGGTAAACATCTTCGAATTTTCGTCCGACATGTCTATAACGTGTTCAGCACTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCga                                                                                                                                                                                                                      >  2:156206/1‑300 (MQ=255)                                                                            ttACCCTTTTCCACGGTAAACATCTTCGAATTTTCGTCCGACATGTCTATAACGTGTTCAGCACTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGCGGAACTCTCCTTGCGCCAGACGTAAGGCGCACCTCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTACACCGGCGCACCGTCGGCCCGAAAGCGCAAGTcac                                                                                                                                                                                                               >  1:519941/1‑299 (MQ=255)                                                                               ccttttgtcACGGGAAACATGTTCGAAGTTTCGTGCGACAGTTCGATAACGGGTGCAGCACTGATCGCGCTGACAAAGACCGGTGATGGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAGTGCGTAAGTCCGCtt                                                                                                                                                                                                            <  1:144693/293‑1 (MQ=255)                                                                                 cTTTTCCACGGTAAACATCTTCGAATTTTCGTCCGACATGTCTATAACGTGTTCAGCACTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGgcagatcggaagagcgtcgtgacg                                                                                                                                                                                                          >  1:311646/1‑279 (MQ=255)                                                                                    ttCCACGGTAAACATCTTCGAATTTTCGTCCGACATGTCTATAACGTGTTCAGCACTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGc                                                                                                                                                                                                       <  2:478910/301‑1 (MQ=255)                                                                                       cACGGTAAACATCTTCGAATTTTCGTCCGACATGTCTATAACGTGTTCAGCACTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCggg                                                                                                                                                                                                    >  2:351204/1‑301 (MQ=255)                                                                                                       gAATTTTCGTCCGACATGTCTATAACGTGTTCAGCACTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAgg                                                                                                                                                                                    >  2:314211/1‑301 (MQ=255)                                                                                                         aTTTTCGTCCGACATGTCTATAACGTGTTCAGCACTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGta                                                                                                                                                                                  >  2:96724/1‑300 (MQ=255)                                                                                                           tttCGTCCGACATGTCTATAACGTGTTCAGCACTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTCAg                                                                                                                                                                                >  1:109664/1‑301 (MQ=255)                                                                                                                 ccGACGTGTCTATACCTTGGTCAGCACTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGAAAGCAGCCCGCGACGCGCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGTGTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCg                                                                                                                                                                          <  1:350741/301‑1 (MQ=255)                                                                                                                  cGACATGTCTATAACGTGTTCAGCACTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGGTAGCTGGCGa                                                                                                                                                                         >  2:347321/1‑301 (MQ=255)                                                                                                                          atgtaatgtgttcggtccgggtgcgtcTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCGATCGCGTTc                                                                                                                                                                 <  1:474314/275‑1 (MQ=255)                                                                                                                           tataACGTGTTCAGCACTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCGATCGCGTTCa                                                                                                                                                                >  2:267813/1‑301 (MQ=255)                                                                                                                            ataACGTGTTCAGCACTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCGATCGCGTTCaa                                                                                                                                                               >  2:35878/1‑301 (MQ=255)                                                                                                                                cGTGTTCAGCACTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAAAGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCGATCGCGTTCAAAAtt                                                                                                                                                           <  2:244954/301‑1 (MQ=255)                                                                                                                                    gtcagcgcgatcagtgcTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTACGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGGGTCAGTTAGCTGGCCATCGCGTTCCAAATTACGt                                                                                                                                                       >  1:129219/16‑301 (MQ=255)                                                                                                                                     tCAGCACTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCGATCGCGTTCAAAATTACGtt                                                                                                                                                      <  2:224256/301‑1 (MQ=255)                                                                                                                                        gCACTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCGAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCGATCGCGTTCAAAATTACGTTCCa                                                                                                                                                   >  2:387676/1‑301 (MQ=255)                                                                                                                                         tgttcagacgtgtgctcttccgatctgATTGCTTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCGATCGCGTTCAAAATTACGTTCCAg                                                                                                                                                  <  2:324579/276‑1 (MQ=255)                                                                                                                                             gATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCGATCGCGTTCAAAATTACGTTCCAGCAcc                                                                                                                                              <  2:189318/301‑1 (MQ=255)                                                                                                                                                        aaaggACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCGATCGCGTTCAAAATTACGTTCCAGCACCTCAGCATATTc                                                                                                                                   <  1:35827/297‑1 (MQ=255)                                                                                                                                                                tgtggtgttgTCGCTTTTAATCGGCTTGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGTCAAAATCATCTAGCAGGTAGAGACACCGCCGCCCGCTTTCACTGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACCTCAACAGCTTCAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCGATCGCGTTCAAAATTACGTTCCAGCACCTCAGCATATTCTGTCTCtt                                                                                                                           <  1:626775/292‑1 (MQ=255)                                                                                                                                                                 gtgATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCGATCGCGTTCAAAATTACGTTCCAGAACCTCAGCATATTCTGTCTCttt                                                                                                                          >  2:523108/1‑301 (MQ=255)                                                                                                                                                                   gATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATACGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGGGCGCGGTGTAGGATAGCCAGCAACCGCGATCCAAATTACGTTTCAgaaatcaggagagcgtcgtggtgggga                                                                                                                        >  1:324175/1‑275 (MQ=255)                                                                                                                                                                      tgctttgCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCGATCGCGTTCAAAATTACGTTCCAGCACCTCAGCATATTCTGTCTCTTTCTccc                                                                                                                     <  2:538221/295‑1 (MQ=255)                                                                                                                                                                             cTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCGATCGCGTTCCAAATTACGTTCCAGCACCTCAGCATATTCTGTCTCTTTCTCCCAGCcgcc                                                                                                              >  1:440148/1‑300 (MQ=255)                                                                                                                                                                                   aGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCGATCGCGTTCAAAATTACGTTCCAGCACCTCAGCATATTCTGTCTCTTTCTCCCAGCCGCGCTGGaa                                                                                                        <  2:520613/301‑1 (MQ=255)                                                                                                                                                                                            cAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGCCCGTGCGCGGTGTAGGTTAGCTGGCGAGCGCGTTCAACATTACGTTCAACCACCTCAGCATATTCTGTCCTTTTCTCCCcgccgcccccgaccgagcaaat                                                                                               >  1:328290/1‑285 (MQ=255)                                                                                                                                                                                                 agcCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCGATCGCGTTCAAAATTACGTTCCAGCACCTCAGCATATTCTGTCTCTTTCTCCCAGCCGCGCTGGAAAGAGAAAGTCagag                                                                                          <  2:109788/301‑1 (MQ=255)                                                                                                                                                                                                      gcgACTCTTATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCACACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCGATCGCGTTCAAAATTACGTTCCAGCACCTCAGCATATTCTGTCTCTTTCTCCCAGCCGCGCTGGAAAGAGAAAGTCAGAGAAAAc                                                                                     <  2:49627/301‑1 (MQ=255)                                                                                                                                                                                                         ccGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCGATCGCGTTCAAAATTACGTTCCAGCACCTCAGCATATTCTGTCTCTTTCTCCCAGCCGCGCTGGAAAGAGAAAGTCAGAGAAAACTCa                                                                                  <  2:142461/300‑1 (MQ=255)                                                                                                                                                                                                           ccttattattcaGCTTAGAGGCAAAATCATTTATTAGGGTGAGACACCGGCGCCCGCTTTCACGGGGGAGGAAATCTTCTTGGGGCAGAAGTAAGGCGGAAATTAACAGCTTAAGCCGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCGATCGCGTTCAAAATTACGTTCCAGCACCTCAGCATATTCTGTCTCTTTCTCCCAGCCGCGCTGGAAAGAGAAAGTCAGAGAAAACTCAgg                                                                                <  2:578441/290‑1 (MQ=255)                                                                                                                                                                                                                tcatcaAGCTTATGGGCAAAAGCATCTAGCCGGTAGAGACTCCGACGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCGATCGCGTTCAAAATTACGTTCCAGCACCTCAGCATATTCTGTCTCTTTCTCCCAGCCGCGCTGGAAAGAGAAAGTCAGAGAAAACTCAGGGAGaa                                                                           <  1:96600/301‑1 (MQ=255)                                                                                                                                                                                                                     ttccgatctggCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCGATCGCGTTCAAAATTACGTTCCAGCACCTCAGCATATTCTGTCTCTTTCTCCCAGCCGCGCTGGAAAGAGAAAGTCAGAGAAAACTCAGGGAGAAATTGc                                                                      <  2:138889/292‑1 (MQ=255)                                                                                                                                                                                                                          cAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCGATCGCGTTCAAAATTACGTTCCAGCACCCCAGCATATTCTGTCTCTTTCTCCCAGCCGCGCTGGAAAGAGAAAGTCAGAGACAACTCAGGGAGACATTGcatcca                                                                 >  1:163019/1‑296 (MQ=255)                                                                                                                                                                                                                              ggCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCGATCGCGTTCAAAATTACGTTCCAGCACCTCAGCATATTCTGTCTCTTTCTCCCAGCCGCGCTGGAAAGAGAAAGTCAGAGAAAACTCAGGGAAAAATTGcagatcggaa                                                             >  1:138764/1‑292 (MQ=255)                                                                                                                                                                                                                                   aaTCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCGATCGCGTTCAAAATTACGTTCCAGCCACCCAGCATATTCTGCCTCTTTCTCCCAGCCGCGCTGGAAAGAGAAAGTCAGACAAAACTCAGGAAGAAATTGCTTACAGATAACacc                                                        >  1:561694/1‑299 (MQ=255)                                                                                                                                                                                                                                   aaTCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCGATCGCGTTCAAAATTACGTTCCAGCACCTCAGCATATTCTGTCTCTTTCTCCCAGCCGCGCTGGAAAGAGAAAGTCAGAGAAAACTCAGGGAGAAATTGCTTACAGGTATCAGc                                                        <  1:263148/301‑1 (MQ=255)                                                                                                                                                                                                                                               ggTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCGATCGCGTTCAAAATTACGTTCCAGCACCTCAGCATATTCTGTCTCTTTCTCCCAGCCGCGCTGGAAAGAGAAAGTCAGAGAAAACTCAGGGAGAAATTGCTTACAGGTATCAGCCATATCGGCcgg                                            >  2:4062/1‑300 (MQ=255)                                                                                                                                                                                                                                                       cacCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCGATCGCGTTCAAAATTACGTTCCAGCACCTCAGCATATTCTGTCTCTTTCTCCCAGCCGCGCTGGAAAGAGAAAGTCAGAGAAAACTCAGGGAGAAATTGCTTACAGGTATCAGCCATATCGGCCGCGATACCgg                                    >  2:116570/1‑301 (MQ=255)                                                                                                                                                                                                                                                             cgcccgcTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCGATCGCGTTCAAAATTACGTTCCAGCACCTCAGCATATTCTGTCTCTTTCTCCCAGCCGCGCTGGAAAGAGAAAGTCAGAGAAAACTCAGGGAGAAATTGCTTACAGGTATCAGCCATATCGGCCGCGATACCGGCGCtat                              <  2:328694/301‑1 (MQ=255)                                                                                                                                                                                                                                                              gcccgcTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTCCGCGGGGTAGGTTAGCTGGCGATCGCGTTCAAAATTACGTTCCAGCACCTCCGCATCTTCTGTCACTTTCTCCCAGCCGCGCTGGACAGGGAAAGTCCGAGAAAACACAGGGAGAAATTGCTAACAGGTATCCTCCATATCGGCCGCGATCACGGCGCtaaa                             >  1:588041/1‑299 (MQ=255)                                                                                                                                                                                                                                                               cccggTTTCACGGGTGCGGAACTTTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCGATCGCGTTCAAAATTACGTTCCAGCACCTCAGCATATTCTGTCTCTTTCTCCCAGCCGCGCTGGAAAGAGAAAGTCAGAGAAAACTCAGGGAGAAATTGCTTACAGGTATCAGCCATATCGGCCGCGATACCGGCGCTATAc                            <  1:4062/301‑1 (MQ=255)                                                                                                                                                                                                                                                               cccgcTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCGATCGCGTTCAAAATTACGTTCCAGCACCTCAGCATATTCTGTCTCTTTCTCCCAGCCGCGCTGGAAAGAGAAAGTCAGAGAAAACTCAGGGAGAAATTGCTTACAGGTATCAGCCATATCGGCCGCGATACCGGCGCTATAc                            >  2:389256/1‑301 (MQ=255)                                                                                                                                                                                                                                                                   cttgtacgggggtgggtctttcctttgcccAGACTTAAGGCGCACTGCAACAGCTTCAGCTGCCCACGCGAGAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCGATCGCGTTCAAAATTACGTTCCAGCACCTCAGCATATTCTGTCTCTTTCTCCCAGCCGCGCTGGAAAGAGAAAGTCAGAGAAAACTCAGGGAGAAATTGCTTACAGGTATCAGCCATATCGGCCGCGATACCGGCGCTATACTCcg                        <  1:313807/273‑1 (MQ=255)                                                                                                                                                                                                                                                                          gggTGGGGAACTCTCCTGGGGCCAGAGGTAAGGCGCACCGCAACAGGTGAAGCTGCCCACGCGGTAAGGTATCTTCCACCGGCGGACCGTCGGGGGGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCGATCGCGTTCAAAATTACGTTCCAGCACCTCAGCATATGCTGTCTCTTTCTCCCAGCCGCGCTGGAAAGAGAAAGTCAGAGAAAACTCAGGGAGAAATTGCTTACAGGTATCAGCCATATCGGCCGCGATACCGGCGCTATACTCCGCGCGCCa                 <  1:522436/301‑1 (MQ=255)                                                                                                                                                                                                                                                                             tGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCGATCGCGTTCAAAATTACGTTCCAGCACCTCAGCATATTCTGTCTCTTTCTCCCAGCCGCGCTGGAAAGAGAAAGTCAGAGAAAACTCAGGGAGAAATTGCCTACAGGTATCAGCCATATCGGCCGCGATACCGGCGCTATACTCCGCGCGCCAggg              >  2:219694/1‑300 (MQ=255)                                                                                                                                                                                                                                                                              gatgAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCGATCGCGTTCAAAATTACGTTCCAGCACCTCAGCATATTCTGTCTCTTTCTCCCAGCCGCGCTGGAAAGAGAAAGTCAGAGAAAACTCAGGGAGAAATTGCTTACAGGTATCAGCCATATCGGCCGCGATACCGGCGCTATACTCCGCGCGCCAGGTg             <  2:624196/298‑1 (MQ=255)                                                                                                                                                                                                                                                                               tttccgatcttgtGCGCCTTACGTCTGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCGATCGCGTTCAAAATTACGTTCCAGCACCTCAGCATATTCTGTCTCTTTCTCCCAGCCGCGCTGGAAAGAGAAAGTCAGAGAAAACTCAGGGAGAAATTGCTTACAGGTATCAGCCATATCGGCCGCGATACCGGCGCTATACTCCGCGCGCCAGGTGc            <  2:362923/289‑1 (MQ=255)                                                                                                                                                                                                                                                                                ggaatcctgtgtGCGCCAGCCGTAAGGTGAACATCATCAGCTTCTGCGGCCCACGCTGTAAGTGTGCTAGCACGGGGGCACCGTCGTGGCGAATGCGTATGTCCGCTTTGTGCGGGCCGTGCGCGTTGAAGGTTAGCTGGCGATCGCGTTCAAAATTACGTTCCAGCACCTCAGCATTTTCTGTCGCTTGCTCCCAGCCGCGCTGGAAAGAGAAAGTCAGAGAAAACTCAGGGAGAAATTGCTTACAGGTATCAGCCATATCGGCCGCGATACCGGCGCTATACTCCACGCGCCAGGTGCt           <  1:267545/290‑1 (MQ=255)                                                                                                                                                                                                                                                                                 gctgttgatgtGCGCCTTACGTCTGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCGATCGCGTTCAAAATTACGTTCCAGCACCTCAGCATATTCTGTCTCTTTCTCCCAGCCGCGCTGGAAAGAGAAAGTCAGAGAAAACTCAGGGAGAAATTGCTTACAGGTATCAGCCATATCGGCCGCGATACCGGCGCTATACTCCGCGCGCCAGGTGCTg          >  2:11317/11‑301 (MQ=255)                                                                                                                                                                                                                                                                                    ctctCCTTGCGCCAGACGTAAGGCGCACATCAACAGCGTAAGCTGCCCACGCGATAAGGTATCTGCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCGATCGCGTTCAAAATTACGTTCCAGCACCTCAGCATATTCTGTCTCTTTCTCCCAGCCGCGCTGGAAAGAGAAAGTCAGAGAAAACTCAGGGAGAAATTGCTTACAGGTATCAGCCATATCGGCCGCGATACCGGCGCTATACTCCGCGCGCCAGGTGCTGAtt       <  1:116446/301‑1 (MQ=255)                                                                                                                                                                                                                                                                                        ccTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCGATCGCGTTCAAAATTACGTTCCAGCACCTCAGCATATTCTGTCTCTTTCTCCCAGCCGCGCTGGAAAGAGAAAGTCAGAGAAAACTCAGGGAGAAATTGCTTACAGGTATCAGCCATATCGGCCGCGATACCGGCGCTATACTCCGCGCGCCAGGTGCTGATTTGCt   >  2:126704/1‑301 (MQ=255)                                                                                                                                                                                                                                                                                         tgtGCGCCTTACGTCTGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCGATCGCGTTCAAAATTACGTTCCAGCACCTCAGCATATTCTGTCTCTTTCTCCCAGCCGCGCTGGAAAGAGAAAGTCAGAGAAAACTCAGGGAGCAATTGCTTACAGGTATCAGCCATATCGGCAGCGATACCGGCCCTATACTCCGCGCACCAGgcgcagcacggaag  >  1:362439/3‑291 (MQ=255)                                                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                                                                    CTTTGATACTGGAGGAGTCATAAGAATTCGACATCAACGTTTCTCGCTCATTTATACTTGGGTTAATCCGTTATTTTACCCTTTTCCACGGTAAACATCTTCGAATTTTCGTCCGACATGTCTATAACGTGTTCAGCACTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTGTAGGTTAGCTGGCGATCGCGTTCAAAATTACGTTCCAGCACCTCAGCATATTCTGTCTCTTTCTCCCAGCCGCGCTGGAAAGAGAAAGTCAGAGAAAACTCAGGGAGAAATTGCTTACAGGTATCAGCCATATCGGCCGCGATACCGGCGCTATACTCCGCGCGCCAGGTGCTGATTTGCTC  >  NC_000913/3880086‑3880666

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