BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,344,036	G→C	6.3%	T104S (ACT→AGT)	yciT ←	global regulator of transcription; DeoR family

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,344,036	0	G	C	6.3%	546.1 / 21.8	159	T104S (ACT→AGT)	yciT	global regulator of transcription; DeoR family
Reads supporting (aligned to +/- strand): ref base G (75/74); new base C (5/5); total (80/79)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.37e-01

Only 100 of 161 total aligned reads displayed.

CTCTTTTTCCAGCACGGCATTGACCACATCGGTACGCATCATATCGCGACCGGTAAATCCAGTTTCAGGTTGCCAGCCATCAATACCAATAAATGCTTTGCTGAAATGCACCTGTTGGATGCACTGGCGTGTCAAAGGGCCAACCATACTTTCGCTTTTTTTCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGACGTCATCACTATCAAGCGAAACAGCAAAGCCATGTGCACGGCGGAGGTAACTCAGTTTTTCGAGTGTGTTGAGATCCTGGCGAATGGTAACTTCAGAAACTCCAGTGGCTTTTGCCAGATCGGTTACGCTAACCTGA  >  NC_000913/1343765‑1344299

ctctTTTTCCAGCACGGCATTGACCACATCGGTACGCATCATATCGCGACCGGTAAATCCAGTTTCAGGTTGCCAGCCATCAATACCAATAAATGCTTTGCTGAAATGCACCTGTTGGATGCACTGGCGTGTCAAAGGGCCAACCATACTTTCGCTTTTTTTCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTc                                                                                                                                                                                                                                                                       <  2:338127/274‑1 (MQ=255)

  ctTTTACCAGAACGGAATTGACCACATCGGTACGCATCATATCGCGACCGGTAAATCCAGTTTCAGGTTGCCAGCCATCAATACCAATAAATGCTTTGCTGAAATGCACCTGTTGGATGCACTGGCGTGTCAAAGGGCCAACCATACTTTCGCTTTTTTTCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCgg                                                                                                                                                                                                                                                                    <  2:43957/275‑1 (MQ=255)

        ccaggccGGCATTGACCACATCGGTACGCATCATATCGCGACCGGTAAACCCCGTTTCAGGTTGCCAGCCATCAATACCAATAAATGCTTTGCTGAAATGCACCTGTTGGATGCACTGGCGTGTCAAAGGGGCAACCATACTTTCGCTTTTTTTCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTGTTCTTCTGCTCGCCCAGAGTCCGGGCGAg                                                                                                                                                                                                                                                               <  1:173363/268‑1 (MQ=255)

         acgaaCGGCATTGACCACATCGGAACGAATCATATCGCGACCGGTAAATCCAGTTTCAGGTTGCCAGCCATCAATACCAATAAATGCTTTGCTGAAATGCACCTGTTGGATGCACTGGCGTGTCAAAGGGCCAACCATACTTTCGCTTTTTTTCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTa                                                                                                                                                                                                                                                             <  4:181656/271‑1 (MQ=255)

          agaaCGGCATTGACCACATCGGTACGCATCATATCGCGACCGGTAAATCCAGTTTCAGGTTGCCAGCCATCAATACCAATAAATGCTTTGCTGAAATGCACCTGTTGGATGCACTGGCGTGTCAAAGGGCCAACCATACTTTCGCTTTTTTTCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAg                                                                                                                                                                                                                                                            <  1:107777/272‑1 (MQ=255)

             aCGGCATTGACCACATCGGTACGCATCATATCGCGACCGGTAAATCCAGTTTCAGGTTGCCAGCCATCAATACCAATAAATGCTTTGCTGAAATGCACCTGTTGGATGCACTGGCGTGTCAAAGGGCCAACCATACTTTCGCTTTTTTTCTGATACACGCCACCGAGCAAAATAACTTCACCAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGCTGATACTGACATTTTTCTTCTGCCCGCCCAGAGTCCGGGCGCGTAGGGc                                                                                                                                                                                                                                                         >  2:164000/1‑275 (MQ=255)

                                                                   ggTTGCCAGCCATCAATACCAATAAATGCTTTGCTGAAATGCACCTGTTGGATGCACTGGCGTGTCAAAGGGCCAACCATACTTTCGCTTTTTTTCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCg                                                                                                                                                                                                                                                                 >  3:175280/1‑213 (MQ=255)

                                                                       gccagccaTCAATACCAATAAATGCTTTGCTGAAATGCACCTGTTGGATGCACTGGCGTGTCAAAGGGCCAACCATACTTTCGCTTTTTTTCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTa                                                                                                                                                                                                                                                             <  1:175281/213‑1 (MQ=255)

                                                                                             tGCTTTGCTGAAATGCACCTGTTGGATGCACTGGCGTGTCAAAGGGCCAACCATACTTTCGCTTTTTTTCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCa                                                                                                                                                                         <  2:525961/275‑1 (MQ=255)

                                                                                                 ttGCTGAAATGCACCTGTTGGATGCACTGGCGTGTCAAAGGGCCAACCATACTTTCGCTTTTTTTCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCAc                                                                                                                                                                               >  2:573189/1‑265 (MQ=255)

                                                                                                      gAAATGCACCTGTTGGATGCACTGGCGTGTCAAAGGGCCAACCATACTTTCGCTTTTTTTCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCAtt                                                                                                                                                                                                                                          <  2:122478/201‑1 (MQ=255)

                                                                                                              ccTGTTGGATGCACTGGCGTGTCAAAGGGCCAACCATACTTTCGCTTTTTTTCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACCTTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGAGGATAGTTGCgcatc                                                                                                                                                        >  1:134068/1‑275 (MQ=255)

                                                                                                                            tGGCGTGTCAAAGGGCCAACCATACTTTCGCTTTTTTTCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCCAACTCGGCGATTTCACGCTTCAGCGTAGGGTTGCtcatcagctgtcccttat                                                                                                                                           >  4:85957/1‑262 (MQ=255)

                                                                                                                              gggtgtCAAAGGGCAAACGAACACTTCGCTGTGTTTGAGAAGCACGACACCGGGAAAAATAACTTGACAAGGCGCGTCTTTCAGCAAATGCGGGTGGTAGCTGCTGACCGGGATGTTAGTGACATTTTTCTTCGGCTCGCCAAGAGTGCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGACGt                                                                                                                                        <  4:163999/273‑1 (MQ=255)

                                                                                                                                 tgtCAAAGGGCCAACCATACTTTCGCTTTTTTTCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCggg                                                                                                                                                                                                                                                                   <  7:703006/149‑1 (MQ=255)

                                                                                                                                 tgtCAAAGGGCCAACCATACTTTCGCTTTTTTTCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGACGtcat                                                                                                                                     >  4:264438/1‑275 (MQ=255)

                                                                                                                                 tgtCAAAGGGCCAACCATACTTTCGCTTTTTTTCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCCCCGCATACGGGTTTCGACGtcat                                                                                                                                     >  2:385164/1‑275 (MQ=255)

                                                                                                                                      aaGGGCCAACCATACTTTCGCTTTTTTTCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAg                                                                                                                                                                                                                                                               >  5:879200/1‑148 (MQ=255)

                                                                                                                                      aaGGGCCAACCATACTTTCGCTTTTTTTCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGt                                                                                                                                                                                                                                                              >  7:438764/1‑149 (MQ=255)

                                                                                                                                       aGGGCCAACCATACTTTCGCTTTTTTTCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGt                                                                                                                                                                                                                                                              >  5:573159/1‑148 (MQ=255)

                                                                                                                                       aGGGCCAACCATACTTTCGCTTTTTTTCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTa                                                                                                                                                                                                                                                             <  7:160866/149‑1 (MQ=255)

                                                                                                                                        gggCCAACCATACTTTCGCTTTTTTTCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAg                                                                                                                                                                                                                                                            <  7:268712/149‑1 (MQ=255)

                                                                                                                                         ggCCAACCATACTTTCGCTTTTTTTCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAgg                                                                                                                                                                                                                                                           >  7:183705/1‑149 (MQ=255)

                                                                                                                                          gCCAACCATACTTTCGCTTTTTTTCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTACTCGCCCGGACTCTGGGCGAGCAgaa                                                                                                                                                                                                                                                          <  8:758208/149‑3 (MQ=255)

                                                                                                                                           ccAACCATACTTTCGCTTTTTTTCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGc                                                                                                                                                                                                                                                         >  8:157908/1‑149 (MQ=255)

                                                                                                                                             aaCCATACTTTCGCTTTTTTTCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGACGTCATCACTATCAAGCg                                                                                                                         >  2:519927/1‑275 (MQ=255)

                                                                                                                                               ccATACTTTCGCTTTTTTTCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATtg                                                                                                                                                                                                                                                     >  7:123506/1‑149 (MQ=255)

                                                                                                                                                 aTACTTTCGCTTTTTTTCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATtgct                                                                                                                                                                                                                                                   >  8:853888/1‑149 (MQ=255)

                                                                                                                                                    cTTTCGCTTTTTTTCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATtgctgc                                                                                                                                                                                                                                                 <  6:680672/148‑1 (MQ=255)

                                                                                                                                                      ttCGCTTTTTTTCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATtgctgctg                                                                                                                                                                                                                                               >  6:250740/1‑148 (MQ=255)

                                                                                                                                                       tCGCTTTTTTTCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATtgctgctgc                                                                                                                                                                                                                                              >  5:294874/1‑148 (MQ=255)

                                                                                                                                                             tttttCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCAtttt                                                                                                                                                                                                                                        <  6:573159/148‑1 (MQ=255)

                                                                                                                                                               tttCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCg                                                                                                                                                                                                                                      <  6:711829/148‑1 (MQ=255)

                                                                                                                                                               tttCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCg                                                                                                                                                                                                                                      <  6:563532/148‑1 (MQ=255)

                                                                                                                                                               tttCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGa                                                                                                                                                                                                                                     <  7:791832/149‑1 (MQ=255)

                                                                                                                                                               tttCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGa                                                                                                                                                                                                                                     <  8:109928/149‑1 (MQ=255)

                                                                                                                                                                ttCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGACGTCATCACTATCAAGCGAAACAGCAAAGCCATGTGc                                                                                                      <  4:519925/275‑1 (MQ=255)

                                                                                                                                                                   ttaTACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCAt                                                                                                                                                                                                                                                       <  3:84972/125‑1 (MQ=255)

                                                                                                                                                                   tGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCgg                                                                                                                                                                                                                                                                    >  8:246786/1‑114 (MQ=255)

                                                                                                                                                                   tGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCgg                                                                                                                                                                                                                                                                    <  7:246786/114‑1 (MQ=255)

                                                                                                                                                                   tGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATaaa                                                                                                                                                                                                                                 <  8:391460/149‑1 (MQ=255)

                                                                                                                                                                     aTACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAg                                                                                                                                                                                                                                >  6:379514/1‑148 (MQ=255)

                                                                                                                                                                         acGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAgg                                                                                                                                                                                                                                                           >  8:59983/1‑117 (MQ=255)

                                                                                                                                                                           gCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGttt                                                                                                                                                                                                                         >  7:673894/1‑149 (MQ=255)

                                                                                                                                                                                 gtgCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCg                                                                                                                                                                                                                   >  7:709230/3‑149 (MQ=255)

                                                                                                                                                                                 gAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCAtt                                                                                                                                                                                                                                          <  7:494870/126‑1 (MQ=255)

                                                                                                                                                                                  aGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCg                                                                                                                                                                                                                   >  5:457879/1‑148 (MQ=255)

                                                                                                                                                                                  aGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGc                                                                                                                                                                       >  4:445791/1‑192 (MQ=255)

                                                                                                                                                                                     aaaaTAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTg                                                                                                                                                                                                               >  7:61005/1‑149 (MQ=255)

                                                                                                                                                                                      aaaTAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGtat                                                                                                                                                                   <  2:445793/192‑1 (MQ=255)

                                                                                                                                                                                            cTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGc                                                                                                                                                                                                        <  7:737680/149‑1 (MQ=255)

                                                                                                                                                                                             ttCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCg                                                                                                                                                                                                                                                                     <  5:374986/87‑1 (MQ=255)

                                                                                                                                                                                             ttCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCg                                                                                                                                                                                                                                                                     >  6:374986/1‑87 (MQ=255)

                                                                                                                                                                                              tCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTACTCGCCCGGACTCTGGGCGAGCAga                                                                                                                                                                                                                                                           >  5:794574/1‑95 (MQ=38)

                                                                                                                                                                                                acaAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGAc                                                                                                                                                                                                     >  6:244732/1‑148 (MQ=255)

                                                                                                                                                                                                  aaGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGAcgc                                                                                                                                                                                                   <  6:309165/148‑1 (MQ=255)

                                                                                                                                                                                                            tttgcgCAATTGGGGGAGGTGGCTGCTGCCCGGTATGATAGGGACATTTTTCGTCTGCTCGCCCAGAGTCCGGGGGAGTAGTGCATTGCGGCTGTCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGAAGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTGGCTCATCATACGGGTGTCGACTTCATCACTATCAAGCGAAACAGCAAAGCCATGTGCACGGCGGAGGTAACTCAGTTTTTCGAGTGTGTTGAGATCCTGGc                                                          <  4:703556/270‑1 (MQ=255)

                                                                                                                                                                                                             ttCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCgg                                                                                                                                                                                        <  5:756556/148‑1 (MQ=255)

                                                                                                                                                                                                             ttCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGc                                                                                                                                                                                       <  8:717815/149‑1 (MQ=255)

                                                                                                                                                                                                             ttCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGc                                                                                                                                                                                       <  8:475716/149‑1 (MQ=255)

                                                                                                                                                                                                                 gCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGt                                                                                                                                                                                   >  7:894042/1‑149 (MQ=255)

                                                                                                                                                                                                                  cAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGt                                                                                                                                                                                   >  6:276902/1‑148 (MQ=255)

                                                                                                                                                                                                                  cAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGt                                                                                                                                                                                   >  5:53849/1‑148 (MQ=255)

                                                                                                                                                                                                                    aaTGCGGGTTGTAGCTGGTGACCGTGATGTTGGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCCGGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGACGTCATCACTATCAAGCGAAACAGCAAAGCCATGTGCACGGCGGAGGTAACTCAGTTTTTCGAGTGTGTTGAGATCCTGGCGAATGGt                                                   <  4:570256/274‑1 (MQ=255)

                                                                                                                                                                                                                    aaTGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGt                                                                                                                                                >  2:213711/1‑181 (MQ=255)

                                                                                                                                                                                                                        cgcgATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCg                                                                                                                                            <  4:213710/181‑1 (MQ=255)

                                                                                                                                                                                                                           gATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCa                                                                                                                                                                         >  8:625209/1‑149 (MQ=255)

                                                                                                                                                                                                                            aTGTAGCTGCTGCCCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGACGTCATCACTATCAAGCGAAACAGCAAAGCCATGTGCACGGCGGAGGTAACTCAGTTTTTCGAGTGTGTTGAGATCCTGGCGAATGGTAACTTCAGa                                          <  2:234246/275‑1 (MQ=255)

                                                                                                                                                                                                                            aTGTAGCTGCTGACCGGGATGTTAGGGACATTTTTTTTCTTCTCGCCCAGGGGCCGGGCGAGTGGGGCATTGCTGCTGCCGTTTTCGATAATGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTAGAGTTGCTCATCATACGGGTTTCGACGTCATCACTATCAAGCGAAACAGCAAAGCCATGTGCACGGCGGAGGTAACTCAGTTTTTCGAGTGTGTTGAGATCCTGGCGAATGGTAACTTCAGa                                          <  2:478166/275‑1 (MQ=255)

                                                                                                                                                                                                                                  ctgctgACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGtata                                                                                                                                                                  <  7:853888/149‑1 (MQ=255)

                                                                                                                                                                                                                                   tgctgACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAg                                                                                                                                                                 >  7:627251/1‑149 (MQ=255)

                                                                                                                                                                                                                                       gACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCTGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGACGTCACCTCTATCAAGCGAAACAGCAAAGCCATGTGCACGGCGGAGGCAACTCAGTTTTTCGAGTGTGTTGAGATCCTGGCGCATGGTAACTTCAGAAACTCCAg                                  >  1:4116/1‑272 (MQ=255)

                                                                                                                                                                                                                                       gACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGc                                                                                                                                                             <  8:421665/149‑1 (MQ=255)

                                                                                                                                                                                                                                           gTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCtca                                                                                                                                                          >  5:631945/1‑148 (MQ=255)

                                                                                                                                                                                                                                           gTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCtca                                                                                                                                                          <  5:250740/148‑1 (MQ=255)

                                                                                                                                                                                                                                             gatgatAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTc                                                                                                                                                                          >  4:366373/1‑130 (MQ=255)

                                                                                                                                                                                                                                                gatAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCa                                                                                                                                                                         <  2:735258/128‑1 (MQ=255)

                                                                                                                                                                                                                                                      gACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGt                                                                                                                                                                     >  5:842038/1‑126 (MQ=255)

                                                                                                                                                                                                                                                      gACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGt                                                                                                                                                                     <  6:842038/126‑1 (MQ=255)

                                                                                                                                                                                                                                                      gACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATa                                                                                                                                                     >  6:30644/1‑142 (MQ=255)

                                                                                                                                                                                                                                                       aCATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGttt                                                                                                                                              >  5:298112/1‑148 (MQ=255)

                                                                                                                                                                                                                                                       aCATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTc                                                                                                                                             <  7:393462/149‑1 (MQ=255)

                                                                                                                                                                                                                                                       aCATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTc                                                                                                                                             >  8:348696/1‑149 (MQ=255)

                                                                                                                                                                                                                                                       aCATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTc                                                                                                                                             >  7:551363/1‑149 (MQ=255)

                                                                                                                                                                                                                                                       aCATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGACGTCATCACTATCAAGCGAAACAGCAAAGCCATGTGCACGGCGGAGGTAACTCAGTTTTTCGAGTGTGTTGAGATCCTGGCGAATGGTAACTTCAGAAACTCCAGTGGCTTTTGCCAGAGCgg                >  2:399247/1‑274 (MQ=255)

                                                                                                                                                                                                                                                       aCATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTCCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGttt                                                                                                                                              <  6:750791/148‑1 (MQ=255)

                                                                                                                                                                                                                                                           ttttCTTCTGCTCGCCCAGGGTCCGGGCGAGTAGTGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGACGTCATCACTATCAATCGAAACAGCAAAGCCATGTGCACGGCGGAGGTAACTCAGTTTTTCGAGTGTGTTGAGATCCTGGCGAATGGTAACTTCAGAAACTCCAGTGGCTTTTGCCAGATCGGTTACg           <  2:573183/275‑1 (MQ=255)

                                                                                                                                                                                                                                                                 tctGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGACGtcatc                                                                                                                                    >  5:784700/1‑148 (MQ=255)

                                                                                                                                                                                                                                                                  ctacTCGCCCGGACTCTGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCgg                                                                                                                                                                                                 >  8:485087/4‑86 (MQ=38)

                                                                                                                                                                                                                                                                  ctacTCGCCCGGACTCTGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCgg                                                                                                                                                                                                 <  7:485087/83‑1 (MQ=38)

                                                                                                                                                                                                                                                                   tacTCGCCCGGACTCTGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGa                                                                                                                                           >  6:595157/3‑139 (MQ=255)

                                                                                                                                                                                                                                                                     cTCGCCCAGAGTCCGGGCGAGTGGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGta                                                                                                                                                                    <  8:132318/112‑1 (MQ=255)

                                                                                                                                                                                                                                                                     cTCGCCCAGAGTCCGGGCGAGTGGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGta                                                                                                                                                                    >  7:132318/1‑112 (MQ=255)

                                                                                                                                                                                                                                                                     cTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGACGTCATCACTATCAAGCGAAACAGCAAAGCCATGTGCACGGCGGAGGTAACTCAGTTTTTCGAGTGTGTTGAGATCCTGGCGAAGGGTAAATTCAGAAACTCAAGGGGCTTTTGCCAGATAGGTTACGCTaaaccta  >  4:707896/1‑269 (MQ=255)

                                                                                                                                                                                                                                                                      tCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGtt                                                                                                                                               >  8:211269/1‑132 (MQ=255)

                                                                                                                                                                                                                                                                      tCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGACGTCATCACTATc                                                                                                                              >  7:724046/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                          ccAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGACGTCATCACTATCAAGc                                                                                                                          >  8:562936/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                          ccAGAGTCCGGGCGAGTAAGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGACGTCATCACTATCAAGc                                                                                                                          <  7:348696/149‑1 (MQ=255)

                                                                                                                                                                                                                                                                            agagTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGACGTCATCACTATCAAGCg                                                                                                                         >  5:211889/1‑148 (MQ=255)

                                                                                                                                                                                                                                                                               gTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGACGTCATCACTATCAAGCGaaa                                                                                                                      <  5:894077/148‑1 (MQ=255)

CTCTTTTTCCAGCACGGCATTGACCACATCGGTACGCATCATATCGCGACCGGTAAATCCAGTTTCAGGTTGCCAGCCATCAATACCAATAAATGCTTTGCTGAAATGCACCTGTTGGATGCACTGGCGTGTCAAAGGGCCAACCATACTTTCGCTTTTTTTCTGATACACGCCACCGAGCAAAATAACTTCACAAGGCGCGTCTTTCAGCAAATGCGCGATGTAGCTGCTGACCGTGATGATAGTGACATTTTTCTTCTGCTCGCCCAGAGTCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGACGTCATCACTATCAAGCGAAACAGCAAAGCCATGTGCACGGCGGAGGTAACTCAGTTTTTCGAGTGTGTTGAGATCCTGGCGAATGGTAACTTCAGAAACTCCAGTGGCTTTTGCCAGATCGGTTACGCTAACCTGA  >  NC_000913/1343765‑1344299

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 11 ≤` `ATCG/ATCG` `< 25 ≤` `ATCG/ATCG` `< 37 ≤` `ATCG/ATCG` `< 40 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`