BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,992,357	A→T	5.8%	D96E (GAT→GAA)	ycaC‑III ←	putative hydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,992,357	0	A	T	5.8%	670.7 / 19.9	200	D96E (GAT→GAA)	ycaC‑III	putative hydrolase
Reads supporting (aligned to +/- strand): ref base A (109/77); new base T (10/2); total (119/81)
Fisher's exact test for biased strand distribution p-value = 1.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 208 total aligned reads displayed.

TCATTGAAGGTGCCGGAAGCATCGGTCACCACGAACACCTCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGTTGAAGAACTTGGCCA  >  NC_002947/5992208‑5992505

     gAAGGTGCCGGAAGCATCGGTCACCACGAACACCTCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTc                                                                                                                                                 <  2:4817108/150‑1 (MQ=255)

              ggAAGCATCGGTCACCACGAACACCTCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCa                                                                                                                                        <  1:3070224/150‑1 (MQ=255)

                 aGCATCGGTCACCACGAACACCTCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGc                                                                                                                                     >  1:4255927/1‑150 (MQ=255)

                      cGGTCACCACGAACACCTCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAg                                                                                                                                       >  2:1800044/1‑143 (MQ=255)

                      cGGTCACCACGAACACCTCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAg                                                                                                                                       <  1:1800044/143‑1 (MQ=255)

                         tCACCACGAACACCTCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCt                                                                                                                             >  1:3756250/1‑150 (MQ=255)

                           accacGAACACCTCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTgg                                                                                                                           >  2:5125950/1‑150 (MQ=255)

                             cacGAACACCTCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGcc                                                                                                                         >  2:2362472/1‑150 (MQ=255)

                              acGAACACCTCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGaaa                                                                                                                                                      <  1:1843274/120‑1 (MQ=255)

                                   cacCTCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAgggcg                                                                                                                   <  1:3954489/150‑1 (MQ=255)

                                    acCTCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAgggcgg                                                                                                                  >  1:353518/1‑150 (MQ=255)

                                     ccTCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAgggcggg                                                                                                                 <  1:1265163/150‑1 (MQ=255)

                                      cTCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAgggcgggc                                                                                                                >  2:4961060/1‑150 (MQ=255)

                                       tCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCa                                                                                                               >  2:425837/1‑150 (MQ=255)

                                       tCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCa                                                                                                               >  2:4475911/1‑150 (MQ=255)

                                       tCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCa                                                                                                               >  2:1190883/1‑150 (MQ=255)

                                       tCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCa                                                                                                               >  1:2792817/1‑150 (MQ=255)

                                        cAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCaa                                                                                                              <  2:2334376/150‑1 (MQ=255)

                                          aaTTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATg                                                                                                            >  2:3564463/1‑150 (MQ=255)

                                               ttcCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACgg                                                                                                       >  2:692892/1‑150 (MQ=255)

                                               ttcCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACgg                                                                                                       >  1:1760549/1‑150 (MQ=255)

                                               gtcctccAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACgg                                                                                                       <  1:1063752/149‑1 (MQ=255)

                                                   tccAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCa                                                                                                   >  2:962365/1‑150 (MQ=255)

                                                   tccAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCa                                                                                                   >  2:514821/1‑150 (MQ=255)

                                                     cAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATc                                                                                                 >  2:1782154/1‑150 (MQ=255)

                                                       aGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACCGCGCATCCg                                                                                               <  2:3067346/150‑1 (MQ=255)

                                                        ggCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCgg                                                                                              <  2:1085672/150‑1 (MQ=255)

                                                        ggCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCgg                                                                                              <  2:4110205/150‑1 (MQ=255)

                                                           ggacaatgcCGGGAAAGCCACGCAGACTTCGGTCACAACACCTGCGATGAGCAACTGCTTCTTGCCGGTGGCCTTCACTGCCTTGACGAATTCCTCATTGTCCCAGGCGTTGATCTGGCCGGGGCGAGCGATATAAGGCGCATCCGGGaa                                                                                           >  1:4130655/8‑150 (MQ=14)

                                                               aatgcCGGGAAAGCCACGCAGACTTCGGTCACAACACCTGCGATGAGCAACTGCTTCTTGCCGGTGGCCTTCACTGCCTTGACGAATTCCTCATTGTCCCAGGCGTTGATCTGGCCGGGGCGAGCGATATAAGGCGCATCCGGGaatgtt                                                                                       >  1:902679/4‑146 (MQ=14)

                                                               aGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAgg                                                                                       <  2:1824554/150‑1 (MQ=255)

                                                                 cgcCGGGAACGCCACGCAAACCGCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCACTGTACGGCGCATCCGGGAACAGGGc                                                                                     >  2:1856161/1‑150 (MQ=255)

                                                                  gcCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTAAGGCGCATCCGGGAACAGGGCt                                                                                    >  2:5450138/1‑150 (MQ=255)

                                                                    cGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCttt                                                                                  <  1:2259333/150‑1 (MQ=255)

                                                                     gggAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTc                                                                                 <  1:1782154/150‑1 (MQ=255)

                                                                     gggAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTc                                                                                 <  2:2063915/150‑1 (MQ=255)

                                                                      ggAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCa                                                                                >  2:403475/1‑150 (MQ=255)

                                                                       gAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGc                                                                                                                          >  2:769314/1‑107 (MQ=255)

                                                                       gAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGc                                                                                                                          <  1:769314/107‑1 (MQ=255)

                                                                         aCGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGtt                                                                             <  1:129293/150‑1 (MQ=255)

                                                                           gCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAg                                                                                        >  1:1658189/1‑137 (MQ=255)

                                                                              aCGCAGACTTCGGTCACAACACCTGCGATGAGCAACTGCTTCTTGCCGGTGGCCTTCACTGCCTTGACGAATTCCTCATTGTCCCAGGCGTTGATCTGGCCGGGGCGAGCGATATAAGGCGCATCCGGGAATGTTTCTTTCAGCTCCGGc                                                                        >  1:2179180/1‑150 (MQ=9)

                                                                              aCGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGc                                                                        >  1:1363502/1‑150 (MQ=255)

                                                                               cGCAGACTTCGGTCACAACACCTGCGATGAGCAACTGCTTCTTGCCGGTGGCCTTCACTGCCTTGACGAATTCCTCATTGTCCCAGGCGTTGATCTGGCCGGGGCGAGCGATATAAGGCGCATCCGGGaat                                                                                          <  2:5183214/131‑2 (MQ=14)

                                                                                gCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCAc                                                                      >  2:2871777/1‑150 (MQ=255)

                                                                                  agacttCGGTCACAACACCTGCGATGAGCAACTGCTTCTTGCCGGTGGCCTTCACTGCCTTGACGAATTCCTCATTGTCCCAGGCGTTGATCTGGCCGGGGCGAGCGATATAAGGCGCATCCGGGAATGTTTCTTTCAGCTCCGGCACCa                                                                    >  1:1007241/6‑150 (MQ=9)

                                                                                      cTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCgg                                                                >  2:2921808/1‑150 (MQ=255)

                                                                                       tCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAAt                                                                                                             <  2:2083784/104‑1 (MQ=255)

                                                                                       tCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAAt                                                                                                             >  1:2083784/1‑104 (MQ=255)

                                                                                            cacaACACCTGCCAGGCTCAGCTGCGTCGTGCCGGTCGCCTTCACCGCCTGGACGAACTCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGTGGGCCGtt                                                          <  2:3936728/150‑1 (MQ=17)

                                                                                             acaacaCCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAg                                                                               >  1:2925092/1‑128 (MQ=255)

                                                                                             acaacaCCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTg                                                         >  1:5327710/1‑150 (MQ=255)

                                                                                             acaacaCCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTg                                                         >  1:1550718/1‑150 (MQ=255)

                                                                                               aacaCCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTggg                                                       <  2:3749129/150‑1 (MQ=255)

                                                                                               aacaCCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTggg                                                       >  2:129017/1‑150 (MQ=255)

                                                                                                  acCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGcc                                                    <  1:4278665/150‑1 (MQ=255)

                                                                                                   ccTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCt                                                   >  2:266979/1‑150 (MQ=255)

                                                                                                    cTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCtt                                                  >  2:3807284/1‑150 (MQ=255)

                                                                                                     tGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTg                                                 >  2:3568527/1‑150 (MQ=255)

                                                                                                     tGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTg                                                 >  1:5194371/1‑150 (MQ=255)

                                                                                                      gCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCCCCAACGGCCCGCTGGTGccggaa                                                <  2:166865/150‑5 (MQ=17)

                                                                                                        aaTGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTc                                              <  2:2675797/150‑1 (MQ=255)

                                                                                                        aaTGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTc                                              >  1:5104914/1‑150 (MQ=255)

                                                                                                                gctgctTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCtgg                                      >  1:4439795/1‑150 (MQ=255)

                                                                                                                 ctgctTCTTGCCGGTGGCCTTCACTGCCTTGACGAATTCCTCATTGTCCCAGGCGTTGATCTGGCCGGGGCGAGCGATATAAGGCGCATCCGGGAATGTTTCTTTCAGCTCCGGCACCAATGGGCCATTCGGGCCGTTTTCAAAACtggt                                     <  2:3653484/150‑1 (MQ=2)

                                                                                                                 ctgctTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCtggt                                     >  2:5363785/1‑150 (MQ=255)

                                                                                                                 ctgctTCTTGCCGGTCGCCTTCACCACCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGaa                                           <  1:916746/144‑1 (MQ=255)

                                                                                                                 ctgctTCTTGCCGGTCGCCTTCACCACCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGaa                                           >  2:916746/1‑144 (MQ=255)

                                                                                                                  tgctTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCtggtg                                    <  1:2521532/150‑1 (MQ=255)

                                                                                                                   gctTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCtggtgg                                   <  1:4091934/150‑1 (MQ=255)

                                                                                                                      tcttGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGa                                >  1:1978067/1‑150 (MQ=255)

                                                                                                                      tcttGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGa                                <  2:2678491/150‑1 (MQ=255)

                                                                                                                        ttGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATcc                                                                                                <  2:1652858/84‑1 (MQ=255)

                                                                                                                           ccGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATg                           <  1:3941579/150‑1 (MQ=255)

                                                                                                                              gTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTc                        >  2:3130161/1‑150 (MQ=255)

                                                                                                                                  ccTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGTCGTTGATCTGGCCAGGGCGGGAAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCa                    <  2:1509872/150‑1 (MQ=255)

                                                                                                                                  ccTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGTGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCa                    <  1:3007022/150‑1 (MQ=255)

                                                                                                                                   cTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAg                   >  1:3923047/1‑150 (MQ=255)

                                                                                                                                   cTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAg                   >  2:2925840/1‑150 (MQ=255)

                                                                                                                                    ttCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAgg                  <  1:5009997/150‑1 (MQ=255)

                                                                                                                                    ttCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAgg                  >  1:946252/1‑150 (MQ=255)

                                                                                                                                     tCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGt                 >  1:1554167/1‑150 (MQ=255)

                                                                                                                                       aCCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGTTg               <  2:4558131/150‑1 (MQ=255)

                                                                                                                                       aCCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGTTg               >  1:2748543/1‑150 (MQ=255)

                                                                                                                                         cGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGTTgaa             <  2:2611004/150‑1 (MQ=255)

                                                                                                                                          gCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCtggtgg                                   <  1:2702953/127‑1 (MQ=255)

                                                                                                                                          gCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCtggtgg                                   >  2:2702953/1‑127 (MQ=255)

                                                                                                                                            cTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGTTgaagaa          <  2:3923047/150‑1 (MQ=255)

                                                                                                                                             ttGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCg                                             >  1:4468837/1‑114 (MQ=255)

                                                                                                                                             ttGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGTTGAAGAAc         >  1:3818604/1‑150 (MQ=255)

                                                                                                                                               gACGAACTCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGTTGAAGAACtt       <  2:3086693/150‑1 (MQ=255)

                                                                                                                                               gACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCgg                                                                                              >  1:490344/1‑63 (MQ=255)

                                                                                                                                               gACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGTTGAAGAACtt       <  2:4655223/150‑1 (MQ=255)

                                                                                                                                               gACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGTTGAAGAACtt       >  1:2461923/1‑150 (MQ=255)

                                                                                                                                               gACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTAGAAGCtgg                                      <  1:101547/119‑1 (MQ=255)

                                                                                                                                               gACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTAGAAGCtgg                                      >  2:101547/1‑119 (MQ=255)

                                                                                                                                                  gAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCa                                                                                >  2:2000834/1‑74 (MQ=255)

                                                                                                                                                  gAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGTTGAAGAACTTGGc    >  2:822708/1‑150 (MQ=255)

                                                                                                                                                   aaaTCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCtggt                                     <  2:3268223/116‑1 (MQ=255)

                                                                                                                                                    aaTCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGTTGAAGAACTTGGCCa  >  1:315907/1‑150 (MQ=255)

TCATTGAAGGTGCCGGAAGCATCGGTCACCACGAACACCTCAAATTCTTCCTCCAAGGCCGACAGCGCCGGGAACGCCACGCAAACCTCGGTCACAACACCTGCAATGATCAGCTGCTTCTTGCCGGTCGCCTTCACCGCCTTGACGAAATCCTCGTTGTCCCAGGCGTTGATCTGGCCAGGGCGGGCAATGTACGGCGCATCCGGGAACAGGGCTTTCAGTTCCGGCACCAGCGGGCCGTTGGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGTTGAAGAACTTGGCCA  >  NC_002947/5992208‑5992505

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