BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,550,520	C→A	5.3%	A293D (GCT→GAT)	groL →	chaperonin GroEL

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,550,520	0	C	A	5.3%	352.3 / 2.5	109	A293D (GCT→GAT)	groL	chaperonin GroEL
Reads supporting (aligned to +/- strand): ref base C (49/53); new base A (4/2); total (53/56)
Fisher's exact test for biased strand distribution p-value = 4.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 109 total aligned reads displayed.

CACTGCTGATCGTTGCCGAAGACGTTGAAGGTGAAGCGCTGGCTACCCTGGTAGTCAACAACATGCGCGGCATCGTCAAGGTTGCTGCAGTCAAGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCTGGAGCACCTGGGTAACGCCAAGCGCGTCATCCTGTCCAAGGAAAACACCACCATCATCGACGGTGCTGGCGCTGACACCGAGAT  >  NC_002947/1550379‑1550667

cACTGCTGATCGTTGCCGAAGACGTTGAAGGTGAAGCGCTGGCTACCCTGGTAGTCAACAACATGCGCGGCATCGTCAAGGTTGCTGCAGTCAAGGCTCCGGGCGGCTGCGAGCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGa                                                                                                                                             >  2:863399/1‑150 (MQ=255)

 aCTGCTGATCGTTGCCGAAGACGTTGAAGGTGAAGCGCTGGCTACCCTGGTAGTCAACAACATGCGCGGCATCGTCAAGGTTGCTGCAGTCAAGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGAc                                                                                                                                            >  2:1849566/1‑150 (MQ=255)

       gATCGTTGCCGAAGACGTTGAAGGTGAAGCGCTGGCTACCCTGGTAGTCAACAACATGCGCGGCATCGTCAAGGTTGCTGCAGTCAAGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTggcgg                                                                                                                                      <  2:731558/150‑1 (MQ=255)

              gCCGAAGACGTTGAAGGTGAAGCGCTGGCTACCCTGGTAGTCAACAACATGCGCGGCATCGTCAAGGTTGCTGCAGTCAAGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTc                                                                                                                               <  2:1831774/150‑1 (MQ=255)

                 gAAGACGTTGAAGGTGAAGCGCTGGCTACCCTGGTAGTCAACAACATGCGCGGCATCGTCAAGGTTGCTGCAGTCAAGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCAtc                                                                                                                            >  1:732948/1‑150 (MQ=255)

                  aaGACGTTGAAGGTGAAGCGCTGGCTACCCTGGTAGTCAACAACATGCGCGGCATCGTCAAGGTTGCTGCAGTCAAGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCAtct                                                                                                                           <  1:320249/150‑1 (MQ=255)

                       gTTGAAGGTGAAGCGCTGGCTACCCTGGTAGTCAACAACATGCGCGGCATCGTCAAGGTTGCTGCAGTCAAGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCgaa                                                                                                                      <  2:1818137/150‑1 (MQ=255)

                       gTTGAAGGTGAAGCGCTGGCTACCCTGGTAGTCAACAACATGCGCGGCATCGTCAAGGTTGCTGCAGTCAAGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCgaa                                                                                                                      >  1:997896/1‑150 (MQ=255)

                         tGAAGGTGAAGCGCTGGCTACCCTGGTAGTCAACAACATGCGCGGCATCGTCAAGGTTGCTGCAGTCAAGGCTCCGGGCTTCGGCGATCCCCGCAAGGCCACGGTGCAGGACACCGCTGTCCTGACTGGCGGCCAGGTCATCTCCgaaga                                                                                                                    >  1:260261/1‑150 (MQ=255)

                          gAAGGTGAAGCGCTGGCTACCCTGGTAGTCAACAACATGCGCGGCATCGTCAAGGTTGCTGCAGTCAAGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCgaagaa                                                                                                                   >  1:497839/1‑150 (MQ=255)

                            aGGTGAAGCGCTGGCTACCCTGGTAGTCAACAACATGCGCGGCATCGTCAAGGTTGCTGCAGTCAAGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAAt                                                                                                                 <  2:1814410/150‑1 (MQ=255)

                             ggTGAAGCGCTGGCTACCCTGGTAGTCAACAACATGCGCGGCATCGTCAAGGTTGCTGCAGTCAAGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATc                                                                                                                <  2:1314580/150‑1 (MQ=255)

                               tGAAGCGCTGGCTACCCTGGTAGTCAACAACATGCGCGGCATCGTCAAGGTTGCTGCAGTCAAGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCgg                                                                                                              <  2:968302/150‑1 (MQ=255)

                                gAAGCGCTGGCTACCCTGGTAGTCAACAACATGCGCGGCATCGTCAAGGTTGCTGCAGTCAAGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGc                                                                                                             <  2:2362395/150‑1 (MQ=255)

                                gAAGCGCTGGCTACCCTGGTAGTCAACAACATGCGCGGCATCGTCAAGGTTGCTGCAGTCAAGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGc                                                                                                             >  2:942852/1‑150 (MQ=255)

                                 aaGCGCTGGCTACCCTGGTAGTCAACAACATGCGCGGCATCGTCAAGGTTGCTGCAGTCAAGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACAGCGATGTCCTg                                                                                                                                              >  2:2298618/1‑116 (MQ=255)

                                 aaGCGCTGGCTACCCTGGTAGTCAACAACATGCGCGGCATCGTCAAGGTTGCTGCAGTCAAGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACAGCGATGTCCTg                                                                                                                                              <  1:2298618/116‑1 (MQ=255)

                                        ggCTACCCTGGTAGTCAACAACATGCGCGGCATCGTCAAGGTTGCTGCAGTCAAGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATc                                                                                                                >  1:386561/1‑139 (MQ=255)

                                        ggCTACCCTGGTAGTCAACAACATGCGCGGCATCGTCAAGGTTGCTGCAGTCAAGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATc                                                                                                                <  2:386561/139‑1 (MQ=255)

                                         gCTACCCTGGTAGTCAACAACATGCGCGGCATCGTCAAGGTTGCTGCAGTCAAGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTcc                                                                                                                                                <  1:1022444/106‑1 (MQ=255)

                                         gCTACCCTGGTAGTCAACAACATGCGCGGCATCGTCAAGGTTGCTGCAGTCAAGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTcc                                                                                                                                                >  2:1022444/1‑106 (MQ=255)

                                         gCTACCCTGGTAGTCAACAACATGCGCGGCATCGTCAAGGTTGCTGCAGTCAAGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTg                                                                                                    <  2:425117/150‑1 (MQ=255)

                                         gCTACCCTGGTAGTCAACAACATGCGCGGCATCGTCAAGGTTGCTGCAGTCAAGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTg                                                                                                    >  1:2057347/1‑150 (MQ=255)

                                          cTACCCTGGTAGTCAACAACATGCGCGGCATCGTCAAGGTTGCTGCAGTCAAGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTgg                                                                                                   <  1:32973/150‑1 (MQ=255)

                                           tACCCTGGTAGTCAACAACATGCGCGGCATCGTCAAGGTTGCTGCAGTCAAGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTc                                                                                                        <  2:2015756/144‑1 (MQ=255)

                                           tACCCTGGTAGTCAACAACATGCGCGGCATCGTCAAGGTTGCTGCAGTCAAGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTc                                                                                                        >  1:2015756/1‑144 (MQ=255)

                                             cccTGGTAGTCAACAACATGCGCGGCATCGTCAAGGTTGCTGCAGTCAAGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGaaa                                                                                                <  1:1925732/150‑1 (MQ=255)

                                              ccTGGTAGTCAACAACATGCGCGGCATCGTCAAGGTTGCTGCAGTCAAGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAg                                                                                                                                  >  1:2149669/1‑115 (MQ=255)

                                              ccTGGTAGTCAACAACATGCGCGGCATCGTCAAGGTTGCTGCAGTCAAGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAg                                                                                                                                  <  2:2149669/115‑1 (MQ=255)

                                              ccTGGTAGTCAACAACATGCGCGGCATCGTCAAGGTTGCTGCAGTCAAGGCTCCGGGCTTCGGCGATCGCCGCAAGGCATGGCTTCTGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAAc                                                                                               >  1:1332929/1‑150 (MQ=255)

                                                tGGTAGTCAACAACATGCGCGGCATCGTCAAGGTTGCTGCAGTCAAGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCg                                                                                             <  1:2174024/150‑1 (MQ=255)

                                                tGGTAGTCAACAACATGCGCGGCATCGTCAAGGTTGCTGCAGTCAAGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCg                                                                                             >  2:1739117/1‑150 (MQ=255)

                                                tGGTAGTCAACAACATGCGCGGCATCGTCAAGGTTGCTGCAGTCAAGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCg                                                                                             >  2:551720/1‑150 (MQ=255)

                                                     gTCAACAACATGCGCGGCATCGTCAAGGTTGCTGCAGTCAAGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCgaa                                                                                                                      <  1:2436483/120‑1 (MQ=255)

                                                     gTCAACAACATGCGCGGCATCGTCAAGGTTGCTGCAGTCAAGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCgaa                                                                                                                      >  2:2436483/1‑120 (MQ=255)

                                                     gTCAACAACATGCGCGGCATCGCCAAGGTTGCTGTTGTCGAGGCTCCGTGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTAcc                                                                                        >  1:1871456/1‑150 (MQ=255)

                                                       caacaaCATGCGCGGCATCGTCAAGGTTGCTGCAGTCAAGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCt                                                                                      >  2:1059538/1‑150 (MQ=255)

                                                        aacaacATGCGCGGCATCGTCAAGGTTGCTGCAGTCAAGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGt                                                                                                                                                  >  2:2211737/1‑89 (MQ=255)

                                                        aacaacATGCGCGGCATCGTCAAGGTTGCTGCAGTCAAGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGt                                                                                                                                                  <  1:2211737/89‑1 (MQ=255)

                                                        aacaacATGCGCGGCATCGTCAAGGTTGCTGCAGTCAAGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCTg                                                                                     >  2:1523245/1‑150 (MQ=255)

                                                        aacaacATGCGCGGCATCGTCAAGGTTGCTGCAGTCAAGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCTg                                                                                     >  2:1481184/1‑150 (MQ=255)

                                                        aacaacATGCGCGGCATCGTCAAGGTTGCTGCAGTCAAGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCTg                                                                                     >  1:902622/1‑150 (MQ=255)

                                                        aacaacATGCGAGGCATCGTCAAGGTTGCTGCAGTCAAGGCTCCGGGCTTCGGCGATCGCCGCACGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCTg                                                                                     <  1:491355/150‑1 (MQ=255)

                                                          caacaCGCGCGGCATCGTCAAGGTTGCTGCAGTCACGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCCGGACCTCGCTGTCCTGACTGGCGGCCAGGTCAACTCCGAGGGAGTCGTTCTGTCCCTGGAAGCCGCTACCCTGGa                                                                                   <  2:1207476/150‑1 (MQ=255)

                                                             caTGCGCGGCATCGTCAAGGTTGCTGCAGTCAAGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCa                                                                                                                                   <  1:2096135/99‑1 (MQ=255)

                                                              aTGCGCGGCATCGTCAAGGTTGCTGCAGTCAAGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCTGGAGCAc                                                                               >  2:2036611/1‑150 (MQ=255)

                                                              aTGCGCGGCATCGTCAAGGTTGCTGCAGTCAAGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCTGGAGCAc                                                                               <  2:1180/150‑1 (MQ=255)

                                                                 cgcgGCATCGTCAAGGTTGCTGCAGTCAAGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCTGGAGCACCTg                                                                            <  1:1626830/150‑1 (MQ=255)

                                                                  gcgGCATCGTCAAGGTTGCTGCAGTCAAGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCTGGAGCACCTgg                                                                           <  2:948723/150‑1 (MQ=255)

                                                                    ggCATCGTCAAGGTTGCTGCAGTCAAGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCTGGAGCACCTGGGt                                                                         <  2:1156321/150‑1 (MQ=255)

                                                                       aTCGTCAAGGTTGCTGCAGTCAAGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCTGGAGCACCTGGGTAAc                                                                      >  1:1918733/1‑150 (MQ=255)

                                                                          gTCAAGGTTGCTGCAGTCAAGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAAACGCTACCCTGGAGCACCTGGGTAACGcc                                                                   >  1:639133/1‑150 (MQ=255)

                                                                          gTCAAGGTTGCTGCAGTCAAGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGACGTTCTGACTGGCGGTCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCTGGGGCACCTGGGTAACGcc                                                                   >  1:766744/1‑150 (MQ=255)

                                                                           tCAAGGTTGCTGCAGTCAAGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCTGGAGCACCTGGGTAACGCCa                                                                  >  1:2402367/1‑150 (MQ=255)

                                                                           tCAAGGTTGCTGCAGTCAAGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCTGGAGCACCTGGGTAACGCCa                                                                  >  2:1569570/1‑150 (MQ=255)

                                                                           tCAAGGTTGCTGCAGTCAAGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCTGGAGCACCTGGGTAACGCCa                                                                  >  2:203005/1‑150 (MQ=255)

                                                                             aaGGTTGCTGCAGTCAAGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCTGGAGCACCTGGGTAACGCCAAg                                                                <  1:551720/150‑1 (MQ=255)

                                                                                gTTGCTGCAGTCAAGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCTGGAGCACCTGGGTAACGcc                                                                   >  1:1619976/1‑144 (MQ=255)

                                                                                gTTGCTGCAGTCAAGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCTGGAGCACCTGGGTAACGcc                                                                   <  2:1619976/144‑1 (MQ=255)

                                                                                gTTGCTGCAGTCAAGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCTGGAGCACCTGGGTAACGCCAAgcgc                                                             <  1:1849566/150‑1 (MQ=255)

                                                                                     tgcAGTCAAAGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCTGGAGCACCTGGGTAACGCCAAGCGCGTCAt                                                        >  2:410952/1‑150 (MQ=255)

                                                                                      gcAGTCAAGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCTGGAGCACCTGGGTAACGCCAAGCGCGTCATc                                                       >  2:1861631/1‑150 (MQ=255)

                                                                                      gcAGTCAAGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCTGGAGCACCTGGGTAACGCCAAGCGCGTCATc                                                       <  1:2402728/150‑1 (MQ=255)

                                                                                       cAGTCAAGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCTGGAGCACCTGGGTAACGCCAAgcgc                                                             >  1:195973/1‑143 (MQ=255)

                                                                                       cAGTCAAGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCTGGAGCACCTGGGTAACGCCAAgcgc                                                             <  2:195973/143‑1 (MQ=255)

                                                                                        aGTCAAGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCTGGAGCACCTGGGTAACGCCAAGCGCGTCATCCt                                                     >  2:1524099/1‑150 (MQ=255)

                                                                                          tCAAGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCa                                                                                                                              <  2:806247/75‑1 (MQ=255)

                                                                                          tCAAGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCa                                                                                                                              >  1:806247/1‑75 (MQ=255)

                                                                                          tCAAGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCAAACTGACGGAGCACGTTGTCCTGACTGGCGGCCAGGTCCAATCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCTGGAGCACCTGGGTAACGCCAAGCGCGTCATCCTGt                                                   <  1:1116536/150‑1 (MQ=255)

                                                                                          tCAAGACTCCGAGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCTGGAGCACCTGGGTAACGCCAAGCGCGTCATCCTGt                                                   >  2:1167181/1‑150 (MQ=255)

                                                                                           cAAGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCTGGAGCACCTGGGTAACGCCAAGCGCGTCATCCTGTc                                                  <  1:942852/150‑1 (MQ=255)

                                                                                           cAAGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCTGGAGCACCTGGGTAACGCCAAGCGCGTCATCCTGTc                                                  <  1:1630006/150‑1 (MQ=255)

                                                                                             aGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCTGGAGCACCTGGGTAACGCCAAGCGCGTCATCCTGTCCa                                                >  1:2036461/1‑150 (MQ=255)

                                                                                               gCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCTGGAGCACCTGGGTAACGCCAAGCGCGTCATCCTGTCCAAg                                              <  2:732948/150‑1 (MQ=255)

                                                                                               gCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCAGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCTGGAGCACCTGGGTAACGCCAAGCGCGTCATCCTGTCCAAg                                              >  1:1342550/1‑150 (MQ=255)

                                                                                                 tCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCAGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCTGGAGCACCTGGGTAACGCCAAGCGCGTCATCCTGGCCAAggg                                            >  2:1237865/1‑149 (MQ=255)

                                                                                                    gggCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCTGGAGCACCTGGGTAACGCCAAGCGCGTCATCCTGTCCAAGGaaaa                                         >  1:1479357/1‑150 (MQ=255)

                                                                                                      gCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTTGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCTGGAGCACCTGGGTAACGCCAAGCGCGTCATCCTGTCCAAGGAAAaca                                       <  1:143510/150‑1 (MQ=255)

                                                                                                       cTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCTGGAGCACCTGGGTAACGCCAAGCGCGTCATCCTGTCCAAGGAAAacac                                      <  1:1059538/150‑1 (MQ=255)

                                                                                                       cTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCTGGAGCACCTGGGTAACGCCAAGCGCGTCATCCTGTCCAAGGAAAacac                                      <  2:2036461/150‑1 (MQ=255)

                                                                                                       cTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCTGGAGCACCTGGGTAACGCCAAGCGCGTCATCCTGTCCAAGGAAAacac                                      <  1:1569570/150‑1 (MQ=255)

                                                                                                        ttCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCCATGCAGCATCACCTGGGTAACGCCAAGCGCGTCATCCTGTCCAAGGAAAAcacc                                     >  1:145956/1‑150 (MQ=255)

                                                                                                         tCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCTGGAGCACCTGGGTAACGCCAAGCGCGTCATCCTGTCCAAGGAAAAcacca                                    <  2:902622/150‑1 (MQ=255)

                                                                                                            gCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCTGGAGCACCTGGGTAACGCCAAGCGCGTAATCCTGTCCAAGGAAAACaccaccc                                 >  2:825646/1‑149 (MQ=255)

                                                                                                                tCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCTGGAGCACCTGGGTAACGCCAAGCGCGTCATCCTGTCCAAGGAAAACACCACcatcat                             >  2:1847426/1‑150 (MQ=255)

                                                                                                                     gcAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCTGGAGCACCTGGGTAACGCCAAGCGCGTCATCCTGTCCAAGGAAAACACCACCATCATCGACg                        <  1:1574806/150‑1 (MQ=255)

                                                                                                                     gcAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCTGGAGCACCTGGGTAACGCCAAGCGCGTCATCCTGTCCAAGGAAAACACCACCATCATCGACg                        >  1:1650453/1‑150 (MQ=255)

                                                                                                                     gcAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCTGGAGCACCTGGGTAACGCCAAGCGCGTCATCCTGTCCAAGGAAAACACCACCATCATCGACg                        <  2:2007723/150‑1 (MQ=255)

                                                                                                                     gcAAGGCCATGCTGCAGGACATCGCTGTCCAGACTGTCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCTGGAGCACCTGGGTAACGCCAAGCGCGTCATCCTGTCCAAGGAAAACACCACCATCATCGACg                        <  2:187946/150‑1 (MQ=255)

                                                                                                                         ggCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCTGGAGCACCTGGGTAACGCCAAGCGCGTCATCCTGTCCAAGGAAAACACCACCATCATCGACGGTGc                    <  2:1332929/150‑1 (MQ=255)

                                                                                                                               gctgcAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCTGGAGCACCTGGGTAACGCCAAGCGCGTCATCCTGTCCAAGGAAAACACCACCATCATCGACGGTGCTGgcgc              <  1:1481184/150‑1 (MQ=255)

                                                                                                                               cagtcAGGACAGCGATGTCCTGAATGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCTGGAGCACCTGGGTAACGCCAAGCGCGTCATCCTGTCCAAGGAAAACACCACCATCATCGACGGTGCTGgcgc              >  2:1183299/5‑150 (MQ=255)

                                                                                                                                   cAGGACAGCGATGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCTGGAGCACCTGGGTAACGCCAAGCGCGTCATCCTGTCCAAGGAAAACACCACCATCATCGACGGTGCTGGCGCt             <  2:2441188/147‑1 (MQ=255)

                                                                                                                                   cAGGACAGCGATGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCTGGAGCACCTGGGTAACGCCAAGCGCGTCATCCTGTCCAAGGAAAACACCACCATCATCGACGGTGCTGGCGCt             >  1:2441188/1‑147 (MQ=255)

                                                                                                                                    aGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCAGCCTGTCCCTGGAAACCGCTACCCTGGAGCACCTGGGTAACGCCAAGCGCGTCATCCTGTCCAAGGAAAACACCACCATCATCGACGGTGCTGGCGCTGaca         >  1:1133408/1‑150 (MQ=255)

                                                                                                                                      gACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCTGGAGCACCTGGGTAACGCCAAGCGCGTCATCCTGTCCAAGGAAAACACCACCATCATCGACGGTGCTGGCGCTGACAcc       <  2:817554/150‑1 (MQ=255)

                                                                                                                                         aTCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCTGGAGCACCTGGGTAACGCCAAGCGCGTCATCCTGTCCAAGGAAAACACCACCATCATCGACGGTGGTGGCGCTGACACCgag    >  2:304529/1‑150 (MQ=255)

                                                                                                                                         aTCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCTGGAGCACCTGGGTAACGCCAAGCGCGTCATCCTGTCCAAGGAAAACACCACCATCATCGACGGTGCTGGCCCGGACACCgag    >  2:61470/1‑150 (MQ=255)

                                                                                                                                          tCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCTGGAGCACCTGGGTAACGCCAAGCGCGTCATCCTGTCCAAGGAAAACACCACCATCATCGACGGTGCTGGCGCTGACACCgaga   <  1:699805/150‑1 (MQ=255)

                                                                                                                                           cGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCTGGAGCACCTGGGTAACGCCAAGCGCGTCATCCTGTCCAAGGCAAACACCCCCATCATCGACGGTGCTGGCGCTGACACCGAGAt  <  2:1594005/150‑1 (MQ=255)

CACTGCTGATCGTTGCCGAAGACGTTGAAGGTGAAGCGCTGGCTACCCTGGTAGTCAACAACATGCGCGGCATCGTCAAGGTTGCTGCAGTCAAGGCTCCGGGCTTCGGCGATCGCCGCAAGGCCATGCTGCAGGACATCGCTGTCCTGACTGGCGGCCAGGTCATCTCCGAAGAAATCGGCCTGTCCCTGGAAACCGCTACCCTGGAGCACCTGGGTAACGCCAAGCGCGTCATCCTGTCCAAGGAAAACACCACCATCATCGACGGTGCTGGCGCTGACACCGAGAT  >  NC_002947/1550379‑1550667

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: `‑`