BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,516,110	C→T	A190T (GCC→ACC)	fecA ←	ferric citrate outer membrane transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,516,110	0	C	T	93.7%	372.8 / NA	130	A190T (GCC→ACC)	fecA	ferric citrate outer membrane transporter
Reads supporting (aligned to +/- strand): ref base C (1/0); major base T (48/70); minor base G (6/2); total (58/72)
Fisher's exact test for biased strand distribution p-value = 7.31e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 130 total aligned reads displayed.

AAGAGGTTGGGCTGAGCTGACCTTCCACGCCCGCCTCGATACCAAAGTCCTGCGGAATGGCACGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGCGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTT‑AAGTACCTCACGCATGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGATCACGTCACGCGC‑GCCAGCATGTTCAAATACGTCGTTTTCACGCGCATCACCCAGCCAGTCGCCGACCACGGTCAGGGCATCTTCTTTTGGT  >  NC_000913/4515867‑4516351

aaGAGGTTGGGCTGAGCTGACCTTCCACGCCCGCCTCGATACCAAAGTCCTGCGGAATGGCACGGGTAACAAAGTTCACCACTCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGGGGGgcgc                                                                                                                                                                                                                                              >  2:440487/1‑251 (MQ=255)

   gggtgggttggtgtggaccacctccgcccgccTCGATACCAAAGTCCTGAGGACGGGCACGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCGCCGCTCCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAAAGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGTGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGTGAGGCGCggg                                                                                                                                                                                                                                           <  2:593852/227‑1 (MQ=255)

             gagctgaCCCTCCCCGCCCGCCTCGATACCAAAGTCCTGCGGAATGGCACGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGcc                                                                                                                                                                                                                                   <  1:563122/249‑1 (MQ=255)

             cgccgacccTTCCCCGCCCGCCCCGATACCAAAGTCCTGCGGAATGGCACGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCAACCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGTGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGcc                                                                                                                                                                                                                                   <  2:348887/242‑1 (MQ=255)

                      ttCCACGCCCGCCTCGATACCAAAGTCCTGCGGAATGGCACGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTACGGGGCGAAGGGGACGGGGCTGCCGTCCATCAGGACGGTCGAGCGGCTGGAGAGGCGCGGGTTCAGGCCCCGGATGCCa                                                                                                                                                                                                                         >  1:620071/1‑250 (MQ=255)

                       cccACGCCCGCCTCGATACCAAAGTCCTGCGGAATGGCACGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGGATGCCaa                                                                                                                                                                                                                        <  1:278821/249‑1 (MQ=255)

                          ccgcccgccTCGATACCAAAGTCCTGCGGAATGGCACGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGtt                                                                                                                                                                                                                    <  1:680650/250‑1 (MQ=255)

                              ccgccTCGATACCAAAGGCCTGCGGAATGGCACGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATc                                                                                                                                                                                                                <  1:634683/251‑1 (MQ=255)

                                    cGATACCAAAGTCCTGCGGAATGGCACGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAg                                                                                                                                                                                                           <  1:610625/250‑1 (MQ=255)

                                      aTACCAAAGTCCAGCGGAATGGCACGGGTAACAAAGTTCACCACGCCGCCCCCGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCACTGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGCCCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCCTCAGAACGTGCGAGCGGCTGGAGAGGCGCGGGTTCAAGCCCCCGATGCCACAGGTCAACGCCAgggg                                                                                                                                                                                                        >  2:422489/1‑249 (MQ=255)

                                       aaCCAAAGTCCTGCGGAATGGCACGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCg                                                                                                                                                                                                       <  1:183314/250‑1 (MQ=255)

                                          taaaGACCTGCGGAATGGCACGGGTAACAAAGTTCACCACGCCGCCAACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTCCAACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGAGGCTGACCGTAAGGGGCGACGGGGACGGGGTTGCCTGCCATCAGGCCGGTCGAGAGGATGGTGGGGCGCGGGTTCGGGCCCCGGATGCCAAAGGTCCACGCCAGGTCGTgg                                                                                                                                                                                                    >  2:25084/2‑251 (MQ=255)

                                               tCCTGCGGAATGGCACGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGCGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCCGGTCGTGGCTGcc                                                                                                                                                                                               >  1:253906/1‑251 (MQ=255)

                                                ccTGCGGTATGGCCCTGTTAACAAAGGTCCCCACGCCGCCCACGCACTTCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAAGGGCATCCATGTTGCCGAGCGAAACGGGATCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGTCTGcc                                                                                                                                                                                               <  2:147097/250‑1 (MQ=255)

                                                   gCGGAATGGCACGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGt                                                                                                                                                                                            <  1:378969/250‑1 (MQ=255)

                                                    cGGAATGGCACGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGc                                                                                                                                                                                          <  1:525728/251‑1 (MQ=255)

                                                     ggAATGGCACGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCCGGTCGTGGCTGCCGGTGcc                                                                                                                                                                                         >  1:573488/1‑251 (MQ=255)

                                                     ggAATGGCACGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGcc                                                                                                                                                                                         <  1:131013/251‑1 (MQ=255)

                                                      gAATGGCACGGGTACCAAAGTTCAACACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATACCGTCCATCAGGACGGTCGAACGGCTGGTGAGGCGCGGGTTCAGGACCCGGATGCCAAAAGTCCTCGCCAGGGCGTGGccggcgggggcg                                                                                                                                                                                        >  2:505773/1‑240 (MQ=255)

                                                       aaTGGCCCGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCgt                                                                                                                                                                                       <  2:163335/251‑1 (MQ=255)

                                                             aCGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTg                                                                                                                                                                                                     >  1:423842/1‑231 (MQ=255)

                                                                ggTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCGTTGTTTTCCg                                                                                                                                                                              <  1:469461/251‑1 (MQ=255)

                                                                   aaCAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCGTTGTTTTCCGgcg                                                                                                                                                                           >  1:73134/1‑251 (MQ=255)

                                                                        aGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTg                                                                                                                                                                      <  1:400607/251‑1 (MQ=255)

                                                                            cacctcGCCGCCTACGCGCGGCGGTCCGTTACGCACCGCACCTCCCCCGCGTACCACGGCAATGGCCTCCAGGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGTTGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACg                                                                                                                                                                   <  2:284648/250‑1 (MQ=255)

                                                                            caccacGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGc                                                                                                                                                                  <  1:394955/251‑1 (MQ=255)

                                                                              caccgccgcCCCCGCTCTGCGGTCCGTTACTCTCCGCACCACCACCGCTTACCACGGTAATGGCAGCCATGTTGCCGAGCGACACGGGAGCCAGTGAATGCTGCGGCTGACGGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGTTGAGGCGCTGGTTCAGGCCCCGGATTCCTAAGTTCAGCGCCAGTTCGTGGCTGGCGGTGCCGTTGTTTTCCGGCGCGCTGACGcga                                                                                                                                                                <  2:434887/248‑3 (MQ=255)

                                                                               cacGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAg                                                                                                                                                               <  1:35404/251‑1 (MQ=255)

                                                                                        cACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGGATGTCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGt                                                                                                                                                      <  1:655427/251‑1 (MQ=255)

                                                                                          cGCTCTGCGGCCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGGCCATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGt                                                                                                                                                      <  2:146203/249‑1 (MQ=255)

                                                                                             tctGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGCCCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGGGAGGCGCGGGTTCCGGCCCCGGGTGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCGTTGTTTTACGGCGCGCTGACGCCCggggcgcgggtaaat                                                                                                                                                 >  2:546163/1‑239 (MQ=255)

                                                                                             tctGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGGAGGCCAACGTTCATCGCCAGGTCGGGGCGGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCCGGGATGCGGTG‑AAgg                                                                                                                                                >  2:497128/1‑250 (MQ=255)

                                                                                             tctGCGGTCCGAAACGCACCGAACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGACACGGGAGCCAGTGAAAGCTGCGGATGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGACGAGCGGCTGGGGAGGCGCGGGTTCAGGACCCGGATTCCAAAGTTCCTCGCCAGGTCGTGGCTGCCGGGGCCGGTGTTTTCCCGCGCGCTGACGCCAGGGATGCGGTT‑AAgg                                                                                                                                                >  2:347742/1‑250 (MQ=255)

                                                                                              ctGCGGTCCGTAACGCACCGCCCCACCACCGCGTACCACGGCAATGGCATCCAGGTGGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGACCATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGTTCGTGGCTGCCGGTGCCGTGGTTTTCCGGCGCGCTGACGACAGGTATGCTGTT‑AAg                                                                                                                                                 <  2:74140/249‑1 (MQ=255)

                                                                                                ggggTCCGCAACGCACCGCTCCCCCTCCGCGTACCACGTCAATGGCACCCATGTTGCCGAGCGAAACGGGTGTCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGGCCAGCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTT‑AAGTAc                                                                                                                                              <  2:544300/248‑1 (MQ=255)

                                                                                                gCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTT‑a                                                                                                                                                   >  1:667171/1‑245 (MQ=255)

                                                                                                gCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTT‑AAGTa                                                                                                                                               <  1:142397/249‑1 (MQ=255)

                                                                                                gCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTT‑AAGTAc                                                                                                                                              <  1:37086/250‑1 (MQ=255)

                                                                                                  ggTCCGTATCGCAACGCACCACCCCCGCGTACCACGGCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGGTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTT‑AAGTACCTc                                                                                                                                           <  2:248041/251‑1 (MQ=255)

                                                                                                     ccGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTT‑AAGTACCTCAc                                                                                                                                         <  1:167167/250‑1 (MQ=255)

                                                                                                         aaCGCACCGCGCCCCCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTT‑AAGTACCTCACGCAt                                                                                                                                     <  1:433901/250‑1 (MQ=255)

                                                                                                         aaCGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTT‑AAGTACCTCCCGCAtg                                                                                                                                    >  2:526241/1‑251 (MQ=255)

                                                                                                              cccgcaccACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCAGCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTT‑TAGTACCTCACGCATggtggt                                                                                                                               <  1:513825/250‑1 (MQ=255)

                                                                                                                cgcaccACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGc                                                                                                                                                                                          >  1:412789/1‑191 (MQ=255)

                                                                                                                  caccaccacCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTT‑AAGTACCTCACGCATGGTGGTTgcg                                                                                                                           >  1:18541/1‑251 (MQ=255)

                                                                                                                  caccaccacCGCGAACCAGGTCAATGGCATCCATGGTGCCGGGCGAAACGGGAGCTAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTT‑AAGTACCTCACGAATGGTGGTTgcg                                                                                                                           <  2:612704/251‑1 (MQ=255)

                                                                                                                              ggaCCACGTCAATGGCATCCATTTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCATGTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTT‑aa                                                                                                                                                  <  2:583023/214‑1 (MQ=255)

                                                                                                                                 ccACGTCACTGACATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCACCAGGACGGTCGGGCGGCTGGTGAGGCGCGGGTTCAGGGCCCGGATGCCAACGTTGTTCGCCAGGTCATGGCAGCTGTGGCCGTTGCTCTCCTGTGCCCGGTCGCGAGGTATGCGGTT‑AAGTACTTCTCCAATGGTGGTTGCGCCGGGTTTAGCGaa                                                                                                             >  2:550427/1‑250 (MQ=255)

                                                                                                                                   ccGTCAATCGCATCCATGGTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCAGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGGTCAGGCCCCGGATGCCAAAATTCATCGCCAGGTCGTGGCAGGCCGCGCCGTTGTTTTCCGGAGCCGTGCCGCCCGGGATGCGGTTCAAG‑ACCTCACGCATGTTGGTTGCGCCGGATTTGGCGAAAt                                                                                                           >  2:438002/2‑250 (MQ=255)

                                                                                                                                   aCGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTT‑AAGTACCTCACGCATGGTGGTTGCGCCGGTTTTGGCGAAATc                                                                                                          >  1:165392/1‑251 (MQ=255)

                                                                                                                                    cGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTT‑AAGTACCTCACGCATGGTGGTTGCGCCGGTTTTGGCGAAATcc                                                                                                         <  1:326364/251‑1 (MQ=255)

                                                                                                                                     gTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGCAAGCTGCGGCTAACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGGCTGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGAGCCGTTGTTTTCCGGCGCGCCGACGCCAGGGAAGCGGTT‑AAGTACCTCACGCATGGTGGTTGCGCCGGTTTTGGCGACATCCt                                                                                                        >  2:625313/1‑251 (MQ=255)

                                                                                                                                       cAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTT‑AAGTACCTCACGCATGGTGGTTGCGCCGGTTTTGGCGAAATCCTc                                                                                                       <  1:299434/250‑1 (MQ=255)

                                                                                                                                       cAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTT‑AAGTACCTCACGCATGGTGGTTGCGCCGGGTTTGGCGAAATCCTc                                                                                                       >  1:377031/1‑250 (MQ=255)

                                                                                                                                        aaTGGCATCCAGGTTGCCGAGCGAAACGGGCGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGAGGCCGTCCATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCGTTGTGTTCCGGCGCGCTGACGCCAGGGATGCGGTT‑AAGTACCTCACGCATGGTGGTTGCGCCGGTTTTGGCGAAATCCTCAc                                                                                                     <  2:261548/251‑1 (MQ=255)

                                                                                                                                           ggCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGAGCTGACGCCCGGGATGCGGTT‑AAGTACCTAACGCATGGTGGTTGCGCCGGTTTTGGCGAAATCCTCCcggc                                                                                                  >  1:429795/1‑251 (MQ=255)

                                                                                                                                               tCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGGGAGGCGCGGGTTCAGGCCCCGGATGCCACAGTTCATCGCCAGGTCGTGGCCGCCGGGGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGCT‑ACGGACCACACGCATGGTGGTTGCGCCCGTTTTGGCGAAATCCTCCcggcgggc                                                                                              >  2:111211/1‑249 (MQ=255)

                                                                                                                                                 cATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTT‑AAGTACCTCACGCATGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGAtc                                                                                             >  1:165907/1‑250 (MQ=255)

                                                                                                                                                      ggCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGTTGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTT‑AAGTACCTCACGCATGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGAtcacgtc                                                                                        <  1:673565/249‑1 (MQ=255)

                                                                                                                                                       gCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCGGGTGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTT‑AAGTACCTCACGCATGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGAtcacgtcac                                                                                      <  2:591602/251‑1 (MQ=255)

                                                                                                                                                        ccGAGCGACACGGAAGCCAGTGAAAGCTGCGGCTGACCGGAAGGGGCGAAGGGGGAGGGGATGCCGTCCATCAAGACGGTCGAGCGGCTGGAGAGGCGCGGGTTCCAGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCCGACGCCAGGGTTGCGGTT‑AAGAACATCACGCATGGTGGTTGAGCCCGTTTTGGCGAAATCCTCACGGCGGagtacgtacgc                                                                                     >  1:95650/1‑241 (MQ=255)

                                                                                                                                                          gaccGAACCGGGAGCCAGTTAAAGCTGTGGTGGACCGTAAGGGGCGAAGGGTACGGGGATGCCGTCCATCAGGACGGTCCGGCGGCTGTTGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCATGGCGGGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGATGCTGTT‑AAGTACCTCACGCTTGGTGGTTGCTCCGGTTTTGGCGAAATCCTCACGGCGGAGCACGTCAcgcg                                                                                   <  2:680906/248‑1 (MQ=255)

                                                                                                                                                            gCGAAGCGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAGGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTGCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTT‑AAGTACCTCACGCATGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGATCACGTCAcgcgc‑g                                                                                <  2:73134/251‑1 (MQ=255)

                                                                                                                                                             cGAAACGGGAGCCAGTGACAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCCGGTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGAAGCCCGGGCTGCGGTT‑AAGTACCTCAAGCATGGTGGTTGCGCCGGTTTTGGGGCAATCCTCAccggggagaccggcacaccatc                                                                                >  2:451447/1‑230 (MQ=255)

                                                                                                                                                                aaCGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTT‑AAGTACCTCACGCATGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGATCACGTCACGCGC‑GCCAg                                                                            <  1:231054/251‑1 (MQ=255)

                                                                                                                                                                   gggAGCTAGTGTAAGCGGCGGATGACCGTAAGGGGCGAAGGGGCAGGGGATGCCGTCCATCAGGACGTTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCTGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGTTGCCGTTGTTTTCCGTCGCGCTGACGCCAGGGATGCGGTT‑AAGTACCTCACGCATGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACTGCGGATCACGTCACGCTC‑GCCAGCAt                                                                         <  2:624064/251‑1 (MQ=255)

                                                                                                                                                                   gggAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTT‑AAGTACCTCACGCATGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGATCACGTCACGCGC‑GCCAGCAt                                                                         <  1:505773/251‑1 (MQ=255)

                                                                                                                                                                    ggAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTT‑AAGTACCTCACGCATGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGATCACGTCACGCGC‑GCCAGCATg                                                                        <  2:314569/251‑1 (MQ=255)

                                                                                                                                                                      agacAGTGAAAGCGGCGGCTGGCGGTAAGGGGCGAGGGGGACGGGGATGCCGGCCATCAGGACGGTCGAGCGGCGGGGGAGGCTCGGGTGTAGGCCCCGGATGACAATGTTCATCGCCAGGTCGTGGCTGCCGGGGCCGTTTGTGTCCGGCGCGCTGACGCCAGGGATGCGGTT‑AAGTACCTCACGCAAGGTGGTTGCGCCGGTTTTGGCGAAAGCCTCACTGCGGATCACGTCACGCGC‑GCCAGCATGtt                                                                      <  2:563344/248‑1 (MQ=255)

                                                                                                                                                                         cAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCGTTGTTTTCCTGCGCGCTGACGCCAGGGATGCGGTT‑AAGTACCTCACGCCTGGTGGTTGCGCCGGTTTTGGCGAAATCCGCACGGCGGATCACGTCAC‑CGCCGCCAGCATGTTCa                                                                    >  1:236835/1‑250 (MQ=255)

                                                                                                                                                                          aGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTT‑AAGTACCTCACGCATGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGATCACGTCACGCGC‑GCCCGCATGGTCaa                                                                   >  1:132593/1‑250 (MQ=255)

                                                                                                                                                                            tGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTT‑AAGTACCTCACGCATGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGATCACGTCACGCGC‑GCCAGCATGTTCAAAt                                                                 >  1:248439/1‑250 (MQ=255)

                                                                                                                                                                               aaGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTAGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTT‑AAGTACCTCACGCATGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGATCACGTCACGCGC‑GCCCGCATGTTCAAATAcgt                                                             >  1:469005/1‑251 (MQ=255)

                                                                                                                                                                                  cTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGGACGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTT‑AAGTACCTCACGCATGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGATCACGTCACGCGC‑GCCAGCATGTTCAACTAcgtcgt                                                          >  1:268468/1‑251 (MQ=255)

                                                                                                                                                                                      ggCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTT‑AAGTACCTCACGCATGGTGGTTGCGCCGGTTTTTGCGAAATCCTCACGGCGGCTACCGTCACACGC‑GCCCGCATGTTCAAATACGTCGGTTTc                                                      >  2:120225/1‑251 (MQ=255)

                                                                                                                                                                                      ggCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTT‑AAGTACCTCACGCATGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGATCACGTCACGCGC‑GCCAGCATGTTCAAATACGTCGTTTTc                                                      <  1:92114/251‑1 (MQ=255)

                                                                                                                                                                                       gCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTT‑AAGTACCTCACGCATGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGATCACGTCACGCGC‑GCCAGCATGTTCAAATACGTCGTTTTCa                                                     <  1:358455/251‑1 (MQ=255)

                                                                                                                                                                                          gACCGTAAGGGGCTAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAGTGATGCGGTT‑AAGTACCTCACGCATGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGATCACGTCACGCGC‑GCCAGCATGTTCAAATACGTCGTTTTCAcgc                                                  <  1:649631/251‑1 (MQ=255)

                                                                                                                                                                                           aCCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTT‑AAGTACCTCACGCATGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGATCACGTCACGCGC‑GCCAGCATGTTCAAATACGTCGGTTTCCcgcg                                                 >  1:577507/1‑251 (MQ=255)

                                                                                                                                                                                            ccGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGGGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCCGGGATGCGGTT‑AAGTACCTCACGCATGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGATCACGTCCCGCCC‑GCCAGCCTGGTCCAATACGTCGTTTTCCCgcgc                                                >  1:416002/1‑251 (MQ=255)

                                                                                                                                                                                              gTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTT‑AAGTACCTCACGCATGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGATCACGTCACGCGC‑GCCAGCCTGTTCAAATACGTCGTTTTCACGCGCAt                                              >  1:543408/1‑251 (MQ=255)

                                                                                                                                                                                               tAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCAGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTT‑AAGTACCTCACGCATGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGATCACGTCACGCGC‑GCCAGCATGTTCAAATACGTCGTTTTCACGCGCATc                                             >  1:592610/1‑251 (MQ=255)

                                                                                                                                                                                                  ggggggAAGGGGACGGGGATGCCGGCCACCAGGGCGGGCGTGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGTATGCCAAAGTTCATCGCCAGGTCGTTGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTT‑AAGTACCTCACGCATGGTGGTTGCGCCGGTTTTGTCGAAATCCTCACGGCGGATCACGTCACGCGC‑GCCAGCATGTTCAAATACGTCGTTTTCACGCGCATCAcc                                          <  2:128186/251‑1 (MQ=255)

                                                                                                                                                                                                        aaGGGGACGGGTATGCCGTCCATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTT‑AAGTACCTCACGCATGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGTCGGATCACGTCACGCGC‑GCCAGCATGTTCAAATACGGCGTTTTCACGCGCATCACccagcca                                    <  1:201009/251‑1 (MQ=255)

                                                                                                                                                                                                         aGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTT‑AAGTACCTCACGCATGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGATCACGTCACGCGC‑GCCAGCCTGTTCAAATACGTCGTTTTCACGCGCATCACccagcc                                     >  1:105231/1‑249 (MQ=255)

                                                                                                                                                                                                         aGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTCCATCGCCAGGTCGTGGCTGCCGGTGCCGTTGTTTTTCGGCGCGCTGACGCCAGGGAt                                                                                                                                                           <  1:402191/133‑1 (MQ=255)

                                                                                                                                                                                                                 gggATGCCGTCCATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTT‑AAGTACCTCACGCATGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGATCACGTCACGCAC‑GCCAGCCAGTTCAAATACGTCGATTTCACGCGCCTCATCCAGCCAGTCGCCGAc                           >  1:289928/1‑251 (MQ=255)

                                                                                                                                                                                                                 gggATGCCGGCCATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCTAGGTCGTGGCTGCCGGTGCCGTTGGTTTCCGGCGTGCTGACGCCAGGGATGCGGTT‑AAGTACCTCACGCATGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGATCACGTCACGGGC‑GTCAGCATGTTCAAATACGTCGTTTTCACGCGCATCACCCAGCCAGTCGCCGAc                           <  1:390656/251‑1 (MQ=255)

                                                                                                                                                                                                                  ggATGCCGTCCATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTT‑AAGTACCTCACGCATGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGATCACGTCACGCGC‑GCCAGCATGTTCAAATACGTCGTTTTCACGCGCATCACCCAGCCAGTCGCCGAcc                          <  1:5002/251‑1 (MQ=255)

                                                                                                                                                                                                                    aTGCCGTCCATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTT‑AAGTACCTCACGCATGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGATCACGTCACGCGC‑GCCAGCATGTTCAAATACGTCGTTTTCACGCGCATCACCCAGCCAGTCGCCGACCAc                        >  1:97235/1‑251 (MQ=255)

                                                                                                                                                                                                                      gCCGTCCATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTT‑AAGTACCTCACGCATGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGATCACGTCACGCGC‑GCCAGCATGTTCAAATACGTCGTTTTCACGCGCATCACCCAGCCCGTCGCCGGCCACgg                      >  1:263652/1‑251 (MQ=255)

                                                                                                                                                                                                                        cGTCCATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTT‑AAGTACCTCACGCATGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGATCACGTCACACGC‑GCCAGCCTGTTCAAATACGTCGTTTTCACGCGCATCACCCAGCCAGACGCCGACCACGGTc                    >  1:534329/1‑251 (MQ=255)

                                                                                                                                                                                                                         gtacctCAGGACGGGAGAGTGGCTGGTGAGGCGCGGGTGCAGTCCCGGGGTGCCAAAGGTCATCGCAAGGTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCGGACGCTAGGGATGCGGTT‑AAGGGCCTCACGCATGGTGGTCGCGCCGGTTTTTGCGAAATCCTCACGGCGGATCACGTCACGCGC‑GCCAGCATGTTCAGATACGTCGTTTTCACGCGCATCACCCAGCCAGTCGCCGACCACGGTCa                   <  2:416002/246‑1 (MQ=255)

                                                                                                                                                                                                                          gccATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTT‑AAGTACCTCACGCATGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGATCACGTCACGCGC‑GCCAGCATGTTCAAATACGTCGTTTTCACGCGCATCACCCAGCCAGTCGCCGACCACGGTc                    <  1:309713/248‑1 (MQ=255)

                                                                                                                                                                                                                           ccATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTT‑AAGTACCTCACGCATGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGATCACGTCACGCGC‑GCCAGCATGTTCAAATACGTCGTTTTCACGCGCCTCACCCAGCCAGTCGCCGACCACGGTCAgg                 >  1:613719/1‑251 (MQ=255)

                                                                                                                                                                                                                           ccATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTT‑AAGTACCTCACGCATGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGATCACGTCACGCGC‑GCCAGCATGTTCAAATACGTCGTTTTCACGCGCATCCCCCAGCCAGTCGCCGACCACGGTCAgg                 >  1:546397/1‑251 (MQ=255)

                                                                                                                                                                                                                           ccATCAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTT‑AAGTACCTCACGCATGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGATCACGTCACGCGC‑GCCAGCATGTTCAAATACGTCGTTTTCACGCGCATCACCCAGCCAGTCGCCGACCACGGTCAgg                 <  1:216500/251‑1 (MQ=255)

                                                                                                                                                                                                                               cAGGACGGTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTT‑AAGTACCTCACGCATGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGATCACGTCACGCGC‑GCCAGCATGTTCAAATACGTCGTTTTCACGCGCATCACCCAGCCAGTCGCCGACCACGGTCAGGGCAt             <  1:37646/251‑1 (MQ=255)

                                                                                                                                                                                                                                     ggTCGAGCGGCTGGTGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTT‑AAGTACCTCACGCATGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGATCACGTCACGCGC‑GCCAGCATGTTTAAATACGTCGTTTTCACGCGCCACCCCCCGCCAGTCGCCGACCACGGGCAGGGCATcttctt       >  1:460074/1‑251 (MQ=255)

                                                                                                                                                                                                                                          ggCGGCTGGGGAGGCGCGGGGTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGTCGTGGTTGTCCGGCGCGCGGACGCCAGGGATGCGGTT‑AAGTACCTCACGCATGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGATCACGTCCCGCGC‑GCCAGCATGTTCAAATACGTCGTTTTCACGCGCATCTCACAGCCAGTCGCCGACCACGGTCAGGGCATCTTCTTTTGGt  <  2:248439/250‑1 (MQ=255)

AAGAGGTTGGGCTGAGCTGACCTTCCACGCCCGCCTCGATACCAAAGTCCTGCGGAATGGCACGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGCGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCTGACGCCAGGGATGCGGTT‑AAGTACCTCACGCATGGTGGTTGCGCCGGTTTTGGCGAAATCCTCACGGCGGATCACGTCACGCGC‑GCCAGCATGTTCAAATACGTCGTTTTCACGCGCATCACCCAGCCAGTCGCCGACCACGGTCAGGGCATCTTCTTTTGGT  >  NC_000913/4515867‑4516351

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