BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,336,923	Δ1 bp	6.5%	coding (604/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,923	0	A	.	6.5%	606.7 / 34.8	185	coding (604/1704 nt)	PP_2932	amidase family protein
Reads supporting (aligned to +/- strand): ref base A (79/93); new base . (5/7); total (85/100)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.62e-01

Only 100 of 190 total aligned reads displayed.

GCGCCACAGGTCGCCGCGCTTGTCGGCGTCATCGGCTACCACCACGTCGAGGATTTCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTGGCGGCGGTGGCGGTACCGGCACCGTTGGACGAGCCCGAGGCGAAGGGCGCGGTCAGGTAGT  >  NC_002947/3336774‑3337072

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

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

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

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

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

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

            gccgcGCTTGTCGGCGTCATCGGCTACCACCACGTCGAGGATTTCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGgc                                                                                                                                           >  2:1033314/1‑150 (MQ=255)

            gccgcGCTTGTCGGCGTCATCGGCTACCACCACGTCGAGGATTTCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGgc                                                                                                                                           <  1:2405900/150‑1 (MQ=255)

                 gcTTGTCGGCGTCATCGGCTACCACCACGTCGAGGATTTCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGcaca                                                                                                                                      <  1:3613741/150‑1 (MQ=255)

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

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

                        ggCGTCATCGGCTACCACCACGTCGAGGATTTCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGGCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCgtt                                                                                                                               <  2:5247766/150‑1 (MQ=255)

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

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

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

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

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

                                 ggCTACCACCACGTCGAGGATTTCAAGGAGGTCGGCCATGGTGCGGGCATCGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGc                                                                                                                      >  1:2058279/1‑150 (MQ=255)

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

                                    tACCACCACGTCGAGGATTTCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCgg                                                                                                                   >  1:1293428/1‑150 (MQ=255)

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

                                           acGTCGAGGATTTCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACgg                                                                                                            <  2:226921/150‑1 (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)

                                                 aGGATTTCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTg                                                                                                      >  1:946591/1‑150 (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)

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

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

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

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

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

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

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

                                                             aggTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTttct                                                                                          <  1:3832593/150‑1 (MQ=255)

                                                               gTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGAGCAGGTttcttc                                                                                        <  2:4967155/150‑1 (MQ=255)

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

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

                                                                    ccATGGTGCGGGCATAGGGCACCACCACGCCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCggcca                                                                                   <  1:3133402/150‑1 (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)

                                                                                  tAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTg                                                                     <  2:3974516/150‑1 (MQ=255)

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

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

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

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

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

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

                                                                                            ccacGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACgg                                                                                                                                                    >  1:2891880/1‑61 (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)

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

                                                                                                    aTGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTGGCggcggtggcg                                                   <  2:3165265/150‑1 (MQ=255)

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

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

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

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

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

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

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

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

                                                                                                                  cGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTGGCGGCGGTGGCGGTACCGGCACCGtt                                     <  2:2932832/150‑1 (MQ=255)

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

                                                                                                                   ggCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTGGCGGCGGTGGCGGTACCGGCACCGTTg                                    <  1:907996/150‑1 (MQ=255)

                                                                                                                   ggCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTGGCGGCGGTGGCGGTACCGGCACCGTTg                                    >  2:2333686/1‑150 (MQ=255)

                                                                                                                    gCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTGGCGGCGGTGGCGGTACCGGCACCGTTgg                                   >  1:319080/1‑150 (MQ=255)

                                                                                                                      cAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCTAAGGCGCTGAAGCTGGCGGCGGTGGCGGTACCGGCACCGTTGGAc                                 <  2:1511462/150‑1 (MQ=255)

                                                                                                                      cAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGACGCTGGCGGCGGTGGCGGTACCGGCACCGTTGGAc                                 <  2:3401155/150‑1 (MQ=255)

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

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

                                                                                                                            ccACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTGGCGGCGGTGGCGGTACCGGCACCGTTGGACGAGccc                           <  2:2369053/150‑1 (MQ=255)

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

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

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

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

                                                                                                                              aCGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTGGCGGCGGTGGCGGTACCGGCACCGTTGGACGAGCCCGa                         <  1:82053/150‑1 (MQ=255)

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

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

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

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

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

                                                                                                                                         tCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCGAAGGCGCTGAAGCTGGCGGCGGTGGCGGTACCGGCACCGTTGGACGAGCCCGAGGCGAAGGgcg              >  1:3341322/1‑150 (MQ=255)

                                                                                                                                         tCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTGGCGGCGGTGGCGGTACCGGCACCGTTGGACGAGCCCGAGGCGAAGGgcg              <  2:2279967/150‑1 (MQ=255)

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

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

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

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

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

                                                                                                                                                  cgcgACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTGGCGGCGGTGGCGGTACCGGCACCGTTGGACGAGCCCGAGGCGAAGGGCGCGGTCAgg      <  1:65301/150‑1 (MQ=255)

                                                                                                                                                    cgACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTGGCGGCGGTGGCGGTACCGGCACCGTTGGACGAGCCCGAGGCGAAGGGCGCGGTCAGGTa    <  2:1293428/150‑1 (MQ=255)

                                                                                                                                                     gACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTGGCGGCGGTGGCGGTACCGGCACCGTTGGACGAGCCCGAGGCGAAGGGCGCGGTCAGGTAg   <  1:1824696/150‑1 (MQ=255)

                                                                                                                                                     gACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTGGCGGCGGTGGCGGTACCGGCACCGTTGGACGAGCCCGAGGCGAAGGGCGCGGTCAGGTAg   >  1:136291/1‑150 (MQ=255)

                                                                                                                                                      aCGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTGGCGGCGGTGGCGGTACCGGCACCGTTGGACGAGCCCGAGGCGAAGGGCGCGGTCAGGTAGt  <  2:4640903/150‑1 (MQ=255)

                                                                                                                                                      aCGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTGGCGGCGGTGGCGGTACCGGCACCGTTGGACGAGCCCGAGGCGAAGGGCGCGGTCAGGTAGt  >  2:4315046/1‑150 (MQ=255)

                                                                                                                                                      aCGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTGGCGGCGGTGGCGGTACCGGCACCGTTGGACGAGCCCGAGGCGAAGGGCGCGGTCAGGTAGt  <  1:1033314/150‑1 (MQ=255)

GCGCCACAGGTCGCCGCGCTTGTCGGCGTCATCGGCTACCACCACGTCGAGGATTTCAAGGAGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTGGCGGCGGTGGCGGTACCGGCACCGTTGGACGAGCCCGAGGCGAAGGGCGCGGTCAGGTAGT  >  NC_002947/3336774‑3337072

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