BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,289,993	T→A	5.9%	L90* (TTG→TAG) ‡	PP_1122 →	OmpA family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,289,993	0	T	A	5.9%	923.8 / 30.2	279	L90* (TTG→TAG) ‡	PP_1122	OmpA family protein
Reads supporting (aligned to +/- strand): ref base T (102/160); new base A (7/10); total (109/170)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

Only 100 of 285 total aligned reads displayed.

CAGCCTGCTCACGGGCTGCGCTACGCACAGCGATGGCAGCGCGCCCCTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCGATCACTGCCCCGAGACCCCGGCCAACACGGCGGTTGACCACATGGGCTGCCCACTGAAACAGTACCCGGCC  >  NC_002947/1289844‑1290141

cAGCCTGCTCACGGGCTGCGCTACGCACAGCGATGGCAGCGCGCCCCTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCtt                                                                                                                                                      <  2:6656620/150‑1 (MQ=255)

cAGCCTGCTCACGGGCTGCGCTACGCACAGCGATGGCAGCGCGCCCCTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCtt                                                                                                                                                      >  2:3215031/1‑150 (MQ=255)

 aGCCTGCTCACGGGCTGCGCTACGCACAGCGATGGCAGCGCGCCCCTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTg                                                                                                                                                     >  1:3406891/1‑150 (MQ=255)

  gCCTGCTCACGGGCTGCGCTACGCACAGCGATGGCAGCGCGCCCCTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCtttg                                                                                                                                                    >  1:5198766/1‑148 (MQ=255)

  gCCTGCTCACGGGCTGCGCTACGCACAGCGATGGCAGCGCGCCCCTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTgg                                                                                                                                                    <  2:5530417/150‑1 (MQ=255)

    cTGCTCACGGGCTGCGCTACGCACAGCGATGGCAGCGCGCCCCTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGgc                                                                                                                                                  <  1:3114507/150‑1 (MQ=255)

     tGCTCACGGGCTGCGCTACGCACAGCGATGGCAGCGCGCCCCTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGgcg                                                                                                                                                 >  1:3788680/1‑150 (MQ=255)

     tGCTCACGGGCTGCGCTACGCACAGCGATGGCAGCGCGCCCCTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGgcg                                                                                                                                                 >  1:5474632/1‑150 (MQ=255)

     tGCTCACGGGCTGCGCTACGCACAGCGATGGCAGCGCGCCCCTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGgcg                                                                                                                                                 >  2:105812/1‑150 (MQ=255)

        tCACGGGCTGCGCTACGCACAGCGATGGCAGCGCGCCCCTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCa                                                                                                                                              >  2:3693700/1‑150 (MQ=255)

        tCACGGGCTGCGCTACGCACAGCGATGGCAGCGCGCCCCTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCa                                                                                                                                              >  1:591602/1‑150 (MQ=255)

        tCACGGGCTGCGCTACGCACAGCGATGGCAGCGCGCCCCTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCa                                                                                                                                              >  1:6627506/1‑150 (MQ=255)

        tCACGGGCTGCGCTACGCACAGCGATGGCAGCGCGCCCCTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCa                                                                                                                                              <  1:981281/150‑1 (MQ=255)

          aCGGGCTGCGCTACGCACAGCGATGGCAGCGCGCCCCTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCAtt                                                                                                                                            <  1:2245610/150‑1 (MQ=255)

          aCGGGCTGCGCTACGCACAGCGATGGCAGCGCGCCCCTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCAtt                                                                                                                                            <  2:6196654/150‑1 (MQ=255)

            gggCTGCGCTACGCACAGCGATGGCAGCGCGCCCCTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGc                                                                                                                                          >  1:2387271/1‑150 (MQ=255)

              gCTGCGCTACGCACAGCGATGGCAGCGCGCCCCTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCAGTGCcg                                                                                                                                        >  1:4858399/1‑150 (MQ=255)

              gCTGCGCTACGCACAGCGATGGCAGCGCGCCCCTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGCAATGGCGCCCAAGGCGCCaccaccg                                                                                                                                        <  1:5696133/150‑7 (MQ=255)

                    cTACGCACAGCGATGGCAGCGCGCCCCTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGc                                                                                                                                  <  2:4823220/150‑1 (MQ=255)

                    cTACGCACAGCGATGGCAGCGCGCCCCTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGc                                                                                                                                  <  2:5289061/150‑1 (MQ=255)

                     tACGCACAGCGATGGCAGCGCGCCCCTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCt                                                                                                                                 <  1:4129364/150‑1 (MQ=255)

                      aCGCACAGCGATGGCAGCGCGCCCCTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGCAATGGCGCCCAAGGCGCCACCACcggcggccca                                                                                                                                <  2:634748/150‑5 (MQ=37)

                       cGCACAGCGATGGCAGCGCGCCCCTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGa                                                                                                                               <  1:2355382/150‑1 (MQ=255)

                       cGCACAGCGATGGCAGCGCGCCCCTCAATCAAAGGACCTGGCCCATCTACAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCAGCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGa                                                                                                                               >  2:180962/1‑150 (MQ=255)

                          acaGCGATGGCAGCGCGCCCCTCAATCAAAGGACCTGGCCCATCTGCAGCGTGCTGGGCGGCTTGGTCGTTGGTTGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGAttt                                                                                                                            <  1:1680936/150‑1 (MQ=255)

                          acaGCGATGGCAGCGCGCCCCTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGAttt                                                                                                                            <  2:1400504/150‑1 (MQ=255)

                          acaGCGATGGCAGCGCGCCCCTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGAttt                                                                                                                            >  2:274251/1‑150 (MQ=255)

                          acaGCGATGGCAGCGCGCCCCTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGAttt                                                                                                                            <  1:5097085/150‑1 (MQ=255)

                          acaGCGATGGCAGCGCGCCCCTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGAttt                                                                                                                            <  1:5002167/150‑1 (MQ=255)

                          acaGCGATGGCAGCGCGCCCCTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGAttt                                                                                                                            <  2:2585938/150‑1 (MQ=255)

                            aGCGATGGCAGCGCGCCCCTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGt                                                                                                                          <  2:232366/150‑1 (MQ=255)

                               gATGGCAGCGCGCCCCTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCGGATTTGTTAc                                                                                                                       >  1:2894723/1‑150 (MQ=255)

                                         cgcCCCTCAATCCAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACg                                                                                                             <  1:2100854/150‑1 (MQ=255)

                                         cgcCCCTCAATCCAAGGACCTGGCCCATCGGCAGCCTGCGGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCGTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGCCg                                                                                                             <  1:4447391/150‑1 (MQ=255)

                                           ccccTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGt                                                                                                           <  2:2400387/150‑1 (MQ=255)

                                            cccTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTg                                                                                                          <  2:4730483/150‑1 (MQ=255)

                                             ccTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGa                                                                                                         <  2:2062008/150‑1 (MQ=255)

                                              cTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGAc                                                                                                        >  2:4925513/1‑150 (MQ=255)

                                              cTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGAc                                                                                                        >  1:1263082/1‑150 (MQ=255)

                                              cTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGAc                                                                                                        >  2:1610731/1‑150 (MQ=255)

                                               tcaatcaaAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACg                                                                                                       >  2:3770230/1‑150 (MQ=255)

                                               tcaatcaaAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACg                                                                                                       >  2:3533482/1‑150 (MQ=255)

                                                 aatcaaAGGACCTGGCCCATCTGCTGCCTGCTGGGCGGATTGGTCGGTGGTGGCCTTGGTGCAATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCTTTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGaa                                                                                                     <  1:4140176/150‑1 (MQ=255)

                                                  atcaaAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCCTTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAg                                                                                                    >  2:4061020/1‑150 (MQ=255)

                                                  atcaaAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAg                                                                                                    <  1:2610320/150‑1 (MQ=255)

                                                   tcaaAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGa                                                                                                   <  2:4252249/150‑1 (MQ=255)

                                                    caaAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGAt                                                                                                  <  2:564716/150‑1 (MQ=255)

                                                    caaAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGAt                                                                                                  <  1:2469416/150‑1 (MQ=255)

                                                     aaaGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATg                                                                                                 <  1:4222882/150‑1 (MQ=255)

                                                      aaGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCCTTGCCGGCGGGCGGATTTGTTACGCCCCGGACGGTGACGAAGATgg                                                                                                >  1:79988/1‑150 (MQ=255)

                                                      aaGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATgg                                                                                                <  2:2167929/150‑1 (MQ=255)

                                                       aGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGt                                                                                               <  2:5474632/150‑1 (MQ=255)

                                                         gACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGa                                                                                             <  2:5051517/150‑1 (MQ=255)

                                                            cTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACgg                                                                                          >  2:511965/1‑150 (MQ=255)

                                                                  cATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCAtttt                                                                                    >  1:1558364/1‑150 (MQ=255)

                                                                  cATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCAtttt                                                                                    >  2:4400552/1‑150 (MQ=255)

                                                                   aTCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCCTTGCCGGCGGGCTGATTTGTTAAGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTc                                                                                   >  1:2273531/1‑150 (MQ=255)

                                                                   aTCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTc                                                                                   <  2:4520632/150‑1 (MQ=255)

                                                                   aTCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTc                                                                                   <  2:4672901/150‑1 (MQ=255)

                                                                      tGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGAc                                                                                <  1:3670905/150‑1 (MQ=255)

                                                                      tGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGAc                                                                                >  1:4788032/1‑150 (MQ=255)

                                                                       gCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGAcc                                                                               >  2:896387/1‑150 (MQ=255)

                                                                          gCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGgc                                                                            <  1:2061194/150‑1 (MQ=255)

                                                                            ctgctgGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGgcgc                                                                          <  2:392603/150‑1 (MQ=255)

                                                                             tgctgGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTgg                                                                                                                                                    >  2:3366687/1‑75 (MQ=255)

                                                                              gcggggCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCGa                                                                        <  1:6099095/147‑1 (MQ=255)

                                                                               ctgGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCGAt                                                                       <  2:5820224/150‑1 (MQ=255)

                                                                                 gggCGGCTTGGTCGGTGGGGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTTCGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCGATCa                                                                     <  2:3322820/150‑1 (MQ=255)

                                                                                  ggcggcTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCGATCAc                                                                    >  1:3915401/1‑150 (MQ=255)

                                                                                  ggcggcTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCGATCAc                                                                    <  1:6137487/150‑1 (MQ=255)

                                                                                   gcggcTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTCCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCGATCACt                                                                   <  1:2159070/150‑1 (MQ=255)

                                                                                        ttGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTTGACCGGCGCGATCACTGcccc                                                              >  1:5412354/1‑150 (MQ=255)

                                                                                        ttGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCGATCACTGcccc                                                              <  2:2434040/150‑1 (MQ=255)

                                                                                        ttGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCGATCACTGcccc                                                              <  1:187083/150‑1 (MQ=255)

                                                                                          ggTCGGTGGTGGCCGTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCGATCACTGCCCCga                                                            <  2:5690247/150‑1 (MQ=255)

                                                                                               gtggtgGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCGATCACTGCCCCGAGAccc                                                       >  1:6274633/1‑150 (MQ=255)

                                                                                               gtggtgGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTAAGCCCAGGACGGTGAAGAAGATGGTGAAGGCATTTTCGACCGGCGCGATCACTGCCCCGAGAccc                                                       >  1:3532554/1‑150 (MQ=255)

                                                                                                  gtggCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCGATCACTGCCCCGAGACCCCgg                                                    <  2:3106955/150‑1 (MQ=255)

                                                                                                       cTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCGATCACTGCCCCGAGACCCCGGCCAac                                               >  1:1125119/1‑150 (MQ=255)

                                                                                                            tGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCGATCACTGCCCCGAGACCCCGGCCAACAcggc                                          <  1:3455581/150‑1 (MQ=255)

                                                                                                                aTCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCGATCACTGCCCCGAGACCCCGGCCAACACGGCGGtt                                      <  2:3868317/150‑1 (MQ=255)

                                                                                                                aTCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCGATCACTGCCCCGAGACCCCGGCCAACACGGCGGtt                                      >  2:829231/1‑150 (MQ=255)

                                                                                                                aTCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCGATCACTGCCCCGAGACCCCGGCCAACACGGCGGtt                                      <  2:485845/150‑1 (MQ=255)

                                                                                                                 tCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGTTTTGTTACGCCCAGGGCGGTGACGTCGATGGTGACGGCATTTTCGACCGGCGCGGTCACTGCCCCGAGACCCCGGCCAACACCGCGGGTg                                     >  1:4723429/1‑150 (MQ=255)

                                                                                                                   gagaGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCGATCACTGCCCCGAGACCCCGGCCAACACGGCGGTTGAc                                   <  1:2692419/150‑1 (MQ=255)

                                                                                                                   gagaGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCGATCACTGCCCCGAGACCCCGGCCAACACGGCGGTTGAc                                   >  2:1227518/1‑150 (MQ=255)

                                                                                                                     gagTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCGATCACTGCCCCGAGACCCCGGCCAACACGGCGGTTGACca                                 <  1:3515260/150‑1 (MQ=255)

                                                                                                                     gagTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTAAGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCGATCACTGCCCCGAGACCCCGGCCAACACGGCGGTTGACca                                 >  1:4236203/1‑150 (MQ=255)

                                                                                                                       gTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCGATCACTGCCCCGAGACCCCGGCCAACACGGCGGTTGACcaca                               >  2:5350175/1‑150 (MQ=255)

                                                                                                                          cttcCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCGATCACTGCCCCGAGACCCCGGCCAACACGGCGGTTGACCACATgg                            <  2:6586339/150‑1 (MQ=255)

                                                                                                                              cTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCGATCACTGCCCCGAGACCCCGGCCAACACGGCGGTTGACCACATGGGCTg                        >  1:6811659/1‑150 (MQ=255)

                                                                                                                                 ggCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCGATCACTGCCCCGAGACCCCGGCCAACACGGCGGTTGACCACATGGGCTGccc                     <  1:1397841/150‑1 (MQ=255)

                                                                                                                                  cccgccgGCAATGGCGCCCAAGGCGCCACCACCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCGATCACTGCCCCGAGACCCCGGCCAACACGGCGGTTGACCACATGGGCTGCCCa                    >  2:3493761/2‑150 (MQ=38)

                                                                                                                                        ggtggtGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCGATCACTGCCCCGAGACCCCGGCCAACACGGCGGTTGACCACATGGGCTGCCCACTGaaa              <  2:5089444/150‑1 (MQ=255)

                                                                                                                                          tggtggCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCGATCACTGCCCCGAGACCCCGGCCAACACGGCGGTTGACCACATGGGCTGCCCACTGAAACa            <  1:880816/150‑1 (MQ=255)

                                                                                                                                            gtggCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCGATCACTGCCCCGAGACCCCGGCCAACACGGCGGTTGACCACATGGGCTGCCCACTGAAACAGt          <  1:5070947/150‑1 (MQ=255)

                                                                                                                                              ggCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCGATCACTGCCCCGAGACCCCGGCCAACACGGCGGTTGACCACATGGGCTGCCCACTGAAACAGTAc        >  2:14487/1‑150 (MQ=255)

                                                                                                                                               gcgcCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCGATCACTGCCCCGAGACCCCGGCCAACACGGCGGTTGACCACATGGGCTGCCCACTGAAACAGTAcc       <  2:6022526/150‑1 (MQ=255)

                                                                                                                                               gcgcCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCGATCACTGCCCCGAGACCCCGGCCAACACGGCGGTTGACCACATGGGCTGCCCACTGAAACAGTAcc       >  1:1020293/1‑150 (MQ=255)

                                                                                                                                                    ttGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCGATCACTGCCCCGAGACCCCGGCCAACACGGCGGTTGACCACATGGGCTGCCCACTGAAACAGTACCCGGcc  >  1:562397/1‑150 (MQ=255)

CAGCCTGCTCACGGGCTGCGCTACGCACAGCGATGGCAGCGCGCCCCTCAATCAAAGGACCTGGCCCATCTGCAGCCTGCTGGGCGGCTTGGTCGGTGGTGGCCTTGGTGCCATCGAGAGTTCTTCCTGGGCCGCCGGTGGTGGCGCCTTGGGCGCCATTGCCGGCGGGCTGATTTGTTACGCCCAGGACGGTGACGAAGATGGTGACGGCATTTTCGACCGGCGCGATCACTGCCCCGAGACCCCGGCCAACACGGCGGTTGACCACATGGGCTGCCCACTGAAACAGTACCCGGCC  >  NC_002947/1289844‑1290141

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