BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	14,850	C→T	6.2%	H507H (CAC→CAT)	gyrB →	DNA gyrase subunit B

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	14,850	0	C	T	6.2%	761.4 / 25.1	243	H507H (CAC→CAT)	gyrB	DNA gyrase subunit B
Reads supporting (aligned to +/- strand): ref base C (109/119); new base T (8/7); total (117/126)
Fisher's exact test for biased strand distribution p-value = 7.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 245 total aligned reads displayed.

GCGCTTCGACAAGATGATTTCGTCCCAGGAAGTGGGCACGCTGATCACTGCGCTGGGCTGTGGCATCGGCCGCGAAGAGTACAACATCGACAAACTGCGTTATCACAACATCATCATCATGACCGATGCTGACGTTGACGGTTCGCACATCCGTACGCTGCTGCTGACCTTCTTCTTCCGTCAGCTGCCGGAGCTGGTCGAGCGTGGCTACATCTATATTGCCCAGCCGCCGCTGTACAAGGTGAAACGAGGCAAGCAGGAGCAGTACATCAAGGACGATGAGGCCATGGAAGAG  >  NC_002947/14703‑14997

gcgcTTCGACAAGATGATTTCGTCCCAGGAAGTGGGCACGCTGATCACTGCGCTGGGCTGTGGCATCGGCCGCGAAGAGTACAACATCGACAAACTGCGTTATCACAACATCATCATCATGACCGATGCTGACGTTGACGGTTCGCACAt                                                                                                                                                   <  1:4762567/150‑1 (MQ=255)

    ttCGACAAGATGATTTCGTCCCAGGAAGTGGGCACGCTGATCACTGCGCTGGGCTGTGGCATCGGCCGCGAAGAGTACAACATCGACAAACTGCGTTATCACAACATCATCATCATGACCGATGCTGACGTTGACGGTTCGCACATCCGt                                                                                                                                               >  2:1414151/1‑150 (MQ=255)

     tCGACAAGATGATTTCGTCCCAGGAAGTGGGCACGCTGATCACTGCGCTGGGCTGTGGCATCGGCCGCGAAGAGTACAACATCGACAAACTGCGTTATCACAACATCATCATCATGACCGATGCTGACGTTGACGGTTCGCACATCCGTa                                                                                                                                              <  2:4748285/150‑1 (MQ=255)

           aGATGATTTCGTCCCAGGAAGTGGGCACGCTGATCACTGCGCTGGGCTGTGGCATCGGCCGCGAAGAGTACAACATCGACAAACTGCGTTATCACAACATCATCATCATGACCGATGCTGACGTTGACGGTTCGCACATCCGTACgctgc                                                                                                                                        >  1:2927543/1‑150 (MQ=255)

               gatTTCGTCCCAGGAAGTGGGCACGCTGATCACTGCGCTGGGCTGTGGCATCGGCCGCGACGAGTACAACATCGACACCCTGCGTTATCACAACATCATCATCAGGACCGATGCTGACGTTGACGGTTCGCACATCCGTACgctgcggct                                                                                                                                    <  2:4263541/150‑1 (MQ=255)

                atTTCGTCCCAGGAAGTGGGCACGCTGATCACTGCGCTGGGCTGTGGCATCGGCCGCGAAGAGTACAACATCGACAAACTGCGTTATCACAACATCATCATCATGACCGATGCTGACGTTGACGGTTCGCACATCCGTACGctgctgctg                                                                                                                                   <  2:9081341/150‑1 (MQ=255)

                atTTCGTCCCAGGAAGTGGGCACGCTGATCACTGCGCTGGGCTGTGGCATCGGCCGCGAAGAGTACAACATCGACAAACTGCGTTATCACAACATCATCATCATGACCGATGCTGACGTTGACGGTTCGCACATCCGTACGctgctgctg                                                                                                                                   <  1:3924239/150‑1 (MQ=255)

                 tttCGTCCCAGGAAGTGGGCACGCTGATCACTGCGCTGGGCTGTGGCATCGGCCGCGAAGAGTACAACATCGACAAACTGCGTTATCACAACATCATCATCATGACCGATGCTGACGTTGACGGTTCGCACATCCGTACGCTGCTGCTGa                                                                                                                                  >  2:1957563/1‑150 (MQ=255)

                  ttCGTCCCAGGAAGTGGGCACGCTGATCACTGCGCTGGGCTGTGGCATCGGCCGCGAAGAGTACAACATCGACAAACTGCGTTATCACAACATCATCATCATGACCGATGCTGACGTTGACGGTTCGCACATCCGTACGCTGCTGCTGAc                                                                                                                                 >  2:5892875/1‑150 (MQ=255)

                      tCCCAGGAAGTGGGCACGCTGATCACTGCGCTGGGCTGTGGCATCGGCCGCGAAGAGTACAACATCGACAAACTGCGTTATCACAACATCATCATCATGACCGATGCTGACGTTGACGGTTCGCACATCCGTACGCTGCTGCTGACcttc                                                                                                                             <  2:8657999/150‑1 (MQ=255)

                        ccAGGAAGTGGGCACGCTGATCACTGCGCTGGGCTGTGGCATCGGCCGCGAAGAGTACAACATCGACAAACTGCGTTATCACAACATCATCATCATGACCGATGCTGACGTTGACGGTTCGCACATCCGTACGCTGCTGCTGACcttctt                                                                                                                           >  1:241018/1‑150 (MQ=255)

                        ccAGGAAGTGGGCACGCTGATCACTGCGCTGGGCTGTGGCATCGGCCGCGAAGAGTACAACATCGACAAACTGCGTTATCACAACATCATCATCATGACCGATGCTGACGTTGACGGTTCGCACATCCGTACGCTGCTGCTGACcttctt                                                                                                                           >  1:2035271/1‑150 (MQ=255)

                         cacgAAGTGCGCACGCTGATCACTGCGCTGCGCTGTGCCATCGCCCGCGAACACTACAACATCGACAAACTGCGTTATCACAACAACAACATCATGACCGATGCTGACGTTGACGGTTCGCACATCCGTACGCTACTGCTGACcttcttc                                                                                                                          <  2:4574647/147‑1 (MQ=255)

                         cAGGAAGTGGGCACGCTGATCACTGCGCTGGGCTGTGGCATCGGCCGCGAAGAGTACAACATCGACAAACTGCGTTATCACAACATCATCATCATGACCGATGCTGACGTTGACGGTTCGCACATCCGTACGCTGCTGCTGACcttcttc                                                                                                                          <  1:3691971/150‑1 (MQ=255)

                          aGGAAGTGGGCACGCTGATCACTGCGCTGGGCTGTGGCATCGGCCGCGAAGAGTACAACATCGACAAACTGCGTTATCACAACATCATCATCATGACCGATGCTGACGTTGACGGTTCGCACATCCGTACGCTGCTGCTGACcttcttct                                                                                                                         <  1:1869106/150‑1 (MQ=255)

                             aaGTGGGCACGCTGATCACTGCGCTGGGCTGTGGCATCGGCCGCGAAGAGTACAACATCGACAAACTGCGTTATCACAACATCATCATCATGACCGATGCTGACGTTGACGGTTCGCACATCCGTACGCTGCTGCTGACCTTCTTCTTcc                                                                                                                      <  2:8061398/150‑1 (MQ=255)

                              aGTGGGCACGCTGATCACTGCGCTGGGCTGTGGCATCGGCCGCGAAGAGTACAACATCGACAAACTGCGTTATCACAACATCATCATCATGACCGATGCTGACGTTGACGGTTCGCACATCCGTACGCTGCTGCTGACctt                                                                                                                              <  1:7416311/141‑1 (MQ=255)

                                   gCACGCTGATCACTGCGCTGGGCTGTGGCATCGGCCGCGAAGAGTACAACATCGACAAACTGCGTTATCACAACATCATCATCATGACCGATGCTGACGTTGACGGTTCGCACATCCGTACGCTGCTGCTGACCTTCTTCTTCCGTCAGc                                                                                                                <  2:3012211/150‑1 (MQ=255)

                                     aCGCTGATCACTGCGCTGGGCTGTGGCATAGGCCGCGAAGAGTACAACATCGACAAACTGCGTTATCACAACATCATCATCATGACCGATGCTGACGTTGAAGGTTCGCACATCCGTACGCTGCTGCTGACCTTCTTCTTCCGTCAGCTg                                                                                                              >  1:8259798/1‑150 (MQ=255)

                                       gCTGATCACTGCGCTGGGCTGTGGCATCGGCCGCGAAGAGTACAACATCGACAAACTGCGTTATCACAACATCATCATCATGACCGATGCTGACGTTGACGGTTCGCACATCCGTACGCTGCTGCTGACCTTCTTCTTCCGTCAGCTGcc                                                                                                            >  1:6417192/1‑150 (MQ=255)

                                       gCTGATCACTGCGCTGGGCTGTGGCATCGGCCGCGAAGAGTACAACATCGACAAACTGCGTTATCACAACATCATCATCATGACCGATGCTGACGTTGACGGTTCGCACATCCGTACGCTGCTGCTGACCTTCTTCTTCCGTCAGCTGcc                                                                                                            >  1:754299/1‑150 (MQ=255)

                                         tGATCACTGCGCTGGGCTGTGGCATCGGCCGCGAAGAGTACAACATCGACAAACTGCGTTATCACAACATCATCATCATGACCGATGCTGACGTTGACGGTTCGCACATCCGTACGCTGCTGCTGACCTTCTTCTTCCGTCAGCTGCCgg                                                                                                          <  1:6171190/150‑1 (MQ=255)

                                          gATCACTGCGCTGGGCTGTGGCATCGGCCGCGAAGAGTACAACATCGACAAACTGCGTTATCACAACATCATCATCATGACCGATGCTGACGTTGACGGTTCGCACATCCGTACGCTGCTGCTGACCTTCTTCTTCCGTCAGCTGCCGGa                                                                                                         <  1:7929875/150‑1 (MQ=255)

                                             cACTGCGCTGGGCTGTGGCATCGGCCGCGAAGAGTACAACATCGACAAACTGCGTTATCACAACATCATCATCATGACCGATGCTGACGTTGACGGTTCGCACATCCGTACGCTGCTGCTGACCTTCTTCTTCCGTCAGCTGCCGGAGCt                                                                                                      >  1:1181980/1‑150 (MQ=255)

                                             cACTGCGCTGGGCTGTGGCATCGGCCGCGAAGAGTACAACATCGACAAACTGCGTTATCACAACATCATCATCATGACCGATGCTGACGTTGACGGTTCGCACATCCGTACGCTGCTGCTGACCTTCTTCTTCCGTCAGCTGCCGGAGCt                                                                                                      <  2:8021318/150‑1 (MQ=255)

                                               cTGCGCTGGGCTGTGGCATCGGCCGCGAAGAGTACAACATCGACAAACTGCGTTATCACAACATCATCATCATGACCGATGCTGACGTTGACGGTTCGCACATCCGTACGCTGCTGCTGACCTTCTTCTTCCGTCAGCTGCCGGAGCTgg                                                                                                    >  2:5313370/1‑150 (MQ=255)

                                                       ggCTGTGGCATCGGCCGCGAAGAGTACAACATCGACAAACTGCGTTATCACAACATCATCATCATGACCGATGCTGACGTTGACGGTTCGCACATCCGTACGCTGCTGCTGACCTTCTTCTTCCGTCAGCTGCCGGAGCTGGTCGAGCGt                                                                                            >  2:1788058/1‑150 (MQ=255)

                                                        gCTGTGGCATCGGCCGCGAAGAGTACAACATCGACAAACTGCGTTATCACAACATCATCATCATGACCGATGCTGACGTTGACGGTTCGCACATCCGTACGCTGCTGCTGACCTTCTTCTTCCGTCAGCTGCCGGAGCTGGTCGAGCGTg                                                                                           <  1:1656926/150‑1 (MQ=255)

                                                        gCTGTGGCATCGGCCGCGAAGAGTACAACATCGACAAACTGCGTTATCACAACATCATCATCATGACCGATGCTGACGTTGACGGTTCGCACATCCGTACGCTGCTGCTGACCTTCTTCTTCCGTCAGCTGCCGGAGCTGGTCGAGCGTg                                                                                           >  2:1447800/1‑150 (MQ=255)

                                                         cTGTGGCATCGGCCGCGAAGAGTACAACATCGACAAACTGCGTTATCACAACATCATCATCATGACCGATGCTGACGTTGACGGTTCGCACATCCGTACGCTGCTGCTGACCTTCTTCTTCCGTCAGCTGCCGGAGCTGGTCGAGCGTgg                                                                                          <  2:4224887/150‑1 (MQ=255)

                                                             ggCATCGGCCGCGAAGAGTACAACATCGACAAACTGCGTTATCACAACATCATCATCATGACCGATGCTGACGTTGACGGTTCGCACATCCGTACGCTGCTGCTGACCTTCTTCTTCCGTCAGCTGCCGGAGCTGGTCGAGCGTGGCTAc                                                                                      <  1:597973/150‑1 (MQ=255)

                                                              gCATCGGCCGCGAGGAGTACAACATCGACAAACTGCGTTATCACAACATCATCATCATGACCGATGCTGACGTTGACGGTTCGCACATCCGTACGCTGCTGCTGACCTTCTTCTTCCGTCAGCTGCCGGAGCTGGTCGAGCGTGGCTACa                                                                                     >  2:6841266/1‑150 (MQ=255)

                                                                aTCGGCCGCGAAGAGTACAACATCGACAAACTGCGTTATCACAACATCATCATCATGACCGATGCTGACGTTGACGGTTCGCACATCCGTACGCTGCTGCTGACCTTCTTCTTCCGTCAGCTGCCGGAGCTGGTCGAGCGTGGCTACATc                                                                                   <  2:1597272/150‑1 (MQ=255)

                                                                           agagTACAACATCGACAAACTGCGTTATCACAACATCATCATCATGACCGATGCTGACGTTGACGGTTCGCACATCCGTACGCTGCTGCTGACCTTCTTCTTCCGTCAGCTGCCGGAGCTGGTCGAGCGTGGCTACATCTATATTGCCCa                                                                        <  2:6335927/150‑1 (MQ=255)

                                                                               tACAACATCGACAAACTGCGTTATCACAACATCATCATCATGACCGATGCTGACGTTGACGGTTCGCACATCCGTACGCTGCTGCTGACCTTCTTCTTCCGTCAGCTGCCGGAGCTGGTCGAGCGTGGCTACATCTATATTGCCCAgccg                                                                    >  2:6495101/1‑150 (MQ=255)

                                                                                  aacaTCGACAAACTGCGTTATCACAACATCATCATCATGACCGATGCTGACGTTGACGGTTCGCACATCCGTACGCTGCTGCTGACCTTCTTCTTCCGTCAGCTGCCGGAGCTGGTCGAGCGTGGCTACATCTATATTGCCCAGccgccg                                                                 >  1:3815789/1‑150 (MQ=255)

                                                                                    tcctgACAAGCTGCGCTACCATAAGATCATCATCATGACGGATGCTGACGTCGACGGTGCGCATATTCGTACGCTGCTGTTGACCTTCTTCTTCCGTCAGATGCCGGAACTGATCGAGCGTGGCCATATCTATATCGCGCAGCCGCCGCt                                                               >  1:2850332/5‑150 (MQ=255)

                                                                                    tcctgACAAGCTGCGCTACCATAAGATCATCATCATGACGGATGCTGACGTCGACGGTGCGCATATTCGTACGCTGCTGTTGACCTTCTTCTTCCGTCAGATGCCGGAACTGATCGAGCGTGGCCATATCTATATCGCGCAGCCGCCGCt                                                               >  2:4986176/5‑150 (MQ=255)

                                                                                     cctgACAAGCTGCGCTACCATAAGATCATCATCATGACGGATGCTGACGTCGACGGTGCGCATATTCGTACGCTGCTGTTGACCTTCTTCTTCCGTCAGATGCCGGAACTGATCGAGCGTGGCCATATCTATATCGCGCAGCCGCCGCtc                                                              <  2:3218730/147‑2 (MQ=255)

                                                                                     aTCGACAAACTGCGTTATCACAACATCATCATCATGACCGATGCTGACGTTGACGGTTCGCACATCCGTACGCTGCTGCTGACCTTCTTCTTCCGTCAGCTGCCGGAGCTGGTCGAGCGTGGCTACATCTATATTGCCCAGCCGCCGCTg                                                              >  2:3142622/1‑150 (MQ=255)

                                                                                     aTCGACAAACTGCGTTATCACAACATCATCATCATGACCGATGCTGACGTTGACGGTTCGCACATCCGTACGCTGCTGCTGACCTTCTTCTTCCGTCAGCTGCCGGAGCTGGTCGAGCGTGGCTACATCTATATTGCCCAGCCGCCGCTg                                                              >  1:7303616/1‑150 (MQ=255)

                                                                                      tCGACAAACTGCGTTATCACAACATCATCATCATGACCGATGCTGACGTTGACGGTTCGCACATCCGTACGCTGCTGCTGACCTTCTTCTTCCGTCAGCTGCCGGAGCTGGTCGAGCGTGGCTACATCTATATTGCCCAGCCGCCGCTGt                                                             >  2:1967393/1‑150 (MQ=255)

                                                                                        gACAAACTGCGTTATCACAACATCATCATCATGACCGATGCTGACGTTGACGGTTCGCACATCCGTACGCTGCTGCTGACCTTCTTCTTCCGTCAGCTGCCGGAGCTGGTCGAGCGTGGCTACATCTATATTGCCCAGCCGCCGCTGTAc                                                           <  1:2856331/150‑1 (MQ=255)

                                                                                          cAAGCTGCGCTACCATAAGATCATCATCATGACGGATGCTGACGTCGACGGTGCGCATATTCGTACGCTGCTGTTGACCTTCTTCTTCCGTCAGATGCCGGAACTGATCGAGCGTGGCCATATCTATATCGCGCAGCCGCCGCTCTATaa                                                         <  1:2656823/150‑1 (MQ=255)

                                                                                          cAAGCTGCGCTACCATAAGATCATCATCATGACGGATGCTGACGTCGACGGTGCGCATATTCGTACGCTGCTGTTGACCTTCTTCTTCCGTCAGATGCCGGAACTGATCGAGCGTGGCCATATCTATATCGCGCAGCCGCCGCTCTATaa                                                         <  2:7299176/150‑1 (MQ=255)

                                                                                          cAAACTGCGTTATCACAACATCATCATCATGACCGATGCTGACGTTGACGGTTCGCACATCCGTACGCTGCTGCTGACCTTCTTCTTCCGTCAGCTGCCGGAGCTGGTCGAGCGTGGCTACATCTATATTGCCCAGCCGCCGCTGTACaa                                                         <  1:144965/150‑1 (MQ=255)

                                                                                          cAAACTGCGTTATCACAACATCATCATCATGACCGATGCTGACGTTGACGGTTCGCACATCCGTACGCTGCTGCTGACCTTCTTCTTCCGTCAGCTGCCGGAGCTGGTCGAGCGTGGCTACATCTATATTGCCCAGCCGCCGCTGTACaa                                                         <  1:5018793/150‑1 (MQ=255)

                                                                                          cAAACTGCGTTATCACAACATCATCATCATGACCGATGCTGACGTTGACGGTTCGCACATCCGTACGCTGCTGCTGACCTACTTCGTCCGTCAGCTGCCGGAGCTGGTCGAGCGTGGCTACATCTATATTGCCCAGCCGCCGCTGTACaa                                                         <  1:3361349/150‑1 (MQ=255)

                                                                                            aaCTGCGTTATCACAACATCATCATCATGACCGATGCTGACGTTGACGGTTCGCACATCCGTACGCTGCTGCTGACCTTCTTCTTCCGTCAGCTGCCGGAGCTGGTCGAGCGTGGCTACATCTATATTGCCCAGCCGCCGCTGTACAAgg                                                       <  1:3259960/150‑1 (MQ=255)

                                                                                            aaCTGCGTTATCACAACATCATCATCATGACCGATGCTGACGTTGACGGTTCGCACATCCGTACGCTGCTGCTGACCTTCTTCTTCCGTCAGCTGCCGGAGCTGGTCGAGCGTGGCTACATCTATATTGCCCAGCCGCCGCTGTACAAgg                                                       >  1:301998/1‑150 (MQ=255)

                                                                                            aaCTGCGTTATCACAACATCATCATCATGACCGATGCTGACGTTGACGGTTCGCACATCCGTACGCTGCTGCTGACCTTCTTCTTCCGTCAGCTGCCGGAGCTGGTCGAGCGTGGCTACATCTATATTGCCCAGCCGCCGCTGTACAAgg                                                       <  2:4048133/150‑1 (MQ=255)

                                                                                             gcTGCGCTACCATAAGATCATCATCATGACGGATGCTGACGTCGACGGTGCGCATATTCGTACGCTGCTGTTGACCTTCTTCTTCCGTCAGATGCCGGAACTGATCGAGCGTGGCCATATCTATATCGCGCAGCCGCCGCTCTATAAGGt                                                      >  2:5704951/2‑150 (MQ=255)

                                                                                             gcTGCGCTACCATAAGATCATCATCATGACGGATGCTGACGTCGACGGTGCGCATATTCGTACGCTGCTGTTGACCTTCTTCTTCCGTCAGATGCCGGAACTGATCGAGCGTGGCCATATCTATATCGCGCAGCCGCCGCTCTATAAGGt                                                      >  2:4212119/2‑150 (MQ=255)

                                                                                               tGCGTTATCACAACATCATCATCATGACCGATGCTGACGTTGACGGTTCGCACATCCGTACGCTGCTGCTGACCTTCTTCTTCCGTCAGCTGCCGGAGCTGGTCGAGCGTGGCTACATCTATATTGCCCAGCCGCCGCTGTACAAGGTGa                                                    >  2:1019022/1‑150 (MQ=255)

                                                                                               tGCGTTATCACAACATCATCATCATGACCGATGCTGACGTTGACGGTTCGCACATCCGTACGCTGCTGCTGACCTTCTTCTTCCGTCAGCTGCCGGAGCTGGTCGAGCGTGGCTACATCTATATTGCCCAGCCGCCGCTGTACAAGGTGa                                                    >  2:7213606/1‑150 (MQ=255)

                                                                                                gcgcTACCATAAGATCATCATCATGACGGATGCTGACGTCGACGGTGCGCATATTCGTACGCTGCTGTTGACCTTCTTCTTCCGTCAGATGCCGGAACTGATCGAGCGTGGCCATATCTATATCGCGCAGCCGCCGCTCTATAAGGTgtc                                                   >  2:3890478/1‑148 (MQ=255)

                                                                                                gCGTTATCACAACATCATCATCATGACCGATGCTGACGTTGACGGTTCGCACATCCGTACGCTGCTGCTGACCTTCTTCTTCCGTCAGCTGCCGGAGCTGGTCGAGCGTGGCTACATCTATATTGCCCAGCCGCCGCTGTACAAGGTGaa                                                   >  2:3470967/1‑150 (MQ=255)

                                                                                                 cGTTATCACAACATCATCATCATGACCGATGCTGACGTTGACGGTTCGCACATCCGTACGCTGCTGCTGACCTTCTTCTTCCGTCAGCTGCCGGAGCTGGTCGAGCGTGGCTACATCTATATTGCCCAGCCGCCGCTGTACAAGGTGaaa                                                  <  2:863001/150‑1 (MQ=255)

                                                                                                 cGTTATCACAACATCATCATCATGACCGATGCTGACGTTGACGGTTCGCACATCCGTACGCTGCTGCTGACCTTCTTCTTCCGTCAGCTGCCGGAGCTGGTCGAGCGTGGCTACATCTATATTGCCCAGCCGCCGCTGTACAAGGTGaaa                                                  <  1:7171441/150‑1 (MQ=255)

                                                                                                  gatATCACAACATCAAAATCATGACCGATGCTGACGTTGACGGTTCGCACATCCGTACGCTGCTGCTGACCTTCTTCTTCCGTCAGCTGCAGGAGCTGGTCGAGCGTGGCTACATCTATATTGCCCAGCCGCCGCTGTACAAGGTGAAAc                                                 >  2:2260637/3‑150 (MQ=255)

                                                                                                  gTTATCACAACATCATCATCATGACCGATGCTGACGTTGACGGTTCGCACATCCGTACGCTGCTGCTGACCTTCTTCTTCCGTCAGCTGCCGGAGCTGGTCGAGCGTGGCTACATCTATATTGCCCAGCCGCCGCTGTACAAGGTGAAAc                                                 <  2:5566960/150‑1 (MQ=255)

                                                                                                  gTTATCACAACATCATCATCATGACCGATGCTGACGTTGACGGTTCGCACATCCGTACGCTGCTGCTGACCTTCTTCTTCCGTCAGCTGCCGGAGCTGGTCGAGCGTGGCTACATCTATATTGCCCAGCCGCCGCTGTACAAGGTGAAAc                                                 >  1:6389589/1‑150 (MQ=255)

                                                                                                    tATCACAACATCATCATCATGACCGATGCTGACGTTGACGGTTCGCACATCCGTACGCTGCTGCTGACCTTCTTCTTCCGTCAGCTGCCGGAGCTGGTCGAGCGTGGCTACATCTATATTGCCCAGccgccg                                                                 >  2:3697209/1‑132 (MQ=255)

                                                                                                    tATCACAACATCATCATCATGACCGATGCTGACGTTGACGGTTCGCACATCCGTACGCTGCTGCTGACCTTCTTCTTCCGTCAGCTGCCGGAGCTGGTCGAGCGTGGCTACATCTATATTGCCCAGCCGCCGCTGTACAAGGTGAAACGa                                               >  2:5928160/1‑150 (MQ=255)

                                                                                                    tATCACAACATCATCATCATGACCGATGCTGACGTTGACGGTTCGCACATCCGTACGCTGCTGCTGACCTTCTTCTTCCGTCAGCTGCCGGAGCTGGTCGAGCGTGGCTACATCTATATTGCCCAGCCGCCGCTGTACAAGGTGAAACGa                                               <  1:8407378/150‑1 (MQ=255)

                                                                                                     aTCACAACCTCATCATCATGACCGATGCTGACGTTGACGGTTCGCACATCCGTACGCTGCTGCTGACCTTCTTCTTCCGTCAGCTGCCGGAGCTGGTCGAGCGTGGCTACATCTATATTGCCCAGCCGCCGCTGTACAAGGTGAAACGAg                                              <  2:8799652/150‑1 (MQ=255)

                                                                                                     aTCACAACATCATCATCATGACCGATGCTGACGTTGACGGTTCGCACATCCGTACGCTGCTGCTGACCTTCTTCTTCCGTCAGCTGCCGGAGCTGGTCGAGCGTGGCTACATCTATATTGCCCAGCCGCCGCTGTACAAGGTGAAACGAg                                              <  2:238976/150‑1 (MQ=255)

                                                                                                     aTCACAACATCATCATCATGACCGATGCTGACGTTGACGGTTCGCACATCCGTACGCTGCTGCTGACCTTCTTCTTCCGTCAGCTGCCGGAGCTGGTCGAGCGTGGCTACATCTATATTGCCCAGCCGCCGCTGTACAAGGTGAAACGAg                                              <  2:849877/150‑1 (MQ=255)

                                                                                                       cacaACATCATCATCATGACCGATGCTGACGTTGACGGTTCGCACATCCGTACGCTGCTGCTGACCTTCTTCTTCCGTCAGCTGCCGGAGCTGGTCGAGCGTGGCTACATCTATATTGCCCAGCCGCCGCTGTACAAGGTGAAACGAGGc                                            <  1:5102739/150‑1 (MQ=255)

                                                                                                        acaacaTCATCATCATGACCGATGCTGACGTTGACGGTTCGCACATCCGTACGCTGCTGCTGACCTTCTTCTTCCGTCAGCTGCCGGAGCTGGTCGAGCGTGGCTACATCTATATTGCCCAGCCGCCGCTGTACAAGGTGAAACGAGGCa                                           >  2:8522073/1‑150 (MQ=255)

                                                                                                        acaacaTCATCATCATGACCGATGCTGACGTTGACGGTTCGCACATCCGTACGCTGCTGCTGACCTTCTTCTTCCGTCAGCTGCCGGAGCTGGTCGAGCGTGGCTACATCTATATTGCCCAGCCGCCGCTGTACAAGGTGAAACGAGGCa                                           >  2:5309581/1‑150 (MQ=255)

                                                                                                        acaacaTCATCATCATGACCGATGCTGACGTTGACGGTTCGCACATCCGTACGCTGCTGCTGACCTTCTTCTTCCGTCAGCTGCCGGAGCTGGTCGAGCGTGGCTACATCTATATTGCCCAGCCGCCGCTGTACAAGGTGAAACGAGGCa                                           >  1:3115567/1‑150 (MQ=255)

                                                                                                           acaTCATCATCATGACCGATGCTGACGTTGACGGTTCGCACATCCGTACGCTGCTGCTGACCTTCTTCTTCCGTCAGCTGCCGGAGCTGGTCGAGCGTGGCTACATCTATATTGCCCAGCCGCCGCTGTACAAGGTGAAACGAGGCAAGc                                        >  1:6810256/1‑150 (MQ=255)

                                                                                                            catcatcatcatGACCGATGCTGACGTTGACGGTTCGCACATCCGTACGCTGCTGCTGACCTTCTTCTTCCGTCAGCTGCCGGAGCTGGTCGAGCGTGGCTACATCTATATTGCCCAGCCGCCGCTGTACAAGGTGAAACGAGGCAAGCa                                       <  2:1107478/150‑1 (MQ=255)

                                                                                                            catcatcatcatGACCGATGCTGACGTTGACGGTTCGCACATCCGTACGCTGCTGCTGACCTTCTTCTTCCGTCAGCTGCCGGAGCTGGTCGAGCGTGGCTACATCTATATTGCCCAGCCGCCGCTGTACAAGGTGAAACGAGGCAAGCa                                       >  2:7826326/1‑150 (MQ=255)

                                                                                                               catcatcatGACCGATGCTGACGTTGACGGTTCGCACATCCGTACGCTGCTGCTGACCTTCTTCTTCCGTCAGCTGCCGGAGCTGGTCGAGCGTGGCTACATCTATATTGCCCAGCCGCCGCTGTACAAGGTGAAACGAGGCAAGCAGGa                                    <  1:4304034/150‑1 (MQ=255)

                                                                                                                 tcatcatGACCGATGCTGACGTTGACGGTTCGCACATCCGTACGCTGCTGCTGACCTTCTTCTTCCGTCAGCTGCCGGAGCTGGTCGAGCGTGGCTACATCTATATTGCCCAGCCGCCGCTGTACAAGGTGAAACGAGGCAAGCAGGAGc                                  >  2:3948737/1‑150 (MQ=255)

                                                                                                                 tcatcatGACCGATGCTGACGTTGACGGTTCGCACATCCGTACGCTGCTGCTGACCTTCTTCTTCCGTCAGCTGCCGGAGCTGGTCGAGCGTGGCTACATCTATATTGCCCAGCCGCCGCTGTACAAGGTGAAACGAGGCAAGCAGGAGc                                  >  2:6253586/1‑150 (MQ=255)

                                                                                                                  catcatGACCGATGCTGACGTTGACGGTTCGCACATCCGTACGCTGCTGCTGACCTTCTTCTTCCGTCAGCTGCCGGAGCTGGTCGAGCGTGGCTACATCTATATTGCCCAGCCGCCGCTGTACAAGGTGAAACGAGGCAAGCAGGAGCa                                 <  1:6200877/150‑1 (MQ=255)

                                                                                                                    tcatGACCGATGCTGACGTTGACGGTTCGCACATCCGTACGCTGCTGCTGACCTTCTTCTTCCGTCAGCTGCCGGAGCTGGTCGAGCGTGGCTACATCTATATTGCCCAGCCGCCGCTGTACAAGGTGAAACGAGGCAAGCAGGAGCAGt                               <  1:7826326/150‑1 (MQ=255)

                                                                                                                      atGGCCGATGCTGACGTTGACGGTTCGCACATCCGTACGCTGCTGCTGACCTTCTTCTTCCGTCAGCTGCCGGAGCTGGTCGAGCGTGGCTACATCTATATTGCCCAGCCGCCGCTGTACAAGGTGAAACGAGGCAAGCAGGAGCAGTAc                             >  1:4446642/1‑150 (MQ=255)

                                                                                                                      atGACCGATGCTGACGTTGACGGTTCGCACATCCGTACGCTGCTGCTGACCTTCTTCTTCCGTCAGCTGCCGGAGCTGGTCGAGCGTGGCTACATCTATATTGccc                                                                         <  2:1571445/106‑1 (MQ=255)

                                                                                                                      atGACCGATGCTGACGTTGACGGTTCGCACATCCGTACGCTGCTGCTGACCTTCTTCTTCCGTCAGCTGCCGGAGCTGGTCGAGCGTGGCTACATCTATATTGCCCAGCCGCCGCTGTACAAGGTGAAACGAGGCAAGCAGGAGCAGTAc                             >  1:4551221/1‑150 (MQ=255)

                                                                                                                         aCCGATGCTGACGTTGACGGTTCGCACATCCGTACGCTGCTGCTGACCTTCTTCTTCCGTCAGCTGCCGGAGCTGGTCGAGCGTGGCTACATCTATATTGCCCAGCCGCCGCTGTACAAGGTGAAACGAGGCAAGCAGGAGCAGTACATc                          >  2:5542971/1‑150 (MQ=255)

                                                                                                                               gCTGACGTTGACGGTTCGCACATCCGTACGCTGCTGCTGACCTTCTTCTTCCGTCAGCTGCCGGAGCTGGTCGAGCGTGGCTACATCTATATTGCCCAGCCGCCGCTGTACAAGGTGAAACGAGGCAAGCAGGAGCAGTACATCAAGGAc                    <  2:8171438/150‑1 (MQ=255)

                                                                                                                                 tGACGTTGACGGTTCGCACATCCGTACGCTGCTGCTGACCTTCTTCTTCCGTCAGCTGCCGGAGCTGGTCGAGCGTGGCTACATCTATATTGCCCAGCCGCCGCTGTACAAGGTGAAACGAGGCAAGCAGGAGCAGTACATCAAGGACGa                  >  1:8006767/1‑150 (MQ=255)

                                                                                                                                 tGACGTTGACGGTTCGCACATCCGTACGCTGCTGCTGACCTTCTTCTTCCGTCAGCTGCCGGAGCTGGTCGAGCGTGGCTACATCTATATTGCCCAGCCGCCGCTGTACAAGGTGAAACGAGGCAAGCAGGAGCAGTACATCAAGGACGa                  >  1:5353113/1‑150 (MQ=255)

                                                                                                                                 tGACGTTGACGGTTCGCACATCCGTACGCTGCTGCTGACCTTCTTCTTCCGTCAGCTGCCGGAGCTGGTCGAGCGTGGCTACATCTATATTGCCCAGCCGCCGCTGTACAAGGTGAAACGAGGCAAGCAGGAGCAGTACATCAAGGACGa                  >  1:2809312/1‑150 (MQ=255)

                                                                                                                                  gACGTTGACGGTTCGCACATCCGTACGCTGCTGCTGACCTTCTTCTTCCGTCAGCTGCCGGAGCTGGTCGAGCGTGGCTACATCTATATTGCCCAGCCGCCGCTGTACAAGGTGAAACGAGGCAAGCAGGAGCAGTACATCAAGGACGAt                 <  2:6455077/150‑1 (MQ=255)

                                                                                                                                  gACGTTGACGGTTCGCACATCCGTACGCTGCTGCTGACCTTCTTCTTCCGTCAGCTGCCGGAGCTGGTCGAGCGTGGCTACATCTATATTGCCCAGCCGCCGCTGTACAAGGTGAAACGAGGCAAGCAGGAGCAGTACATCAAGGACGAt                 >  2:5351613/1‑150 (MQ=255)

                                                                                                                                    cGTTGACGGTTCGCACATCCGTACGCTGCTGCTGACCTTCTTCTTCCGTCAGCTGCCGGAGCTGGTCGAGCGTGGCTACATCTATATTGCCCAGCCGCCGCTGTACAAGGTGAAACGAGGCAAGCAGGAGCAGTACATCAAGGACGATGa               <  2:4348980/150‑1 (MQ=255)

                                                                                                                                          cGGTTCGCACATCCGTACGCTGCTGCTGACCTTCTTCTTCCGTCAGCTGCCGGAGCTGGTCGAGCGTGGCTACATCTATATTGCCCAGCCGCCGCTGTACAAGGTGAAACGAGGCAAGCAGGAGCAGTACATCAAGGACGATGAGGCCAt         <  2:1306398/150‑1 (MQ=255)

                                                                                                                                           gggtCGCACATCCGTACGCTGCTGCTGACCTTCTTCTTCCGTCAGCTGCCGGAGCTGGTCGAGCGTGGCTACATCTATATTGCCCAGCCGCCGCTGTACAAGGTGAAACGAGGCAAGCAGGAGCAGTACATCAAGGACGATGAGGCCATg        <  2:2628063/147‑1 (MQ=255)

                                                                                                                                           ggTTCGCACATCCGTACGCTGCTGCTGACCTTCTTCTTCCGTCAGCTGCCGGAGCTGGTCGAGCGTGGCTACATCTATATTGCCCAGCCGCCGCTGTACAAGGTGAAACGAGGCAAGCAGGAGCAGTACATCAAGGACGATGAGGCCATg        <  2:7871589/150‑1 (MQ=255)

                                                                                                                                           ggTTCGCACATCCGTACGCTGCTGCTGACCTTCTTCTTCCGTCAGCTGCCGGAGCTGGTCGAGCGTGGCTACATCTATATTGCCCAGCCGCCGCTGTACAAGGTGAAACGAGGCAAGCAGGAGCAGTACATCAAGGACGATGAGGCCATg        >  1:424849/1‑150 (MQ=255)

                                                                                                                                              tCGCACATCCGTACGCTGCTGCTGACCTTCTTCTTCCGTCAGCTGCCGGAGCTGGTCGAGCGTGGCTACATCTATATTGCCCAGCCGCCGCTGTACAAGGTGAAACGAGGCAAGCAGGAGCAGTACATCAAGGACGATGAGGCCATGGaa     <  1:5892875/150‑1 (MQ=255)

                                                                                                                                              tCGCACATCCGTACGCTGCTGCTGACCTTCTTCTTCCGTCAGCTGCCGGAGCTGGTCGAGCGTGGCTACATCTATATTGCCCAGCCGCCGCTGTACAAGGTGAAACGAGGCAAGCAGGAGCAGTACATCAAGGACGATGAGGCCATGGaa     <  1:549596/150‑1 (MQ=255)

                                                                                                                                              tCGCACATCCGTACGCTGCTGCTGACCTTCTTCTTCCGTCAGCTGCCGGAGCTGGTCGAGCGTGGCTACATCTATATTGCCCAGCCGCCGCTGTACAAGGTGAAACGAGGCAAGCAGGAGCAGTACATCAAGGACGATGAGGCCATGGaa     <  1:8815819/150‑1 (MQ=255)

                                                                                                                                                gCACATCCGTACGCTGCTGCTGACCTTCTTCTTCCGTCAGCTGCCGGAGCTGGTCGAGCGTGGCTACATCTATATTGCCCAGCCGCCGCTGTACAAGGTGAAACGAGGCAAGCAGGAGCAGTACATCAAGGACGATGAGGCCATGGAaga   >  2:5372593/1‑150 (MQ=255)

                                                                                                                                                 cacaTCCGTACGCTGCTGCTGACCTTCTTCTTCCGTCAGCTGCCGGAGCTGGTCGAGCGTGGCTACATCTATATTGCCCAGCCGCCGCTGTACAAGGTGAAACGAGGCAAGCAGGAGCAGTACATCAAGGACGATGAGGCCATGGAagag  >  2:5857244/1‑150 (MQ=255)

GCGCTTCGACAAGATGATTTCGTCCCAGGAAGTGGGCACGCTGATCACTGCGCTGGGCTGTGGCATCGGCCGCGAAGAGTACAACATCGACAAACTGCGTTATCACAACATCATCATCATGACCGATGCTGACGTTGACGGTTCGCACATCCGTACGCTGCTGCTGACCTTCTTCTTCCGTCAGCTGCCGGAGCTGGTCGAGCGTGGCTACATCTATATTGCCCAGCCGCCGCTGTACAAGGTGAAACGAGGCAAGCAGGAGCAGTACATCAAGGACGATGAGGCCATGGAAGAG  >  NC_002947/14703‑14997

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