BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,378,234	C→G	6.0%	S434T (AGC→ACC) ‡	dnaK ←	chaperone protein DnaK

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,378,234	0	C	G	6.0%	1153.2 / 44.2	343	S434T (AGC→ACC) ‡	dnaK	chaperone protein DnaK
Reads supporting (aligned to +/- strand): ref base C (126/196); new base G (8/13); total (134/209)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 349 total aligned reads displayed.

GGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCA  >  NC_002947/5378085‑5378383

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

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

    gCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTg                                                                                                                                                   <  2:7073360/150‑1 (MQ=255)

      cGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGt                                                                                                                                                 >  1:3513816/1‑150 (MQ=255)

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

       gTTGGCGTCGATGTCAAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGtt                                                                                                                                                <  2:4254919/150‑1 (MQ=255)

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

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

          ggCGTCGATGTCGAAGGTCACTTCGATCTGTGGAACGCCACGCGGTGCCGGCGGAATGCCAACGAGGTCGAACTGGCCGAGGATCTTGTTGTCGGCAGCCATTTCACGCTCGCCCTGGAACACGCGGATCGTCACGGCCGACTGATTgtc                                                                                                                                             <  1:59778/150‑1 (MQ=255)

             gTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTgtcgtc                                                                                                                                          <  1:1294043/150‑1 (MQ=255)

                gATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGc                                                                                                                                       <  2:6701924/150‑1 (MQ=255)

                 aTGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCa                                                                                                                                      <  2:4345821/150‑1 (MQ=255)

                 aTGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCa                                                                                                                                      >  1:4076678/1‑150 (MQ=255)

                 aTGTCGAAGGTCACTTCGATCTGTGGAACGCCACGCGGTGCCGGCGGAATGCCAACGAGGTCGAACTGGCCGAGGATCTTGTTGTCGGCAGCCATTTCACGCTCGCCCTGGAACACGCGGATCGTCACGGCCGACTGATTGTCTTCAGcg                                                                                                                                      >  1:6823426/1‑149 (MQ=255)

                   gTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTg                                                                                                                                                   <  1:2830043/135‑1 (MQ=255)

                   gTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTg                                                                                                                                                   >  2:2830043/1‑135 (MQ=255)

                      gAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCga                                                                                                                                 <  2:1402539/150‑1 (MQ=255)

                             aCTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACAcc                                                                                                                          >  1:1884304/1‑150 (MQ=255)

                              cTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCt                                                                                                                         >  2:4032261/1‑150 (MQ=255)

                               ttCGATTTGCGGCACGCCGCGTGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTg                                                                                                                        >  1:3640013/1‑150 (MQ=255)

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

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

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

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

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

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

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

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

                                 cGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCg                                                                                                                      <  1:5573790/150‑1 (MQ=255)

                                    tttGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGAct                                                                                                                   >  2:2993297/1‑150 (MQ=255)

                                       gCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGActtct                                                                                                                >  1:1015648/1‑150 (MQ=255)

                                         ggCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTg                                                                                                              >  1:7810700/1‑150 (MQ=255)

                                         ggCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTg                                                                                                              <  2:3078909/150‑1 (MQ=255)

                                             cgccgcGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCg                                                                                                          <  2:4832786/150‑1 (MQ=255)

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

                                               ccgcGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCggg                                                                                                        <  1:984795/150‑1 (MQ=255)

                                                cgcgcgGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGa                                                                                                       <  2:8005347/150‑1 (MQ=255)

                                                cgcgcgGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGa                                                                                                       <  2:1119050/150‑1 (MQ=255)

                                                cgcgcgGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGa                                                                                                       <  1:5272409/150‑1 (MQ=255)

                                                   gcgGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGAtgg                                                                                                    <  2:3739937/150‑1 (MQ=255)

                                                   gcgGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGAtgg                                                                                                    <  1:1591035/150‑1 (MQ=255)

                                                    cgGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGAtggt                                                                                                   <  2:6768281/150‑1 (MQ=255)

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

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

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

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

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

                                                           ggCGGAATGCCAACGAGGTCGAACTGGCCGAGGATCTTGTTGTCGGCAGCCATTTCACGCTCGCCCTGGAACGCGCGGATCGTCACGGCCGACTGATTGTCTTCAGCGGTCGAGAAGGTCTGCGACTTCTTGGTCGGGATCGTCGTGttg                                                                                            >  2:8219784/1‑149 (MQ=255)

                                                            gCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGttct                                                                                           >  1:790320/1‑150 (MQ=255)

                                                              ggAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGttcttc                                                                                         <  2:7650110/150‑1 (MQ=255)

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

                                                                    tCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATc                                                                                   >  2:4885082/1‑150 (MQ=255)

                                                                       gCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAgc                                                                                <  2:414669/150‑1 (MQ=255)

                                                                       gCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAgc                                                                                <  2:7928718/150‑1 (MQ=255)

                                                                        ccAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGATCTTCTCGATCAgcg                                                                               >  2:1496352/1‑150 (MQ=255)

                                                                         cAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTTCTTACGCTCACCTTGCATCACGTGAATTGTCACGGCGCTCTGGTTGTCGTGGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTATCGATCAgcgc                                                                              <  1:83208/150‑1 (MQ=255)

                                                                         cAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAgcgc                                                                              <  2:3242563/150‑1 (MQ=255)

                                                                            gTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGt                                                                           >  1:4189784/1‑150 (MQ=255)

                                                                            gTCGAACTGGCCGAGGATCTTGTTGTCGGCAGCCATTTCACGCTCGCCCTGGAACACGCGGATCGTCACGGCCGACTGATTGTCTTCAGCGGTCGAGAAGGTCTGCGACTTCTTGGTCGGGATCGTCGTGTTGCGTTCGATCAggcgggt                                                                           <  1:2387832/150‑7 (MQ=255)

                                                                                aaCTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTCCGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGtcatc                                                                       <  2:3448910/150‑1 (MQ=255)

                                                                                aaCTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGtcatc                                                                       >  2:6311324/1‑150 (MQ=255)

                                                                                 aCTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGtcatca                                                                      <  1:4313655/150‑1 (MQ=255)

                                                                                  cTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCAc                                                                     <  1:5149213/150‑1 (MQ=255)

                                                                                   ttGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCAcg                                                                    >  2:1220706/1‑150 (MQ=255)

                                                                                     gCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCga                                                                                                                                 >  2:7793185/1‑87 (MQ=255)

                                                                                        cAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACgccgcc                                                               >  1:3396691/1‑150 (MQ=255)

                                                                                         aGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACgccgcct                                                              >  2:3544617/1‑149 (MQ=255)

                                                                                            gACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATg                                                           <  1:4920696/150‑1 (MQ=255)

                                                                                                tGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGttt                                                       <  2:8946374/150‑1 (MQ=255)

                                                                                                tGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGttt                                                       <  1:6306296/150‑1 (MQ=255)

                                                                                                 gTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTc                                                      <  2:319136/150‑1 (MQ=255)

                                                                                                   tCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGa                                                    <  1:2108455/150‑1 (MQ=255)

                                                                                                   tCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGa                                                    <  1:2205409/150‑1 (MQ=255)

                                                                                                   tCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGa                                                    <  1:7331589/150‑1 (MQ=255)

                                                                                                    cTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGAt                                                   <  1:9070121/150‑1 (MQ=255)

                                                                                                     tGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATa                                                  >  2:3846985/1‑150 (MQ=255)

                                                                                                          aGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAg                                             <  1:8460319/150‑1 (MQ=255)

                                                                                                           gCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAgg                                            >  2:8866229/1‑150 (MQ=255)

                                                                                                           gCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAgg                                            <  2:2381398/150‑1 (MQ=255)

                                                                                                           gCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAgg                                            <  2:2335131/150‑1 (MQ=255)

                                                                                                           gCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCAATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAgg                                            <  1:6523355/150‑1 (MQ=255)

                                                                                                                cTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGttctt                                                                                          <  2:1905778/99‑1 (MQ=255)

                                                                                                                 ttACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGc                                      <  1:1782984/150‑1 (MQ=255)

                                                                                                                 ttACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGc                                      <  1:4489895/150‑1 (MQ=255)

                                                                                                                     gCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGc                                  >  2:254861/1‑150 (MQ=255)

                                                                                                                       tCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTg                                >  2:6932651/1‑150 (MQ=255)

                                                                                                                       tCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTg                                <  1:3064695/150‑1 (MQ=255)

                                                                                                                         aCCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGAc                              >  1:805546/1‑150 (MQ=255)

                                                                                                                         aCCTTGCAGCACGAGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTACTTGGTCGGGATGGTGGTGTACTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGAc                              <  2:5917111/150‑1 (MQ=255)

                                                                                                                           cTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGt                            >  2:8432592/1‑150 (MQ=255)

                                                                                                                             tGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTcc                          <  1:6152599/150‑1 (MQ=255)

                                                                                                                             tGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCCCGCTGAccctg                          <  2:5131239/150‑5 (MQ=255)

                                                                                                                              gcagcaCGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCCCGCTGAccctgg                         <  2:6507178/150‑6 (MQ=255)

                                                                                                                               cagcaCGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCcag                        <  1:198416/150‑1 (MQ=255)

                                                                                                                                  caCGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCcagcag                     >  2:395330/1‑150 (MQ=255)

                                                                                                                                          tcgTCACGGCCGACTGATTGTCTTCAGCGGTCGAGAAGGTCTGCGACTTCTTGGTCGGGATCGTCGTGTTGCGTTCGATCAGGCGGGTGAACACGCCACCGAGCGTTTCAATGCCGAGCGAAAGCGGGGTCACGTCGAGCAGCAGAACGt             >  1:4374251/3‑150 (MQ=255)

                                                                                                                                           cgTCACGGCCGACTGATTGTCTTCAGCGGTCGAGAAGGTCTGCGACTTCTTGGTCGGGATCGTCGTGTTGCGTTCGATCAGGCGGGTGAACACGCCACCGAGCGTTTCAATGCCGAGCGAAAGCGGGGTCACGTCGAGCAGCTGAACGTc            <  2:2131293/149‑1 (MQ=255)

                                                                                                                                            gTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCTGCTGCTGGACGTca           <  2:2226383/150‑2 (MQ=255)

                                                                                                                                             tCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCTGCTGCTGGACGTcag          <  1:3805746/150‑3 (MQ=255)

                                                                                                                                               aCGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTc        <  2:6067863/150‑1 (MQ=255)

GGATGCCGTTGGCGTCGATGTCGAAGGTCACTTCGATTTGCGGCACGCCGCGCGGTGCTGGCGGAATGTCAGCCAGGTCGAACTTGCCCAGCGACTTGTTCTGCGCAGCCTGCTTACGCTCACCTTGCAGCACGTGAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCA  >  NC_002947/5378085‑5378383

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