BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,378,158	C→G	7.4%	L459L (CTG→CTC)	dnaK ←	chaperone protein DnaK

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,378,158	0	C	G	7.4%	709.0 / 35.0	216	L459L (CTG→CTC)	dnaK	chaperone protein DnaK
Reads supporting (aligned to +/- strand): ref base C (110/90); new base G (7/9); total (117/99)
Fisher's exact test for biased strand distribution p-value = 4.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.29e-01

Only 100 of 218 total aligned reads displayed.

CCGACAGGCCGGAGTTGGCCTTGATCACGATCGACTGAGTCTTGCCGGTGGCCTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCG  >  NC_002947/5378010‑5378306

  gACAGGCCGGAGTTGGCCTTGATCACGATCGACTGAGTCTTGCCGGTGGCCTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAgg                                                                                                                                                   >  2:383848/1‑150 (MQ=255)

    cAGGCCGGAGTTGGCCTTGATCACGATCGACTGAGTCTTGCCGGTGGCCTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTc                                                                                                                                                 <  1:1511875/150‑1 (MQ=255)

       gCCGGAGTTGGCCTTGATCACGATCGACTGAGTCTTGCCGGTGGCCTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGaa                                                                                                                                              >  2:5030180/1‑150 (MQ=255)

        ccGGAGTTGGCCTTGATCACGATCGACTGAGTCTTGCCGGTGGCCTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAAc                                                                                                                                             >  1:1458800/1‑150 (MQ=255)

              ttGGCCTTGATCACGATCGACTGAGTCTTGCCGGTGGCCTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGCACTTGccc                                                                                                                                       >  2:221667/1‑150 (MQ=255)

              ttGGCCTTGATCACGATCGACTGAGTCTTGCCGGTGGCCTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGccc                                                                                                                                       >  1:636455/1‑150 (MQ=255)

                 gCCTTGATCACGATCGACTGAGTCTTGCCGGTGGCCTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGc                                                                                                                                    >  2:960691/1‑150 (MQ=255)

                 gCCTTGATCACGATCGACTGAGTCTTGCCGGTGGCCTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGc                                                                                                                                    >  2:351911/1‑150 (MQ=255)

                  ccTTGATCACGATCGACTGAGTCTTGCCGGTGGCCTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCg                                                                                                                                   >  2:1022832/1‑150 (MQ=255)

                   cTTGATCACGATCGACTGAGTCTTGCCGGTGGCCTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGa                                                                                                                                  >  2:3090941/1‑150 (MQ=255)

                   cTTGATCACGATCGACTGAGTCTTGCCGGTGGCCTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGa                                                                                                                                  >  1:2329558/1‑150 (MQ=255)

                     tGATCACGATCGACTGAGTCTTGCCGGTGGCCTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACt                                                                                                                                <  2:127694/150‑1 (MQ=255)

                       aTCACGATCGACTGAGTCTTGCCGGTGGCCTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTg                                                                                                                              <  2:623984/150‑1 (MQ=255)

                       aTCACGATCGACTGAGTCTTGCCGGTGGCCTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTg                                                                                                                              <  1:2465300/150‑1 (MQ=255)

                       aTCACGATCGACTGAGTCTTGCCGGTGGCCTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTg                                                                                                                              <  2:1519438/150‑1 (MQ=255)

                         cACGATCGACTGAGTCTTGCCGGTGGCCTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGtt                                                                                                                            >  1:951746/1‑150 (MQ=255)

                          aCGATCGACTGAGTCTTGCCGGTGGCCTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTc                                                                                                                           >  2:1227248/1‑150 (MQ=255)

                           cGATCGACTGAGTCTTGCCGGTGGCCTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCt                                                                                                                          <  1:4639916/150‑1 (MQ=255)

                            gATCGACTGAGTCTTGCCGGTGGCCTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTg                                                                                                                         >  2:3924795/1‑150 (MQ=255)

                            gATCGACTGAGTCTTGCCGGTGGCCTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTg                                                                                                                         >  1:5448859/1‑150 (MQ=255)

                             ctCGACTGAGTCTTGCCGGTGGCCTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTgc                                                                                                                        <  1:3148498/149‑1 (MQ=255)

                              tCGACTGAGTCTTGCCGGTGGCCTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTgcg                                                                                                                       <  1:4081090/150‑1 (MQ=255)

                                  gtgttccTTGCCTGTGCCCTTGTCCTTTGCGGACACGTTCACGATACCGTTGGCGTCGATGTCGAAGGTCACTTCGATCTGTGGAACGCCACGCGGTGCCGGCGGAATGCCAACGAGGTCGAACTGGCCGAGGATCTTGTTGTCGGCAGc                                                                                                                   >  1:4976651/7‑150 (MQ=255)

                                   tGAGTCTTGCCGGTGGCCTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGcc                                                                                                                  >  1:2767857/1‑150 (MQ=255)

                                   tGAGTCTTGCCGGTGGCCTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGcc                                                                                                                  >  1:4053634/1‑150 (MQ=255)

                                      gTCTTGCCGGTGGCCTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGc                                                                                                               <  1:3973712/150‑1 (MQ=255)

                                       tCTTGCCGGTGGCCTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCt                                                                                                              >  2:4078634/1‑150 (MQ=255)

                                         ttGCCTGTGCCCTTGTCCTTTGCGGACACGTTCACGATACCGTTGGCGTCGATGTCGAAGGTCACTTCGATCTGTGGAACGCCACGCGGTGCCGGCGGAATGCCAACGAGGTCGAACTGGCCGAGGATCTTGTTGTCGGCAGccatttca                                                                                                            >  2:2424572/1‑144 (MQ=255)

                                         ttGCCGGTGGCCTTGCCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTa                                                                                                            >  2:4462899/1‑150 (MQ=255)

                                          tGCCGGTGGCCTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTAc                                                                                                           <  1:4349243/150‑1 (MQ=255)

                                                  gCCTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCt                                                                                                   <  1:383848/150‑1 (MQ=255)

                                                  gCCTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCt                                                                                                   >  2:2441368/1‑150 (MQ=255)

                                                    cTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTg                                                                                                 >  1:100934/1‑150 (MQ=255)

                                                    cTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTg                                                                                                 >  1:4275496/1‑150 (MQ=255)

                                                    cTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTg                                                                                                 >  1:4867909/1‑150 (MQ=255)

                                                     ttGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTgc                                                                                                <  2:2163465/150‑1 (MQ=255)

                                                     ttGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTgc                                                                                                >  1:98178/1‑150 (MQ=255)

                                                      tGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTgca                                                                                               <  1:2956362/150‑1 (MQ=255)

                                                      tGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTgca                                                                                               <  2:4981605/150‑1 (MQ=255)

                                                        tCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGTTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTgcagc                                                                                             >  1:1062923/1‑150 (MQ=255)

                                                          tttCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCAc                                                                                           >  1:191831/1‑150 (MQ=255)

                                                               cACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGaa                                                                                      <  2:1823170/150‑1 (MQ=255)

                                                                aCCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAAt                                                                                     <  2:2329558/150‑1 (MQ=255)

                                                                      gTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCAc                                                                               >  1:3796675/1‑150 (MQ=255)

                                                                       tGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTgc                                                                                                >  1:198669/1‑132 (MQ=255)

                                                                       tGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTgc                                                                                                <  2:198669/132‑1 (MQ=255)

                                                                        gCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGACTGGTCACgg                                                                             >  2:4199764/1‑150 (MQ=255)

                                                                        gCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACgg                                                                             >  1:676058/1‑150 (MQ=255)

                                                                            gATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGct                                                                         >  1:3248575/1‑150 (MQ=255)

                                                                             aTGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGaa                                                                                                                                              <  1:4117737/80‑1 (MQ=255)

                                                                             aTGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGaa                                                                                                                                              >  2:4117737/1‑80 (MQ=255)

                                                                               gCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTgc                                                                                                                        >  2:5303098/1‑100 (MQ=255)

                                                                                ccGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTgg                                                                     <  1:3742465/150‑1 (MQ=255)

                                                                                ccGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTgg                                                                     >  1:5200124/1‑150 (MQ=255)

                                                                                  gTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGtt                                                                   >  1:5332562/1‑150 (MQ=255)

                                                                                  gTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGtt                                                                   <  2:5406617/150‑1 (MQ=255)

                                                                                     ggCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTgtc                                                                <  2:3611848/150‑1 (MQ=255)

                                                                                      gCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTgtcg                                                               >  1:754823/1‑150 (MQ=255)

                                                                                      gCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTgtcg                                                               >  1:4152048/1‑150 (MQ=255)

                                                                                         tCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGtcgtcg                                                            <  2:5240399/150‑1 (MQ=255)

                                                                                         tCGATGTCGAAGGTCACTTCGATCTGTGGAACGCCACGCGGTGCCGGCGGAATGCCAACGAGGTCGAACTGGCCGAGGATCTTGTTGTCGGCAGCCATTTCACGCTCGCCCTGGAACACGCGGATCGTCACGGCCGACTGATTgtcttca                                                            <  2:1653256/150‑2 (MQ=255)

                                                                                               tCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTc                                                      >  2:4015637/1‑150 (MQ=255)

                                                                                               tCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTc                                                      <  1:3534648/150‑1 (MQ=255)

                                                                                               tCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTc                                                      <  2:609203/150‑1 (MQ=255)

                                                                                                cGAAGGTCACTTCGATGTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCg                                                     <  2:4352931/150‑1 (MQ=255)

                                                                                                 gAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTc                                                                                                                           >  2:1052796/1‑79 (MQ=255)

                                                                                                 gAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTgca                                                                                               >  2:4874392/1‑107 (MQ=255)

                                                                                                 gAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTgca                                                                                               <  1:4874392/107‑1 (MQ=255)

                                                                                                 gAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCga                                                    >  2:1192294/1‑150 (MQ=255)

                                                                                                 gAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCga                                                    >  2:4431275/1‑150 (MQ=255)

                                                                                                      tCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAaca                                               >  2:2829555/1‑150 (MQ=255)

                                                                                                      tCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAaca                                               >  1:4489759/1‑150 (MQ=255)

                                                                                                       cACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAacac                                              >  1:4694574/1‑150 (MQ=255)

                                                                                                          ttCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCTGGTCGAACTTGCCCAGCGAATAGTTCTCCGCCGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTg                                           >  2:2763355/1‑150 (MQ=255)

                                                                                                          ttCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTg                                           >  1:3357684/1‑150 (MQ=255)

                                                                                                          ttCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGACCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACATTGCAGCACGTGAATGGTCACGGGGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTg                                           >  2:3672717/1‑150 (MQ=255)

                                                                                                           tCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGc                                          <  1:1031645/150‑1 (MQ=255)

                                                                                                             gATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGa                                        <  2:3948635/150‑1 (MQ=255)

                                                                                                                 tGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGActtc                                    <  2:3414053/150‑1 (MQ=255)

                                                                                                                   cGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGActtctt                                  <  1:2019787/150‑1 (MQ=255)

                                                                                                                   cGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGActtctt                                  <  2:1611958/150‑1 (MQ=255)

                                                                                                                     gCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTc                                                      >  1:11832/1‑128 (MQ=255)

                                                                                                                       aCGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTc                              >  1:2201014/1‑150 (MQ=255)

                                                                                                                         gccgcGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCgg                            >  2:518669/1‑150 (MQ=255)

                                                                                                                           cgcgcgGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGGCTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGa                          <  2:2959103/150‑1 (MQ=255)

                                                                                                                                ggTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGAtggtg                     >  2:3763025/1‑150 (MQ=255)

                                                                                                                                ggTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGAtggtg                     <  1:3593595/150‑1 (MQ=255)

                                                                                                                                 gTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGAtggtgg                    <  2:2890598/150‑1 (MQ=255)

                                                                                                                                   gcgggCGGACTGTCAGCCAGGTCGACCGTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCGTGCAGCACGTGACTGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGgtggtg                  <  2:4767440/147‑1 (MQ=255)

                                                                                                                                   gCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGgtggtg                  <  1:130684/150‑1 (MQ=255)

                                                                                                                                     tGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTgca                                                                                               >  2:2905534/1‑71 (MQ=255)

                                                                                                                                     tGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTgca                                                                                               <  1:2905534/71‑1 (MQ=255)

                                                                                                                                     tGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGtt                >  2:2702662/1‑150 (MQ=255)

                                                                                                                                      ggCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGttc               <  1:5030180/150‑1 (MQ=255)

                                                                                                                                       gCGGAATGCCAACGAGGTCGAACTGGCCGAGGATCTTGTTGTCGGCAGCCATTTCACGCTCGCCCTGGAACACGCGGATCGTCACGGCCGACTGATTGTCTTCAGCGGTCGAGAAGGTCTGCGACTTCTTGGTCGGGATCGTCGTGttgc              >  1:2789669/1‑148 (MQ=255)

                                                                                                                                       gCGGAATGCCAACGAGGTCGAACTGGCCGAGGATCTTGTTGTCGGCAGCCATTTCACGCTCGCCCTGGAACACGCGGATCGTCACGGCCGACTGATTGTCTTCAGCGGTCGAGAAGGTCTGCGACTTCTTGGTCGGGATCGTCGTGttgc              <  1:1080889/150‑3 (MQ=255)

                                                                                                                                          gAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTtcttct           <  1:4766858/150‑1 (MQ=255)

                                                                                                                                            aTGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCg         <  2:1385888/150‑1 (MQ=255)

                                                                                                                                             tGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGa        <  1:2401906/150‑1 (MQ=255)

                                                                                                                                             tGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGa        >  1:1722147/1‑150 (MQ=255)

CCGACAGGCCGGAGTTGGCCTTGATCACGATCGACTGAGTCTTGCCGGTGGCCTTGTCTTTCGCACCTACGTGCAGGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCG  >  NC_002947/5378010‑5378306

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