BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	6,169,837	T→A	5.2%	S172C (AGC→TGC)	PP_5408 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	6,169,837	0	T	A	5.2%	473.9 / 12.7	154	S172C (AGC→TGC)	PP_5408	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (65/80); new base A (2/6); total (67/87)
Fisher's exact test for biased strand distribution p-value = 4.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.73e-01

Only 100 of 160 total aligned reads displayed.

AACAATTTCATTGTCTATGTTTTTTAGTTTGCTTTCGAGCACCCGTAGAACCGTTTCAAGTGCAAGGTGGTGCTCAAGGAGGCAGATCCTTCCGTAGAATTGTGCCAGCCGCAATGCCAATGGGTCCCTGTTGGCCGTAGGCCTTGCGCTGCCTAGGCTGCGCCCGCAACATCCGCTCCACATTACCTGGACAGAGTGGCCATCGCGAGAAAATGGCTGCTCACAGAACTGGCAGTACGAGAGCAAATCGACGTTGTGCTTCGCGCACACAGTTACATGTGGGTGAAGACGTCGCCA  >  NC_002947/6169688‑6169984

aacaaTTTCATTGTCTATGTTTTTTAGTTTGCTTTCGAGCACCCGTAGAACCGTTTCAAGTGCAAGGTGGTGCTCAAGGAGGCAGATCCTTCCGTAGAATTGTGCCAGCCGCAATGCCAATGGGTCGCTGTTGGCCGTAGGCCTTGCGCt                                                                                                                                                     <  2:1182554/150‑1 (MQ=255)

aacaaTTTCATTGTCTATGTTTTTTAGTTTGCTTTCGAGCACCCGTAGAACCGTTTCAAGTGCAAGGTGGTGCTCAAGGAGGCAGATCCTTCCGTAGAATTGTGCCAGCCGCAATGCCAATGGGTCCCTGTTGGCCGTAGGCCTTGCGCt                                                                                                                                                     >  2:1533277/1‑150 (MQ=255)

