BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,629,794	A→C	8.0%	L229R (CTT→CGT)	PP_3195 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,629,794	0	A	C	8.0%	563.9 / 5.8	215	L229R (CTT→CGT)	PP_3195	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (112/78); new base C (2/21); total (114/101)
Fisher's exact test for biased strand distribution p-value = 2.93e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 217 total aligned reads displayed.

CCACTGCGCGCGGACAGTAGTAGCCGCCCAGGGAAACGCCTTCCAGGCCGATGCGCCTGGCATCGACATCAGCACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTGGTTGAATCCAGCCCGTTGACTTGCACCAGCAGCGGTGCCGGGCCACTCACGC  >  NC_002947/3629645‑3629941

ccACTGCGCGCGGACAGTAGTAGCCGCCCAGGGAAACGCCTTCCAGGCCGATGCGCCTGGCATCGACATCAGCACGGGGTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGaa                                                                                                                                                     <  2:5382688/150‑1 (MQ=255)

       gcgcgGACAGTAGTAGCCGCCCAGGGAAACGCCTTCCAGGCCGATGCGCCTGGCATCGACATCAGCACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTc                                                                                                                                              <  2:4149928/150‑1 (MQ=255)

       gcgcgGACAGTAGTAGCCGCCCAGGGAAACGCCTTCCAGGCCGATGCGCCTGGCATCGACATCAGCACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTc                                                                                                                                              <  2:4576446/150‑1 (MQ=255)

         gcgGACAGTAGTAGCCGCCCAGGGAAACGCCTTCCAGGCCGATGCGCCTGGCATCGACATCAGCACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGc                                                                                                                                            <  2:4320988/150‑1 (MQ=255)

          cgGACAGTAGTAGCCGCCCAGGGAAACGCCTTCCAGGCCGATGCGCCTGGCATCGACATCAGCACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGcc                                                                                                                                           <  2:3190155/150‑1 (MQ=255)

          cgGACAGTAGTAGCCGCCCAGGGAAACGCCTTCCAGGCCGATGCGCCTGGCATCGACATCAGCACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGcc                                                                                                                                           <  2:2244216/150‑1 (MQ=255)

          cgGACAGTAGTAGCCGCCCAGGGAAACGCCTTCCAGACCGATGCGCCTGGCATCGACATCAGCACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGcc                                                                                                                                           <  2:1350531/150‑1 (MQ=255)

           ggACAGTAGTAGCCGCCCAGGGAAACGCCTACCAGGCCGATGCGCCTGGCATCGACATCAGCACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCTCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCg                                                                                                                                          >  2:515166/1‑150 (MQ=255)

            gACAGTAGTAGCCGCCCAGGGAAACGCCTTCCAGGCCGATGCGCCTGGCATCGACATCAGCACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCgg                                                                                                                                         >  2:5103137/1‑150 (MQ=255)

               agtagtAGCCGCCCAGGGAAACGCCTTCCAGGCCGATGCGCCTGGCATCGACATCAGCACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGc                                                                                                                                      >  2:3144039/1‑150 (MQ=255)

                 tagtagCCGCCCAGGGAAACGCCTTCCAGGCCGATGCGCCTGGCATCGACATCAGCACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCg                                                                                                                                    >  2:5446523/1‑150 (MQ=255)

                        cgccCAGGGAAACGCCTTCCAGGCCGATGCGCCTGGCATCGACATCAGCACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTggg                                                                                                                             >  2:4224764/1‑149 (MQ=255)

                         gccCAGGGAAACGCCTTCCAGGCCGATGCGCCTGGCATCGACATCAGCACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTc                                                                                                                            >  2:3871201/1‑150 (MQ=255)

                         gccCAGGGAAACGCCTTCCAGGCCGATGCGCCTGGCATCGACATCAGCACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCGTCGCCGGTGGcgggcggggc                                                                                                                            >  1:2420916/1‑148 (MQ=255)

                          cccAGGGAAACGCCTTCCAGGCCGATGCGCCTGGCATCGACATCAGCACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCg                                                                                                                           <  2:4353079/150‑1 (MQ=255)

                             aGGGAAACGCCTTCCAGGCCGATGCGCCTGGCATCGACATCAGCACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGCAGCGCGTCGCCCGTGCCGGGCGGGTCGACa                                                                                                                        >  1:2861784/1‑150 (MQ=255)

                             aGGGAAACGCCTTCCAGGCCGATGCGCCTGGCATCGACATCAGCACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACa                                                                                                                        >  2:1288250/1‑150 (MQ=255)

                               ggAAACGCCTTCCAGGCCGATGCGCCTGGCATCGACATCAGCACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAAt                                                                                                                      >  2:4988025/1‑150 (MQ=255)

                                 aaaCGCCTTCCAGGCCGATGCGCCTGGCATCGACATCAGCACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGCAGCGCGTCGCCGGTGCCGGGCTGGTCGGCAATCa                                                                                                                    >  1:201029/1‑150 (MQ=255)

                                 aaaCGCCTTCCAGGCCGATGCGCCTGGCATCGACATCAGCACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCa                                                                                                                    <  2:3265457/150‑1 (MQ=255)

                                   ccGCCTTCCTGGCCGATGCGCCTGGAATCGACATCAGCCGGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCTGTGTTGGCTGTCATGGCGTGCCGTCATGCCGTGGTGGCGAAGCGCTCCCCCGGTGCCGGGCTGGTCGACAATCAGc                                                                                                                  <  1:2637300/149‑1 (MQ=255)

                                     gCCTTCCAGGCCGATGCGCCTGGCATCGACATCAGCACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGCAGCGCTTCGCCGGTGCCGGGCGGGTCGACAATAAGcga                                                                                                                >  1:658195/1‑150 (MQ=255)

                                     gCCTTCCAGGCCGATGCGCCTGGCATCGACATCAGCACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATAAGcga                                                                                                                >  1:4129968/1‑150 (MQ=255)

                                            aGGCCGATGCGCCTGGCATCGACATCAGCACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACAcg                                                                                                         >  1:4273187/1‑150 (MQ=255)

                                              gCCGATGCGCCTGGCATCGACATCAGCACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACAcgcc                                                                                                       >  2:4066500/1‑150 (MQ=255)

                                              gCCGATGCGCCTGGCATCGACATCAGCACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACAcgcc                                                                                                       >  1:4336098/1‑150 (MQ=255)

                                              gCCGATGCGCCTGGCATCGACATCAGCACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAAGCAGCGACGACAcgcc                                                                                                       >  2:4045556/1‑150 (MQ=255)

                                                     cgcCTGGCATCGACATCCGCTGGGGTTTCCAGCCAGTCCAGTTTGCGGCGGGCCCAGTGTTGGGTGTCTTGGGGTGCCGTCATGCGGTGGTGGCGTAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTg                                                                                                <  1:2559620/150‑1 (MQ=255)

                                                       ccTGGCATCGACATCAGCACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGgc                                                                                              <  2:1673276/150‑1 (MQ=255)

                                                        cTGGCATCGACATCAGCACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGGAGCGCTTAGCCGGTGCCGGGCGGGTCGACAATCAGCGACCACACGCCACGCCTGgcc                                                                                             >  1:4729412/1‑150 (MQ=255)

                                                         tGGCATCGACATCAGCACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGgcca                                                                                            <  2:2953201/150‑1 (MQ=255)

                                                             aTCGACATCAGCACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCACGCCTGgccagcca                                                                                        >  1:5520662/1‑150 (MQ=255)

                                                              tCGACATCCGCAGGGGTTTCCAGCCAGTCCAGTATGCGGCTGGCCAGGTTTGCGCTGTAATGGCGTGCCGTCATGACGAGGAGGCGAAGCTCTGCGCCGGTGCCGGGCTGGTCGCCAATCATCGACGACACGCATCGCCTGGccagccag                                                                                       <  2:3333576/150‑1 (MQ=255)

                                                                 aCATCAGCACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCGGGTCGACAATCAGCGACGACACGCCACGCCTGGCCACCcaggca                                                                                    >  1:4903278/1‑150 (MQ=255)

                                                                   acccgCACGGTTTTCCAGCCACGCCACTCGGCGGGGGGCCCTGTGGGCGGGGTCTTGGCGTGCCGACAGGCCGGGGCGGCGCCGCGCTGCCCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCcaggcagg                                                                                  <  1:531605/146‑1 (MQ=255)

                                                                   aTCAGCAGGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCAGGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCTGGCcaggcgag                                                                                  <  2:1388389/150‑4 (MQ=255)

                                                                    tCAGCACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGACGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCT‑GCCTGGCTGGCcaggcgaggc                                                                                <  2:3279717/150‑6 (MQ=255)

                                                                        cACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCCGTGTTCGCTGTCTTGGCGTGCCGTCATGCCGTGGAGGCGATGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCaggcaggcaggc                                                                             <  2:2659327/150‑1 (MQ=255)

                                                                         aCGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAAAAAGCCACGACACGCCACGCCAGGCCAGCCAGGCAGGCAGGcc                                                                            >  1:217691/1‑150 (MQ=255)

                                                                           gggTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGa                                                                          >  2:714131/1‑150 (MQ=255)

                                                                           gggTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGa                                                                          >  2:102921/1‑150 (MQ=255)

                                                                           gggTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGa                                                                          >  2:4547099/1‑150 (MQ=255)

                                                                           gggTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCGGGTCGACAAACAGCGACGACACCCCACGCCTGGCCAGCCACGCAGGCCGGCCGa                                                                          >  1:4832637/1‑150 (MQ=255)

                                                                               ttCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGACGCGCTTCGCCCGTGCCGGGCGGGGCGACCATCAGCGACGACACCCCTCGCCTGGCCCCCCCGGCCGGCAGGCCGAcgcg                                                                      >  1:3478414/1‑150 (MQ=255)

                                                                               ttCCACCCACTCCACTATGCGGCGGGCCCAGTGTTCGCTCTCACGGCGTGCCGTCATGCCGGGCAGGCTAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGAcgcg                                                                      <  1:3717359/150‑1 (MQ=255)

                                                                                  cagccagTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGCCGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTTGCCAGCCAGGCAGGCAGGCCGACGCGGTa                                                                   <  2:3446016/150‑1 (MQ=255)

                                                                                  cagccagTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTa                                                                   >  2:2876259/1‑150 (MQ=255)

                                                                                  cagccagTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTa                                                                   >  2:2791821/1‑150 (MQ=255)

                                                                                  cagccagTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCGTCGCCGGTGCCGGGCGGGTCGACAATCAGCGACGACACCCCGCGCCTGGCCCGCCCCGCAGGCAGGCCGCCGCGGTa                                                                   >  1:4817501/1‑150 (MQ=255)

                                                                                    gccagTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTatt                                                                 >  2:2519639/1‑150 (MQ=255)

                                                                                     ccagTCCACTTTGCGGCTGGCCCTGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGACGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTattt                                                                <  1:2329462/150‑1 (MQ=255)

                                                                                     ccagTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTattt                                                                <  2:2692887/150‑1 (MQ=255)

                                                                                      ccgTCCACTATGCGGCGGGCCCCGTGTCCGCGGACATGGCGGGCCGTCCGGCCGTGGCGGCCACGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTatttc                                                               <  1:4513191/148‑1 (MQ=255)

                                                                                      cagTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTctttc                                                               >  2:3268710/1‑150 (MQ=255)

                                                                                      cagTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGGatttc                                                               >  2:599856/1‑150 (MQ=255)

                                                                                      cagTCAACTAGGCGGCTGGCCCCGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGCGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTatttc                                                               <  1:5139359/150‑1 (MQ=255)

                                                                                       agTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTatttca                                                              <  2:4587034/150‑1 (MQ=255)

                                                                                        gTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTatttcat                                                             >  2:4739633/1‑150 (MQ=255)

                                                                                         tCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGCAGCGCGTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTatttcatt                                                            >  1:833351/1‑150 (MQ=255)

                                                                                         cccACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCGGCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTatttcatt                                                            <  1:2193274/149‑1 (MQ=255)

                                                                                            aCTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAATCGCTTCGCCGGTGCCGGGCGGGTCGACAAACACCGAAGACACGCCACGCCAGGCCAGCCAGGCAGGCAGGCCCACCCGGTatttaaattgc                                                         >  1:781207/1‑143 (MQ=255)

                                                                                            aCTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTTCCGGGCTGGTCGACAATCAGCGCCGACACCACTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTcc                                                         >  1:1527599/1‑150 (MQ=255)

                                                                                            aCTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCAGGCCAGCCAGGCAGGCAGGCCGACGCGGTCTTTCATTTcc                                                         >  1:2159910/1‑150 (MQ=255)

                                                                                            aCTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAAACAGCGACGACACGTCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTAATTCATTTcc                                                         >  1:3746286/1‑150 (MQ=255)

                                                                                              tttGCGGCGGGCCCTGTGTTCGCTGTCTGGGCGTGCCGTCATGCGGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCtt                                                       <  1:5103137/148‑1 (MQ=255)

                                                                                               aTGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTg                                                      <  2:1736154/150‑1 (MQ=255)

                                                                                                tGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGGGTTTCATGTCCTTgg                                                     >  2:2578956/1‑150 (MQ=255)

                                                                                                tGCGGCGGGCCCCGTGTTCGCTGTCAGGGCGTGCCGTCATGCCGTGGCGGCGACGCGCTCCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTgg                                                     <  1:1730883/150‑1 (MQ=255)

                                                                                                tGCGGCGGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCGTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTgg                                                     <  2:3111331/150‑1 (MQ=255)

                                                                                                  cggcggGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCGGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTGGtt                                                   <  2:3539869/150‑1 (MQ=255)

                                                                                                  cGGCTGGCCCACTGTTCGCTGGCATGGCGTCCCGTCAAGCCGCGGAGGCGAAGCGCTTCGGCGGTGCCGGGCTGGTCGACACTCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTGGtt                                                   <  2:5067264/150‑1 (MQ=255)

                                                                                                     cTGGCCCCGTGTTCGCTGTCTTGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTGGTTGaa                                                <  2:654282/150‑1 (MQ=255)

                                                                                                        gCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTGGTTGAATcc                                             >  1:1267978/1‑150 (MQ=255)

                                                                                                         cccAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTGGTTGAATCCa                                            <  2:5548614/150‑1 (MQ=255)

                                                                                                         cccAGTGTTCGCTGTCATGGCGAGCCGTCATGCCGGGGAGGCGAACCGCTTCGCCCCTGCCGGGCTCGTCGACAATCAGCCACCACACCCCACGCCTGGCCACCCAGGCAGGCAGCCCGACGCCGTAATTCATTACCAGGCTTGAAACCa                                            >  2:41848/1‑150 (MQ=255)

                                                                                                          cctgtgtTCGGTGTCAGGGCGTCCGGACAGGCGGGGGCGGCGCCGCGCGGCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTGGTTGAATCCAg                                           <  1:4036794/147‑1 (MQ=255)

                                                                                                           cAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCGGGGCGACAATCAGCGACGACACGCCTCGCCAGGCCAGCCACGCAGGCCGGCCGACGCGGTATTTCATTTCCTTGGTTGGATCCacc                                          >  1:5377045/1‑148 (MQ=255)

                                                                                                                ttCGCTGTCAGGGCGTGCCGTCAGGCCGTGGCGGCGACGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTGGTTGAATCCAGCCCGtt                                     <  2:184426/150‑1 (MQ=255)

                                                                                                                  ggcggTCTTGGGGTGCCGTCTTGCGGTGGCGGCGACGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTGGTTGAATCCAGCCCGTTGa                                   <  1:3451958/146‑1 (MQ=255)

                                                                                                                  cGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTGGTTGAATCCAGCCCGTTGa                                   <  2:4336098/150‑1 (MQ=255)

                                                                                                                  cGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTGGTTGAATCCAGCCCGTTGa                                   <  2:4697529/150‑1 (MQ=255)

                                                                                                                   gCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCGGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTGGTTGAATCCAGCCCGTTGAc                                  >  1:1457040/1‑150 (MQ=255)

                                                                                                                   gCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTAGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACCCCACGCCAGGCCAGCCACGCAGGCAGGCCGACGCGGTATTTCATTTCCTTGGTTGAATCCAGCCCGTTGAc                                  >  1:5035066/1‑150 (MQ=255)

                                                                                                                    cTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTGGTTGAATCCAGCCCGTTGACt                                 >  2:3817599/1‑150 (MQ=255)

                                                                                                                     tgccttGGGGTGCCGTCTGGCGGTGGAGGCGCCGCGCTCCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTGGTTGAATCCAGCCCGTTGACtt                                <  1:3144039/145‑1 (MQ=255)

                                                                                                                      gTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCCGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCGTGGTTGAATCCAGCCCGTTGACTTg                               >  2:1718770/1‑150 (MQ=255)

                                                                                                                       tCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTGGTTGAATCCAGCCCGTTGACTTGc                              <  2:1342157/150‑1 (MQ=255)

                                                                                                                          tGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCAGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTGGTTGAATCCAGCCCGTTGACTTGCAcc                           >  2:123997/1‑150 (MQ=255)

                                                                                                                                  cgtccggCCGGGGCGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTGGTTGAATCCAGCCCGTTGACTTGCACCAGCAGCgg                   <  1:2769570/144‑1 (MQ=255)

                                                                                                                                  cGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGGGCCGGGCTGGTCGACAAACAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTGGTTGAATCCAGCCCGTTGACTTGCACCAGCAGCgg                   >  1:1235260/1‑150 (MQ=255)

                                                                                                                                   gTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTGGTTGAATCCAGCCCGTTGACTTGCACCAGCAGCGGt                  >  2:1831238/1‑150 (MQ=255)

                                                                                                                                        gCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTGGTTGAATCCAGCCCGTTGACTTGCACCAGCAGCGGTGCCgg             <  2:4630119/150‑1 (MQ=255)

                                                                                                                                        gCCGGGGCGGCGCCGCGCGGCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTGGTTGAATCCAGCCCGTTGACTTGCACCAGCAGCGGTGCCgg             <  1:3035130/150‑1 (MQ=255)

                                                                                                                                              gAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTGGTTGAATCCAGCCCGTTGACTTGCACCAGCAGCGGTGCCGGGCCACt       >  2:2048624/1‑150 (MQ=255)

                                                                                                                                              gAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTGGTTGAATCCAGCCCGTTGACTTGCACCAGCAGCGGTGCCGGGCCACt       >  1:2503548/1‑150 (MQ=255)

                                                                                                                                              gAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTGGTTGAATCCAGCCCGTTGACTTGCACCAGCAGCGGTGCCGGGCCACt       >  2:3621077/1‑150 (MQ=255)

                                                                                                                                                 gCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTGGTTGAATCCAGCCCGTTGACTTGCACCAGCAGCGGTGCCGGGCCACTCAc    >  1:125850/1‑150 (MQ=255)

                                                                                                                                                   gAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTGGTTGAATCCAGCCCGTTGACTTGCACCAGCAGCGGTGCCGGGCCACTCACGc  <  2:4980252/150‑1 (MQ=255)

                                                                                                                                                   gAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTGGTTGAATCCAGCCCGTTGACTTGCACCAGCAGCGGTGCCGGGCCACTCACGc  >  2:211611/1‑150 (MQ=255)

                                                                                                                                                   gAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTGGTTGAATCCAGCCCGTTGACTTGCACCAGCAGCGGTGCCGGGCCACTCACGc  >  1:3118489/1‑150 (MQ=255)

CCACTGCGCGCGGACAGTAGTAGCCGCCCAGGGAAACGCCTTCCAGGCCGATGCGCCTGGCATCGACATCAGCACGGGTTTCCAGCCAGTCCACTATGCGGCTGGCCCAGTGTTCGCTGTCATGGCGTGCCGTCATGCCGTGGAGGCGAAGCGCTTCGCCGGTGCCGGGCTGGTCGACAATCAGCGACGACACGCCTCGCCTGGCCAGCCAGGCAGGCAGGCCGACGCGGTATTTCATTTCCTTGGTTGAATCCAGCCCGTTGACTTGCACCAGCAGCGGTGCCGGGCCACTCACGC  >  NC_002947/3629645‑3629941

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