BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,336,918:1	+T	6.9%	coding (609/1704 nt)	PP_2932 ←	amidase family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,336,918	1	.	T	6.9%	728.7 / 46.8	175	coding (609/1704 nt)	PP_2932	amidase family protein
Reads supporting (aligned to +/- strand): ref base . (78/85); new base T (5/7); total (83/92)
Fisher's exact test for biased strand distribution p-value = 7.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

Only 100 of 182 total aligned reads displayed.

TGCAGGCGCCACAGGTCGCCGCGCTTGTCGGCGTCATCGGCTACCACCACGTCGAGGATTTCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTGGCGGCGGTGGCGGTACCGGCACCGTTGGACGAGCCCGAGGCGAAGGGCGCGGTCAG  >  NC_002947/3336769‑3337067

tGCAGGCGCCACAGGTCGCCGCGCTTGTCGGCGTCATCGGCTACCACCACGTCGAGGATTTCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCcg                                                                                                                                                        <  1:1361769/150‑1 (MQ=255)

     gcgcCACAGGTCGCCGCGCTTGTCGGCGTCATCGGCTACCACCACGTCGAGGATTTCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGa                                                                                                                                                  >  1:618681/1‑150 (MQ=255)

         cacaGGTCGCCGCGCTTGTCGGCGTCATCGGCTACCACCACGTCGAGGATTTCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACggg                                                                                                                                              >  1:4251622/1‑150 (MQ=255)

         cacaGGTCGCCGCGCGTGTCGGCGTCATCGGCTACCACCACGTCGAGGATTTCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACggg                                                                                                                                              >  2:878280/1‑150 (MQ=255)

           caGGTCGCCGCGCTTGTCGGCGTCATCGGCTACCACCACGTCGAGGATTTCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGgt                                                                                                                                            <  2:723163/150‑1 (MQ=255)

           caGGTCGCCGCGCTTGTCGGCGTCATCGGCTACCACCACGTCGAGGATTTCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGgt                                                                                                                                            >  1:3401155/1‑150 (MQ=255)

              gTCGCCGCGCTTGTCGGCGTCATCGGCTACCACCACGTCGAGGATTTCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCGTCGCG‑CGGGgtgat                                                                                                                                         <  2:6390/150‑3 (MQ=255)

                cgccgcGCTTGTCGGCGTCATCGGCTACCACCACGTCGAGGATTTCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCGTCGCG‑CGGGgtgattt                                                                                                                                       <  2:2067606/150‑5 (MQ=255)

                cgccgcGCTTGTCGGCGTCATCGGCTACCACCACGTCGAGGATTTCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCGTCGCG‑CGGGgtgattt                                                                                                                                       <  2:1677691/150‑5 (MQ=255)

                 gccgcGCTTGTCGGCGTCATCGGCTACCACCACGTCGAGGATTTCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCGTCGCG‑CGGGgtgatttc                                                                                                                                      <  2:2326375/150‑6 (MQ=255)

                  ccgcGCTTGTCGGCGTCATCGGCTACCACCACGTCGAGGATTTCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCGTCGCG‑CGGGgtgatttcg                                                                                                                                     <  1:4112474/150‑7 (MQ=255)

                        ttGTCGGCGTCATCGGCTACCACCACGTCGAGGATTTCAAGTAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAgg                                                                                                                               <  1:2895913/150‑1 (MQ=255)

                             ggCGTCATCGGCTACCACCACGTCGAGGATTTCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCgtt                                                                                                                          <  2:394932/150‑1 (MQ=255)

                                gTCATCGGCTACCACCACGTCGAGGATTTCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGgc                                                                                                                                      >  2:2363770/1‑135 (MQ=255)

                                gTCATCGGCTACCACCACGTCGAGGATTTCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGgc                                                                                                                                      <  1:2363770/135‑1 (MQ=255)

                                gTCATCGGCTACCACCACGTCGAGGATTTCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCgttgtt                                                                                                                       <  2:4998453/150‑1 (MQ=255)

                                gTCATCGGCTACCACCACGTCGAGGATTTCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCgttgtt                                                                                                                       <  2:3194125/150‑1 (MQ=255)

                                gTCATCGGCGACCACCACGTCGAGGATTTCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCgttgtt                                                                                                                       <  1:3457486/150‑1 (MQ=255)

                                 tCATCGGCTACCACCACGTCGAGGATTTCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGttgttg                                                                                                                      <  1:130501/150‑1 (MQ=255)

                                  cATCGGCTACCACCACGTCGAGGATTTCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTgg                                                                                                                     >  1:4640903/1‑150 (MQ=255)

                                        cTACCACCACGTCGAGGATTTCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCg                                                                                                               <  1:3991440/150‑1 (MQ=255)

                                          accaccACGTCGAGGATTTCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCAACCCG‑CCCGACCGGGTGTAGGCGCACAGGCCGGTGTTTGACCCCggg                                                                                                             >  1:2868721/1‑150 (MQ=255)

                                              ccacGTCGAGGATTTCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCAc                                                                                                         >  1:982817/1‑150 (MQ=255)

                                              ccacGTCGAGGATTTCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCAc                                                                                                         >  1:2801147/1‑150 (MQ=255)

                                                   tCGAGGATTTCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACgggg                                                                                                                                             <  2:4835656/109‑1 (MQ=255)

                                                     gAGGATTTCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCt                                                                                                  <  1:4826709/150‑1 (MQ=255)

                                                      aGGATTTCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTg                                                                                                 <  2:3642899/150‑1 (MQ=255)

                                                        gATTTCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGgt                                                                                                                                            >  1:4919091/1‑105 (MQ=255)

                                                        gATTTCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGgt                                                                                                                                            <  2:4919091/105‑1 (MQ=255)

                                                         aTTTCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCTCGGCCGCTGGAc                                                                                              <  1:486869/150‑1 (MQ=255)

                                                         aTTTCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGAc                                                                                              <  1:938631/150‑1 (MQ=255)

                                                         aTTTCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGAc                                                                                              <  2:2272260/150‑1 (MQ=255)

                                                          tttCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGAcc                                                                                             >  1:2706006/1‑150 (MQ=255)

                                                          tttCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGAcc                                                                                             >  2:1164874/1‑150 (MQ=255)

                                                          tttCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGAcc                                                                                             >  1:662401/1‑150 (MQ=255)

                                                          tttCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGAcc                                                                                             <  1:80413/150‑1 (MQ=255)

                                                           ttCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCa                                                                                            <  2:1139358/150‑1 (MQ=255)

                                                            tCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAg                                                                                           <  2:3559103/150‑1 (MQ=255)

                                                             cAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAgg                                                                                          <  1:4536755/150‑1 (MQ=255)

                                                              atggaggTCGGCCATGGTGCGGGCATAGGGCACCTCCTCGTACATGGTCGGCGTTAGCGGCCTGTTGCCACGCACCGAATTCACCCCG‑CGAGACGGGGTGTAGGCGCACAGGCCGTTGTTTGACGCCGGGCCACGGCCGCTGGACCAGGt                                                                                         <  1:816635/148‑1 (MQ=255)

                                                              aaGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGCCGGGGTGTAGGCGCACAGGCCGTTGGTGGACGCCGGGCCACGGCCGCTGGACCAGGt                                                                                         >  2:5275581/1‑150 (MQ=255)

                                                              aaGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCATGt                                                                                         >  1:1274407/1‑150 (MQ=255)

                                                              aaGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGt                                                                                         >  1:1972515/1‑150 (MQ=255)

                                                              aaGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGt                                                                                         <  2:650562/150‑1 (MQ=255)

                                                                ggaggTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGttt                                                                                       >  2:4822613/1‑150 (MQ=255)

                                                                    gTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTttcttc                                                                                   <  1:1003388/150‑1 (MQ=255)

                                                                     tCGGCCATGGTGCGGGCATAGGGTACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCg                                                                                  >  2:4586955/1‑150 (MQ=255)

                                                                     tCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGGTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCg                                                                                  >  2:1349593/1‑150 (MQ=255)

                                                                          cATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCggccag                                                                             >  2:3303106/1‑150 (MQ=255)

                                                                              gTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCggccaggcca                                                                         >  1:4460962/1‑150 (MQ=255)

                                                                                   ggCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGgc                                                                    <  1:4768727/150‑1 (MQ=255)

                                                                                     cATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGgcgc                                                                  >  1:3493018/1‑150 (MQ=255)

                                                                                     cATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGgcgc                                                                  <  2:4118349/150‑1 (MQ=255)

                                                                                        aGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCggg                                                                                                             >  1:4099872/1‑104 (MQ=255)

                                                                                            caccaccacGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCAGGCCACGGCCGCTGGACCAg                                                                                           <  1:2683446/118‑1 (MQ=255)

                                                                                            caccaccacGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCAGGCCACGGCCGCTGGACCAg                                                                                           >  2:2683739/1‑118 (MQ=255)

                                                                                             accaccacGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTg                                                          >  1:3165265/1‑150 (MQ=255)

                                                                                              ccaccacGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTgg                                                         <  2:4456898/150‑1 (MQ=255)

                                                                                              ccaccacGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTgg                                                         <  1:3245364/150‑1 (MQ=255)

                                                                                               caccacGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTggc                                                        >  1:1650066/1‑150 (MQ=255)

                                                                                                 ccacGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACgg                                                                                                                                               <  2:2891880/61‑1 (MQ=255)

                                                                                                  cacGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTggcggc                                                     <  1:2923405/150‑1 (MQ=255)

                                                                                                     gTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTGGCggcggt                                                  <  2:969578/150‑1 (MQ=255)

                                                                                                     gTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCACTGAAGCTGGCggcggt                                                  <  2:2588908/150‑1 (MQ=255)

                                                                                                      tCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTGGCggcggtg                                                 <  2:2462780/150‑1 (MQ=255)

                                                                                                        cATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTGGCggcggtggc                                               >  1:2614805/1‑150 (MQ=255)

                                                                                                            gTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTGGCGGCGGTGGCGGTa                                           >  1:4732920/1‑150 (MQ=255)

                                                                                                             tCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTGGCGGCGGTGGCGGTAc                                          >  1:3908828/1‑150 (MQ=255)

                                                                                                             tCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTGGCGGCGGTGGCGGTAc                                          >  1:5105701/1‑150 (MQ=255)

                                                                                                              cGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTGGCGGCGGTGGCGGTAcc                                         <  2:3754560/150‑1 (MQ=255)

                                                                                                               ggCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTGGCGGCGGTGGCGGTACCg                                        <  2:4571046/150‑1 (MQ=255)

                                                                                                               ggCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTGGCGGCGGTGGCGGTACCg                                        <  2:4696572/150‑1 (MQ=255)

                                                                                                                  gTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTGGCGGCGGTGGCGGTACCGGCa                                     >  1:1935045/1‑150 (MQ=255)

                                                                                                                    tAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTGGCGGCGGTGGCGGTACCGGCAcc                                   <  1:2822447/150‑1 (MQ=255)

                                                                                                                      gCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTGGCGGCGGTGGCGGTACCGGCACCGt                                 >  2:1975387/1‑150 (MQ=255)

                                                                                                                       cGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTGGCGGCGGTGGCGGTACCGGCACCGtc                                <  1:2091007/150‑2 (MQ=255)

                                                                                                                             gTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTGGCGGCGGTGGCGGTACCGGCACCGTTGGACGa                          >  1:3118732/1‑150 (MQ=255)

                                                                                                                                gCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTGGCggcggtggcg                                              <  1:3103537/127‑1 (MQ=255)

                                                                                                                                 ccACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTGGCGGCGGTGGCGGTACCGGCACCGTTGGACGAGccc                      >  1:4270624/1‑150 (MQ=255)

                                                                                                                                   aCGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTGGCGGCGGTGGCGGTACCGGCACCGTTGGACGAGCCCGa                    >  1:4489648/1‑150 (MQ=255)

                                                                                                                                   aCGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTGGCGGCGGTGGCGGTACCGGCACCGTTGGACGAGCCCGa                    >  2:789515/1‑150 (MQ=255)

                                                                                                                                    cGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTGGCGGCGGTGGCGGTACCGGCACCGTTGGACGAGCCCGAg                   <  2:4460962/150‑1 (MQ=255)

                                                                                                                                     gCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTGGCGGCGGTGGCGGTACCGGCACCGTTGGACGAGCCCGAgg                  <  1:2399451/150‑1 (MQ=255)

                                                                                                                                     gCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTGGCGGCGGTGGCGGTACCGGCACCGTTGGACGAGCCCGAgg                  <  1:761572/150‑1 (MQ=255)

                                                                                                                                       aCCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTGGCGGCGGTGGCGGTACCGGCACCGTTGGACGAGCCCGAGGCg                >  2:1575481/1‑150 (MQ=255)

                                                                                                                                         cGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTGGCGGCGGTGGCGGTACCGGCACCGTTGGACGAGCCCGAGGCGaa              >  1:1037084/1‑150 (MQ=255)

                                                                                                                                            aaTCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTGGCGGCGGTGGCGGTACCGGCACCGTTGGACGAGCCCGAGGCGAAggg           <  1:3790595/150‑1 (MQ=255)

                                                                                                                                             tacACCCCGTCGCG‑CGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTGGCGGCGGTGGCGGTACCGGCACCGTTGGACGAGCCCGAGGCGAAGGgc          >  2:5416434/3‑150 (MQ=255)

                                                                                                                                             aTCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTGGCGGCGGTGGCGGTACCGGCACCGTTGGACGAGCCCGAGGCGAAGGgc          >  2:4019376/1‑150 (MQ=255)

                                                                                                                                               cACCCCGTCGCG‑CGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTGGCGGCGGTGGCGGTACCgg                                       >  1:3050333/1‑119 (MQ=255)

                                                                                                                                               cACCCCGTCGCG‑CGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTGGCGGCGGTGGCGGTACCgg                                       <  2:3050333/119‑1 (MQ=255)

                                                                                                                                               cACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAATGGCGCTGAAGCTGGCGGCGGTGGCGGTACCGGCACCGTTGGACGAGCCCGAGGCGAAGGgcgc        >  2:106126/1‑150 (MQ=255)

                                                                                                                                                aCCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTGGCGGCGGTGGCGGTACCGGCACCGTTGGACGAGCCCGAGGCGAAGGGcgcg       <  2:3593355/150‑1 (MQ=255)

                                                                                                                                                  cccG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTGGCGGCGGTGGCGGTACCGGCACCGTTGGACGAGCCCGa                    >  2:3374069/1‑135 (MQ=255)

                                                                                                                                                  cccG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTGGCGGCGGTGGCGGTACCGGCACCGTTGGACGAGCCCGAGGCGTAGGGCGCGGt     >  1:1817044/1‑150 (MQ=255)

                                                                                                                                                  cccG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTGGCGGCGGTGGCGGTACCGGCACCGTTGGACGAGCCCGAGGCGAAGGGCGCGGt     >  2:3230903/1‑150 (MQ=255)

                                                                                                                                                  cccG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTGGCGGCGGTGGCGGTACCGGCACCGTTGGACGAGCCCGAGGCGAAGGGCGCGGt     >  1:2200652/1‑150 (MQ=255)

                                                                                                                                                  cccG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTGGCGGCGGTGGCGGTACCGGCACCGTTGGACGAGCCCGAGGCGAAGGGCGCGGt     <  2:3118732/150‑1 (MQ=255)

                                                                                                                                                     g‑cgcgACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTGGCGGCGGTGGCGGTACCGGCACCGTTGGACGAGCCCGAGGCGAAGGGCGCGGTCAg  <  1:33466/150‑1 (MQ=255)

                                                                                                                                                     g‑cgcgACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTGGCGGCGGTGGCGGTACCGGCACCGTTGGACGAGCCCGAGGCGAAGGGCGCGGTCAg  <  1:878280/150‑1 (MQ=255)

TGCAGGCGCCACAGGTCGCCGCGCTTGTCGGCGTCATCGGCTACCACCACGTCGAGGATTTCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTGGCGGCGGTGGCGGTACCGGCACCGTTGGACGAGCCCGAGGCGAAGGGCGCGGTCAG  >  NC_002947/3336769‑3337067

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