aacaaTTTCATTGTCTATGTTTTTTAGTTTGCTTTCGAGCACCCGTAGAACCGTTTCAAGTGCAAGGTGGTGCTCAAGGAGGCAGATCCTTCCGTAGAATTGTGCCAGCCGCAATGCCAATGGGTCCCTGTTGGCCGTAGGCCTTGCGCt                                                                                                                                                     <  1:1408990/150‑1 (MQ=255)

     tttCATTGTCTATGTTTTTTAGTTTGCTTTCGAGCACCCGTAGAACCGTTTCAAGTGCAAGGTGGTGCTCAAGGAGGCAGATCCTTCCGTAGAATTGTGCCAGCCGCAATGCCAATGGGTCCCTGTTGGCCGTAGGCCTTGCGCTGCCTa                                                                                                                                                <  1:859714/150‑1 (MQ=255)

     tttCATTGTCTATGTTTTTTAGTTTGCTTTCGAGCACCCGTAGAACCGTTTCAAGTGCAAGGTGGTGCTCAAGGAGGCAGATCCTTCCGTAGAATTGTGCCAGCCGCAATGCCAATGGGTCCCTGTTGGCCGTAGGCCTTGCGCTGCCTa                                                                                                                                                >  1:183881/1‑150 (MQ=255)

         aTTGTCTATGTTTTTTCGTTTGCTTTCGAGCACCCGTAGAACCGTTTCAAGTGCAAGGTGGTGCTCAAGGAGGCAGATCCTTCCGTAGAATTGTGCCAGCCGCAATGCCAATGGGTCCCTGTTGGCCGTAGGCCTTGCGCTGCCTAGGCt                                                                                                                                            <  2:1839288/150‑1 (MQ=255)

           tGTCTATGTTTTTTAGTTTGCTTTCGAGCACCCGTAGAACCGTTTCAAGTGCAAGGTGGTGCTCAAGGAGGCAGATCCTTCCGTAGAATTGTGCCAGCCGCAATGCCAATGGGTCCCTGTTGGCCGTAGGCCTTGCGCTGCCTAGGCTgc                                                                                                                                          <  1:2337625/150‑1 (MQ=255)

               tATGTTTTTTAGTTTGCTTTCGAGCACCCGTAGAACCGTTTCAAGTGCAAGGTGGTGCTCAAGGAGGCAGATCCTTCCGTAGAATTGTGCCAGCCGCAATGCCAATGGGTCCCTGTTGGCCGTAGGCCTTGCGCTGCCTAGGCTGCGccc                                                                                                                                      >  2:814833/1‑150 (MQ=255)

               tATGTTTTTTAGTTTGATTTCGAGAGCACGTAGCGTCGTGTCAAGGGCAAGGTGTTGCTCAAGGAGGCAGATCCTTCCGTAGAATTGTGCCAGCCGCAATGCCAATGGGTCCCTGTTGGCCGTAGGCCTTGCGCAGCCTAGGCAgcgcaa                                                                                                                                      <  2:857389/150‑3 (MQ=255)

                aTGTTTTTTAGTTTGCTTTCGAGCACCCGTAGAACCGTTTCAAGTGCAAGGTGGTGCTCAAGGAGGCAGATCCTTCCGTAGAATTGTGCCAGCCGCAATGCCAATGGGTCCCTGTTGGCCGTAGGCCTTGCGCTGcc                                                                                                                                                  <  1:105108/137‑1 (MQ=255)

                aTGTTTTTTAGTTTGCTTTCGAGCACCCGTAGAACCGTTTCAAGTGCAAGGTGGTGCTCAAGGAGGCAGATCCTTCCGTAGAATTGTGCCAGCCGCAATGCCAATGGGTCCCTGTTGGCCGTAGGCCTTGCGCAGCCTAGGCAgcgcaag                                                                                                                                     <  1:781154/150‑4 (MQ=255)

                     ttttAGTTTGCTTTCGAGCACCCGTAGAACCGTTTCAAGTGCAAGGTGGTGCTCAAGGAGGCAGATCCTTCCGTAGAATTGTGCCAGCCGCAATGCCAATGGGTCCCTGTTGGCCGTAGGCCTTGCGCTGCCTAGGCTGCGCCCGCAACa                                                                                                                                >  1:328760/1‑150 (MQ=255)

                       ttAGTTTGCTTTCGAGCACCCGTAGAACCGTTTCAAGTGCAAGGTGGTGCTCAAGGAGGCAGATCCTTCCGTAGAATTGTGCCAGCCGCAATGCCAATGGGTCCCTGTTGGCCGTAGGCCTTGCGCTGCCTAGGCTGCGCCCGCAACATc                                                                                                                              <  1:551297/150‑1 (MQ=255)

                            ttGCTTTCGAGCACCCGTAGAACCGTTTCAAGTGCAAGGTGGTGCTCAAGGAGGCAGATCCTTCCGTAGAATTGTGCCAGCCGCAATGCCAATGGGTCCCTGTTGGCCGTAGGCCTTGCGCTGCCTAGGCTGCGCCCGCAACATCCGCTc                                                                                                                         <  1:2420049/150‑1 (MQ=255)

                                   cGAGCACCCGTAGAACCGTTTCAAGTGCACGGTGGTCCCCAAGGCCCCACAGCCGTCCGTAGAATTGTGCCAGACGCAATGCCAATGGGTCCCTGTTGGCCGTAGGCCTTGCGCTGCCTAGGCTGCGCCCGCAACATCCGCTCCACATTa                                                                                                                  <  1:2377733/150‑1 (MQ=255)

                                    gAGCACCCGTAGAACCGTTTCAAGTGCAAGGTGGTGCTCAAGGAGGCAGATCCTTCCGTAGAATTGTGCCAGCCGCAATGCCAATGGGTCCCTGTTGGCCGTAGGCCTTGCGCTGCCTAGGCTGCGCCCGCAACATCCGCTCCACATTAc                                                                                                                 >  1:833807/1‑150 (MQ=255)

                                          ccGTAGAACCGTTTCAAGTGCAAGGTGGTGCTCAAGGAGGCAGATCCTTCCGTAGAATTGTGCCAGCCGCAATGCCAATGGGTCCCTGTTGGCCGTAGGCCTTGCGCTGCCTAGGCTGCGCCCGCAACATCCGCTCCACATTACCTg                                                                                                              >  2:525465/1‑147 (MQ=255)

                                          ccGTAGAACCGTTTCAAGTGCAAGGTGGTGCTCAAGGAGGCAGATCCTTCCGTAGAATTGTGCCAGCCGCAATGCCAATGGGTCCCTGTTGGCCGTAGGCCTTGCGCTGCCTAGGCTGCGCCCGCAACATCCGCTCCACATTACCTg                                                                                                              <  1:525465/147‑1 (MQ=255)

                                           cGTAGAACCGTTTCAAGTGCAAGGTGGTGCTCAAGGAGGCAGATCCTTCCGTAGAATTGTGCCAGCCGCAATGCCAATGGGTCCCTGTTGGCCGTAGGCCTTGCGCTGCCTAGGCTGCGCCCGCAACATCCGCTCCACATTACCTGGACa                                                                                                          <  1:1533277/150‑1 (MQ=255)

                                           cGTAGAACCGTTTCAAGTGCAAGGTGGTGCTCAAGGAGGCAGATCCTTCCGTAGAATTGTGCCAGCCGCAATGCCAATGGGTCCCTGTTGGCCGTAGGCCTTGCGCTGCCTAGGCTGCGCCCGCAACATCCGCTCCACATTACCTGGACa                                                                                                          <  2:779088/150‑1 (MQ=255)

                                           cGTAGAACCGTTTCAAGTGCAAGGTGGTGCTCAAGGAGGCAGATCCTTCCGTAGAATTGTGCCAGCCGCAATGCCAATGGGTCCCTGTTGGCCGTAGGCCTTGCGCTGCCTAGGCTGCGCCCGCAACATCCGCTCCACATTACCTGGACa                                                                                                          <  2:96547/150‑1 (MQ=255)

                                              aGAACCGTTTCAAGTGCAAGGTGGTGCTCAAGGAGGCAGATCCTTCCGTAGAATTGTGCCAGCCGCAATGCCAATGGGTCCCTGTTGGCCGTAGGCCTTGCGCTGCCTAGGCTGCGCCCGCAACATCCGCTCCACATTACCTGGACagag                                                                                                       <  1:1786594/150‑1 (MQ=255)

                                                aaCCGTTTCAAGTTCAAGGTGGTGCGCAAGGACGCAGATCCTTCCGTAGAATTGTGCCAGCCGCAATTCCAATGGGTCCCTGTTTGCCGTAGGCCTTGCGCTGCCTAGGCTGCGTCCGTAACATCCGCTCCACATTACCTGGACAGAGTg                                                                                                     <  2:1076486/150‑1 (MQ=255)

                                                aaCCGTTTCAAGTGCAAGGTGGTGCTCAAGGAGGCAGATCCTTCCGTAGAATTGTGCCAGCCGCAATGCCAATGGGTCCCTGTTGGCCGTAGGCCTTGCGCTGCCTAGGCTGCGCCCGCAACATCCGCTCCACATTACCt                                                                                                               <  2:1144828/140‑1 (MQ=255)

                                                aaCCGTTTCAAGTGCAAGGTGGTGCTCAAGGAGGCAGATCCTTCCGTAGAATTGTGCCAGCCGCAATGCCAATGGGTCCCTGTTGGCCGTAGGCCTTGCGCTGCCTAGGCTGCGCCCGCAACATCCGCTCCACATTACCTGGACAGAGTg                                                                                                     >  2:1861802/1‑150 (MQ=255)

                                                aaCCGTTTCAAGTGCAAGGTGGTGCTCAAGGAGGCAGATCCTTCCGTAGAATTGTGCCAGCCGCAATGCCAATGGGTCCCTGTTGGCCGTAGGCCTTGCGCTGCCTAGGCTGCGCCCGCAACATCCGCTCCACATTACCTGGACAGAGTg                                                                                                     <  1:1169897/150‑1 (MQ=255)

                                                  ccGTTTCAAGTGCAAGGTGGTGCTCAAGGAGGCAGATCCTTCCGTAGAATTGTGCCAGCCGCAATGCCAATGGGTCCCTGTTGGCCGTAGGCCTTGCGCTGCCTAGGCTGCGCCCGCAACATCCGCTCCACATTACCTGGACAGAGTGGc                                                                                                   >  2:1131817/1‑150 (MQ=255)

                                                  ccGTTTCAAGTGCAAGGTGGTGCTCAAGGAGGCAGATCCTTCCGTAGAATTGTGCCAGCCGCAATGCCAATGGGTCCCTGTTGGCCGTAGGCCTTGCGCTGCCTAGGCTGCGCCCGCAACATCCGCTCCACATTACCTGGACAGAGTGGc                                                                                                   >  2:249120/1‑150 (MQ=255)

                                                    gTTTCAAGTGCAAGGTGGTGCTCAAGGAGGCAGATCCTTCCGTAGAATTGTGCCAGCCGCAATGCCAATGGGTCCCTGTTGGCCGTAGGCCTTGCGCTGCCTAGGCTGCGCCCGCAACATCCGCTCCACATTACCTGGACAGAGTGGCCa                                                                                                 >  1:966184/1‑150 (MQ=255)

                                                     tttCAAGTGCAAGGTGGTGCTCAAGGAGGCAGATCCTTCCGTAGAATTGTGCCAGCCGCAATGCCAATGGGTCCCTGTTGGCCGTAGGCCTTGCGCTGCCTAGGCTGCGCCCGCAACATCCGCTCCACATTACCTGGACAGAGTGGCCAt                                                                                                >  2:1356139/1‑150 (MQ=255)

                                                     tttCAAGTGCAAGGTGGTGCTCAAGGAGGCAGATCCTTCCGTAGAATTGTGCCAGCCGCAATGCCAATGGGTCCCTGTTGGCCGTAGGCCTTGCGCTGCCTAGGCTGCGCCCGCAACATCCGCTCCACATTACCTGGACAGAGTGGCCAt                                                                                                >  1:200867/1‑150 (MQ=255)

                                                     tttCAAGTGCAAGGTGGTGCTCAAGGAGGCAGATCCTTCCGTAGAATTGTGCCAGCCGCAATGCCAATGGGTCCCTGTTGGCCGTAGGCCTTGCGCTGCCTAGGCTGCGCCCGCAACATCCGCTCCACATTACCTGCACAGAGTGGCCAt                                                                                                <  2:433621/150‑1 (MQ=255)

                                                        cAAGTGCAAGGTGGTGCTCAAGGAGGCAGATCCTTCCGTAGAATTGTGCCAGCCGCAATGCCAATGGGTCCCTGTTGGCCGTAGGCCTTGCGCTGCCTAGGCTGCGCCCGCAACATCCGCTCCACATTACCTGGACAGAGTGGCCATcgc                                                                                             <  1:20438/150‑1 (MQ=255)

                                                          aGTGCAAGGTGGTGCTCAAGGAGGCAGATCCTTCCGTAGAATTGTGCCAGCCGCAATGCCAATGGGTCCCTGTTGGCCGTAGGCCTTGCGCTGCCTAGGCTGCGCCCGCAACATCCGCTCCACATTACCTGGACAGAGTGGCCATCGCga                                                                                           <  1:2090624/150‑1 (MQ=255)

                                                           gTGCAAGGTGGTGCTCAAGGAGGCAGATCCTTCCGTAGAATTGTGCCAGCCGCAATGCCAATGGGTCCCTGTTGGCCGTAGGCCTTGCGCTGCCTAGGCTGCGCCCGCAACATCCGCTCCACATTACCTGGACAGAGTGGCCATCGCgag                                                                                          >  2:1571241/1‑150 (MQ=255)

                                                                aGGTGGTGCTCAAGGAGGCAGATCCTTCCGTAGAATTGTGCCAGCCGCAATGCCAATGGGTCCCTGTTGGCCGTAGGCCTTGCGCTGCCTAGGCTGCGCCCGCAACATCCGCTCCACATTACCTGGACAGAGTGGCCATCGCGAGAAAAt                                                                                     >  1:806983/1‑150 (MQ=255)

                                                                aGGTGGTGCTCAAGGAGGCAGATCCTTCCGTAGAATTGTGCCAGCCGCAATGCCAATGGGTCCCTGTTGGCCGTAGGCCTTGCGCTGCCTAGGCTGCGCCCGCAACATCCGCTCCACATTACCTGGACAGAGTGGCCATCGCGAGAAAAt                                                                                     >  1:776407/1‑150 (MQ=255)

                                                                   tggtgCTCAAGGAGGCAGATCCTTCCGTAGAATTGTGCCAGCCGCAATGCCAATGGGTCCCTGTTGGCCGTAGGCCTTGCGCTGCCTAGGCTGCGCCCGCAACATCCGCTCCACATTACCTGGACAGAGTGGCCATCGCGAGAAAATGgc                                                                                  >  1:1725835/1‑150 (MQ=255)

                                                                   tggtgCTCAAGGAGGCAGATCCTTCCGTAGAATTGTGCCAGCCGCAATGCCAATGGGTCCCTGTTGGCCGTAGGCCTTGCGCTGCCTAGGCTGCGCCCGCAACATCCGCTCCACATTACCTGGACAGAGTGGCCATCGCGAGAAAATGgc                                                                                  >  1:1726035/1‑150 (MQ=255)

                                                                    ggtgCTCAAGGAGGCAGATCCTTCCGTAGAATTGTGCCAGCCGCAATGCCAATGGGTCCCTGTTGGCCGTAGGCCTTGCGCTGCCTAGGCTGCGCCCGCAACATCCGCTCCACATTACCTGGACAGAGTGGCCATCGCGAGAAAATGgct                                                                                 <  2:886903/150‑1 (MQ=255)

                                                                     gtgCTCAAGGAGGCAGATCCTTCCGTAGAATTGTGCCAGCCGCAATGCCAATGGGTCCCTGTTGGCCGTAGGCCTTGCGCTGCCTAGGCTGCGCCCGCAACATCCGCTCCACATTACCTGGACAGAGTGGCCATCGCGAGAAAATGgctg                                                                                >  2:1417992/1‑150 (MQ=255)

                                                                      tgCTCAAGGAGGCAGATCCTTCCGTAGAATTGTGCCAGCCGCAATGCCAATGGGTCCCTGTTGGCCGTAGGCCTTGCGCTGCCTAGGCTGCGCCCGCAACATCCGCTCCACATTACCTGGACAGAGTGGCCATCGCGAGAAAATGgctgc                                                                               >  1:713117/1‑150 (MQ=255)

                                                                        cTCAAGGAGGCAGATCCTTCCGTAGAATTGTGCCAGCCGCAATGCCAATGGGTCCCTGTTGGCCGTAGGCCTTGCGCTGCCTAGGCTGCGCCCGCAACATCCGCTCCACATTACCTGGACAGAGTGGCCATCGCGAGAAAATGGCTGCTc                                                                             <  2:1725835/150‑1 (MQ=255)

                                                                         tCAAGGAGGCAGATCCTTCCGTAGAATTGTGCCAGCCGCAATGCCAATGGGTCCCTGTTGGCCGTAGGCCTTGCGCTGCCTAGGCTGCGCCCGCAACATCCGCTCCACATTACCTGGACAGAGTGGCCATCGCGAGAAAATGGCTGCTca                                                                            >  2:505535/1‑150 (MQ=255)

                                                                         tCAAGGAGGCAGATCCTTCCGTAGAATTGTGCCAGCCGCAATGCCAATGGGTCCCTGTTGGCCGTAGGCCTTGCGCTGCCTAGGCTGCGCCCGCAACATCCGCTCCACATTACCTGGACAGAGTGGCCATCGCGAGAAAATGGCTGCTca                                                                            >  2:848234/1‑150 (MQ=255)

                                                                         tCAAGGAGGCAGATCCTTCCGTAGAATTGTGCCAGCCGCAATGCCAATGGGTCCCTGTTGGCCGTAGGCCTTGCGCTGCCTAGGCTGCGCCCGCAACATCCGCTCCACATTACCTGGACAATTTGTCCATCGCGAGAAAATGGCTGCTca                                                                            <  1:1146304/150‑1 (MQ=255)

                                                                          cAAGGAGGCAGATCCTTCCGTAGAATTGTGCCAGCCGCAATGCCAATGGGTCCCTGTTGGCCGTAGGCCTTGCGCTGCCTAGGCTGCGCCCGCAACATc                                                                                                                              >  1:1601966/1‑99 (MQ=255)

                                                                          cAAGGAGGCAGATCCTTCCGTAGAATTGTGCCAGCCGCAATGCCAATGGGTCCCTGTTGGCCGTAGGCCTTGCGCTGCCTAGGCTGCGCCCGCAACATCCGCTCCACATTACCTGGACAGAGTGGCCATCGCGAGAAAATGGCTGCTca                                                                            <  2:919166/149‑1 (MQ=255)

                                                                          cAAGGAGGCAGATCCTTCCGTAGAATTGTGCCAGCCGCAATGCCAATGGGTCCCTGTTGGCCGTAGGCCTTGCGCTGCCTAGGCTGCGCCCGCAACATCCGCTCCACATTACCTGGACAGAGTGGCCATCGCGAGAAAATGGCTGCTca                                                                            >  1:919166/1‑149 (MQ=255)

                                                                          cAAGGAGGCAGATCCTTCCGTAGAATTGTGCCAGCCGCAATGCCAATGGGTCCCTGTTGGCCGTAGGCCTTGCGCTGCCTAGGCTGCGCCCGACACATc                                                                                                                              <  2:1601966/99‑1 (MQ=255)

                                                                           aaGGAGGCAGATCCTTCCGTAGAATTGTGCCAGCCGCAATGCCAATGGGTCCCTGTTGGCCGTAGGCCTTGCGCTGCCTAGGCTGCGCCCGCAACATCCGCTCCACATTACCTGGACAGAGTGGCCATCGCGAGAAAATGGCTGCTcaca                                                                          >  2:2124238/1‑150 (MQ=255)

                                                                           aaGGAGGCAGATCCTTCCGTAGAATTGTGCCAGCCGCAATGCCAATGGGTCCCTGTTGGCCGTAGGCCTTGCGCAGCCTAGGCAgcgcaagg                                                                                                                                    <  2:527979/92‑5 (MQ=255)

                                                                              gaggCAGATCCTTCCGTAGAATTGTGCCAGCCGCAATGCCAATGGGTCCCTGTTGGCCGTAGGCCTTGCGCTTCCTAGGCTGCGCCCGCAACATCCGCTCCACATTACCTGGACAGAGGGGCCATCGCGAGAAAATGgctgct                                                                              >  2:362145/1‑143 (MQ=255)

                                                                              gaggCAGATCCTTCCGTAGAATTGTGCCAGCCGCAATGCCAATGGGTCCCTGTTGGCCGTAGGCCTTGCGCTGCCTAGGCTGCGCCCGCAACATCCGCTCCACATTACCTGGACAGAGTGGCCATCGCGAGAAAATGGCTGCTCACAGaa                                                                       >  2:282494/1‑150 (MQ=255)

                                                                               aggCAGATCCTTCCGTAGAATTGTGCCAGCCGCAATGCCAATGGGTCCCTGTTGGCCGTAGGCCTTGCGCTGCCTAGGCTGCGCCCGCAACATCCGCTCCACATTACCTGGACAGAGTGGCCATCGCGAGAAAATgtcgg                                                                                >  1:1074968/1‑136 (MQ=255)

                                                                               aggCAGATCCTTCCGTAGAATTGTGCCAGCCGCAATGCCAATGGGTCCCTGTTGGCCGTAGGCCTTGCGCTGCCTAGGCTGCGCCCGCAACATCCGCTCCACATTACCTGGACAGAGTGGCCATCGCGAGAAAATGgctg                                                                                <  2:1074968/140‑1 (MQ=255)

                                                                                ggCAGATCCTTCCGTAGAATTGTGCCAGCCGCAATGCCAATGGGTCCCTGTTGGCCGTAGGCCTTGCGCTGCCTAGGCTGCGCCCGCAACATCCGCTCCACATTACCTGGACAGAGTGGCCATCGCGAGAAAATGGCTGCTCACAGAACt                                                                     <  1:2005694/150‑1 (MQ=255)

                                                                                  cAGATCCTTCCGTATAATTGTGCCAGCCGCAATGCCAATGGGTCCCTGTTGGCCGTAGGCCTTGCGCTGCCTAGGCTGTGCCCGCAACATCCGCTCCACATTACCTGGACAGAGTGGCCATCGCGAGAAAATGGCTGCTCACAGAACTgg                                                                   <  1:988477/150‑1 (MQ=255)

                                                                                  cAGATCCTTCCGTATAATTGTGCCAGCCGCAATGCCAATGGGTCCCTGTTGGCCGTAGGCCTTGCGCTGCCTAGGCTGTGCCCGCAACATCCGCTCCACATTACCTGGACAGAGTGGCCATCGCGAGAAAATGGCTGCTCACAGAACTgg                                                                   <  1:988426/150‑1 (MQ=255)

                                                                                  cAGATCCTTCCGTAGAATTGTGCCAGCCGCAATGCCAATGGGTCCCTGTTGGCCGTAGGCCTTGCGCTGCCTAGGCTGCGCCCGCAACATCCGCTCCACATTACCTGGACAGAGTGGCCATCGCGAGAAAATGGCTGCTCACAGAACTgg                                                                   <  1:1131817/150‑1 (MQ=255)

                                                                                  cAGATCCTTCCGTAGAATTGTGCCAGCCGCAATGCCAATGGGTCCCTGTTGGCCGTAGGCCTTGCGCTGCCTAGGCTGCGCCCGCAACATCCGCTCCACATTACCTGGACAGAGTGGCCATCGCGAGAAAATGGCTGCTCACAGAACTgg                                                                   <  1:2168450/150‑1 (MQ=255)

                                                                                  cAGATCCTTCCGTAGAATTGTGCCAGCCGCAATGCCAATGGGTCCCTGTTGGCCGTAGGCCTTGCGCTGCCTAGGCTGCGCCCGCAACATCCGCTCCACATTACCTGGACAGAGTGGCCATCGCGAGAAAATGGCTGCTCACAGAACTgg                                                                   <  1:1649554/150‑1 (MQ=255)

                                                                                     aTCCTTCCGTAGAATTGTGCCAGCCGCAATGCCAATTGGTCCCTGTTGGCCGTAGGCCTTGCGCAGCCTAGGCAgcgcaagg                                                                                                                                    >  2:1642364/1‑78 (MQ=38)

                                                                                     aTCCTTCCGTAGAATTGTGCCAGCCGCAATGCCAATTGGTCCCTGTTGGCCGTAGGCCTTGCGCAGCCTAGGCAgcgcaagg                                                                                                                                    <  1:1642364/82‑5 (MQ=38)

                                                                                      tCCTTCCGTAGAATTGTGCCAGCCGCAATGCCAATGGGTCCCTGTTGGCCGTAGGCCTTGCGCTGCCTAGGCTGCGCCCGCAACATCCGCTCCACATTACCTGGACAGAGTGGCCATCGCGAGAAAATGGCTGCTCACAGAACTGGCAGt                                                               >  2:232094/1‑150 (MQ=255)

                                                                                       ccTTCCGTAGAATTGTGCCAGCCGCAATGCCAATGGGTCCCTGTTGGCCGTAGGCCTTGCGCTGCCTAGGCTGCGCCCGCAACATCCGCTCCACATTACCTGGACAGAGTGGCCATCGCGAGAAAATGGCTGCTCACAGAACTGGCAGTa                                                              >  1:2344411/1‑150 (MQ=255)

                                                                                        cTTCCGTAGAATTGTGCCAGCCGCAATGCCAATGGGTCCCTGTTGGCCGTAGGCCTTGCGCTGCCTAGGCTGCGCCCGCAACATCCGCTCCACATTACCTGGACAGAGTGGCCATCGCGAGAAAATGGCTGCTCACAGAACTGGCAGTAc                                                             >  1:480776/1‑150 (MQ=255)

                                                                                          tCCGTAGAATTGTGCCAGCCGCAATGCCAATGGGTCCCTGTTGGCCGTAGGCCTTGCGCTGCCTAGGCTGCGCCCGCAACATCCGCTCCACATTACCTGGACAGAGTGGCCATCGCGAGAAAATGGCTGCTCACAGAACTGGCAGTACga                                                           <  2:776407/150‑1 (MQ=255)

                                                                                            cGTAGAATTGTGCCGGCCGCAATGCCAATGGGTCCCTGTTGGCCGTAGGCCTTGCGCTGCCTAGGCTGCGCCCGCAACATCCGCTCCACATTACCTGTACAGAGTGGCCATCGCGAGAAAATGGCTGCTCACAGAACTGGCAGTACgaga                                                         <  2:345166/150‑1 (MQ=255)

                                                                                            cGTAGAATTGTGCCAGCCGCAATGCCAATGGGTCCCTGTTGGCCGTAGGCCTTGCGCTGCCTAGGCTGCGCCCGCAACATCCGCTCCACATTACCTGGACAGAGTGGCCATCGCGAGAAAATGGCTGCTCACAGAACTGGCAGTACgaga                                                         <  2:1766831/150‑1 (MQ=255)

                                                                                            cGTAGAATTGTGCCAGCCGCAATGCCAATGGGTCCCTGTTGGCCGTAGGCCTTGCGCTGCCTAGGCTGCGCCCGCAACATCCGCTCCACATTACCTGGACAGAGTGGCCATCGCGAGAAAATGGCTGCTCACAGAACTGGCAGTACgaga                                                         >  2:57216/1‑150 (MQ=255)

                                                                                            cGTAGAATTGTGCCAGCCGCAATGCCAATGGGTCCCTGTTGGCCGTAGGCCTTGCGCTGCCTAGGCTGCGCCCGCAACATCCGCTCCACATTACCTGGACAGAGTGGCCATCGCGAGAAAATGGCTGCTCACAGAACTGGCAGTACgaga                                                         <  2:2176608/150‑1 (MQ=255)

                                                                                              tAGAATTGTGCCAGCCGCAATGCCAATGGGTCCCTGTTGGCCGTAGGCCTTGCGCTGCCTAGGCTGCGCCCGCAACATCCGCTCCACATTACCTGGACAGAGTGGCCATCGCGAGAAAATGGCTGCTCACAGAACTGGCAGTACGAGAGc                                                       <  1:1417992/150‑1 (MQ=255)

                                                                                                      tgCNAGCCGCAATGCCAATGGGTCCCTGTTGGCCGTAGGCCTTGCGCTGCCTAGGCTGCGCCCGCAACATCCGCTCCACATTACCTGGACAGAGTGGCCATCGCGAGAAAATGGCTGCTCACAGAACTGGCAGTACGAGAGCAAATCGAc                                               >  2:1500276/1‑150 (MQ=255)

                                                                                                         cAGCCGCAATGCCAATGGGTCCCTGTTGGCCGTAGGCCTTGCGCTGCCTAGGCTGCGCCCGCAACATCCGCTCCACATTACCTGGACAGAGTGGCCATCGCGAGAAAATGGCTGCTCACAGAACTGGCAGTACGAGAGCAAATCGACGtt                                            >  1:596441/1‑150 (MQ=255)

                                                                                                           gCCGCAATGCCAATGGGTCCCTGTTGGCCGTTCGCCTGGCGCTGCCTAGGCTGCGCCCGCAACATCCTAGTCACATTACCTGGACAGAGTGGCCATCGCGAGAAAATGGCTGCTCACAGAACTGGCAGTACGAGAGCAAATCGACGTtgt                                          <  1:1861802/150‑1 (MQ=255)

                                                                                                           gCCGCAATGCCAATGGGTCCCTGTTGGCCGTAGGCCTTGCGCTGCCTAGGCTGCGCCCGCAACATCCGCTCCACATTACCTGGACAGAGTGGCCATCGCGAGAAAATGGCTGCTCACAGAACTGGCAGTACGAGAGCAAATCGACGTtgt                                          <  2:713117/150‑1 (MQ=255)

                                                                                                            ccGCAATGCCAATGGGTCCCTGTTGGCCGTAGGCCTTGCGCTGCCTAGGCTGCGCCCGCAACATCCGCTCCACATTACCTGGACAGAGTGGCCATCGCga                                                                                           <  2:711390/100‑1 (MQ=255)

                                                                                                            ccGCAATGCCAATGGGTCCCTGTTGGCCGTAGGCCTTGCGCTGCCTAGGCTGCGCCCGCAACATCCGCTCCACATTACCTGGACAGAGTGGCCATCGCga                                                                                           >  1:711390/1‑100 (MQ=255)

                                                                                                            ccGCAATGCCAATGGGTCCCTGTTGGCCGTAGGCCTTGCGCTGCCTAGGCTGCGCCCGCAACATCCGCTCCACATTACCTGGACAGAGTGGCCATCGCGAGAAAATGGCTGCTCACAGAACTGGCAGTACGAGAGCAAATCGACGTtgtg                                         <  1:571909/150‑1 (MQ=255)

                                                                                                             ccgaaTGCCAATGGGTCCCTGTTGGCCGTAGGCCTTGCGCTGCCTAGGCTGCGCCCGCAACATCCGCTCCACATTACCTGGACAGAGTGGCCATCGCGAGAAAATGGCTGATCACAGAACTGGCAGTACGAGAGCAAATCGACGTTGTGc                                        >  1:82673/4‑150 (MQ=255)

                                                                                                                    ccAATGGGTCCCTGTTGGCCGTAGGCCTTGCGCTGCCTAGGCTGCGCCCGCAACATCCGCTCCACATTACCTGGACAGAGTGGCCATCGCGAGAAAATGGCTGCTCACAGAACTGGCAGTACGAGAGCAAATCGACGTTGTGCTTCgcgc                                 >  2:2257620/1‑150 (MQ=255)

                                                                                                                     cAATGGGTCCCTGTTGGCCGTAGGCCTTGCGCTGCCTAGGCTGCGCCCGCAACATCCGCTCCACATTACCTGGACAGAGTGGCCATCGCGAGAAAATGGCTGCTCACAGAACTGGCAGTACGAGAGCAAATCGACGTTGTGCTTCGCGca                                >  1:2069996/1‑150 (MQ=255)

                                                                                                                      aaTGGGTCCCTGTTGGCCGTAGGCCTTGCGCTGCCTAGGCTGCGCCCGCAACATCCGCTCCACATTACCTGGACAGAGTGGCCATCGCGAGAAAATGGCTGCTCACAGAACTgg                                                                   >  1:2051343/1‑114 (MQ=255)

                                                                                                                      aaTGGGTCCCTGTTGGCCGTAGGCCTTGCGCTGCCTAGGCTGCGCCCGCAACATCCGCTCCACATTACCTGGACAGAGTGGCCATCGCGAGAAAATGGCTGCTCACAGAACTgg                                                                   <  2:2051343/114‑1 (MQ=255)

                                                                                                                         gggTCCCTGTTGGCCGTAGGCCTTGCGCTGCCTAGGCTGCGCCCGCAACATCCGCTCCACATTACCTGGACAGAGTGGCCATCGCGAGAAAATGGCTGCTCACAGAACTGGCAGTACGAGAGCAAATCGACGTTGTGCTTCGCGcacaca                            <  1:282494/150‑1 (MQ=255)

                                                                                                                             cccTGTTGGCCGTAGGCCTTGCGCTGCCTAGGCTGCGCCCGCAACATCCGCTCCACATTACCTGGACAGAGTGGCCATCGCGAGAAAATGGCTGCTCACAGAACTGGCAGTACGAGAGCAAATCGACGTTGTGCTTCGCGCACACAGTTa                        >  2:2162182/1‑150 (MQ=255)

                                                                                                                              ccTGTTGGCCGTAGGCCTTGCGCTGCCTAGGCTGCGCCCGCAACATCCGCTCCACATTACCTGGACAGAGTGGCCATCGCGAGAAAATGGCTGCTCACAGAACTGGCAGTACGAGAGCAAATCGACGTTGTGCTTCGCGCACACAGTTAc                       <  2:1339175/150‑1 (MQ=255)

                                                                                                                                 gTTGGCCGTAGGCCTTGCGCTGCCTAGGCTGCGCCCGCAACATCCGCTCCACATTACCTGGACAGAGTGGCCATCGCGAGAAAATGGCTGCTCACAGAACTGGCAGTACGAGAGCAAATCGACGTTGTGCTTCGCGCACACAGTTACAtg                    <  1:16810/150‑1 (MQ=255)

                                                                                                                                  ttGGCGGTTGTCCTTGCGCTGCCTAGGCTGCGCCCGCAACATCCGCTCCACATTACCTGGACAGAGTGGCCATCGCGAGAAAATGGCTGCTCACAGAACTGGCAGTACGAGAGCAAATCGACGTTGTGCTTCGCGCACACAGTTACAtgt                   >  2:1912967/1‑150 (MQ=255)

                                                                                                                                    ggCCGTAGGCCTTGCGCTGCCTAGGCTGCGCCCGCAACATCCGCTCCACATTACCTGGACAGAGTGGCCATCGCGAGAAAATGGCTGCTCACAGAACTGGCAGTACGAGAGCAAATCGACGTTGTGCTTCGCGCACACAGTTACATGTgg                 <  1:1123683/150‑1 (MQ=255)

                                                                                                                                    ggCCGTAGGCCTTGCGCTGCCTAGGCTGCGCCCGCAACATCCGCTCCACATTACCTGGACAGAGTGGCCATCGCGAGAAAATGGCTGCTCACAGAACTGGCAGTACGAGAGCAAATCGACGTTGTGCTTCGCGCACACAGTTACATGTgg                 <  1:1236044/150‑1 (MQ=255)

                                                                                                                                    ggCCGTAGGCCTTGCGCTGCCTAGGCTGCGCCCGCAACATCCGCGCCACATTACCTGGACAGAGTGGCCATCGCGAGAAAATGGCTGCTCACAGAACTGGCAGTACGAGAGCAAATCGACGTTGTGCTTCGCGCACACAGTTACATGTgg                 <  1:509465/150‑1 (MQ=255)

                                                                                                                                        gTAGGCCTTGCGCTGCCTAGGCTGCGCCCGCAACATACGCTCAACGTTGTGTGGACAGAGTGGCCATCGCGAGAAAATGGCTGCTCACAGAACTGGCAGTACGAGAGCAAATCGACGTTGTGCTTCGCGCACACAGTTACATGTGGGTGa             >  1:1616369/1‑150 (MQ=255)

                                                                                                                                         tAGGCCTTGCGCTGCCTAGGCTGCGCCCGCAACATCCGCTCCACATTACCTGGACAGAGTGGCCATCGCGAGAAAATGGCTGCTCACAGAACTGGCAGTACGAGAGCAAATCGACGTTGTGCTTCGCGCACACAGTTACATGTGGGTGaa            <  2:297779/150‑1 (MQ=255)

                                                                                                                                            gCCTTGCGCTGCCTAGGCTGCGCCCGCAACATCCGCTCCACATTACCTGGACAGAGTGGCCATCGCGAGAAAATGGCTGCTCACAGAACTGGCAGTACGAGAGCAAATCGACGTTGTGCTTCGCGCACACAGTTACATGTgg                 <  1:344778/142‑1 (MQ=255)

                                                                                                                                            gCCTTGCGCTGCCTAGGCTGCGCCCGCAACATCCGCTCCACATTACCTGGACAGAGTGGCCATCGCGAGAAAATGGCTGCTCACAGAACTGGCAGTACGAGAGCAAATCGACGTTGTGCTTCGCGCACACAGTTACATGTgg                 >  2:344778/1‑142 (MQ=255)

                                                                                                                                            gCCTTGCGCTGCCTAGGCTGCGCCCGCAACATCCGCTCCACATTACCTGGACAGAGTGGCCATCGCGAGAAAATGGCTGCTCACAGAACTGGCAGTACGAGAGCAAATCGACGTTGTGCTTCGCGCACACAGTTACATGTGGGTGAAGAc         <  2:1679213/150‑1 (MQ=255)

                                                                                                                                                 gcgcTGCCTAGGCTGCGCCCGCAACATCCGCTCCACATTACCTGGACAGAGTGGCCATCGCGAGAAAATGGCTGCTCACAGAACTGGCAGTACGAGAGCAAATCGACGTTGTGCTTCGCGCACACAGTTACATGTgg                 >  2:1733648/1‑137 (MQ=255)

                                                                                                                                                 gcgcTGCCTAGGCTGCGCCCGCAACATCCGCTCCACATTACCTGGACAGAGTGGCCATCGCGAGAAAATGGCTGCTCACAGAACTGGCAGTACGAGAGCAAATCGACGTTGTGCTTCGCGCACACAGTTACATGTGGGTGAAGACGTCGc    <  2:200867/150‑1 (MQ=255)

AACAATTTCATTGTCTATGTTTTTTAGTTTGCTTTCGAGCACCCGTAGAACCGTTTCAAGTGCAAGGTGGTGCTCAAGGAGGCAGATCCTTCCGTAGAATTGTGCCAGCCGCAATGCCAATGGGTCCCTGTTGGCCGTAGGCCTTGCGCTGCCTAGGCTGCGCCCGCAACATCCGCTCCACATTACCTGGACAGAGTGGCCATCGCGAGAAAATGGCTGCTCACAGAACTGGCAGTACGAGAGCAAATCGACGTTGTGCTTCGCGCACACAGTTACATGTGGGTGAAGACGTCGCCA  >  NC_002947/6169688‑6169984

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