BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,947,936	G→A	5.1%	A317V (GCG→GTG)	PP_3478 ←	secretion protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,947,936	0	G	A	5.1%	1016.9 / 26.2	315	A317V (GCG→GTG)	PP_3478	secretion protein
Reads supporting (aligned to +/- strand): ref base G (155/144); new base A (13/3); total (168/147)
Fisher's exact test for biased strand distribution p-value = 3.65e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

Only 100 of 337 total aligned reads displayed.

CATCGTCGCTCAGCACGCCGAAGGTGTTGGGCCATTGCAGGCCAAGGTCGAGGATGCGCGAGGTGGCCACCTCCATCACCTCCACTTCCAGCACCACCTCGGGGTTGGACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCATGACCAGTTCCTGGTAGTCC  >  NC_002947/3947787‑3948085

cATCGTCGCTCAGCACGCCGAAGGTGTTGGGCCATTGCAGGCCAAGGTCGAGGATGCGCGAGGTGGCCACCTCCATCACCTCCACTTCCAGCACCACCTCGGGGTTGGACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACcg                                                                                                                                                       >  1:5807761/1‑150 (MQ=255)

 aTCGTCGCTCAGCACGCCGAAGGTGTTGGGCCATTGCAGGCCAAGGTCGAGGATGCGCGAGGTGGCCACCTCCATCACCTCCACTTCCAGCACCACCTCGGGGTTGGACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACcgc                                                                                                                                                      >  2:1953846/1‑150 (MQ=255)

  tcgtcgCTCAGCACGCCGAAGGTGTTGGGCCATTGCAGGCCAAGGTCGAGGATGCGCGAGGTGGCCACCTCCATCACCTCCACTTCCAGCACCACCTCGGGGTTGGACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACcgcg                                                                                                                                                     >  2:3879820/1‑150 (MQ=255)

  tcgtcgCTCAGCACGCCGAAGGTGTTGGGCCATTGCAGGCCAAGGTCGAGGATGCGCGAGGTGGCCACCTCCATCACCTCCACTTCCAGCACCACCTCCGGGTTGGACTCGTCCTGCTACACCACCAGCTTCTCGGCCATGCGCACcgcg                                                                                                                                                     >  2:5731267/1‑150 (MQ=255)

   cgtcgCTCAGCACGCCGAAGGTGTTGGGCCATTGCAGGCCAAGGTCGAGGATGCGCGAGGTGGCCACCTCCATCACCTCCACTTCCAGCACCACCTCGGGGTTGGACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGt                                                                                                                                                    <  1:339237/150‑1 (MQ=255)

        cTCAGCACGCCGAAGGTGTTGGGCCATTGCAGGCCAAGGTCGAGGATGCGCGAGGTGGCCACCTCCATCACCTCCACTTCCAGCACCACCTCGGGGTTGGACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGAc                                                                                                                                               <  2:4760968/150‑1 (MQ=255)

                gCCGAAGGTGTTGGGCCATTGCAGGCCAAGGTCGAGGATGCGCGAGGTGGCCACCTCCATCACCTCCACTTCCAGCACCACCTCGGGGTTGGACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTcgcg                                                                                                                                       <  2:2374387/150‑1 (MQ=255)

                   gAAGGTGTTGGGCCATTGCAGGCCAAGGTCGAGGATGCGCGAGGTGGCCACCTCCATCACCTCCACTTCCAGCACCACCTCGGGGTTGGACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCAt                                                                                                                                    >  2:6871598/1‑150 (MQ=255)

                   gAAGGTGTTGGGCCATTGCAGGCCAAGGTCGAGGATGCGCGAGGTGGCCACCTCCATCACCTCCACTTCCAGCACCACCTCGGGGTTGGACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCAt                                                                                                                                    >  1:281294/1‑150 (MQ=255)

                    aaGGTGTTGGGCCATTGCAGGCCAAGGTCGAGGATGCGCGAGGTGGCCACCTCCATCACCTCCACTTCCAGCACCACCTCGGGGTTGGACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATg                                                                                                                                   >  2:4326059/1‑150 (MQ=255)

                      ggTGTTGGGCCATTGCAGGCCAAGGTCGAGGATGCGCGAGGTGGCCACCTCCATCACCTCCACTTCCAGCCCCACCTCGGGGTTGGACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGt                                                                                                                                 <  1:2365634/150‑1 (MQ=255)

                      ggTGTTGGGCCATTGCAGGCCAAGGTCGAGGATGCGCGAGGTGGCCACCTCCATCACCTCCACTTCCAGCACCACCTCGGGGTTGGACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCGACACGTCGGACGCGGTGCGCATggc                                                                                                                                 <  1:3875272/150‑2 (MQ=255)

                      ggTGTTGGGCCATTGCAGGCCAAGGTCGAGGATGCGCGAGGTGGCCACCTCCATCACCTCCACTTCCAGCACCACCTCGGGGTTGGACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGt                                                                                                                                 >  1:2947678/1‑150 (MQ=255)

                       gtgtTGGGCCATTGCAGGCCAAGGTCGAGGATGCGCGAGGTGGCCACCTCCATCACCTCCACTTCCAGCACCACCTCGGGGTTGGACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTc                                                                                                                                <  1:4844341/150‑1 (MQ=255)

                         gtTGGGCCATTGCAGGCCAAGGTCGAGGATGCGCGAGGTGGCCACCTCCATCACCTCCACTTCCAGCACCACCTCGGGGTTGGACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAg                                                                                                                              <  2:2295188/150‑1 (MQ=255)

                          ttGGGCCATTGCAGGCCAAGGTCGAGGATGCGCGAGGTGGCCACCTCCATCACCTCCACTTCCAGCACCACCTCGGGGTTGGACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAgg                                                                                                                             <  1:6118702/150‑1 (MQ=255)

                            gggCCATTGCAGGCCAAGGTCGAGGATGCGCGAGGTGGCCACCTCCATCACCTCCACTTCCAGCACCACCTCGGGGTTGGACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGgt                                                                                                                           <  1:352847/150‑1 (MQ=255)

                            gggCCATTGCAGGCCAAGGTCGAGGATGCGCGAGGTGGCCACCTCCATCACCTCCACTTCCAGCACCACCTCGGGGTTGGACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGgt                                                                                                                           >  1:5710282/1‑150 (MQ=255)

                             ggCCATTGCAGGCCAAGGTCGAGGATGCGCGAGGTGGCCACCTCCATCACCTCCACTTCCAGCACCACCTCGGGGTTGGACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGgtg                                                                                                                          <  2:2056010/150‑1 (MQ=255)

                              gCCATTGCAGGCCAAGGTCGAGGATGCGCGAGGTGGCCACCTCCATCACCTCCACTTCCAGCACCACCTCGGGGTTGGACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGgtgt                                                                                                                         >  2:6309271/1‑150 (MQ=255)

                               ccATTGCAGGCCAAGGTCGAGGATGCGCGAGGTGGCCACCTCCATCACCTCCACTTCCAGCACCACCTCGGGGTTGGACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGtt                                                                                                                        <  2:1421345/150‑1 (MQ=255)

                                 aTTGCAGGCCAAGGTCGAGGATGCGCGAGGTGGCCACCTCCATCACCTCCACTTCCAGCACCACCTCGGGGTTGGACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGa                                                                                                                      >  2:5129418/1‑150 (MQ=255)

                                 aTTGCAGGCCAAGGTCGAGGATGCGCGAGGTGGCCACCTCCATCACCTCCACTTCCAGCACCACCTCGGGGTTGGACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGa                                                                                                                      <  2:1029932/150‑1 (MQ=255)

                                  ttGCAGGCCAAGGTCGAGGATGCGCGAGGTGGCCACCTCCATCACCTCCACTTCCAGCACCACCTCGGGGTTGGACTGGTCCTTCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAg                                                                                                                     <  1:3442761/150‑1 (MQ=255)

                                  ttGCAGGCCAAGGTCGAGGATGCGCGAGGTGGCCACCTCCATCACCTCCACTTCCAGCACCACCTCGGGGTTGGACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAg                                                                                                                     <  1:817683/150‑1 (MQ=255)

                                  ttGCAGGCCAAGGTCGAGGATGCGCGAGGTGGCCACCTCCATCACCTCCACTTCCAGCACCACCTCGGGGTTGGACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAg                                                                                                                     >  2:5805491/1‑150 (MQ=255)

                                   tGCAGGCCAAGGTCGAGGATGCGCGAGGTGGCCACCTCCATCACCTCCCCTTCCAGCACCACCTCGGGGTTGGACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAgg                                                                                                                    <  2:1210194/150‑1 (MQ=255)

                                   tGCAGGCCAAGGTCGAGGATGCGCGAGGTGGCCACCTCCATCACCTCCACTTCCAGCACCACCTCGGGGTTGGACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAgg                                                                                                                    <  1:2548111/150‑1 (MQ=255)

                                      aGGCCAAGGTCGAGGATGCGCGAGGTGGCCACCTCCATCACCTCCACTTCCAGCACCACCTCGGGGTTGGACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGgcgc                                                                                                                 >  1:4595079/1‑150 (MQ=255)

                                        gCCAAGGTCGAGGATGCGCGAGGTGGCCACCTCCATCACCTCCACTTCCAGCACCACCTCGGGGTTGGACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTc                                                                                                               <  1:4944730/150‑1 (MQ=255)

                                         ccAAGGTCGAGGATGCGCGAGGTGGCCACCTCCATCACCTCCACTTCCAGCACCACCTCGGGGTTGGACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCg                                                                                                              >  1:5647828/1‑150 (MQ=255)

                                          cAAGGTCGAGGATGCGCGAGGTGGCCACCTCCGTCACCTCCACTTCCAGCACCACCTCGGGGTTGGACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGt                                                                                                             >  2:6766728/1‑150 (MQ=255)

                                          cAAGGTCGAGGATGCGCGAGGTGGCCACCTCCATCACCTCCACTTCCAGCACCACCTCGGGGTTGGACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGt                                                                                                             >  1:1457517/1‑150 (MQ=255)

                                           aaGGTCGAGGATGCGCGAGGTGGCCACCTCCATCACCTCCACTTCCAGCACCACCTCGGGGTTGGACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTc                                                                                                            <  1:2743786/150‑1 (MQ=255)

                                              gTCGAGGATGCGCGAGGTGGCCACCTCCATCACCTCCACTTCCAGCACCACCTCGGGGTTGGACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAAc                                                                                                         <  2:2865875/150‑1 (MQ=255)

                                                    gATGCGCGAGGTGGCCACCTCCATCACCTCCACTTCCAGCACCACCTCGGGGTTGGACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAacac                                                                                                   >  2:867075/1‑150 (MQ=255)

                                                        cgcgAGGTGGCCACCTCCATCACCTCCACTTCCAGCACCACCTCGGGGTTGGACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTcg                                                                                               <  2:4639467/150‑1 (MQ=255)

                                                         gcgAGGTGGCCACCTCCATCACCTCCACTTCCAGCACCACCTCGGGGTTGGACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTcgc                                                                                              >  1:1135728/1‑150 (MQ=255)

                                                         gcgAGGTGGCCACCTCCATCACCTCCACTTCCAGCACCACCTCGGGGTTGGACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTcgc                                                                                              >  1:229438/1‑150 (MQ=255)

                                                         gcgAGGTGGCCACCTCCATCACCTCCACTTCCAGCACCACCTCGGGGTTGGACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTcgc                                                                                              >  2:4053059/1‑150 (MQ=255)

                                                           gAGGTGGCCACCTCCATCACCTCCACTTCCAGCACCACCTCGGGGTTGGACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCgg                                                                                            <  1:5960003/150‑1 (MQ=255)

                                                           gAGGTGGCCACCTCCATCACCTCCACTTCCAGCACCACCTCGGGGTTGGACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCgg                                                                                            <  1:2862694/150‑1 (MQ=255)

                                                               tGGCCACCTCCATCACCTCCACTTCCAGCACCACCTCGGGGTTGGACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCt                                                                                        >  1:1607743/1‑150 (MQ=255)

                                                                ggCCACCTCCATCACCTCCACTTCCAGCACCACCTCGGGGTTGGACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCtt                                                                                       <  1:757237/150‑1 (MQ=255)

                                                                  ccACCTCCATCACCTCCACTTCCAGCACCACCTCGGGGTTGGACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGa                                                                                     <  1:6257794/150‑1 (MQ=255)

                                                                         cATCACCTCCACTTCCAGCACCACCTCGGGGTTGGACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGt                                                                              >  1:4330571/1‑150 (MQ=255)

                                                                          aTCACCTCCACTTCCAGCACCACCTCGGGGTTGGACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTc                                                                             <  1:6351960/150‑1 (MQ=255)

                                                                              ccTCCACTTCCAGCACCACCTCGGGGTTGGACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGa                                                                         <  2:5578587/150‑1 (MQ=255)

                                                                              ccTCCACTTCCAGCACCACCTCGGGGTTGGACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGa                                                                         <  2:5316183/150‑1 (MQ=255)

                                                                               cTCCACGTCCAGCACCACCTCGGGGTTGGACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGAc                                                                        <  1:6451171/150‑1 (MQ=255)

                                                                                tCCACTTCCAGCACCACCTCGGGGTTGGACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGAcc                                                                       >  1:4246974/1‑150 (MQ=255)

                                                                                tCCACTTCCAGCACCACCTCGGGGTTGGACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGAcc                                                                       <  2:4249206/150‑1 (MQ=255)

                                                                                   aCTTCCAGCACCACCTCGGGGTTGGACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATg                                                                    >  2:5691237/1‑150 (MQ=255)

                                                                                    cTTCCAGCACCACCTCGGGGTTGGACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGt                                                                   >  1:1048882/1‑150 (MQ=255)

                                                                                    cTTCCAGCACCACCTCGGGGTTGGACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGt                                                                   >  2:4493871/1‑150 (MQ=255)

                                                                                     ttCCAGCACCACCTCGGGGTTGGACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGtt                                                                  >  2:3765844/1‑150 (MQ=255)

                                                                                         aGCACCACCTCGGGGTTGGACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAAc                                                              <  2:3446650/150‑1 (MQ=255)

                                                                                          gCACCACCTCGGGGTTGGACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAAcg                                                             >  2:3731222/1‑150 (MQ=255)

                                                                                             ccaccTCGGGGTTGGACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCgg                                                          <  2:3197508/150‑1 (MQ=255)

                                                                                                 cTCGGGGTTGGACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGtt                                                      <  2:6476653/150‑1 (MQ=255)

                                                                                                  tCGGGGTTGGACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTg                                                     <  2:6166302/150‑1 (MQ=255)

                                                                                                      ggTTGGACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGt                                                 <  1:5259132/150‑1 (MQ=255)

                                                                                                       gTTGGACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTc                                                <  2:1497223/150‑1 (MQ=255)

                                                                                                       gTTGGACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTc                                                <  1:5085354/150‑1 (MQ=255)

                                                                                                          ggACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGAt                                             <  1:6853505/150‑1 (MQ=255)

                                                                                                           gACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATg                                            >  1:1175228/1‑150 (MQ=255)

                                                                                                            aCTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGc                                           >  1:4774064/1‑150 (MQ=255)

                                                                                                            aCTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTAGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGc                                           >  1:4669902/1‑150 (MQ=255)

                                                                                                              tGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTg                                         <  2:6073576/150‑1 (MQ=255)

                                                                                                                gTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGt                                       <  2:1153757/150‑1 (MQ=255)

                                                                                                                 tCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTc                                      <  2:3325838/150‑1 (MQ=255)

                                                                                                                  ccTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCa                                     <  2:6584516/150‑1 (MQ=255)

                                                                                                                   cTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAg                                    <  2:5879091/150‑1 (MQ=255)

                                                                                                                     gCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGt                                  >  1:3761287/1‑150 (MQ=255)

                                                                                                                          tGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAAc                             >  2:6214617/1‑150 (MQ=255)

                                                                                                                              gcagcTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTgc                         <  2:6004282/150‑1 (MQ=255)

                                                                                                                              gcagcTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTgc                         <  1:3279958/150‑1 (MQ=255)

                                                                                                                                 gcTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCa                      >  2:1675610/1‑150 (MQ=255)

                                                                                                                                 gcTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCa                      <  2:5016957/150‑1 (MQ=255)

                                                                                                                                  cTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCAt                     <  1:6309271/150‑1 (MQ=255)

                                                                                                                                  cTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCAt                     >  1:5358786/1‑150 (MQ=255)

                                                                                                                                  cTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCAt                     <  1:1778240/150‑1 (MQ=255)

                                                                                                                                     ctcGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCATGAc                  >  2:454991/1‑150 (MQ=255)

                                                                                                                                     ctcGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCATGAc                  <  2:281294/150‑1 (MQ=255)

                                                                                                                                     ctcGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCATGAc                  >  2:4407395/1‑150 (MQ=255)

                                                                                                                                     cctgaccATGCGCGACACGTCGGACGCGGTGCGCATGGCCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCATGAc                  >  1:4363348/6‑150 (MQ=38)

                                                                                                                                      tcGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCATGAcc                 >  1:4026651/1‑150 (MQ=255)

                                                                                                                                      ctgaccATGCGCGACACGTCGGACGCGGTGCGCATGGCCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCATGAcc                 >  2:772705/5‑150 (MQ=38)

                                                                                                                                       tgaccATGCGCGACACGTCGGACGCGGTGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCATGACCa                >  1:2153026/4‑150 (MQ=38)

                                                                                                                                       tgaccATGCGCGACACGTCGGACGCGGTGCGCATGGCCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCCTGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCATGACCa                >  1:2105655/4‑150 (MQ=38)

                                                                                                                                       cGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCATGACCa                <  1:3661932/150‑1 (MQ=255)

                                                                                                                                       cGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCATGACCa                <  2:909231/150‑1 (MQ=255)

                                                                                                                                          ccATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCATGACCAGtt             >  2:3648613/1‑150 (MQ=255)

                                                                                                                                            aTGCGCGACACGTCGGACGCGGTGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCATGACCAGTTcc           >  1:417695/1‑150 (MQ=255)

                                                                                                                                              gcgcACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCATGACCAGTTCCTg         <  1:1953846/150‑1 (MQ=255)

                                                                                                                                                 cACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCATGACCAGTTCCTGGTa      <  2:357577/150‑1 (MQ=255)

                                                                                                                                                 cACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCAGGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCATGACCAGTTCCTGGTa      <  1:6529696/150‑1 (MQ=255)

                                                                                                                                                    cgcgTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCATGACCAGTTCCTGGTAGTc   >  2:4892778/1‑150 (MQ=255)

                                                                                                                                                     gcgTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCATGACCAGTTCCTGGTAGTcc  <  2:1846245/150‑1 (MQ=255)

                                                                                                                                                     gcgTCCGACGTGCCGCGCAGGGTCAGGGTGTTGCGGCGCTCGTCCACGAACACGTCGCGGGACTGGCGCATGGTCTGGACCAGGGTCAACGCGGTGTTGGCGTCGATGCTGGGCAGGTAGAACGTGCGCATGACCAGTTCCTGGTAGTcc  <  2:3158471/150‑1 (MQ=255)

CATCGTCGCTCAGCACGCCGAAGGTGTTGGGCCATTGCAGGCCAAGGTCGAGGATGCGCGAGGTGGCCACCTCCATCACCTCCACTTCCAGCACCACCTCGGGGTTGGACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCATGACCAGTTCCTGGTAGTCC  >  NC_002947/3947787‑3948085

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