BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,503,184	T→G	8.7%	L29R (CTG→CGG) ‡	oprC →	copper receptor OprC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,503,184	0	T	G	8.7%	354.8 / 23.5	126	L29R (CTG→CGG) ‡	oprC	copper receptor OprC
Reads supporting (aligned to +/- strand): ref base T (63/52); new base G (4/7); total (67/59)
Fisher's exact test for biased strand distribution p-value = 3.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.93e-01

Only 100 of 129 total aligned reads displayed.

TTTTGCGCACGAAGCCAGGAGGCTTCGCGCTAACTCATGACTGATCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTC‑CCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGC‑CGTGGCCCCGAGTTCGCCGCTGACTGTGGTCACCAACCCGAAAGACCCGCGCCAGCCAGTGCCGGCCAGCGACGGCGCCG  >  NC_002947/5503040‑5503333

ttttGCGCACGAAGCCAGGAGGCTTCGCGCTAACTCATGACTGATCGACTGGAATCATTATGTCCGGCTGCACCCCTGTGTTTGCCTTGTC‑CCTGCGTGGGACCATCGCGGCCCGGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGgccg                                                                                                                                                   >  1:1331400/1‑150 (MQ=255)

 tttGCGCACGAAGCCAGGAGGCTTCGCGCTAACTCATGACTGATCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTC‑CCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGgccgc                                                                                                                                                  <  1:1438566/150‑1 (MQ=255)

      gcACGAAGCCAGGAGGCTTCGCGCTAACTCATGACTGATCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTC‑CCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGcc                                                                                                                                             >  2:2305708/1‑150 (MQ=255)

        aCGAAGCCAGGAGGCTTCGCGCTAACTCATGACTGATCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTC‑CCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCt                                                                                                                                           >  2:1723374/1‑150 (MQ=255)

        aCGAAGCCAGGAGGCTTCGCGCTAACTCATGACTGATCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTC‑CCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCt                                                                                                                                           >  1:1887408/1‑150 (MQ=255)

          ggaGCCAGGAGGCTTCGCGCTAACTCATGACTGATCGACTGGAATCATTATGTCCGGCTGGGGGGGGGTGGGGGGGGTGTC‑CCTGCGTGGGACCATCGCCGCCCTGTGCTGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTgg                                                                                                                                         <  1:359988/148‑1 (MQ=255)

          gAAGCCAGGAGGCTTCGCGCTAACTCATGACTGATCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTC‑CCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCAGGGCCGGCGGccaggg                                                                                                                                         <  2:2206923/150‑5 (MQ=255)

             gCCAGGAGGCTTCGCGCTAACTCATGACTGATCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTC‑CCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCa                                                                                                                                      >  2:1580964/1‑150 (MQ=255)

              ccAGGAGGCTTCGCGCTAACTCATGACTGATCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTC‑CCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCAt                                                                                                                                     <  1:1944111/150‑1 (MQ=255)

                 ggaggCTTCGCGCTAACTCATGACTGATCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTC‑CCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCAGGGCCGGCGGCCAGGGCCATggg                                                                                                                                  <  1:35828/150‑3 (MQ=255)

                  gaggCTTCGCGCTAACTCATGACTGATCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTC‑CCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAg                                                                                                                                 >  2:1578153/1‑150 (MQ=255)

                  gaggCTTCGCGCTAACTCATGACTGATCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTC‑CCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCAGGGCCGGCGGCCAGGGCCATgggc                                                                                                                                 <  2:1386248/150‑4 (MQ=255)

                   aggCTTCGCGCTAACTCATGACTGATCGACTGGAATCATTATGTCCGGCTGCACCACTGTTTTTGCCTTGTC‑CCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGa                                                                                                                                <  1:1000784/150‑1 (MQ=255)

                      cTTCGCGCTAACTCATGGCTGATCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTC‑CCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCAGGGCCGGCGGCCAGGGCCATgggcgcga                                                                                                                             <  2:1388322/150‑8 (MQ=38)

                       ttCGCGCTAACTCATGACTGATCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTC‑CCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGAACAt                                                                                                                            <  2:409208/150‑1 (MQ=255)

                         cgcgCTAACTCATGACTGATCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTC‑CCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCAtgt                                                                                                                          >  2:166656/1‑150 (MQ=255)

                         cgcgCTAACTCATGACTGATCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTC‑CCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCACCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCAtgt                                                                                                                          <  2:1198101/150‑1 (MQ=255)

                          gcgcTAACTCATGACTGATCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTC‑CCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCa                                                                                                                                      <  2:582805/137‑1 (MQ=255)

                          gcgcTAACTCATGACTGATCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTC‑CCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCAtgtg                                                                                                                         >  2:1410157/1‑150 (MQ=255)

                          gcgcTAACTCATGACTGATCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTC‑CCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCAtgtg                                                                                                                         >  2:1849033/1‑150 (MQ=255)

                          gcgcTAACTCATGACTGATCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTC‑CCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCAtgtg                                                                                                                         >  2:232487/1‑150 (MQ=255)

                           cgcTAACTCATGACTGATCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTC‑CCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGc                                                                                                                        <  1:166656/150‑1 (MQ=255)

                                aCTCATGACTGATCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTC‑CCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTg                                                                                                                                                       <  1:2163396/114‑1 (MQ=255)

                                aCTCATGACTGATCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTC‑CCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTg                                                                                                                                                       >  2:2163396/1‑114 (MQ=255)

                                  tCATGACTGATCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTC‑CCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGc                                                                                                                 >  1:185258/1‑150 (MQ=255)

                                      gACTGATCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGAC‑CCTGCGTGGGGCCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCGCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCAccc                                                                                                             >  1:1267048/1‑150 (MQ=255)

                                          gATCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTC‑CCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGt                                                                                                         <  1:1849033/150‑1 (MQ=255)

                                           aTCGACTGGAATCATTATGTCCGGCTTCACCTCTGTTTTTGCATTGCCAGCGGCAACGG‑CCATCGTCGCCCTGTGCGGCTCGCTGCTCGCGCCCCTGGCCCTGGTCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGtt                                                                                                        <  2:675599/150‑1 (MQ=255)

                                              gACTGGAGTCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTC‑CCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAg                                                                                                     >  2:706839/1‑150 (MQ=255)

                                              gACTGGAATCATTATGACCGGCTGCACCCCTGTTTTTGCCTTGTC‑CCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAg                                                                                                     <  2:848123/150‑1 (MQ=255)

                                                  ggAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTC‑CCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTg                                                                                                                                          <  1:2384682/109‑1 (MQ=255)

                                                  ggAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTC‑CCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTg                                                                                                                                          >  2:2384682/1‑109 (MQ=255)

                                                   gAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTC‑CCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTgg                                                                                                                                                      >  1:354476/1‑96 (MQ=255)

                                                   gAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTC‑CCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGa                                                                                                >  1:2113967/1‑150 (MQ=255)

                                                   gAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTC‑CCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGa                                                                                                >  2:1957140/1‑150 (MQ=255)

                                                     aTCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTC‑CCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACg                                                                                              >  2:334324/1‑150 (MQ=255)

                                                     aTCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTC‑CCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACg                                                                                              <  1:1957140/150‑1 (MQ=255)

                                                      tCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTC‑CCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGc                                                                                                                                        >  1:976871/1‑107 (MQ=255)

                                                      tCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTC‑CCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACgg                                                                                             >  2:256628/1‑150 (MQ=255)

                                                      tCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTC‑CCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACgg                                                                                             >  2:1313470/1‑150 (MQ=255)

                                                        aTTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTC‑CCTGCGTGGTACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTg                                                                                           <  1:135478/150‑1 (MQ=255)

                                                        aTTATGTACGGCTGCACCCCTGTTTTTGCCTTGTC‑CCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTg                                                                                           >  1:1656829/1‑150 (MQ=255)

                                                         ttATGTCCGGCTGCACCCCTGTTTTTGCCTTGTC‑CCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGa                                                                                          >  2:779799/1‑150 (MQ=255)

                                                              tCCGGCTGCACCCATGTTTTTGCCTTGTC‑CCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCAcc                                                                                     <  2:1331400/150‑1 (MQ=255)

                                                               ccGGCTGCACCCCTGTTTTTGCCTTGTC‑CCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTg                                                                                                                                          >  1:1799761/1‑96 (MQ=255)

                                                               ccGGCTGCACCCCTGTTTTTGCCTTGTC‑CCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTg                                                                                                                                          <  2:1799761/96‑1 (MQ=255)

                                                               ccGGCTGCACCCCTGTTTTTGCCTTGTC‑CCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCAccg                                                                                    >  1:705716/1‑150 (MQ=255)

                                                               ccGGCTGCACCCCTGTTTTTGCCTTGTC‑CCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCAccg                                                                                    <  2:705716/150‑1 (MQ=255)

                                                                cGGCTGCACCCCTGTTTTTGCCTTGTC‑CCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCAGGGCCGGCGGCCAGGGCCATgg                                                                                                                                   <  1:517453/102‑2 (MQ=37)

                                                                cGGCTGCACCCCTGTTTTTGCCTTGTC‑CCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCAGGGCCGGCGGCCAGGGCCATgg                                                                                                                                   >  2:517453/1‑101 (MQ=37)

                                                                  gCTGCACCCCTGTTTTTGCCTTGTC‑CCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGccc                                                                                                  >  2:1886251/1‑133 (MQ=255)

                                                                     gCACCCCTGTTTTTGCCTTGTC‑CCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGC‑CGTgg                                                                             <  1:1807022/150‑1 (MQ=255)

                                                                     gCACCCCTGTTTTTGCCTTGTC‑CCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGC‑CGTgg                                                                             <  2:2157813/150‑1 (MQ=255)

                                                                      cACCCCTGTTTTTGCCTTGTC‑CCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGACCCCGAGTTGAGCCCGACGGTGATCACCGC‑CGTGGc                                                                            >  1:1820130/1‑150 (MQ=255)

                                                                       aCCCCTGTTTTTGCCTTGTC‑CCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCAGGGCCGGCGGCCAGGGCCATg                                                                                                                                    >  2:2027859/1‑94 (MQ=37)

                                                                        ccccTGTTTTTGCCTTGTC‑CCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCAtgtg                                                                                                                         >  2:1347144/1‑104 (MQ=255)

                                                                         cccTGTTTTTGCCTTGTC‑CCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGaa                                                                                                                                  <  2:170460/94‑1 (MQ=255)

                                                                         cccTGTTTTTGCCTTGTC‑CCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCAGGGCCGGCGGCCAGGGCCAt                                                                                                                                     >  1:2327232/1‑91 (MQ=37)

                                                                           cTGTTTTTGCCTTGTC‑CCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGC‑CGTGGCCCCGa                                                                       >  2:2090581/1‑150 (MQ=255)

                                                                           cTGTTTTTGCCTTGTC‑CCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGC‑CGTGGCCCCGa                                                                       >  1:1091630/1‑150 (MQ=255)

                                                                             gTTTTTGCCTTGTC‑CCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGC‑CGTGGCCCCGAGt                                                                     <  1:614301/150‑1 (MQ=255)

                                                                                     cTTGTC‑CCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGC‑CGTGGCCCCGAGTTCGCCGCt                                                             >  2:125767/1‑150 (MQ=255)

                                                                                     cTTGTC‑CCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGC‑CGTGGCCCCGAGTTCGCCGCt                                                             >  1:2260841/1‑150 (MQ=255)

                                                                                      ttGTC‑CCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGC‑CGt                                                                               <  2:1741352/131‑1 (MQ=255)

                                                                                      ttGTC‑CCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGC‑CGt                                                                               >  1:1741352/1‑131 (MQ=255)

                                                                                       tGTC‑CCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGC‑CGTGGCCCCGAGTTCGCCGCTGa                                                           <  1:1580964/150‑1 (MQ=255)

                                                                                             cTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGC‑CGTGGCCCCGAGTTCGCCGCTGACtgtg                                                      >  1:1658848/1‑150 (MQ=255)

                                                                                             cTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGC‑CGTGGCCCCGAGTTCGCCGCTGACtgtg                                                      <  1:256628/150‑1 (MQ=255)

                                                                                             cTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGC‑CATGGCCCCGAGTTCGCCGCTGACtgtg                                                      <  1:2075849/150‑1 (MQ=255)

                                                                                             cTGCGTGGGACCATCGCCGCCCTGTCCGCCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGATCATCTGCACGACGCTCCCGAGTTGAGCCCGACGGTGATCACCGC‑CGTGGCCCCGAGTTCGCCGCTGACtgtg                                                      <  2:1331736/150‑1 (MQ=255)

                                                                                              tGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGC‑CGTGGCCCCGAGTTCGCCGCTGACTGTgg                                                     >  1:2314650/1‑150 (MQ=255)

                                                                                               gCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGC‑CGTGGCCCCGAGTTCGCCGCTGACTGTGGt                                                    >  1:683826/1‑150 (MQ=255)

                                                                                                  tGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGC‑CGTGGCCCCGAGTTCGCCGCTGACtgtg                                                      >  1:117821/1‑145 (MQ=255)

                                                                                                  tGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGC‑CGTGGCCCCGAGTTCGCCGCTGACtgtg                                                      <  2:117821/145‑1 (MQ=255)

                                                                                                  tGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGC‑CGTGGCCCCGAGTTCGCCGCTGACTGTGGTCAc                                                 <  2:2314650/150‑1 (MQ=255)

                                                                                                  tGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGC‑CGTGGCCCCGAGTTCGCCGCTGACTGTGGTCAc                                                 <  1:2198083/150‑1 (MQ=255)

                                                                                                      aCCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGC‑CGTGGCCCCGAGTTCGCCGCTGACtgtg                                                      <  2:1644526/141‑1 (MQ=255)

                                                                                                       ccATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGC‑CGTGGCCCCGAGTTCGCCGCTGACTGTGGTCACCAAcc                                            >  2:2076654/1‑150 (MQ=255)

                                                                                                        cATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGC‑CGTGGCCCCGAGTTCGCCGCTGACTGTGGTCACCAAccc                                           >  1:1688500/1‑150 (MQ=255)

                                                                                                            gccgccCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGC‑CGTGGCCCCGAGTTCGCCGCTGACTGTGGTCACCAACCCGaaa                                       >  1:2151736/1‑150 (MQ=255)

                                                                                                             ccgccCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGC‑CGTGGCCCCGAGTTCGCCGCTGACTGTGGTCACCAACCCGAAAg                                      >  2:1534791/1‑150 (MQ=255)

                                                                                                              cgccCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGC‑CGTGGCCCCGAGTTCGCCGCTGACTGTGGTCACCAACCCGAAAGa                                     >  2:901452/1‑150 (MQ=255)

                                                                                                               gccCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGC‑CGTGGCCCCGAGTTCGCCGCTGACTGTGGTCACCAACCCGAAAGAc                                    <  2:683826/150‑1 (MQ=255)

                                                                                                                 ccTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGC‑CGTGGCCCCGAGTTCGCCGCTGACTGTGGTCACCAACCCGAAAGAccc                                  <  2:404752/150‑1 (MQ=255)

                                                                                                                         gCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCACTGGCCATGAAGACAATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGC‑CGTGGCCCCGAGTTCGCCGCTGACTGTGGTCACCAACCCGAAAGACCCGCgccagc                          >  2:41048/1‑150 (MQ=255)

                                                                                                                             gctgctCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGC‑CGTGGCCCCGAGTTCGCCGCTGa                                                           >  1:1455013/1‑113 (MQ=255)

                                                                                                                             gctgctCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGC‑CGTGGCCCCGAGTTCGCCGCTGa                                                           <  2:1455013/113‑1 (MQ=255)

                                                                                                                             gctgctCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGC‑CGTGGCCCCGAGTTCGCCGCTGACTGTGGTCACCAACCCGAAAGACCCGCGCCAGCCAGt                      <  1:879775/150‑1 (MQ=255)

                                                                                                                               tgctCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGC‑CGTGGCCCCGAGTTCGCCGCTGACTGTGGTCACCAACCCGAAAGACCCGCGCCAGCCAGTGc                    >  1:921794/1‑150 (MQ=255)

                                                                                                                                  tCGCGCCCATGGCCCTGGCCGCCGGCCCTGTCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGC‑CGTGGCCCCGAGTTCGCCGCTGACTGTGGTCACCAACCCGAAAGACCCGCGCCAGCCAGTGCCgg                 <  1:1578153/150‑1 (MQ=255)

                                                                                                                                  tCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGC‑CGTGGCCCCGAGTTCGCCGCTGACTGTGGTCACCAACCCGAAAGACCCGCGCCAGCCAGTGCCgg                 >  2:1361192/1‑150 (MQ=255)

                                                                                                                                    gcgcCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGC‑CGTGGCCCCGAGTTCGCCGCTGACTGTGGTCACCAACCCGAAAGAc                                    <  2:80081/129‑1 (MQ=255)

                                                                                                                                    gcgcCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGC‑CGTGGCCCCGAGTTCGCCGCTGACTGTGGTCACCAACCCGAAAGACCCGCGCCAGCCAGTGCCGGcc               >  1:1054689/1‑150 (MQ=255)

                                                                                                                                     cgcCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCCCTTTCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGC‑CGTGGCCCCGAGTTCGCCGCTGACTGTGGTCACCAACCCGAAAGACCCGCGCCAGCCAGTGCCGGCCa              <  1:1000400/150‑1 (MQ=255)

                                                                                                                                        ccATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGC‑CGTGGCCCCGAGTTCGCCGCTGACTGTGGTCACCAACCCGAAAGACCCGCGCCGGCCAGTGCCGGCCAGCg           >  2:2363400/1‑150 (MQ=255)

                                                                                                                                              cccTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGC‑CGTGGCCCCGAGTTCGCCGCTGACTGTGGTCACCAACCCGAAAGACCCGCGCCAGCCAGCGCCGGCCAGcccccgcc     >  2:1942012/1‑143 (MQ=255)

                                                                                                                                               ccTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGC‑CGTGGCCCCGAGTTCGCCGCTGACTGTGGTCACCAACCCGAAAGACCCGCGCAAGCCAGTGCCGGCCAGCGACGgcgc    >  1:1502630/1‑150 (MQ=255)

                                                                                                                                                cTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCAc                                                                                      <  1:1637948/68‑1 (MQ=255)

                                                                                                                                                cTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGAT‑ACGGCTTGTGGCCCCGAGTTCGCCGCTGACTGTGGTCACCAACCCGAAAGACCCGCGCCAGCCAGTGCCGGCCAGCGACGGCGcc   <  1:1731258/150‑1 (MQ=255)

                                                                                                                                                 tGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGC‑CGTGGCCCCGAGTTCGCCGCTGACTGTGGTCAc                                                 >  2:2191718/1‑103 (MQ=255)

TTTTGCGCACGAAGCCAGGAGGCTTCGCGCTAACTCATGACTGATCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTC‑CCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGC‑CGTGGCCCCGAGTTCGCCGCTGACTGTGGTCACCAACCCGAAAGACCCGCGCCAGCCAGTGCCGGCCAGCGACGGCGCCG  >  NC_002947/5503040‑5503333

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