BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,335,001	Δ1 bp	6.0%	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,335,001	0	A	.	6.0%	364.3 / 16.6	117	coding (604/1704 nt)	PP_2932	amidase family protein
Reads supporting (aligned to +/- strand): ref base A (60/50); new base . (3/4); total (63/54)
Fisher's exact test for biased strand distribution p-value = 7.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.11e-01

Only 100 of 117 total aligned reads displayed.

CACAGGTCGCCGCGCTTGTCGGCGTCATCGGCTACCACCACGTCGAGGATTTCAAGG‑AGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTGGCGGCGGTGGCGGTACCGGCACCGTTGGACGAGCCCGAGGCGAAGGGCGCGGTCAGG  >  NC_002947/3334856‑3335146

cctgtggcgccgcGCTTGTCGGCGTCATCGGCTACCACCACGTCGAGGATTTCAAGG‑AGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACggg                                                                                                                                               <  1:1463167/143‑1 (MQ=255)

 acaGGTCGCCGCGCTTGTCGGCGTCATCGGCTACCACCACGTCGAGGATTTCAAGG‑AGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACgg                                                                                                                                                >  1:57048/1‑148 (MQ=255)

 acaGGTCGCCGCGCTTGTCGGCGTCATCGGCTACCACCACGTCGAGGATTTCAAGG‑AGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACgg                                                                                                                                                <  2:57048/148‑1 (MQ=255)

    ggTCGCCGCGCTTGTCGGCGTCATCGGCTACCACCACGTCGAGGATTTCAAGG‑AGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCGTCGCG‑CGGGgtga                                                                                                                                           <  1:2242882/150‑2 (MQ=255)

    ggTCGCCGCGCTTGTCGGCGTCATCGGCTACCACCACGTCGAGGATTTCAAGG‑AGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCGTCGCG‑CGGGgtga                                                                                                                                           <  2:1013758/150‑2 (MQ=255)

       cgccgcGCTTTTCGGCGTCATCGGCTACCACCACGTCGAGGATTTCAAGG‑AGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAgg                                                                                                                                        <  1:1583605/150‑1 (MQ=255)

       cgccgcGCTTGTCGGCGTCATCGGCTCCCACCCACTCGCGG‑TTTCGAGGTATGTCGGCCATGGTGCGGGAATAGGGCACCACCACGTCCATGGCCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCGTCGCG‑CGGGgtgattt                                                                                                                                        <  1:540460/150‑5 (MQ=255)

        gccgcGCTTGTCGGCGTCATCGGCTACCACCACGTCGAGGATTTCAAGG‑AGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACcgggcgccggc                                                                                                                                       >  2:1859122/1‑143 (MQ=255)

        gccgcGCTTGTCGGCGTCATCGGCTACCACCACGTCGAGGATTTCAAGG‑AGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGgc                                                                                                                                       >  1:320202/1‑150 (MQ=255)

        gccgcGCTTGTCGGCGTCATCGGCTACCACCACGTCGAGGATTTCAAGG‑AGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGgc                                                                                                                                       >  1:1557880/1‑150 (MQ=255)

               ttGTCGGCGTCATCGGCTACCACCACGTCGAGGATTTCAAGG‑AGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAgg                                                                                                                                <  2:608510/150‑1 (MQ=255)

               ttGTCGGCGTCATCGGCTACCACCACGTCGAGGATTTCAAGG‑AGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAgg                                                                                                                                >  1:140606/1‑150 (MQ=255)

                  tCGGCGTCATCGGCTACCACCACGTCGAGGATTTCAAGG‑AGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCg                                                                                                                             >  1:1079132/1‑150 (MQ=255)

                    ggCGTCATCGGCTACCACCACGTCGAGGATTTCAAGG‑AGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCgtt                                                                                                                           <  1:2357910/150‑1 (MQ=255)

                       gTCATCGGCTACCACCACGTCGAGGATTTCAAGG‑AGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCgttgtt                                                                                                                        <  2:320202/150‑1 (MQ=255)

                        tCATCGGCTACCACCACGTCGAGGATTTCAAGG‑AGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGGACAGGCCGttgttg                                                                                                                       <  2:1353898/150‑1 (MQ=255)

                          aTCGGCTACCACCACGTCGAGGATTTCAAGG‑AGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGa                                                                                                                     >  1:475774/1‑150 (MQ=255)

                           gcGGCTACCACCACGTCGAGGATTTCAAGG‑AGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGGAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGAc                                                                                                                    >  2:2211434/2‑150 (MQ=255)

                             ggCTACCACCACGTCGAGGATTTCAAGG‑AGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGc                                                                                                                  <  2:2384179/150‑1 (MQ=255)

                               cTACCACCACGTCGAGGATTTCAAGG‑AGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGANANCACCCCN‑CGCGACGGGGTGNAGGCGCACAGGCCGTTGTTGGACGNCg                                                                                                                <  1:1798634/150‑1 (MQ=255)

                               cTACCACCACGGCGAGGATTTCAAGG‑AGGTCGGCCATGGTGCGGGCCTAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGTCAGTTTCCACGCACCGAAATCACCCCG‑CGCGACGGGGAGTAGGCGGACAGGCCGTTGTTTGACGCCg                                                                                                                >  2:427317/1‑150 (MQ=255)

                                  ccaccaCGTCGAGGATTTCAAGG‑AGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGc                                                                                                             >  1:1418065/1‑150 (MQ=255)

                                   caccacGTCGAGGATTTCAAGG‑AGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGcc                                                                                                            >  2:1497219/1‑150 (MQ=255)

                                    accacGTCGAGGATTTCAAGG‑AGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCa                                                                                                           >  2:1416091/1‑150 (MQ=255)

                                     ccacGTCGAGGATTTCAAGG‑AGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCAc                                                                                                          >  2:2184836/1‑150 (MQ=255)

                                            gAGGATTTCAAGG‑AGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCt                                                                                                   >  1:1890139/1‑150 (MQ=255)

                                             aGGATTTCAAGG‑AGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAg                                                                                                                                         >  1:952433/1‑111 (MQ=255)

                                                 tttCAAGG‑AGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCgttgtt                                                                                                                        >  2:152945/1‑124 (MQ=255)

                                                   tCAAGG‑AGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAg                                                                                            <  2:864078/150‑1 (MQ=255)

                                                   tCAAGG‑AGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAg                                                                                            <  1:973859/150‑1 (MQ=255)

                                                    cAAGG‑AGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGTGCCACGGCCGCTGGACCAgg                                                                                           <  2:348663/150‑1 (MQ=255)

                                                    cAAGG‑AGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAgg                                                                                           <  1:1055800/150‑1 (MQ=255)

                                                     aaGG‑AGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGt                                                                                          >  1:1790909/1‑150 (MQ=255)

                                                            gTCGGCCATGGTGCGGGCATTGGGCACCACCACGTCCATGGTCGGCGTTCGCGGGCAGTTGCCACGCACCGCAGTCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTttcttc                                                                                    <  2:442553/150‑1 (MQ=255)

                                                             tCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGGAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCg                                                                                   <  1:2211434/150‑1 (MQ=255)

                                                               ggCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGGGGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCggc                                                                                 >  1:1468820/1‑150 (MQ=255)

                                                               ggCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCggg                                                                                                              >  2:950758/1‑121 (MQ=255)

                                                                 ccATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCggcca                                                                               >  1:2211394/1‑150 (MQ=255)

                                                                    tGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCggccaggc                                                                            <  1:996414/150‑1 (MQ=255)

                                                                        gCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGCCCAGGTTTCTTCGGCCAGGCCaaa                                                                        >  2:814041/1‑150 (MQ=255)

                                                                          gggCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAgg                                                                      >  1:138656/1‑150 (MQ=255)

                                                                            gCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTATTCGGCCAGGCCAAAGgcg                                                                    <  1:1588064/150‑1 (MQ=255)

                                                                             cATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCGCGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGgcgc                                                                   >  1:2362191/1‑150 (MQ=255)

                                                                              aTAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCt                                                                  <  2:1079132/150‑1 (MQ=255)

                                                                              aTAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTCTTGGCGCACAGGCCGTTGTTGGACGCCGAGCCACGGCCGCTTTACCAGGTTTCTTCGGCCAGGCCAAAGGCGCt                                                                  >  1:1966137/1‑150 (MQ=255)

                                                                                aGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGa                                                                >  1:1087024/1‑150 (MQ=255)

                                                                                aGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGa                                                                >  2:557227/1‑150 (MQ=255)

                                                                                  ggCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGgtgt                                                                                                                                           <  1:1057904/73‑1 (MQ=255)

                                                                                   gCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGttgttg                                                                                                                       >  2:1676244/1‑92 (MQ=255)

                                                                                   gCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGttgttg                                                                                                                       <  1:1676244/92‑1 (MQ=255)

                                                                                      ccaccacGTCCATGGTCGGCGTTAGCGGCCGGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTgg                                                          <  1:769921/150‑1 (MQ=255)

                                                                                       caccacGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCCGGCCACAGGCGCTGAAGCTggc                                                         >  2:1489088/1‑150 (MQ=255)

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

                                                                                        accacGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTggcg                                                        >  1:614615/1‑150 (MQ=255)

                                                                                         ccacGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCggccaggc                                                                            >  2:1811497/1‑129 (MQ=255)

                                                                                         ccacGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCggccaggc                                                                            <  1:1811497/129‑1 (MQ=255)

                                                                                         ccacGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTggcgg                                                       >  1:1177829/1‑150 (MQ=255)

                                                                                          cacGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGtt                                                                                         >  2:1109722/1‑115 (MQ=255)

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

                                                                                           acGTCCATGGTCGGCGTTAGCGGTCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGGAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCCGCCAGGCCCAGATCGCTGAAGCTGgcggcg                                                     >  1:1329763/1‑150 (MQ=255)

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

                                                                                               ccATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGa                                                                                                <  1:1013254/103‑1 (MQ=255)

                                                                                               ccATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGa                                                                                                >  2:1013254/1‑103 (MQ=255)

                                                                                                     tCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTGGCGGCGGTGGCGGTAc                                           <  1:230672/150‑1 (MQ=255)

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

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

                                                                                                                ggCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCTAAAGGCGCTGAAGCTGGCGGCGGTGGCGGTACCGGCACCGTTg                                <  2:1790909/150‑1 (MQ=255)

                                                                                                                  ccAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTGGCGGCGGTGGCGGTACCgg                                        >  2:1351332/1‑140 (MQ=255)

                                                                                                                  ccAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTGGCGGCGGTGGCGGTACCgg                                        <  1:1351332/140‑1 (MQ=255)

                                                                                                                  ccAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTGGCGGCGGTGGCGGTACCGGCACCGTTGGa                              >  1:1804127/1‑150 (MQ=255)

                                                                                                                  ccAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTGGCGGCGGTGGCGGTACCAGCACCGTTGGa                              <  2:925656/150‑1 (MQ=255)

                                                                                                                   cAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTGGCGGCGGTGGCGGTACCGGCACCGTTGGAc                             >  1:827805/1‑150 (MQ=255)

                                                                                                                     gTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTttcttc                                                                                    <  1:1032688/93‑1 (MQ=255)

                                                                                                                     gTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTttcttc                                                                                    >  2:1032688/1‑93 (MQ=255)

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

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

                                                                                                                             gCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGACGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTGGCGGCGGTGGCGGTACCGGCACCGTTGGACGAGCCCGAgg                   >  2:1733406/1‑150 (MQ=255)

                                                                                                                              cACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGTACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTGGCGGCGGTGGCGGTACCGGCACCGTTGGACGAGCCCGAGGc                  >  1:306136/1‑150 (MQ=255)

                                                                                                                              cACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTGGCGGCGGTGGCGGTACCGGCACCGTTGGACGAGCCCGAGGc                  <  1:799930/150‑1 (MQ=255)

                                                                                                                              cACCAAAATCTCCCGA‑CGCTAGTGGTTATAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTGGCGGCGGTGGCGGTACCGGCACCGTTGGACGAGCCCGAGGc                  <  2:140606/150‑1 (MQ=255)

                                                                                                                               aCCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTGGCggcggtggc                                                <  1:711788/119‑1 (MQ=255)

                                                                                                                               aCCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTGGCggcggtggc                                                >  2:711788/1‑119 (MQ=255)

                                                                                                                                   aaaTCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTGGCGGCGGTGGCGGTACCGGCACCGTTGGACGAGCCCGAGGCGAAgg             <  1:904717/150‑1 (MQ=255)

                                                                                                                                    ctacACCCCGTCGCG‑CGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTGGCGGCGGTGGCGGTACCGGCACCGTTGGACGAGCCCGAGGCGAAggg            >  2:1967580/4‑150 (MQ=255)

                                                                                                                                    ctacACCCCGTCGCG‑CGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTGGCGGCGGTGGCGGTACCGGCACCGTTGGACGAGCCCGAGGCGAAggg            >  2:1952023/4‑150 (MQ=255)

                                                                                                                                    ctacACCCCGTCGCG‑CGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTGGCGGCGGTGGCGGTACCGGCACCGTTGGACGAGCCCGAGGCGAAggg            >  2:9483/4‑150 (MQ=255)

                                                                                                                                    aaTCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTGGCGGCGGTGGCGGTACCGGCACCGTTGGACGAGCCCGAGGCGAAggg            >  2:1182164/1‑150 (MQ=255)

                                                                                                                                     aTCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTGGCGGCGGTGGCGGTACCgg                                        <  2:660322/121‑1 (MQ=255)

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

                                                                                                                                      tCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTGGCGGCGGTGGCGGTACCGGCACCGTTGGACGAGCCCGAGGCGAAGGgcg          >  1:1704303/1‑150 (MQ=255)

                                                                                                                                        aCCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCCAAGGCGCTGAAGCTGGCGGCGGTGCCGGTACCGCCAcc                                    <  2:1391746/122‑1 (MQ=255)

                                                                                                                                        aCCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTGGCGGCGGTGCCGGTACCGCCAcc                                    >  1:1391746/1‑122 (MQ=255)

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

                                                                                                                                           ccG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCggcc                                                                                <  2:1739959/75‑1 (MQ=255)

                                                                                                                                           ccG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCggcc                                                                                >  1:1739959/1‑75 (MQ=255)

                                                                                                                                            cg‑cgcgACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTGGCGGCGGTGGCGGTACCGGCACCGTTGGACGAGCCCGAg                    <  2:233534/134‑1 (MQ=255)

                                                                                                                                            cg‑cgcgACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTGGCGGCGGTGGCGGTACCGGCACCGTTGGACGAGCCCGAg                    >  1:233534/1‑134 (MQ=255)

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

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

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

CACAGGTCGCCGCGCTTGTCGGCGTCATCGGCTACCACCACGTCGAGGATTTCAAGG‑AGGTCGGCCATGGTGCGGGCATAGGGCACCACCACGTCCATGGTCGGCGTTAGCGGCCAGTTGCCACGCACCGAAATCACCCCG‑CGCGACGGGGTGTAGGCGCACAGGCCGTTGTTGGACGCCGGGCCACGGCCGCTGGACCAGGTTTCTTCGGCCAGGCCAAAGGCGCTGAAGCTGGCGGCGGTGGCGGTACCGGCACCGTTGGACGAGCCCGAGGCGAAGGGCGCGGTCAGG  >  NC_002947/3334856‑3335146

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