BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,623,232	T→A	6.6%	D56V (GAC→GTC)	PP_1422 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,623,232	0	T	A	6.6%	460.9 / 16.6	151	D56V (GAC→GTC)	PP_1422	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (69/71); new base A (5/5); total (74/77)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.07e-01

Only 100 of 151 total aligned reads displayed.

TCTGCTCGATCACCAGGCCGATGGCCAGGTTGCAGCTGTAGTTGATGCCTAGCCGCGTGATGGCGGTATGCAGGTCGGTCACCTCGACCGCAGCACGCAGCAATGGGCTGTTGACGACCTTGATCAGGCGCGCAGACAGGGCCGCGTCGCGGCCGATCACCTTGCTCAGGGCCGCTACGCTGATCTCGCTGTCTTC‑CGCTGCCTCGCGGATGCTCAAGGCAACCTCGGGCAGGGTTGGCAGGACCAGGTCATCCTTTTCGATGGCGTCGAGCAATTGTGCTTGAACCAACTCGGC  >  NC_002947/1623086‑1623380

tCTGCTCCATCACCAGGCCGATGGCCAGGTTGCAGCTGTAGTTGATGCCTAGCCGCGTGATGGCGGTATGCAGGTCGGTCACCTCGACCGCAGCACGCAGCAATGGGCTGTTGACGACCTTGATCAGGCGCGCAGACAGGGCCGCGTcgc                                                                                                                                                    >  1:1842807/1‑150 (MQ=255)

   gCTCGATCACCAGGCCGATGGCCAGGTTGCAGCTGTAGTTGATGCCTAGCCGCGTGATGGCGGTATGCAGGTCGGTCACCTCGACCGCAGCACGCAGCAATGGGCTGTTGACGACCTTGATCAGGCGCGCAGACAGGGCCGCGTCGCGGc                                                                                                                                                 >  2:239526/1‑150 (MQ=255)

   gCTCGATCACCAGGCCGATGGCCAGGTTGCAGCTGTAGTTGATGCCTAGCCGCGTGATGGCGGTATGCAGGTCGGTCACCTCGACCGCAGCACGCAGCAATGGGCTGTTGACGACCTTGATCAGGCGCGCAGACAGGGCCGCGTCGCGGc                                                                                                                                                 <  2:2399905/150‑1 (MQ=255)

      cGATCACCAGGCCGATGGCCAGGTTGCAGCTGTAGTTGATGCCTAGCCGCGTGATGGCGGTATGCAGGTCGGTCACCTCGACCGCAGCACGCAGCAATGGGCTGTTGACGACCTTGATCAGGCGCGCAGACAGGGCCGCGACGCGGccct                                                                                                                                              <  1:1018814/150‑3 (MQ=255)

       gATCACCAGGCCGATGGCCAGGTTGCAGCTGTAGTTGATGCCTAGCCGCGTGATGGCGGTATGCAGGTCGGTCACCTCGACCGCAGCACGCAGCAATGGGCTGTTGACGACCTTGATCAGGCGCGCAGACAGGGCCGCGTCGCGGCCGAt                                                                                                                                             >  1:660549/1‑150 (MQ=255)

       gATCACCAGGCCGATGGCCAGGTTGCAGCTGTAGTTGATGCCTAGCCGCGTGATGGCGGTATGCAGGTCGGTCACCTCGACCGCAGCACGCAGCAATGGGCTGTTGACGACCTTGATCAGGCGCGCAGACAGGGCCGCGACGCGGccct                                                                                                                                              <  1:1557840/149‑3 (MQ=255)

              aGGCCGATGGCCAGGTTGCAGCTGTAGTTGATGCCTAGCCGCGTGATGGCGGTATGCAGGTCGGTCACCTCGACCGCAGCACGCAGCAATGGGCTGTTGACGACCTTGATCAGGCGCGCAGACAGGGCCGCGTCGCGGCCGATCACCTTg                                                                                                                                      <  1:1614047/150‑1 (MQ=255)

               ggCCGATGGCCAGGTTGCAGCTGTAGTTGATGCCTAGCCGCGTGATGGCGGTATGCAGGTCGGTCACCTCGACCGCAGCACGCAGCAATGGGCTGTTGACGACCTTGATCAGGCGCGCAGACAGGGCCGCGTCGCGGCCGATCACCTTGc                                                                                                                                     <  1:1447992/150‑1 (MQ=255)

                  cGATGGCCAGGTTGCAGCTGTAGTTGATGCCTAGCCGCGTGCTGGCGCTATGCAGGTCGGTAACCTCGACCGTAGCACGCAGCAATGGGCTGTTGACGAGATTGATCAGGCGCGCAGACAGGGCCGCGTGGCGGCCGATCACGTTGCTCa                                                                                                                                  >  1:1000334/1‑150 (MQ=255)

                  cGATGGCCAGGTTGCAGCTGTAGTTGATGCCTAGCCGCGTGATGGCGGTATGCAGGTCGGTCACCTCGACCGCAGCACGCAGCAATGGGCTGTTGACGAGCTTGATCAGGCGCGCAGACAGGGCCGCGTCGCGGCCGATCACCTTGCTCa                                                                                                                                  >  2:1270530/1‑150 (MQ=255)

                   gATGGCCAGGTTGCAGCTGTAGTTGATGCCTAGCCGCGTGATGGCGGTATGCAGGTCGGTCACCTCGACCGCAGCACGCAGCAATGGGCTGTTGACGACCTTGATCAGGCGCGCAGACAGGGCCGCGTCGCGGCCGATCACCTTGCTCAg                                                                                                                                 >  2:1984554/1‑150 (MQ=255)

                    aTGGCCAGGTTGCAGCTGTAGTTGATGCCTAGCCGCGTGATGGCGGTATGCAGGTCGGTCACCTCGACCGCAGCACGCAGCAATGGGCTGTTGACGACCTTGATCAGGCGCGCAGACAGGGCCGCGTCGCGGCCGATCACCTTGCTCAgg                                                                                                                                <  2:616053/150‑1 (MQ=255)

                         cAGGTTGCAGCTGTAGTTGATGCCTAGCCGCGTGATGGCGGTATGCAGGTCGGTCACCTCGACCGCAGCACGCAGCAATGGGCTGTTGACGACCTTGATCAGGCGCGCAGACAGGGCCGCGTCGCGGCCGATCACCTTGCTCAGGGCCGc                                                                                                                           >  2:281566/1‑150 (MQ=255)

                              tGCAGCTGTAGTTGATGCCTAGCCGCGTGATGGCGGTATGCAGGTCGGTCACCTCGACCGCAGCACGCAGCAATGGGCTGTTGACGACCTTGATCAGGCGCGCAGACAGGGCCGCGTCGCGGCCGATCACCTTGCTCAGGGCCGCTACGc                                                                                                                      <  1:1270530/150‑1 (MQ=255)

                                      tAGTTGATGCCTAGCCGCGTGATGGCGGTATGCAGGTCGGTCACCTCGACCGCAGCACGCAGCAATGGGCTGTTGACGACCTTGATCAGGCGCGCAGACAGGGCCGCGTCGCGGCCGATCACCTTGCTCa                                                                                                                                  >  2:449999/1‑130 (MQ=255)

                                      tAGTTGATGCCTAGCCGCGTGATGGCGGTATGCAGGTCGGTCACCTCGACCGCAGCACGCAGCAATGGGCTGTTGACGACCTTGATCAGGCGCGCAGACAGGGCCGCGTCGCGGCCGATCACCTTGCTCAGGGcc                                                                                                                             <  1:879373/135‑1 (MQ=255)

                                              gCCTAGCCGCGTGATGGCGGTATGCAGGTCGGTCACCTCGACCGCAGCACGCAGCAATGGGCTGTTGACGACCTTGATCAGGCGCGCAGACAGGGCCGCGTCGCGGCCGATCACCTTGCTCAGGGCCGCTACGCTGATCTCGCTGTCTTc                                                                                                      >  2:293636/1‑150 (MQ=255)

                                              gCCTAGCCGCGTGATGGCGGTATGCAGGTCGGTCACCTCGACCGCAGCACGCAGCAATGGGCTGTTGACGACCTTGATCAGGCGCGCAGACAGGGCCGCGTCGCGGCCGATCACCTTGCTCAGGGCCGCTACGCTGATCTCGCTGTCTTc                                                                                                      <  2:857065/150‑1 (MQ=255)

                                                cTAGCCGCGTGATGGCGGTATGCAGGTCGGTCACCTCGACCGCAGCACGCAGCAATGGGCTGTTGACGACCTTGATCAGGCGCGCAGACAGGGCCGCGTCGCGGCCGATCACCTTGCTCAGGGCCGCTACGCTGATCTCGCTGTCTTC‑Cg                                                                                                   >  2:446195/1‑150 (MQ=255)

                                                cTAGCCGCGTGATGGCGGTATGCAGGTCGGTCACCTCGACCGCAGCACGCAGCAATGGGCTGTTGACGACCTTGATCAGGCGCGCAGACAGGGCCGCGTCGCGGCCGATCACCTTGCTCAGGGCCGCTACGCTGATCTCGCTGTCTTC‑Cg                                                                                                   >  2:347667/1‑150 (MQ=255)

                                                    ccGCGTGATGGCGGTATGCAGGTCGGTCACCTCGACCGCAGCACGCAGCAATGGGCTGTTGACGACCTTGATCAGGCGCGCAGACAGGGCCGCGTCGCGGCCGATCACCTTGCTCAGGGCCGCTACGCTGATCTCGCTGTCTTC‑CGCTGc                                                                                               <  2:233031/150‑1 (MQ=255)

                                                     cgcgTGATGGCGGTATGCAGGTCGGTCACCTCGACCGCAGCACGCAGCAATGGGCTGTTGACGACCTTGATCAGGCGCGCAGACAGGGCCGCGTCGCGGCCGATCACCTTGCTCAGGGCCGCTACGCTGATCTCGCt                                                                                                            >  2:251272/1‑137 (MQ=255)

                                                     cgcgTGATGGCGGTATGCAGGTCGGTCACCTCGACCGCAGCACGCAGCAATGGGCTGTTGACGACCTTGATCAGGCGCGCAGACAGGGCCGCGTCGCGGCCGATCACCTTGCTCAGGGCCGCTACGCTGATCTCGCt                                                                                                            <  1:251272/137‑1 (MQ=255)

                                                     cgcgTGATGGCGGTATGCAGGTCGGTCACCTCGACCGCAGCACGCAGCAATGGGCTGTTGACGACCTTGATCAGGCGCGCAGACAAGCCCGCGTCGCGGCCGATCACCTTGCTCAGGGCCGCTACGCTGCTCTCCATGTCTTC‑CGCTGcc                                                                                              >  2:473387/1‑150 (MQ=255)

                                                      gcgTGATGGCGGTATGCAGGTCGGTCACCTCGACCGCAGCACGCAGCAATGGGCTGTTGACGACCTTGATCAGGCGCGCAGACAGGGCCGCGTCGCGGCCGATCACCTTGCTCATGGCCGCTACGCTGATCTCGCTGTCTTC‑CGCTGCCt                                                                                             >  1:639240/1‑150 (MQ=255)

                                                        gTGATGGCGGTATGCAGGTCGGTCACCTCGACCGCAGCACACAGCAATGGGCTGGTGACGACCTTGATCAGGCGCTCAGACAGGGCCGCGTCGCGGCCGATCACCTTGCTCAGGGCCGCTACGCTGATCTCGCTGTCTTC‑CGCTGCCTcg                                                                                           >  2:2081446/1‑150 (MQ=255)

                                                          gATGGCGGTATGCAGGTCGGTCACCTCGACCGCAGCACCAAGCACTGGGCTGGTTACTACCTGGATCAGGCGCGCAGACAGGGCCGCGTCGCGGCCGATCACCTTGCTCAGGGCCGCTACGCTGATCTCGCTGTCTTC‑CGCTGCCTcgcg                                                                                         >  1:658597/1‑150 (MQ=255)

                                                           aTGGCGGTATGCAGGTCGGTCACCTCGACCGCAGCACGCAGCAATGGGCTGTTGACGACCTTGATCAGGCGCGCAGACAGGGCCGCGTCGCGGCCGATCACCTTGCTCAGGACCGCTACGCTGATCTCGCTGTCTTC‑CGCTGCCTCGCgg                                                                                        >  1:213968/1‑150 (MQ=255)

                                                               cGGTATGCAGGTCGGTCACCTCGACCGCAGCACGCAGCAATGGGCTGTTGACGACCTTGATCAGGCGCGCAGACAGGGCCGCGTCGCGGCCGATCACCTTGCTCAGGGCCGCTACGCTGATCTCGCTGTCTTC‑CGCTGCCTCGCGGATGc                                                                                    <  1:1492444/150‑1 (MQ=255)

                                                               cGGTATGCAGGTCGGTCACCTCGACCGCAGCACGCAGCAATGGGCTGTTGACGACCTTGATCAGGCGCGCAGACAGGGCCGCGTCGCGGCCGATCACCTTGCTCAGGGCCGCTACGCTGATCTCGCTGTCTTC‑CGCTGCCTCGCGGATGc                                                                                    >  1:1176664/1‑150 (MQ=255)

                                                                  tATGCAGGTCGGTCACCTCGACCGCAGCACGCAGCAATGGGCTGTTGACGACCTTGATCAGGCGCGCAGACAGGGCCGCGTCGCGGCCGATCACCTTGCTCAGGGCCGCTACGCTGATCTCGCTGTCTTC‑CGCTGCCTCGCGGATGCTCa                                                                                 >  2:1603264/1‑150 (MQ=255)

                                                                   aTGCAGGTCGGTCACCTCGACCGCAGCACGCAGCAATGGGCTGTTGACGACCTTGATCAGGCGCGCAGACAGGGCCGCGTCGCGGCCGATCACCTTGCTCAGGGCCGCTACGCTGATCTCGCTGTCTTC‑CGCTGCCTCGCGGATGCTCaa                                                                                <  1:473387/150‑1 (MQ=255)

                                                                    tGCAGGTCGGTCACCTCGACCGCAGCACGCAGCAATGGGCTGTTGACGACCTTGATCAGGCGCACAGACAGGGCCGCGTCGCGGCCGATCACCTTGCTCAGGGCCGCTACGCTGATCTCGCTGTCTTC‑CGCTGCCTCGCGGATGCTCAAg                                                                               <  1:1228262/150‑1 (MQ=255)

                                                                     gCAGGTCGGTCACCTCGACCGCAGCACGCAGCAATGGGCTGTTGACGACCTTGATCAGGCGCGCAGACAGGGCCGCGTCGCGGCCGATCACCTTGCTCAGGGCCGCTACGCTGATCTCGCTGTCTTc                                                                                                      >  1:939046/1‑127 (MQ=255)

                                                                     gCAGGTCGGTCACCTCGACCGCAGCACGCAGCAATGGGCTGTTGACGACCTTGATCAGGCGCGCAGACAGGGCCGCGTCGCGGCCGATCACCTTGCTCAGGGCCGCTACGCTGATCTCGCTGTCTTc                                                                                                      <  2:939046/127‑1 (MQ=255)

                                                                      cAGGTCGGTCACCTCGACCGCAGCACGCAGCAATGGGCTGTTGACGACCTTGATCAGGCGCGCAGACAGGGCCGCGTCGCGGCCGATCACCTTGCTCAGGGCCGCTACGCTGATCTCGCTGTCTTC‑CGCTGCCTcgc                                                                                          >  1:2113619/1‑137 (MQ=255)

                                                                       aGGTCGGTCACCTCGACCGCAGCACGCAGCAATGGGCTGTTGACGACCTTGATCAGGCGCGCAGACAGGGCCGCGTCGCGGCCGATCACCTTGc                                                                                                                                     >  2:2212910/1‑94 (MQ=255)

                                                                            ggtcACCTCGACCGCAGCACGCAGCAATGGGCTGTTGACGACCTTGATCAGGCGCGCAGACAGGGCCGCGTCGCGGCCGATCACCTTGCTCAGGGCCGCTACGCTGATCTCGCTGTCTTC‑CGCTGCCTCGCGGATGCTCAAGGCAACCTc                                                                       >  1:72467/1‑150 (MQ=255)

                                                                            ggtcACCTCGACCGCAGCACGCAGCAATGGGCTGTTGACGACCTTGATCAGGCGCGCAGACAGGGCCGCGTCGCGGCCGATCACCTTGCTCAGGGCCGCTACGCTGATCTCGCTGTCTTC‑CGCTGCCTCGCGGATGCTCAAGGCAACCTc                                                                       >  1:1587157/1‑150 (MQ=255)

                                                                              tcACCTCGACCGCAGCACGCAGCAATGGGCTGTTGACGACCTTGATCAGGCGCGCAGACAGGGCCGCGACGCGGccct                                                                                                                                              >  1:1119817/1‑76 (MQ=255)

                                                                              tcACCTCGACCGCAGCACGCAGCAATGGGCTGTTGACGACCTTGATCAGGCGCGCAGACAGGGCCGCGACGCGGccct                                                                                                                                              <  2:1119817/78‑3 (MQ=255)

                                                                              tacgccagACTGCACCACGCAGCAATGGTTTGTTGACGACCTTGATCAGGCGCGCAGACAGGGCCGCGTCGCGGCCGATCACCTTGCTCAGGGCCGCTACGCTGATCTCGCTGTCTTC‑CGCTGCCTCGCGGATGCTCAAGGCAACCTCgg                                                                     <  2:2316996/143‑1 (MQ=255)

                                                                                gccgggACGGCTGCACGCAGCCATGGGCTGTTGACGACCTGGATGAGGCGCGCAAACAGGGCCGCGTCGCGCACGATCACCTTGCTCAGGGCCACTACACTCAACTCGCTGTCTTA‑CGCTTACTCGCGGATGCTCAAGGCAACCTCGGGc                                                                   >  2:127572/6‑150 (MQ=255)

                                                                                 ccTCGACCGCAGCACGCAGCAATGGGCTGTTGACGACCTTGATCAGGCGCGCAGACAGGGCCGCGTCGCGGCCGATCACCTTGCTCAGGGCCGCTACGCTGATCTCGCTGTCTTC‑CGCTGCCTCGCGGATGCTCAAGGCAACCTCGGGCa                                                                  <  1:942654/150‑1 (MQ=255)

                                                                                  cTCGACCGCAGCACGCAGCAATGGGCTGTTGACGACCTTGATCAGGCGCGCAGACAGGGCCGCGTCGCGGCCGATCACCTTGCTCAGGGCCGCTACGCTGATCTCGCTGTCTTC‑CGCTGCCTCGCGGATGCTCAAGGCAACCTCGGGCAg                                                                 <  2:660549/150‑1 (MQ=255)

                                                                                  cTCGACCGCAGCACGCAGCAATGGGCTGTTGACGACCTTGATCAGGCGCGCAGACAGGGCCGCGTCGCGGCCGATCACCTTGCTCAGGGCCGCTACGCTGATCTCGCTGTCTTC‑CGCTGCCTCGCGGATGCTCAAGGCAACCTCGGGCAg                                                                 >  1:1160532/1‑150 (MQ=255)

                                                                                  cTCGACCGCAGCACGCAGCAATGGGCTGTTGACGACCTTGATCAGGCGCGCAGACAGGGCCGCGTCGCGGCCGATCACCTTGCTCAGGGCCGCTACGCTGATCTCGCTGTCTTC‑CGCTGCCTCGCGGATGCTCAAGGCAACCTCGGGCAg                                                                 >  2:1160802/1‑150 (MQ=255)

                                                                                  cTCGACCGCAGCACGCAGCAATGGGCTGTTGACGACCTTGATCAGGCGCGCAGACAGGGCCGCGTCGCGGCCGATCACCTTGCTCAGGGCCGCTACGCTGATCTCGCTGTCTTC‑CGCTGCCTCGCGGATGCTCAAGGCAACCTCGGGCAg                                                                 >  2:1164016/1‑150 (MQ=255)

                                                                                   tCGACCGCAGCACGCAGCAATGGGCTGTTGACGACCTTGATCAGGCGCGCAGACAGGGCCGCGTCGCGGCCGATCACCTTGCTCAGGGCCGCTACGc                                                                                                                      <  2:2221933/97‑1 (MQ=255)

                                                                                   tCGACCGCAGCACGCAGCAATGGGCTGTTGACGACCTTGATCAGGCGCGCAGACAGGGCCGCGTCGCGGCCGATCACCTTGCTCAGGGCCGCTACGc                                                                                                                      >  1:2221933/1‑97 (MQ=255)

                                                                                   tCGACCGCAGCACGCAGCAATGGGCTGTTGACGACCTTGATCAGGCGCGCAGACAGGGCCGCGTCGCGGCCGATCACCTTGCTCAGGGCCGCTACGCTGATCTCGCTGTCTTC‑CGCTGCCTCGCGGATGCTCAAg                                                                               <  2:1187417/135‑1 (MQ=255)

                                                                                   tCGACCGCAGCACGCAGCAATGGGCTGTTGACGACCTTGATCAGGCGCGCAGACAGGGCCGCGTCGCGGCCGATCACCTTGCTCAGGGCCGCTACGCTGATCTCGCTGTCTTC‑CGCTGCCTCGCGGATGCTCAAGGCAACCTCggg                                                                    >  1:466313/1‑146 (MQ=255)

                                                                                   tCGACCGCAGCACGCAGCAATGGGCTGTTGACGACCTTGATCAGGCGCGCAGACAGGGCCGCGTCGCGGCCGATCACCTTGCTCAGGGCCGCTACGCTGATCTCGCTGTCTTC‑CGCTGCCTCGCGGATGCTCAAGGCAACCTCGGGCAgg                                                                >  2:2118222/1‑150 (MQ=255)

                                                                                        cgcagcacgcagcaATGGGCTGTTGACGACCTTGATCAGGCGCGCAGACAGGGCCGCGTCGCGGCCGATCACCTTGCTCAGGGCCGCTACGCTGATCTCGCTGTCTTC‑CGCTGCCTCGCGGATGCTCAAGGCAACCTCGGGCAGGGTTgg                                                           <  2:458953/150‑1 (MQ=255)

                                                                                         gcagcacgcagcaATGGGCTGTTGACGACCTTGATCAGGCGCGCAGACAGGGCCGCGTCGCGGCCGATCACCTTGCTCAGGGCCGCTACGCTGATCTCGCTGTCTTC‑CGCTGCCTCGCGGATGCTCAAGGCAACCTCGGGCAGGGTTGGc                                                          <  2:72467/150‑1 (MQ=255)

                                                                                         gcagcacgcagcaATGGGCTGTTGACGACCTTGATCAGGCGCGCAGACAGGGCCGCGTCGCGGCCGATCACCTTGCTCAGGGCCGCTACGCTGATCTCGCTGTCTTC‑CGCTGCCTCGCGGATGCTCAAGGCAACCTCGGGCAGGGTTGGc                                                          <  1:1164016/150‑1 (MQ=255)

                                                                                             cacgcagcaATGGGCTGTTGACGACCTTGATCAGGCGCGCAGACAGGGCCGCGTCGCGGCCGATCACCTTGCTCAGGGCCGCTACGCTGATCTCGCTGTCTTC‑CGCTGCCTCGCGGATGCTCAAGGCAACCTCGGGCAGGGTTGGCAGGa                                                      >  2:1812810/1‑150 (MQ=255)

                                                                                              acgcagcaATGGGCTGTTGACGACCTTGATCAGGCGCGCAGACAGGGCCGCGTCGCGGCCGATCACCTTGCTCAGGGCCGCTACGCTGATCTCGCTGTCTTC‑CGCTGCCTCGCGGATGATCAAGGCAACCTCGGGCAGGGTTGGCAGGAc                                                     <  2:885478/150‑1 (MQ=255)

                                                                                                gcagcaATGGGCTGTTGACGACTTTGATCAGGCGCGCAGACAGGGCCGCGTCGCGGCCGATCACCTTGCTCAGGGCCGCTACGCTGATCTCGCTGTCTTC‑CGCTGCCTCGCGGATGCTCAAGGGAACCTCGGGCAGGGGGGGCAAGatcg                                                   >  2:1178120/1‑147 (MQ=255)

                                                                                                gcagcaATGGGCTGTTGACGACCTTGATCAGGCGCGCAGACAGGGCCGCGTCGCGGCCGATCACCTTGCTCAGGGCCGCTACGCTGATCTCGCTGTCTTC‑CGCTGCCTCGCGGATGCTCAAGGCAACCTCGGGCAGGGTTGGCAGGACCa                                                   <  1:1301624/150‑1 (MQ=255)

                                                                                                gcagcaATGGGCTGTTGACGACCTTGATCAGGCGCGCAGACAGGGCCGCGTCGCGGCCGATCACCTTGCTCAGGGCCGCTACGCTGATCTCGCTGTCTTC‑CGCTGCCTCGCGGATGCTCAAGGCAACCTCGGGCAGGGTTGGCAGGACCa                                                   >  1:723790/1‑150 (MQ=255)

                                                                                                 cagcaATGGGCTGTTGACGACCTTGATCAGGCGCGCAGACAGGGCCGCGTCGCGGCCGATCACCTTGCTCAGGGCCGCTACGCTGATCTCGCTGTCTTC‑CGCTGCCTCGCGGATGCTCAAGGCAACCTCGGGCAGGGTTGGCAGGACCAg                                                  <  2:653657/150‑1 (MQ=255)

                                                                                                 cagcaATGGGCTGTTGACGACCTTGATCAGGCGCGCAGACAGGGCCGCGTCGCGGCCGATCACCTTGCTCAGGGCCGCTACGCTGATCTCGCTGTCTTC‑CGCTGCCTCGCGGATGCTCAAGGCAACCTCGGGCAGGGTTGGCAGGACCAg                                                  <  2:723790/150‑1 (MQ=255)

                                                                                                 cagcaATGGGCTGTTGACGACCTTGATCAGGCGCGCAGACAGGGCCGCGTCGCGGCCGATAACCTTACCCAGGGCCACTACGCTGATCTCGCTGTCTTC‑CGCTGCCTCGCGGATGCTCAAGGCAACCTCGGGCAGGGTTGGCAGGACCAg                                                  <  1:293636/150‑1 (MQ=255)

                                                                                                  agcaATGGGCTGTTGACGACCTTGATCAGGCGCGCAGACAGGGCCGCGTCGCGGCCGATCACCTTGCTCAGGGCCGCTACGCTGATCTCGCTGTCTTC‑CGCTGCCTCGCGGATGCTCAAGGCAACCTCGGGCAGGGTTGGCAGGACCAgg                                                 <  2:671058/150‑1 (MQ=255)

                                                                                                    caATGGGCTGTTGACGACCTTGATCAGGCGCGCAGACAGGGCCGCGTCGCGGCCGATCACCTTGCTCAGGGCCGCTACGCTGATCTCGCTGTCTTC‑CGCTGCCTCGCGGATGCTCAAGGCAACCTCGGGCAgg                                                                >  1:2386153/1‑133 (MQ=255)

                                                                                                    caATGGGCTGTTGACGACCTTGATCAGGCGCGCAGACAGGGCCGCGTCGCGGCCGATCACCTTGCTCAGGGCCGCTACGCTGATCTCGCTGTCTTC‑CGCTGCCTCGCGGATGCTCAAGGCAACCTCGGGCAgg                                                                <  2:2386153/133‑1 (MQ=255)

                                                                                                    caATGGGCTGTTGACGACCTTGATCAGGCGCGCAGACAGGGCCGCGTCGCGGCCGATCACCTTGCTCAGGGCCGCTACGCTGATCTCGCTGTCTTC‑CGCTGCCTCGCGGATGCTCAAGGCAACCTCGGGCAGGGTTGGCAGGACCAGGTc                                               <  1:1922685/150‑1 (MQ=255)

                                                                                                      aTGGGCTGTTGACGACCTTGATCAGGCGCGCAGACAGGGCCGCGTCGCGGCCGATCACCTTGCTCAGGGCCGCTACGCTGATCTCGCTGTCTTC‑CGCTGCCTCGCGGATGCTCAAGGCAACCTCGGGCAGGGTTGGCAGGACCAGGTCAt                                             <  1:1715346/150‑1 (MQ=255)

                                                                                                       tGGGCTGTTGACGACCTTGATCAGGCGCGCAGACAGGGGCGCGTCGCGGCCGATCACCTTGCTCAGGGCCGCTACGCTGATCTCGCTGTCTTC‑CGCTGCCTCGCGGATGCTCAAGGCAACCTCGGGCAGGGTTGGCAGGACCAGGTCATc                                            >  2:1395172/1‑150 (MQ=255)

                                                                                                         ggCTGTTGACGACCTTGATCAGGCGCGCAGACAGGGCCGCGTCGCGGCCGATCACCTTGCTCAGGGCCGCTACGCTGATCTCGCTGTCTTC‑CGCTGCCTCGCGGATGCTCaa                                                                                <  2:555578/112‑1 (MQ=255)

                                                                                                         ggCTGTTGACGACCTTGATCAGGCGCGCAGACAGGGCCGCGTCGCGGCCGATCACCTTGCTCAGGGCCGCTACGCTGATCTCGCTGTCTTC‑CGCTGCCTCGCGGATGCTCaa                                                                                >  1:555578/1‑112 (MQ=255)

                                                                                                               tgacagccTTGATCAGGCGCGCAGACAGGGCCGGGACGCGGCCGATCACTTTGCTCAGGGCCGCTACTCTGATCTCTCGGTCTTC‑CGCTGCCTCGCGGTTGCTCAAGGCAACCTCGGGCAGGGTTGGCAGGACCAGGGCATCCTTTTCGa                                    <  2:468166/144‑1 (MQ=255)

                                                                                                                 acgacCTTGATCAGGCGCGCAGACAGGGCCGCGTCGCGGCCGATCACCTTGCTCAGGGCCGCTACGCTGATCTCGCTGTCTTC‑CGCTGCCTCGCGGATGCTCAAGGCAACCTCGGGCAGGGTTGGCAGGACCAGGTCATCCTTTTCGATg                                  <  1:281566/150‑1 (MQ=255)

                                                                                                                     ccTTGATCAGGCGCGCATACAGGGCCGCGTCGCGGCCGATCACCTTGCTCAGGGCCGCTACGCTGATCTCGCTGTCTTC‑CGCTGCCTCGCGGATGCTCAAGGCAACCTCGGGCAGGGTTGGCAGGACCAGGTCATCCTTTTCGATGGCGt                              >  2:1845573/1‑150 (MQ=255)

                                                                                                                      cTTGATCAGGCGCGCAGACAGGGCCGCGTCGCGGCCGATCACCTTGCTCAGGGCCGCTACGCTGATCTNGCTGTCTTC‑CGCTGCCTCGCGGATGCTCAAGGCAACCTCGGGCa                                                                  <  1:1314996/113‑1 (MQ=255)

                                                                                                                      cTTGATCAGGCGCGCAGACAGGGCCGCGTCGCGGCCGATCACCTTGCTCAGGGCCGCTACGCTGATCTCGCTGTCTTC‑CGCTGCCTCGCGGATGCTCAAGGCAACCTCGGGCa                                                                  >  2:1314996/1‑113 (MQ=255)

                                                                                                                      cTTGATCAGGCGCGCAGACAGGGCCGCGTCGCGGCCGATCACCTTGCTCAGGGCCGCTACGCTGATCTCGCTGTCTTC‑CGCTGCCTCGCGGATGCTCAAGGCAACCTCGGGCAGGGTTg                                                            >  2:1788598/1‑119 (MQ=255)

                                                                                                                      cTTGATCAGGCGCGCAGACAGGGCCGCGTCGCGGCCGATCACCTTGCTCAGGGCCGCTACGCTGATCTCGCTGTCTTC‑CGCTGCCTCGCGGATGCTCAAGGCAACCTCGGGCAGGGTTg                                                            <  1:1788598/119‑1 (MQ=255)

                                                                                                                      cTTGATCAGGCGCGCAGACAGGGCCGCGTCGCGGCCGATCACCTTGCTCAGGGCCGCTACGCTGATCTCGCTGTCTTC‑CGCTGCCTCGCGGATGCTCAAGGCAACCTCGGGCAGGGTTGGCAGGACCAGGTCATCCTTTTCGATGGCGTc                             >  1:388806/1‑150 (MQ=255)

                                                                                                                      cTTGATCAGGAGCGCAGACAGGGCCGCGTCGCGGCCGATCACCTTGCTCAGGGCCGCTACGCTGATCTCGCTGTCTTC‑CGCTGCCTCGCGGATGCTCAAGGCAACCTCGGGCAGGGTTGGCAGGACCAGGTCATCCTTTTCGATGGCGTc                             >  2:2395900/1‑150 (MQ=255)

                                                                                                                              ggCGCGCAGACAGGGCCGCGTCGCGGCCGATCACCTTGCTCAGGGCCGCTACGCTGATCTCGCTGTCTTC‑CGCTGCCTCGCGGATGCTCAAGGCAACCTCGGGCAGGGTTGGCAGGACCAGGTCATCCTTTTCGATGGCGTCGAGCAAtt                     <  1:448988/150‑1 (MQ=255)

                                                                                                                              ggCGCGCAGACAGGGCCGCGTCGCGGCCGATCACCTTGCTCAGGGCCGCTACGCTGATCTCGCTGTCTTC‑CGCTGCCTCGCGGATGCTCAAGGCAACCTCGGGCAGGGTTGGCAGGACCAGGTCATCCTTTTCGATGGCGTCGAGCAAtt                     <  1:446195/150‑1 (MQ=255)

                                                                                                                                   gcAGACAGGGCCGCGTCGCGGCCGATCACCTTGCTCAGGGCCGCTACGCTGATCTCGCTGTCTTC‑CGCTGCCTCGCGGATGCTCAAGGCAACCTCGGGCAGGGTTGGCAGGACCAGGTCATCCTTTTCGATGGCGTCGAGCAATTGTGCt                <  1:1663107/150‑1 (MQ=255)

                                                                                                                                   gcAGACAGGGCCGCGTCGCGGCCGATCACCTTGCTCAGGGCCGCTACGCTGATCTCGCTGTCTTC‑CGCTGCCTCGCGGATGCTCAAGGCAACCTCGGGCAGGGTTGGCAGGACCAGGTCATCCTTTTCGATGGCGTCGAGCAATTGTGCt                <  1:1812810/150‑1 (MQ=255)

                                                                                                                                      gACAGGGCCGCGTCGCGGCCGATCACCTTGCTCAGGGCCGCTACGCTGATCTCGCTGTCTTC‑CGCTGCCTCGCGGATGCTCAAGGCAACCTCGGGCAGGGTTGGCAGGACCAGGTCATCCTTTTCGATGGCGTCGAGCAATTGTGCTTGa             >  1:2151240/1‑150 (MQ=255)

                                                                                                                                        cAGGGCCGCGTCGCGGCCGATCACCTTGCTCAGGGCCGCTACGCTGATCTCGCTGTCTTC‑CGCTGCCTCGCGGATGCTCAAGGCAACCTCGGGCAGGGTTGGCAGGACCAGGTCATCCTTTTCGATGGCGTCGAGCAATTGTGCTTGAAc           <  1:2081446/150‑1 (MQ=255)

                                                                                                                                        cAGGGCCGCGCCGCGGCCGATCACCTTGCTCAGGGCCGCTACGCTGATCTCGCTGTCTTC‑CGCTGCCTCGCGGATGCTCAAGGCAACCTCGGGCAGGGTTGGCAGGACCAGGTCATCCTTTTCGATGGCGTCGAGCAATTGTGCTTGAAc           <  1:239526/150‑1 (MQ=255)

                                                                                                                                         tcggCCGCGACGCGGCCGATCACCTTGCTCAGGGCCGCTACGCTGATCTCGCTGTCTTC‑CGCTGCCTCGCGGATGCTCAAGGCAACCTCGGGCAGGGTTGGCAGGACCa                                                   <  2:1689292/107‑1 (MQ=255)

                                                                                                                                           ggCCGCGTCGCGGCCGATCACCTTGCTCAGGGCCGCTACGCTGATNTCGCTGTCTTC‑CGCTGCCTCGCGGATGCTCAAGGCAACCTCGGGCAg                                                                 <  1:810416/93‑1 (MQ=255)

                                                                                                                                           ggCCGCGTCGCGGCCGATCACCTTGCTCAGGGCCGCTACGCTGATGTCGCTGTCTTC‑CGCTGCCTCGCGGATGCTCAAGGCAACCTCGGGCAg                                                                 >  2:810416/1‑93 (MQ=255)

                                                                                                                                           ggCCGCGTCGCGGCCGATCACCTTGCTCAGGGCCGCTACGCTGATCTCGCTGTCTTC‑CGCTGCCTCGCGGATGCTCAAGGCAACCTCGGGCAGGGTTGGCAGGACCAGGTCATCCTTTTc                                      <  2:213795/120‑1 (MQ=255)

                                                                                                                                             ccGCGTCGCGGCCGATCACCTTGCTCAGGGCCGCTACGCTGATCTCGCTGTCTTC‑CGCTGCCTCGCGGATGCTCAAGGCAACCTCGGGCAGGGTTGGCAGGACCAGGTCa                                              >  1:1876802/1‑110 (MQ=255)

                                                                                                                                             ccGCGTCGCGGCCGATCACCTTGCTCAGGGCCGCTACGCTGATCTCGCTGTCTTC‑CGCTGCCTCGCGGATGCTCAAGGCAACCTCGGGCAGGGTTGGCAGGACCAGGTCa                                              <  2:1876802/110‑1 (MQ=255)

                                                                                                                                             ccGCGTCGCGGCCGATCACCTTGCTCAGGGCCGCTACGCTGATCTCGCTGTCTTC‑CGCTGCCTCGCGGATGCTCAAGGCAACCTCGGGCAGGGTTGGCAGGACCAGGTCATCCTTTTCGATGGCGTCGAGCAATTGTGCTTGAACCAACt      >  1:2402386/1‑150 (MQ=255)

                                                                                                                                              cgcgTCGCGGCCGATCACCTTGCTCAGGGCCGCTACGCTGATCTCGCTGTCTTC‑CGCTGCCTCGCGGATGCTCAAGGCAACCTCGGGCa                                                                  <  1:382801/89‑1 (MQ=255)

                                                                                                                                              cgcgTCGCGGCCGATCACCTTGCTCAGGGCCGCTACGCTGATCTCGCTGTCTTC‑CGCTGCCTCGCGGATGCTCAAGGCAACCTCGGGCa                                                                  >  2:382801/1‑89 (MQ=255)

                                                                                                                                               gcgTCGCGGCCGATCACCTTGCTCAGGGCCGCTACGCTGATCTCGCTGTCTTC‑CGCTGCCTCGCGGATGCTCAAGGCAACCTCGGGCAGGGTTGGCAGGACCAGGTCATCCTTTTCGATGGCGTCGAGCAATTGTGCTTGAACCAACTCg    >  2:912652/1‑150 (MQ=255)

                                                                                                                                               gcgTCGCGGCCGATCACCTTGCTCAGGGCCGCTACGCTGATCTCGCTGTCTTC‑CGCTGCCTCGCGGATGCTCAAGGCAACCTCGGGCAGGGTTGGCAGGACCAGGTCATCCTTTTCGATGGCGTCGAGCAATTGTGCTTGAACCAACTCg    >  1:2229833/1‑150 (MQ=255)

                                                                                                                                                 gTCGCGGCCGATCAGCTTGCTCAGGGCCGCTACGCTGTTCTCGCTGTCTCCACGC‑GCCTCGCAGCTGCCCAAGGCAACCTCGGGCAGCGTCGGCAGTACCAGGTCCGCCTTGACGATGGAGTCGAGATCGTCTGTTGGAGCCAACTCGGc  <  2:639240/150‑1 (MQ=255)

TCTGCTCGATCACCAGGCCGATGGCCAGGTTGCAGCTGTAGTTGATGCCTAGCCGCGTGATGGCGGTATGCAGGTCGGTCACCTCGACCGCAGCACGCAGCAATGGGCTGTTGACGACCTTGATCAGGCGCGCAGACAGGGCCGCGTCGCGGCCGATCACCTTGCTCAGGGCCGCTACGCTGATCTCGCTGTCTTC‑CGCTGCCTCGCGGATGCTCAAGGCAACCTCGGGCAGGGTTGGCAGGACCAGGTCATCCTTTTCGATGGCGTCGAGCAATTGTGCTTGAACCAACTCGGC  >  NC_002947/1623086‑1623380

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