BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,005,756	G→T	5.3%	P78T (CCG→ACG)	PP_2624 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,005,756	0	G	T	5.3%	1244.2 / 41.0	361	P78T (CCG→ACG)	PP_2624	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (139/203); new base T (10/9); total (149/212)
Fisher's exact test for biased strand distribution p-value = 3.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.33e-01

Only 100 of 363 total aligned reads displayed.

TTTCAGGGCGACCAGTGCTTCGCGCAGCGCCAGCAGCTTGTTCAGCTCGGGCACCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTCGCCCAGCACCAGCATCTTGTGCGGCAGCTCCACC  >  NC_002947/3005607‑3005904

tttCAGGGCGACCAGTGCTTCGCGCAGCGCCAGCAGCTTGTTCAGCTCGGGCACCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACgg                                                                                                                                                      >  2:5963455/1‑150 (MQ=255)

           ccAGTGCTTCGCGCAGCGCCAGCAGCTTGTTCAGCTCGGGCACCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGt                                                                                                                                           >  1:1812505/1‑150 (MQ=255)

               tGCTTCGCGCAGCGCCAGCAGCTTGTTCAGCTCGGGCACCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGTGACGGTACCGTCCAg                                                                                                                                       <  1:920898/150‑1 (MQ=37)

               tGCTTCGCGCAGCGCCAGCAGCTTGTTCAGCTCGGGCACCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGTGACGGTACCGTCCAg                                                                                                                                       <  1:2106693/150‑1 (MQ=37)

                gCTTCGCGCAGCGCCAGCAGCTTGTTCAGCTCGGGCACCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAgg                                                                                                                                      <  1:4873616/150‑1 (MQ=255)

                 cTTCGCGCAGCGCCAGCAGCTTGTTCAGCTCGGGCACCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGTGACGGTACCGTCCAGGc                                                                                                                                     <  1:2936467/150‑1 (MQ=38)

                  ttCGCGCAGCGCCAGCAGCTTGTTCAGCTCGGGCACCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCt                                                                                                                                    >  2:4181498/1‑150 (MQ=255)

                  ttCGCGCAGCGCCAGCAGCTTGTTCAGCTCGGGCACCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCt                                                                                                                                    >  1:4038610/1‑150 (MQ=255)

                  ttCGCGCAGCGCCAGCAGCTTGTTCAGCTCGGGCACCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCt                                                                                                                                    >  1:3331277/1‑150 (MQ=255)

                   tCGCGCAGCGCCAGCAGCTTGTTCAGCTCGGGCACCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTc                                                                                                                                   <  1:725377/150‑1 (MQ=255)

                       gcAGCGCCAGCAGCTTGTTCAGCTCGGGCACCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACAt                                                                                                                               >  2:7014494/1‑150 (MQ=255)

                          ccgcCAGCAGCTGGGTCAGCTCGGCCACCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTGGATGGACCTGACCTGCAGGGGGACGGCCAGCTCGGGATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTgg                                                                                                                            <  2:2287834/149‑1 (MQ=255)

                             ccAGCAGCTTGTTCAGCTCGGGCACCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGcgt                                                                                                                         <  2:3161313/150‑1 (MQ=255)

                              cagctgCTTGTTCAGCTCGGGCTCCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCCCGTCCAGGCTCACATTGGcgtc                                                                                                                        <  1:521459/150‑1 (MQ=255)

                              cagcagCTTGTTCAGCTCGGGCACCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGcgtc                                                                                                                        <  2:816679/150‑1 (MQ=255)

                              cagcagCTTGTTCAGCTCGGGCACCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGcgtc                                                                                                                        <  2:1347471/150‑1 (MQ=255)

                                gcagcTTGTTCAGCTCGGGCACCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGcgtcgt                                                                                                                      <  1:4448087/150‑1 (MQ=255)

                                 cagcTTGTTCAGCTCGGGCACCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGtt                                                                                                                     <  1:6055735/150‑1 (MQ=255)

                                 cagcTTGTTCAGCTCGGGCACCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGtt                                                                                                                     <  1:4708822/150‑1 (MQ=255)

                                  agcTTGTTCAGCTCGGGCACCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTc                                                                                                                    <  1:1046322/150‑1 (MQ=255)

                                   gcTTGTTCAGCTCGGGCACCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCa                                                                                                                   <  1:7440146/150‑1 (MQ=255)

                                    cTTGTTCAGCTCGGGCACCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCaa                                                                                                                  >  2:1614416/1‑150 (MQ=255)

                                     ttGTTCAGCTCGGGCACCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAac                                                                                                                 <  1:222811/150‑1 (MQ=255)

                                       gTTCAGCTCGGGCACCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAacac                                                                                                               >  1:1705116/1‑150 (MQ=255)

                                           aGCTCGGGCACCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTg                                                                                                           >  2:6663077/1‑150 (MQ=255)

                                            gCTCGGGCACCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGt                                                                                                          <  1:4181498/150‑1 (MQ=255)

                                             cTCGGGCACCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGtt                                                                                                         >  2:1144630/1‑150 (MQ=255)

                                             cTCGGGCACCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGtt                                                                                                         <  1:2165029/150‑1 (MQ=255)

                                             cTCGGGCACCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGtt                                                                                                         <  2:1961142/150‑1 (MQ=255)

                                              tCGGGCACCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTg                                                                                                        <  1:2558236/150‑1 (MQ=255)

                                              tCGGGCACCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGGCGTTCAACACGCTGTTg                                                                                                        >  2:5294139/1‑150 (MQ=255)

                                                  gCACCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAatg                                                                                                    >  2:1459051/1‑150 (MQ=255)

                                                   cACCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACAGCGTGTTGaacga                                                                                                   <  1:5334483/150‑4 (MQ=255)

                                                    aCCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAatgta                                                                                                  <  1:5471413/150‑1 (MQ=255)

                                                     ccTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAatgtat                                                                                                 <  1:4627080/150‑1 (MQ=255)

                                                     ccTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAatgtat                                                                                                 <  2:2543785/150‑1 (MQ=255)

                                                     ccTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAatgtat                                                                                                 <  1:1167875/150‑1 (MQ=255)

                                                      cTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAatgtctg                                                                                                >  2:6059772/1‑150 (MQ=255)

                                                           gcgcgATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTAt                                                                                           <  1:553918/150‑1 (MQ=255)

                                                            cgcgATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATc                                                                                          <  2:4331541/150‑1 (MQ=255)

                                                               gATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGAt                                                                                       <  1:1448805/150‑1 (MQ=255)

                                                               gATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGAt                                                                                       >  1:4544351/1‑150 (MQ=255)

                                                               gATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGAt                                                                                       >  2:7215139/1‑150 (MQ=255)

                                                               gATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGAt                                                                                       <  2:4844511/150‑1 (MQ=255)

                                                                aTCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATg                                                                                      <  1:343303/150‑1 (MQ=255)

                                                                   cggtccggCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATgcgc                                                                                   <  1:3718847/150‑1 (MQ=255)

                                                                   cggtccggCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATgcgc                                                                                   <  2:443384/150‑1 (MQ=255)

                                                                   cggtccggCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATgcgc                                                                                   >  1:97866/1‑150 (MQ=255)

                                                                   cggtccggCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATgcgc                                                                                   <  1:2012170/150‑1 (MQ=255)

                                                                    ggtccggCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCa                                                                                  >  1:131846/1‑150 (MQ=255)

                                                                       ccggCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGa                                                                               <  2:1868961/150‑1 (MQ=255)

                                                                        cggCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGAc                                                                              <  2:2168514/150‑1 (MQ=255)

                                                                         ggCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGAcc                                                                             <  2:469753/150‑1 (MQ=255)

                                                                          gCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCt                                                                            >  1:7365749/1‑150 (MQ=255)

                                                                                 aTCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTc                                                                     <  1:2979582/150‑1 (MQ=255)

                                                                                 aTCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTc                                                                     <  2:3700794/150‑1 (MQ=255)

                                                                                 aTCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTc                                                                     <  2:5570928/150‑1 (MQ=255)

                                                                                 aTCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTc                                                                     <  1:2332163/150‑1 (MQ=255)

                                                                                   cGTTGATGGACCTGAACTGCAGGGGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCt                                                                   <  2:2898438/150‑1 (MQ=255)

                                                                                    gTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTc                                                                  >  1:7030155/1‑150 (MQ=255)

                                                                                     ttGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCa                                                                 >  1:5807778/1‑150 (MQ=255)

                                                                                           gACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCa                                                           <  2:6278504/150‑1 (MQ=255)

                                                                                             ccTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCaa                                                                <  1:4294231/143‑1 (MQ=255)

                                                                                               tGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGc                                                       >  2:3236356/1‑150 (MQ=255)

                                                                                               tGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGc                                                       >  2:2007437/1‑150 (MQ=255)

                                                                                                gAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCt                                                      <  1:359539/150‑1 (MQ=255)

                                                                                                   cTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTc                                                   <  1:6578161/150‑1 (MQ=255)

                                                                                                    tGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCa                                                  <  2:7051359/150‑1 (MQ=255)

                                                                                                      cAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATc                                                <  1:7110517/150‑1 (MQ=255)

                                                                                                       aGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATcc                                               <  2:7300179/150‑1 (MQ=255)

                                                                                                         gTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAg                                             >  2:3402961/1‑150 (MQ=255)

                                                                                                           ggACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGc                                           <  2:5121765/150‑1 (MQ=255)

                                                                                                            gACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGcc                                          >  2:5915844/1‑150 (MQ=255)

                                                                                                               ggCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACCTTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGaa                                       <  2:551151/150‑1 (MQ=255)

                                                                                                               ggCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGaa                                       <  1:1519536/150‑1 (MQ=255)

                                                                                                               ggCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGaa                                       <  1:3100672/150‑1 (MQ=255)

                                                                                                               ggCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGaa                                       <  1:6243303/150‑1 (MQ=255)

                                                                                                               ggCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGaa                                       <  2:4784585/150‑1 (MQ=255)

                                                                                                                gCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAg                                      <  1:161647/150‑1 (MQ=255)

                                                                                                                    gCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGGTGAATGTATGTTTATAGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTCGc                                  >  1:4997588/1‑150 (MQ=255)

                                                                                                                       cGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTCGCcca                               >  2:2472977/1‑150 (MQ=255)

                                                                                                                            tCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTCGCccagcacc                          >  1:4079894/1‑150 (MQ=255)

                                                                                                                               ggggCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTCGCccagcaccagc                       >  1:6964110/1‑150 (MQ=255)

                                                                                                                                gggCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGa                                                                               >  1:3880505/1‑93 (MQ=255)

                                                                                                                                gggCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGa                                                                               <  2:3880505/93‑1 (MQ=255)

                                                                                                                                  gCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTCGCCCAGCACCAGCATc                    <  1:3724976/150‑1 (MQ=255)

                                                                                                                                   cGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAatgtat                                                                                                 <  1:114333/72‑1 (MQ=255)

                                                                                                                                   cGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAatgtat                                                                                                 >  2:114333/1‑72 (MQ=255)

                                                                                                                                    ggTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTCGCCCAGCACCAGCATCtt                  <  2:6969366/150‑1 (MQ=255)

                                                                                                                                     gTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGATTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTCGCCCAGCACCAGCATCTtg                 <  1:6459743/150‑1 (MQ=255)

                                                                                                                                           aGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTCGCCCAGCACCAGCATCTTGTGCGGc           <  1:5906309/150‑1 (MQ=255)

                                                                                                                                            gCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTCGCCCAGCACCAGCATCTTGTGCGGCa          <  1:6734153/150‑1 (MQ=255)

                                                                                                                                             ccTGGACGTGACGGTACCGTCCAGGCTCGCATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAgc                                                              <  1:2671653/97‑1 (MQ=21)

                                                                                                                                             ccTGGACGTGACGGTACCGTCCAGGCTCGCATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAgc                                                              >  2:2671653/1‑97 (MQ=21)

                                                                                                                                             ccTGGACGTGACGGTACCGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTCGCCCAGCACCAGCATCTTGTGCGGCAg         >  2:5470548/1‑150 (MQ=18)

                                                                                                                                             ccTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTCGCCCAGCACCAGCATCTTGTGCGGCAg         >  2:3809070/1‑150 (MQ=255)

                                                                                                                                              cTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTCGCCCAGCACCAGCATCTTGTGCGGCAGc        >  1:3956074/1‑150 (MQ=255)

                                                                                                                                               tGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTCGCCCAGCACCAGCATCTTGTGCGGCAGCt       >  2:6179691/1‑150 (MQ=255)

                                                                                                                                               tGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTCGCCCAGCACCAGCATCTTGTGCGGCAGCt       >  2:3498213/1‑150 (MQ=255)

                                                                                                                                                   cGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTCGCCCAGCACCAGCATCTTGTGCGGCAGCTCCAc   <  1:5963455/150‑1 (MQ=255)

TTTCAGGGCGACCAGTGCTTCGCGCAGCGCCAGCAGCTTGTTCAGCTCGGGCACCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTCGCCCAGCACCAGCATCTTGTGCGGCAGCTCCACC  >  NC_002947/3005607‑3005904

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