BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,503,190	C→G	9.7%	A31G (GCC→GGC)	oprC →	copper receptor OprC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,503,190	0	C	G	9.7%	382.9 / 27.3	124	A31G (GCC→GGC)	oprC	copper receptor OprC
Reads supporting (aligned to +/- strand): ref base C (62/50); new base G (5/7); total (67/57)
Fisher's exact test for biased strand distribution p-value = 3.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.48e-01

Only 100 of 128 total aligned reads displayed.

TTTGCGCACGAAGCCAGGAGGCTTCGCGCTAACTCATGACTGATCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTC‑CCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGC‑CGTGGCCCCGAGTTCGCCGCTGACTGTGGTCACCAACCCGAAAGACCCGCGCCAGCCAGTGCCGGCCAGCGACGGCGCCGACTAC  >  NC_002947/5503041‑5503338

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

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

       aCGAAGCCAGGAGGCTTCGCGCTAACTCATGACTGATCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTC‑CCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCt                                                                                                                                                >  1:1887408/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)

                      ttCGCGCTAACTCATGACTGATCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTTTC‑CCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCAt                                                                                                                                 >  2:49274/1‑150 (MQ=255)

                      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‑CCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCa                                                                                                                                           >  1:582805/1‑137 (MQ=255)

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

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

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

                          cgcTAACTCATGACTGATCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTC‑CCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGc                                                                                                                             <  1:166656/150‑1 (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                                                                                                                                               >  2:2384682/1‑109 (MQ=255)

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

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

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

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

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

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

                                                       aTTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTC‑CCTGCGTGGTACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTg                                                                                                <  1:135478/150‑1 (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‑CCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCAccg                                                                                         >  1:705716/1‑150 (MQ=255)

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

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

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

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

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

                                                                      aCCCCTGTTTTTGCCTTGTC‑CCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGC‑CGTGGcc                                                                                >  2:1705649/1‑150 (MQ=255)

                                                                      aCCCCTGTTTTTGCCTTGTC‑CCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCAGGGCCGGCGGCCAGGGCCATg                                                                                                                                         <  1:2027859/94‑1 (MQ=38)

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

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

                                                                          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                                                                  >  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)

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

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

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

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

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

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

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

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

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

                                                                                                          cgccgcCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGC‑CGTGGCCCCGAGTTCGCCGCTGACTGTGGTCACCAACCCGaa                                             <  1:232487/150‑1 (MQ=255)

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

                                                                                                           gccgccCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGC‑CGTGGCCCCGAGTTCGCCGCTGACTGTGGTCACCAACCCGaaa                                            <  1:334324/150‑1 (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                                                                <  2:1455013/113‑1 (MQ=255)

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

                                                                                                                            gctgctCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGC‑CGTGGCCCCGAGTTCGCCGCTGACTGTGGTCACCAACCCGAAAGACCCGCGCCAGCCAGt                           >  2:1327958/1‑150 (MQ=255)

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

                                                                                                                            gctgctAGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGC‑CGTGGCCCCGAGTTCGCCGCTGACTGTGGTCACCAACCCGAAAGACCCGCGCCAGCCAGt                           <  2:783940/150‑1 (MQ=255)

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

                                                                                                                                ctCGCGCCCATTGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGC‑CGTGGCCCCGAGTTCGCCGCTGACTGTGGTCACCAACCCGAAAGACCCGCGCCAGCCAGTGCCg                       >  2:2208965/1‑150 (MQ=255)

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

                                                                                                                                 tCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGC‑CGTGGCCCCGAGTTCGCCGCTGACTGTGGTCACCAACCCGAAAGACCCGCGCCAGCCAGTGCCgg                      <  1:1801555/150‑1 (MQ=255)

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

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

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

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

                                                                                                                                      ttcATGGCCAGGGCCGGCGGCCAGGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGC‑CGTGGCCCCGAGTTCGCCGCTGACTGTGGTCACCAACCCGAAAGACCCGCGCCAGCCAGTGCCGGCCAGc                 >  2:1194798/3‑150 (MQ=255)

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

                                                                                                                                        cATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGC‑CGTGGCCCCGAGTTCGCCGCTGACTGTGGTCACCAACCCGAAAGACCCGCGCCAGCCAGTGCCGGCCAGCGa               <  2:2113967/150‑1 (MQ=255)

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

                                                                                                                                              cagggCCGGCGGCCAGGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGC‑CGTGGCCCCGAGTTCGCCGCTGACTGTGGTCACCAACCCGAAAGACCCGCGCCAGCCAGTGCCGGCCAGCGACGgcgc         >  2:831791/4‑150 (MQ=255)

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

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

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

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

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

                                                                                                                                                tGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGC‑CGTGGCCCCGAGTTCGCCGCTGACTGTGGTCACCAACCCGAAAGACCCGCGCCAGCCAGTGCCGGCCAGCGACGGCGCCg       >  2:1988671/1‑150 (MQ=255)

                                                                                                                                                 ggCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGC‑CGTGGCCCCGAGTTCGCCGCTGACTGTGGTCACCAACCCGAAAGACCCGCGCCAGCCAGTGCCGGCCAGCGACGGCGCCGa      <  2:1658848/150‑1 (MQ=255)

                                                                                                                                                    cgccgGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGC‑CGTGGCCCCGAGTTCGCCGCTGACTGTGGTCACCAACCCGAAAGACCCGCGCCAGCCAGTGCCGGCCAGCGACGGCGCCGACTa   >  1:1468246/1‑150 (MQ=255)

                                                                                                                                                     gccgGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGC‑CGTGGCCCCGAGTTCGCCGCTGACTGTGGTCACCAACCCTAAAGACCCGCGCCAGCCAGTGCCGGCCAGCGACGGCGCCGACTAc  <  2:185258/150‑1 (MQ=255)

                                                                                                                                                     gccgGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGC‑CGTGGCCCCGAGTTCGCCGCTGACTGTGGTCACCAACCCGAAAGACCCGCGCCAGCCAGTGCCGGCCAGCGACGGCGCCGACTAc  >  1:2156505/1‑150 (MQ=255)

TTTGCGCACGAAGCCAGGAGGCTTCGCGCTAACTCATGACTGATCGACTGGAATCATTATGTCCGGCTGCACCCCTGTTTTTGCCTTGTC‑CCTGCGTGGGACCATCGCCGCCCTGTGCGGCTCGCTGCTCGCGCCCATGGCCCTGGCCGCCGGCCCTGGCCATGAAGACCATGTGCACGACGCACCCGAGTTGAGCCCGACGGTGATCACCGC‑CGTGGCCCCGAGTTCGCCGCTGACTGTGGTCACCAACCCGAAAGACCCGCGCCAGCCAGTGCCGGCCAGCGACGGCGCCGACTAC  >  NC_002947/5503041‑5503338

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