BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,638,543	C→G	5.3%	G167G (GGC→GGG)	lepB →	signal peptidase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,638,543	0	C	G	5.3%	1042.0 / 33.9	303	G167G (GGC→GGG)	lepB	signal peptidase I
Reads supporting (aligned to +/- strand): ref base C (164/123); new base G (3/13); total (167/136)
Fisher's exact test for biased strand distribution p-value = 3.48e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

Only 100 of 309 total aligned reads displayed.

AAGGTCATCGAGGTCGGTGACCCGCAACGCGGGGATGTAATGGTGTTCCGCTACCCGAGTGACCCGAATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGATCCGCAAGGAAATGAGCCGTTATCGCATGCCGCCGGACCAGCAGTGGACCGTGCC  >  NC_002947/1638394‑1638692

  ggTCATCGAGGTCGGTGACCCGCAACGCGGGGATGTAATGGTGTTCCGCTACCCGAGTGACCCGAATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCa                                                                                                                                                     >  1:599162/1‑150 (MQ=255)

      aTCGAGGTCGGTGACCCGCAACGCGGGGATGTAATGGTGTTCCGCTACCCGAGTGACCCGAATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGccc                                                                                                                                                 >  1:5868090/1‑150 (MQ=255)

      aTCGAGGTCGGTGACCCGCAACGCGGGGATGTAATGGTGTTCCGCTACCCGAGTGACCCGAATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGccc                                                                                                                                                 >  1:6752989/1‑150 (MQ=255)

      aTCGAGGTCGGTGACCCGCAACGCGGGGATGTAATGGTGTTCCGCTACCCGAGTGACCCGAATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGccc                                                                                                                                                 >  2:6050935/1‑150 (MQ=255)

      aTCGAGGTCGGTGACCCGCAACGCGGGGATGTAATGGTGTTCCGCTACCCGAGTGACCCGAATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGccc                                                                                                                                                 >  1:1617756/1‑150 (MQ=255)

         gAGGTCGGTGACCCGCAACGCGGGGATGTAATGGTGTTCCGCTACCCGAGTGACCCGAATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGtt                                                                                                                                              <  1:3407717/150‑1 (MQ=255)

           ggTCGGTGACCCGCAACGCGGGGATGTAATGGTGTTCCGCTACCCGAGTGACCCGAATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGc                                                                                                                                            <  1:2336161/150‑1 (MQ=255)

            gTCGGTGACCCGCAACGCGGGGATGTAATGGTGTTCCGCTACCCGAGTGACCCGAATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGcc                                                                                                                                           >  1:1448977/1‑150 (MQ=255)

            gTCGGTGACCCGCAACGCGGGGATGTAATGGTGTTCCGCTACCCGAGTGACCCGAATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGcc                                                                                                                                           <  1:7267014/150‑1 (MQ=255)

            gTCGGTGACCCGCAACGCGGGGATGTAATGGTGTTCCGCTACCCGAGTGACCCGAATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGcc                                                                                                                                           <  2:477106/150‑1 (MQ=255)

            gTCGGTGACCCGCAACGCGGGGATGTAATGGTGTTCCGCTACCCGAGTGACCCGAATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGcc                                                                                                                                           >  2:1236380/1‑150 (MQ=255)

              cGGTGACCCGCAACGCGGGGATGTAATGGTGTTCCGCTACCCGAGTGACCCGAATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGTTGTTCGGCAACGGGCTGGCCGTTGACGa                                                                                                                                         <  2:7376095/150‑1 (MQ=255)

              cGGTGACCCGCAACGCGGGGATGTAATGGTGTTCCGCTACCCGAGTGACCCGAATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGGCAACGGGCTGGCCGTTGACGa                                                                                                                                         <  1:1701568/150‑1 (MQ=255)

               ggTGACCCGCAACGCGGGGATGTAATGGTGTTCCGCTACCCGAGTGACCCGAATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGaa                                                                                                                                        >  1:3608194/1‑150 (MQ=255)

               ggTGACCCGCAACGCGGGGATGTAATGGTGTTCCGCTACCCGAGTGACCCGAATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGaa                                                                                                                                        <  1:3042926/150‑1 (MQ=255)

                     ccGCAACGCGGGGATGTAATGGTGTTCCGCTACCCGAGTGACCCGAATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGc                                                                                                                                            <  2:2025681/140‑1 (MQ=255)

                     ccGCAACGCGGGGATGTAATGGTGTTCCGCTACCCGAGTGACCCGAATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGc                                                                                                                                            >  1:2025681/1‑140 (MQ=255)

                     ccGCAACGCGGGGATGTAATGGTGTTCCGCTACCCGAGTGACCCGAATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACtg                                                                                                                                  >  2:6754274/1‑150 (MQ=255)

                      cGCAACGCGGGGATGTAATGGTGTTCCGCTACCCGAGTGACCCGAATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACtgg                                                                                                                                 >  1:6869952/1‑150 (MQ=255)

                        cAACGCGGGGATGTAATGGTGTTCCGCTACCCGAGTGACCCGAATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACtggtg                                                                                                                               <  1:5348242/150‑1 (MQ=255)

                        cAACGCGGGGATGTAATGGTGTTCCGCTACCCGAGTGACCCGAATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACtggtg                                                                                                                               <  1:6559860/150‑1 (MQ=255)

                        cAACGCGGGGATGTAATGGTGTTCCGCTACCCGAGTGACCCGAATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGGCAACGGGCTGGCCGTTGACGaacagccgctt                                                                                                                               <  2:1205324/150‑8 (MQ=38)

                         aaCGCGGGGATGTAATGGTGTTCCGCGCCCCGAGTGACCCGCATGTCAACTCCCGCCAGCGAGGGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCGAGCCCGTTTCCGAACTACtggtgg                                                                                                                              <  2:5956267/150‑1 (MQ=255)

                         aaCGCGGGGATGTAATGGGGTTCCGCTACCCGAGTGACCCGAATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACtggtgg                                                                                                                              <  1:850732/150‑1 (MQ=255)

                                ggATGTAATGGTGTTCCGCTACCCGAGTGACCCGAATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGa                                                                                                                       <  2:6973892/150‑1 (MQ=255)

                                    gTAATGGTGTTCCGCTACCCGAGTGACCCGAATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCg                                                                                                                   >  2:4428018/1‑150 (MQ=255)

                                    gTAATGGTGTTCCGCTACCCGAGTGACCCGAATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCg                                                                                                                   <  2:6357646/150‑1 (MQ=255)

                                     tAATGGTGTTCCGCTACCCGAGTGACCCGAATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCgg                                                                                                                  <  2:3802382/150‑1 (MQ=255)

                                     tAATGGTGTTCCGCTACCCGAGTGACCCGAATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCgg                                                                                                                  <  1:1413527/150‑1 (MQ=255)

                                      aaTGGTGTTCCGCTACCCGAGTGACCCGAATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCCGCGCGCCggg                                                                                                                 >  1:1029200/1‑150 (MQ=255)

                                       aTGGTGTTCCGCTACCCGAGTGACCCGAATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACTGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGt                                                                                                                <  2:3250787/150‑1 (MQ=255)

                                           tgtTCCGCTACCCGAGTGACCCGAATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAAAAACTGGTGGGCAGCGAGCCGGGTAccc                                                                                                            >  2:399397/1‑150 (MQ=255)

                                             ttCCGCTACCCGAGTGACCCGAATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTg                                                                                                          <  1:5891157/150‑1 (MQ=255)

                                               ccGCTACCCGAGTGACCCGAATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTggg                                                                                                        <  1:3508006/150‑1 (MQ=255)

                                               ccGCTACCCGAGTGACCCGAATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTggg                                                                                                        <  2:4082630/150‑1 (MQ=255)

                                               ccGCTACCCGAGTGACCCGAATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTg                                                                                                          >  1:4101761/1‑148 (MQ=255)

                                                 gCTACCCGAGTGACCCGAATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCa                                                                                                      <  1:108256/150‑1 (MQ=255)

                                                 gCTACCCGAGTGACCCGAATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCa                                                                                                      <  2:3621100/150‑1 (MQ=255)

                                                 gCTACCCGAGTGACCCGAATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCa                                                                                                      <  1:3280990/150‑1 (MQ=255)

                                                   tACCCGAGTGACCCGAATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGt                                                                                                    <  2:1448977/150‑1 (MQ=255)

                                                      ccGAGTGACCCGAATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCCGCCCGTTGCCGAACAACAGGTGGGCAGCGAGCCGGGTACCCAGGGCAGTGCt                                                                                                 >  1:7204493/1‑150 (MQ=255)

                                                       cGAGTGACCCGAATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTg                                                                                                          <  1:5724115/140‑1 (MQ=255)

                                                       cGAGTGACCCGAATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTg                                                                                                          >  2:5724115/1‑140 (MQ=255)

                                                        gAGTGACCCGAATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCa                                                                                               <  2:6875186/150‑1 (MQ=255)

                                                        gAGTGACCCGAATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGACCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCa                                                                                               <  2:4912305/150‑1 (MQ=255)

                                                         aGTGACCCGAATGTCAACTACATCAAGCGTGAGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAg                                                                                              <  2:2766778/150‑1 (MQ=255)

                                                         aGTGACCCGAATGTCAACTACAGCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAg                                                                                              <  2:7146754/150‑1 (MQ=255)

                                                              cccGAATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTa                                                                                         <  2:599162/150‑1 (MQ=255)

                                                              cccGAATGTCAACGACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTa                                                                                         <  2:5027195/150‑1 (MQ=255)

                                                               ccGAATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCGGGGCAGTGCTCAGCTGTAc                                                                                        >  1:987406/1‑150 (MQ=255)

                                                                cGAATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACa                                                                                       >  1:6195256/1‑150 (MQ=255)

                                                                cGAATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACAGGTGGGCAGCGAGCCGGGTACCCAGGGCAGTGCTCAGCAGTACa                                                                                       >  1:5052910/1‑150 (MQ=255)

                                                                 gAATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACaa                                                                                      <  1:5498958/150‑1 (MQ=255)

                                                                 gAATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACaa                                                                                      >  2:7052096/1‑150 (MQ=255)

                                                                 gAATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACaa                                                                                      <  2:4223858/150‑1 (MQ=255)

                                                                      tCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGaaa                                                                                 <  2:3270783/150‑1 (MQ=255)

                                                                        aaCTACATCACGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAg                                                                               <  2:5508916/150‑1 (MQ=255)

                                                                        aaCTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTAcc                                                                                                             <  1:6223580/120‑1 (MQ=255)

                                                                        aaCTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAg                                                                               <  1:1339232/150‑1 (MQ=255)

                                                                        aaCTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAg                                                                               >  1:294206/1‑150 (MQ=255)

                                                                              aTCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAACCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGgc                                                                         >  1:1088509/1‑150 (MQ=255)

                                                                                 aaGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGcgcg                                                                      >  2:1548638/1‑150 (MQ=255)

                                                                                  aGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCgg                                                                     <  1:7052096/150‑1 (MQ=255)

                                                                                     gtgtGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCCGCGAGCCGGGTACCCGGGGCAGTGCGCAGCGGTACAAGGAAAAGCTCGGCGCGGCCg                                                                  >  1:33390/1‑150 (MQ=255)

                                                                                       gtgGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTg                                                                                                          >  2:2697341/1‑108 (MQ=255)

                                                                                       gtgGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTg                                                                                                          <  1:2697341/108‑1 (MQ=255)

                                                                                         ggTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAaca                                                              >  2:5415583/1‑150 (MQ=255)

                                                                                          gTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAacac                                                             >  2:6418108/1‑150 (MQ=255)

                                                                                              gCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGa                                                         <  1:1572878/150‑1 (MQ=255)

                                                                                              gCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGa                                                         >  2:1086929/1‑150 (MQ=255)

                                                                                               ccTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGAt                                                        <  1:2274719/150‑1 (MQ=255)

                                                                                                cTGCCGGGCGACGTGGTCCGCTACACCAACGACTAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGATc                                                       >  2:3411137/1‑150 (MQ=255)

                                                                                                   ccGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCCGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGATCCGc                                                    >  1:3848126/1‑150 (MQ=255)

                                                                                                   ccGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGATCCGc                                                    >  2:2546786/1‑150 (MQ=255)

                                                                                                     gggCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCCACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGACCACCTGCTCCGCac                                                  >  2:2493478/1‑149 (MQ=255)

                                                                                                     gggCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGATCCGCaa                                                  <  2:2213519/150‑1 (MQ=255)

                                                                                                     gggCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGATCCGCaa                                                  <  1:6931596/150‑1 (MQ=255)

                                                                                                     gggCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACCGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGATCCGCac                                                  >  2:6242647/1‑149 (MQ=255)

                                                                                                             tGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGATCCGCAAGGAAATGa                                          <  2:1840881/150‑1 (MQ=255)

                                                                                                               gTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGATCCGCAAGGAAATGAGc                                        >  2:809881/1‑150 (MQ=255)

                                                                                                                tCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGATCCGCAAGGAAATGAGcc                                       >  2:3751465/1‑150 (MQ=255)

                                                                                                                tCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGATCCGCAAGGAAATGAGcc                                       >  2:5114218/1‑150 (MQ=255)

                                                                                                                 ccGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGATCCGCAAGGAAATGAGCCg                                      >  1:2846762/1‑150 (MQ=255)

                                                                                                                       cacCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACTACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGATCCGCAAGGAAATGAGCCGTTATCg                                <  1:6449866/150‑1 (MQ=255)

                                                                                                                         ccAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGATCCGCAAGGAAATGAGCCGTTATCGCa                              >  1:4584043/1‑150 (MQ=255)

                                                                                                                          cAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGATCCGCAAGGAAATGAGCCGTTATCGCAt                             >  1:5204996/1‑150 (MQ=255)

                                                                                                                           aaCGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGATCCGCAAGGAAATGAGCCGTTATCGCATg                            >  2:336626/1‑150 (MQ=255)

                                                                                                                               aCAAGCGGCTGTTCGTCAACGGCCAGCCCCTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCAGATCCGCAAGGAAATGAGCCGTTATCGCATgccgc                        >  1:754135/1‑150 (MQ=255)

                                                                                                                                 aaGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGATCCGCAAGGAAATGAGCCGTTATCGCATGccgccg                      <  1:4448200/150‑1 (MQ=255)

                                                                                                                                 aaGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACAGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCACAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGAACCGCAAGGAAATGAGCCGTTATCGCATGcccccg                      >  1:4479121/1‑150 (MQ=255)

                                                                                                                                        tGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGATCCGCAAGGAAATGAGCCGTTATCGCATGCCGCCGGACcagc               >  2:535560/1‑150 (MQ=255)

                                                                                                                                        tGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGATCCGCAAGGAAATGAGCCGTTATCGCATGCCGCCGGACcagc               >  2:3146371/1‑150 (MQ=255)

                                                                                                                                         gTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGATCCGCAAGGAAATGAGCCGTTATCGCATGCCGCCGGACcagca              <  2:3608194/150‑1 (MQ=255)

                                                                                                                                            cGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGATCCGCAAGGAAATGAGCCGTTATCGCATGCCGCCGGACCAGCAGTg           >  1:4717573/1‑150 (MQ=255)

                                                                                                                                                aaCGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGATCCGCAAGGAAATGAGCCGTTATCGCATGCCGCCGGACCAGCAGTGGAcc       >  2:872754/1‑150 (MQ=255)

                                                                                                                                                aaCGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGATCCGCAAGGAAATGAGCCGTTATCGCATGCCGCCGGACCAGCAGTGGAcc       >  2:4630175/1‑150 (MQ=255)

                                                                                                                                                 aCGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGATCCGCAAGGAAATGAGCCGTTATCGCATGCCGCCGGACCAGCAGTGGACCg      <  1:6708683/150‑1 (MQ=255)

                                                                                                                                                  cGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGATCCGCAAGGAAATGAGCCGTTATCGCATGCCGCCGGACCAGCAGTGGACCGt     >  1:5176701/1‑150 (MQ=255)

                                                                                                                                                    gCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGATCCGCAAGGAAATGAGCCGTTATCGCATGCCGCCGGACCAGCAGTGGACCGTGc   >  1:4303813/1‑150 (MQ=255)

                                                                                                                                                    gCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGATCCGCAAGGAAATGAGCCGTTATCGCATGCCGCCGGACCAGCAGTGGACCGTGc   <  2:3018406/150‑1 (MQ=255)

AAGGTCATCGAGGTCGGTGACCCGCAACGCGGGGATGTAATGGTGTTCCGCTACCCGAGTGACCCGAATGTCAACTACATCAAGCGTGTGGTTGGCCTGCCGGGCGACGTGGTCCGCTACACCAACGACAAGCGGCTGTTCGTCAACGGCCAGCCCGTTGCCGAACAACTGGTGGGCAGCGAGCCGGGTACCCTGGGCAGTGCTCAGCTGTACAAGGAAAAGCTCGGCGCGGCCGAACACCTGATCCGCAAGGAAATGAGCCGTTATCGCATGCCGCCGGACCAGCAGTGGACCGTGCC  >  NC_002947/1638394‑1638692

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